Mercurial > repos > artbio > lumpy_smoove

comparison vcf2hrdetect.py @ 6:0dc7ea74bf61 draft default tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/lumpy_smoove commit 5f62b58c65917afbb018296568c81ad132106012
author artbio
date Wed, 05 Apr 2023 15:55:36 +0000
parents 41704ce96c7f
children
comparison
equal deleted inserted replaced
5:d06cb88d6849 6:0dc7ea74bf61
24 else: 24 else:
25 tabdict[newid] = [vcfdict[id][0], vcfdict[id][1], 25 tabdict[newid] = [vcfdict[id][0], vcfdict[id][1],
26 chrom, coordbreak, "TRA"] 26 chrom, coordbreak, "TRA"]
27 for id in list(vcfdict): 27 for id in list(vcfdict):
28 if "_" in id: 28 if "_" in id:
29 del(vcfdict[id]) 29 del vcfdict[id]
30 for id in vcfdict.keys(): # only sv that are not of type TRA or INV 30 for id in vcfdict.keys(): # only sv that are not of type TRA or INV
31 chr1 = vcfdict[id][0] 31 chr1 = vcfdict[id][0]
32 chr2 = vcfdict[id][0] 32 chr2 = vcfdict[id][0]
33 pos1 = vcfdict[id][1] 33 pos1 = vcfdict[id][1]
34 pos2 = vcfdict[id][7].split("END=")[1].split(";")[0] 34 pos2 = vcfdict[id][7].split("END=")[1].split(";")[0]