13
|
1 ##fileformat=VCFv4.1
|
|
2 ##fileDate=2014-05-05 16:18:46
|
|
3 ##source=SoftSearch.pl
|
|
4 ##reference=chr1.small.fa
|
|
5 ##Usage= SoftSearch.pl -l 5 -q 20 -r 5 -d 300 -m 5 -u 0 -s 6 -b Small_test.bam -f chr1.small.fa -o s_Test.vcf
|
|
6 ##INFO=<ID=SVTYPE,Number=1,Type=String,Description="Type of structural variant">
|
|
7 ##INFO=<ID=EVENT,Number=1,Type=String,Description="ID of event associated to breakend">
|
|
8 ##INFO=<ID=MATEID,Number=.,Type=String,Description="ID of mate breakends">
|
|
9 ##INFO=<ID=ISIZE,Number=.,Type=String,Description="Size of the SV">
|
|
10 ##INFO=<ID=END,Number=1,Type=Integer,Description="End position of the variant described in this record">
|
|
11 ##FORMAT=<ID=lSC,Number=1,Type=Integer,Description="Length of the longest soft clips supporting the BND">
|
|
12 ##FORMAT=<ID=nSC,Number=1,Type=Integer,Description="Number of supporting soft-clips">
|
|
13 ##FORMAT=<ID=uRP,Number=1,Type=Integer,Description="Number of unmated read pairs nearby Soft-Clips">
|
|
14 ##FORMAT=<ID=levD_local,Number=1,Type=Float,Description="Levenstein distance between soft-clipped bases and the area around the original soft-clipped site">
|
|
15 ##FORMAT=<ID=levD_distl,Number=1,Type=Float,Description="Levenstein distance between the soft-clipped bases and mate location">
|
|
16 ##FORMAT=<ID=CTX,Number=1,Type=Integer,Description="Number of chromosomal translocations">
|
|
17 ##FORMAT=<ID=DEL,Number=1,Type=Integer,Description="Number of reads supporting Large Deletions">
|
|
18 ##FORMAT=<ID=INS,Number=1,Type=Integer,Description="Number of reads supporting Large insertions">
|
|
19 ##FORMAT=<ID=NOV_INS,Number=1,Type=Integer,Description="Number of reads supporting novel sequence insertion">
|
|
20 ##FORMAT=<ID=TDUP,Number=1,Type=Integer,Description="Number of reads supporting a tandem duplication">
|
|
21 ##FORMAT=<ID=INV,Number=1,Type=Integer,Description="Number of reads supporting inversions">
|
|
22 ##FORMAT=<ID=sDEL,Number=1,Type=Integer,Description="Number of reads supporting novel sequence insertion">
|
|
23 ##INFO=<ID=NO_MATE_SC,Number=1,Type=Flag,Description="When there is no softclipping of the mate read location, an appromiate position is used">
|
|
24 ##FORMAT=<ID=GT,Number=1,Type=String,Description="Dummy value for maintaining VCF-Spec">
|
|
25 #CHROM POS ID REF ALT QUAL FILTER INFO FORMAT s_Test
|
|
26 chr1 524004 BND_cvfiwqny T T[chr1:525966[ 2.56 PASS SVTYPE=BND;EVENT=DEL;END=525966;ISIZE=1962;MATEID=BND_xfnhcrbo GT:CTX:DEL:INS:INV:NOV_INS:TDUP:lSC:nSC:uRP:distl_levD 0/1:2:16:0:0:0:0:0:0:0:0
|
|
27 chr1 525966 BND_xfnhcrbo G ]chr1:524004]G 2.56 PASS SVTYPE=BND;EVENT=DEL;END=524004;ISIZE=1962;MATEID=BND_cvfiwqny GT:CTX:DEL:INS:INV:NOV_INS:TDUP:lSC:nSC:uRP:distl_levD 0/1:2:16:0:0:0:0:31:7:0:0.39
|
|
28 chr1 773527 BND_aznzkzaf T T[chr1:773527[ 1000.00 PASS SVTYPE=BND;EVENT=INV;END=781624;ISIZE=8097;MATEID=BND_cpmhcuwq GT:CTX:DEL:INS:INV:NOV_INS:TDUP:lSC:nSC:uRP:distl_levD 0/1:1:0:0:60:4:0:0:0:0:0
|
|
29 chr1 773527 BND_jybykdcg T A[chr1:773527[ 1000.00 PASS SVTYPE=BND;EVENT=INV;END=781565;ISIZE=8038;MATEID=BND_wverdlsg GT:CTX:DEL:INS:INV:NOV_INS:TDUP:lSC:nSC:uRP:distl_levD 0/1:1:0:0:67:6:0:0:0:0:0
|
|
30 chr1 773904 BND_rycfkqeh G ]chr1:773904]G 2.43 PASS SVTYPE=BND;EVENT=INV;END=781564;ISIZE=7660;MATEID=BND_dasorwlz GT:CTX:DEL:INS:INV:NOV_INS:TDUP:lSC:nSC:uRP:distl_levD 0/1:19:0:0:76:6:0:32:11:6:0.41
|
|
31 chr1 773907 BND_yudcuvey C ]chr1:773907]T 4.74 PASS SVTYPE=BND;EVENT=INV;END=781107;ISIZE=7200;MATEID=BND_zwtirrcp GT:CTX:DEL:INS:INV:NOV_INS:TDUP:lSC:nSC:uRP:distl_levD 0/1:19:0:0:100:12:0:33:7:12:0.21
|
|
32 chr1 773928 BND_dmqfcmpy C A[chr1:773928[ 16.39 PASS SVTYPE=BND;EVENT=INV;END=781562;ISIZE=7634;MATEID=BND_wsczqaaz GT:CTX:DEL:INS:INV:NOV_INS:TDUP:lSC:nSC:uRP:distl_levD 0/1:3:0:0:46:5:0:0:0:0:0
|
|
33 chr1 773963 BND_lcdckfeq T ]chr1:773963]G 2.77 PASS SVTYPE=BND;EVENT=INV;END=781564;ISIZE=7601;MATEID=BND_befaawsf GT:CTX:DEL:INS:INV:NOV_INS:TDUP:lSC:nSC:uRP:distl_levD 0/1:19:0:0:76:6:0:36:10:6:0.36
|
|
34 chr1 773967 BND_wggxcsyf C ]chr1:773967]T 1000.00 PASS SVTYPE=BND;EVENT=INV;END=781107;ISIZE=7140;MATEID=BND_uzdmoyno GT:CTX:DEL:INS:INV:NOV_INS:TDUP:lSC:nSC:uRP:distl_levD 0/1:18:0:0:95:8:0:18:5:8:0.00
|
|
35 chr1 781107 BND_uzdmoyno T C[chr1:781107[ 1000.00 PASS SVTYPE=BND;EVENT=INV;END=773967;ISIZE=7140;MATEID=BND_wggxcsyf GT:CTX:DEL:INS:INV:NOV_INS:TDUP:lSC:nSC:uRP:distl_levD 0/1:18:0:0:95:8:0:0:0:0:0
|
|
36 chr1 781107 BND_zwtirrcp T C[chr1:781107[ 4.74 PASS SVTYPE=BND;EVENT=INV;END=773907;ISIZE=7200;MATEID=BND_yudcuvey GT:CTX:DEL:INS:INV:NOV_INS:TDUP:lSC:nSC:uRP:distl_levD 0/1:19:0:0:100:12:0:0:0:0:0
|
|
37 chr1 781562 BND_wsczqaaz A ]chr1:781562]C 16.39 PASS SVTYPE=BND;EVENT=INV;END=773928;ISIZE=7634;MATEID=BND_dmqfcmpy GT:CTX:DEL:INS:INV:NOV_INS:TDUP:lSC:nSC:uRP:distl_levD 0/1:3:0:0:46:5:0:17:5:5:0.06
|
|
38 chr1 781564 BND_befaawsf G T[chr1:781564[ 2.77 PASS SVTYPE=BND;EVENT=INV;END=773963;ISIZE=7601;MATEID=BND_lcdckfeq GT:CTX:DEL:INS:INV:NOV_INS:TDUP:lSC:nSC:uRP:distl_levD 0/1:19:0:0:76:6:0:0:0:0:0
|
|
39 chr1 781564 BND_dasorwlz G G[chr1:781564[ 2.43 PASS SVTYPE=BND;EVENT=INV;END=773904;ISIZE=7660;MATEID=BND_rycfkqeh GT:CTX:DEL:INS:INV:NOV_INS:TDUP:lSC:nSC:uRP:distl_levD 0/1:19:0:0:76:6:0:0:0:0:0
|
|
40 chr1 781565 BND_wverdlsg A ]chr1:781565]T 1000.00 PASS SVTYPE=BND;EVENT=INV;END=773527;ISIZE=8038;MATEID=BND_jybykdcg GT:CTX:DEL:INS:INV:NOV_INS:TDUP:lSC:nSC:uRP:distl_levD 0/1:1:0:0:67:6:0:34:12:6:0.00
|
|
41 chr1 781624 BND_cpmhcuwq T ]chr1:781624]T 1000.00 PASS SVTYPE=BND;EVENT=INV;END=773527;ISIZE=8097;MATEID=BND_aznzkzaf GT:CTX:DEL:INS:INV:NOV_INS:TDUP:lSC:nSC:uRP:distl_levD 0/1:1:0:0:60:4:0:24:12:4:0.00
|
|
42 chr1 813168 BND_lpicqgph G G[chr1:813168[ 1000.00 PASS SVTYPE=BND;EVENT=INV;END=816792;ISIZE=3624;MATEID=BND_ecpjudes GT:CTX:DEL:INS:INV:NOV_INS:TDUP:lSC:nSC:uRP:distl_levD 0/1:1:0:0:49:5:0:0:0:0:0
|
|
43 chr1 813648 BND_hbwudayx G ]chr1:813648]T 1000.00 PASS SVTYPE=BND;EVENT=INV;END=816343;ISIZE=2695;MATEID=BND_uwtkigle GT:CTX:DEL:INS:INV:NOV_INS:TDUP:lSC:nSC:uRP:distl_levD 0/1:0:0:0:44:3:0:18:10:3:0.00
|
|
44 chr1 813648 BND_wulykabc G G[chr1:813648[ 2.49 PASS SVTYPE=BND;EVENT=INV;END=816790;ISIZE=3142;MATEID=BND_cpobomjn GT:CTX:DEL:INS:INV:NOV_INS:TDUP:lSC:nSC:uRP:distl_levD 0/1:0:0:0:56:7:0:0:0:0:0
|
|
45 chr1 816343 BND_uwtkigle T G[chr1:816343[ 1000.00 PASS SVTYPE=BND;EVENT=INV;END=813648;ISIZE=2695;MATEID=BND_hbwudayx GT:CTX:DEL:INS:INV:NOV_INS:TDUP:lSC:nSC:uRP:distl_levD 0/1:0:0:0:44:3:0:0:0:0:0
|
|
46 chr1 816790 BND_cpobomjn G ]chr1:816790]G 2.49 PASS SVTYPE=BND;EVENT=INV;END=813648;ISIZE=3142;MATEID=BND_wulykabc GT:CTX:DEL:INS:INV:NOV_INS:TDUP:lSC:nSC:uRP:distl_levD 0/1:0:0:0:56:7:0:40:10:7:0.40
|
|
47 chr1 816792 BND_ecpjudes G ]chr1:816792]G 1000.00 PASS SVTYPE=BND;EVENT=INV;END=813168;ISIZE=3624;MATEID=BND_lpicqgph GT:CTX:DEL:INS:INV:NOV_INS:TDUP:lSC:nSC:uRP:distl_levD 0/1:1:0:0:49:5:0:33:16:5:0.00
|
|
48 chr1 854784 BND_migefcff G G[chr1:863642[ 4.52 PASS SVTYPE=BND;EVENT=DEL;END=863642;ISIZE=8858;MATEID=BND_csfcfrly GT:CTX:DEL:INS:INV:NOV_INS:TDUP:lSC:nSC:uRP:distl_levD 0/1:0:99:0:0:9:0:0:0:0:0
|
|
49 chr1 854784 BND_ramkiwhh G G[chr1:863640[ 2.93 PASS SVTYPE=BND;EVENT=DEL;END=863640;ISIZE=8856;MATEID=BND_gcuiwknc GT:CTX:DEL:INS:INV:NOV_INS:TDUP:lSC:nSC:uRP:distl_levD 0/1:0:98:0:0:9:0:0:0:0:0
|
|
50 chr1 855228 BND_xaptufoh C ]chr1:863644]C 2.49 PASS SVTYPE=BND;EVENT=DEL;END=863644;ISIZE=8416;MATEID=BND_ipasaubx GT:CTX:DEL:INS:INV:NOV_INS:TDUP:lSC:nSC:uRP:distl_levD 0/1:0:104:0:0:8:0:30:18:8:0.40
|
|
51 chr1 863640 BND_gcuiwknc C ]chr1:854784]C 2.93 PASS SVTYPE=BND;EVENT=DEL;END=854784;ISIZE=8856;MATEID=BND_ramkiwhh GT:CTX:DEL:INS:INV:NOV_INS:TDUP:lSC:nSC:uRP:distl_levD 0/1:0:98:0:0:9:0:35:10:9:0.34
|
|
52 chr1 863642 BND_csfcfrly G ]chr1:854784]G 4.52 PASS SVTYPE=BND;EVENT=DEL;END=854784;ISIZE=8858;MATEID=BND_migefcff GT:CTX:DEL:INS:INV:NOV_INS:TDUP:lSC:nSC:uRP:distl_levD 0/1:0:99:0:0:9:0:18:6:9:0.22
|
|
53 chr1 863644 BND_ipasaubx T T[chr1:855228[ 2.49 PASS SVTYPE=BND;EVENT=DEL;END=855228;ISIZE=8416;MATEID=BND_xaptufoh GT:CTX:DEL:INS:INV:NOV_INS:TDUP:lSC:nSC:uRP:distl_levD 0/1:0:104:0:0:8:0:0:0:0:0
|
|
54 chr1 903148 BND_xbqyvnjg G C[chr1:903148[ 1000.00 PASS SVTYPE=BND;EVENT=INV;END=915991;ISIZE=12843;MATEID=BND_puxffwzz GT:CTX:DEL:INS:INV:NOV_INS:TDUP:lSC:nSC:uRP:distl_levD 0/1:0:0:0:64:3:0:0:0:0:0
|
|
55 chr1 903615 BND_wpzdzdca T ]chr1:903615]C 3.12 PASS SVTYPE=BND;EVENT=INV;END=915991;ISIZE=12376;MATEID=BND_akjluuhn GT:CTX:DEL:INS:INV:NOV_INS:TDUP:lSC:nSC:uRP:distl_levD 0/1:0:0:0:60:11:0:19:7:11:0.32
|
|
56 chr1 903617 BND_jtvtjzwt C ]chr1:903617]C 1000.00 PASS SVTYPE=BND;EVENT=INV;END=915564;ISIZE=11947;MATEID=BND_oonwkzas GT:CTX:DEL:INS:INV:NOV_INS:TDUP:lSC:nSC:uRP:distl_levD 0/1:0:0:0:64:4:0:29:14:4:0.00
|
|
57 chr1 903617 BND_ysqifxup C G[chr1:903617[ 3.83 PASS SVTYPE=BND;EVENT=INV;END=915989;ISIZE=12372;MATEID=BND_atzsjndp GT:CTX:DEL:INS:INV:NOV_INS:TDUP:lSC:nSC:uRP:distl_levD 0/1:0:0:0:61:7:0:0:0:0:0
|
|
58 chr1 915564 BND_oonwkzas C C[chr1:915564[ 1000.00 PASS SVTYPE=BND;EVENT=INV;END=903617;ISIZE=11947;MATEID=BND_jtvtjzwt GT:CTX:DEL:INS:INV:NOV_INS:TDUP:lSC:nSC:uRP:distl_levD 0/1:0:0:0:64:4:0:0:0:0:0
|
|
59 chr1 915989 BND_atzsjndp G ]chr1:915989]C 3.83 PASS SVTYPE=BND;EVENT=INV;END=903617;ISIZE=12372;MATEID=BND_ysqifxup GT:CTX:DEL:INS:INV:NOV_INS:TDUP:lSC:nSC:uRP:distl_levD 0/1:0:0:0:61:7:0:23:7:7:0.26
|
|
60 chr1 915991 BND_akjluuhn C T[chr1:915991[ 3.12 PASS SVTYPE=BND;EVENT=INV;END=903615;ISIZE=12376;MATEID=BND_wpzdzdca GT:CTX:DEL:INS:INV:NOV_INS:TDUP:lSC:nSC:uRP:distl_levD 0/1:0:0:0:60:11:0:0:0:0:0
|
|
61 chr1 915991 BND_puxffwzz C ]chr1:915991]G 1000.00 PASS SVTYPE=BND;EVENT=INV;END=903148;ISIZE=12843;MATEID=BND_xbqyvnjg GT:CTX:DEL:INS:INV:NOV_INS:TDUP:lSC:nSC:uRP:distl_levD 0/1:0:0:0:64:3:0:42:14:3:0.00
|
|
62 chr1 937624 BND_vzxstibz A T[chr1:937624[ 1000.00 PASS SVTYPE=BND;EVENT=INV;END=958311;ISIZE=20687;MATEID=BND_tftqvjfy GT:CTX:DEL:INS:INV:NOV_INS:TDUP:lSC:nSC:uRP:distl_levD 0/1:0:0:0:52:3:0:0:0:0:0
|
|
63 chr1 938004 BND_gexxmgys T ]chr1:938004]A 4.98 PASS SVTYPE=BND;EVENT=INV;END=958308;ISIZE=20304;MATEID=BND_yxhungrm GT:CTX:DEL:INS:INV:NOV_INS:TDUP:lSC:nSC:uRP:distl_levD 0/1:0:0:0:59:7:0:10:5:7:0.20
|
|
64 chr1 938006 BND_trhijbux A G[chr1:938006[ 3.98 PASS SVTYPE=BND;EVENT=INV;END=958306;ISIZE=20300;MATEID=BND_leqnxqht GT:CTX:DEL:INS:INV:NOV_INS:TDUP:lSC:nSC:uRP:distl_levD 0/1:0:0:0:59:4:0:0:0:0:0
|
|
65 chr1 958306 BND_leqnxqht G ]chr1:958306]A 3.98 PASS SVTYPE=BND;EVENT=INV;END=938006;ISIZE=20300;MATEID=BND_trhijbux GT:CTX:DEL:INS:INV:NOV_INS:TDUP:lSC:nSC:uRP:distl_levD 0/1:0:0:0:59:4:0:40:10:4:0.25
|
|
66 chr1 958308 BND_yxhungrm A T[chr1:958308[ 4.98 PASS SVTYPE=BND;EVENT=INV;END=938004;ISIZE=20304;MATEID=BND_gexxmgys GT:CTX:DEL:INS:INV:NOV_INS:TDUP:lSC:nSC:uRP:distl_levD 0/1:0:0:0:59:7:0:0:0:0:0
|
|
67 chr1 958311 BND_tftqvjfy T ]chr1:958311]A 1000.00 PASS SVTYPE=BND;EVENT=INV;END=937624;ISIZE=20687;MATEID=BND_vzxstibz GT:CTX:DEL:INS:INV:NOV_INS:TDUP:lSC:nSC:uRP:distl_levD 0/1:0:0:0:52:3:0:35:11:3:0.00
|