Mercurial > repos > devteam > vcfgeno2haplo
changeset 0:2bdb74ab3270 draft
Uploaded
author | devteam |
---|---|
date | Thu, 19 Mar 2015 14:32:31 -0400 |
parents | |
children | 401f29f278f0 |
files | macros.xml test-data/vcfgeno2haplo-test1.vcf test-data/vcflib-phix.vcf test-data/vcflib-test-genome-phix.fa tool-data/fasta_indexes.loc.sample tool-data/tool_data_table_conf.xml.sample tool_dependencies.xml vcfgeno2haplo.xml |
diffstat | 8 files changed, 379 insertions(+), 0 deletions(-) [+] |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/macros.xml Thu Mar 19 14:32:31 2015 -0400 @@ -0,0 +1,28 @@ +<macros> + <xml name="requirements"> + <requirements> + <requirement type="package" version="8a5602bf07">vcflib</requirement> + <yield/> + </requirements> + </xml> + <xml name="stdio"> + <stdio> + <exit_code range="1:" level="fatal" /> + </stdio> + </xml> + <xml name="citations"> + <citations> + <citation type="bibtex"> +@misc{Garrison2015, + author = {Garrison, Erik}, + year = {2015}, + title = {vcflib}, + publisher = {GitHub}, + journal = {GitHub repository}, + url = {https://github.com/ekg/vcflib}, +} + </citation> + </citations> + </xml> + <token name="@IS_PART_OF_VCFLIB@">is a part of VCFlib toolkit developed by Erik Garrison (https://github.com/ekg/vcflib).</token> +</macros>
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/vcfgeno2haplo-test1.vcf Thu Mar 19 14:32:31 2015 -0400 @@ -0,0 +1,12 @@ +##fileformat=VCFv4.1 +##fileDate=20140610 +##source=Dan +##reference=file:///galaxy/data/phiX/sam_index/phiX.fa +##filter="AF > .001" +##INFO=<ID=AC,Number=A,Type=Integer,Description="Allele count in genotypes, for each ALT allele, in the same order as listed"> +##INFO=<ID=AF,Number=A,Type=Float,Description="Allele Frequency, for each ALT allele, in the same order as listed"> +##FORMAT=<ID=GT,Number=1,Type=String,Description="Genotype"> +##FORMAT=<ID=AC,Number=.,Type=Integer,Description="Allele count in genotypes, for each ALT allele, in the same order as listed"> +##FORMAT=<ID=AF,Number=.,Type=Float,Description="Allele Frequency, for each ALT allele, in the same order as listed"> +##FORMAT=<ID=NC,Number=.,Type=String,Description="Nucleotide and indel counts"> +#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT groupA groupB
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/vcflib-phix.vcf Thu Mar 19 14:32:31 2015 -0400 @@ -0,0 +1,219 @@ +##fileformat=VCFv4.1 +##fileDate=20140610 +##source=Dan +##reference=file:///galaxy/data/phiX/sam_index/phiX.fa +##filter="AF > .001" +##INFO=<ID=AC,Number=A,Type=Integer,Description="Allele count in genotypes, for each ALT allele, in the same order as listed"> +##INFO=<ID=AF,Number=A,Type=Float,Description="Allele Frequency, for each ALT allele, in the same order as listed"> +##FORMAT=<ID=GT,Number=1,Type=String,Description="Genotype"> +##FORMAT=<ID=AC,Number=.,Type=Integer,Description="Allele count in genotypes, for each ALT allele, in the same order as listed"> +##FORMAT=<ID=AF,Number=.,Type=Float,Description="Allele Frequency, for each ALT allele, in the same order as listed"> +##FORMAT=<ID=NC,Number=.,Type=String,Description="Nucleotide and indel counts"> +#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT groupA groupB +phiX174 1015 . A C 0 . AC=2;AF=0.0555555555556 GT:AC:AF:NC 0:0,2:0,0.2:A=8,C=2,G=1, 0:0,0:0,0:A=26,G=1, +phiX174 1016 . A T 0 . AC=2;AF=0.0555555555556 GT:AC:AF:NC 0:0:0:A=10, 0:2:0.0769230769231:A=24,T=2, +phiX174 1038 . A C 0 . AC=2;AF=0.046511627907 GT:AC:AF:NC 0:2:0.153846153846:A=11,C=2, 0:0:0:A=30, +phiX174 1045 . T G 0 . AC=2;AF=0.0425531914894 GT:AC:AF:NC 0:0:0:T=15,G=1, 0:2:0.0625:T=30,G=2, +phiX174 1048 . C A 0 . AC=2;AF=0.0434782608696 GT:AC:AF:NC 0:2:0.133333333333:A=2,C=13, 0:0:0:C=31,G=1, +phiX174 1050 . A C 0 . AC=3;AF=0.0714285714286 GT:AC:AF:NC 0:0:0:A=12, 0:3:0.1:A=27,C=3, +phiX174 1051 . T G 0 . AC=2;AF=0.0487804878049 GT:AC:AF:NC 0:2:0.166666666667:C=1,T=10,G=2, 0:0:0:T=29,G=1, +phiX174 1054 . T G 0 . AC=2;AF=0.0454545454545 GT:AC:AF:NC 0:2:0.153846153846:T=11,G=2, 0:0:0:T=31,G=1, +phiX174 1076 . G C 0 . AC=2;AF=0.0363636363636 GT:AC:AF:NC 0:2:0.105263157895:A=1,C=2,G=17, 0:0:0:G=36, +phiX174 1080 . G C 0 . AC=2;AF=0.0344827586207 GT:AC:AF:NC 0:2:0.105263157895:A=1,C=2,T=1,G=17, 0:0:0:C=1,G=39, +phiX174 1092 . T G 0 . AC=2;AF=0.0263157894737 GT:AC:AF:NC 0:2:0.0833333333333:C=1,T=22,G=2, 0:0:0:T=52,G=1, +phiX174 1093 . T G 0 . AC=2;AF=0.025974025974 GT:AC:AF:NC 0:2:0.08:T=23,G=2, 0:0:0:T=52, +phiX174 1102 . T C 0 . AC=6;AF=0.0779220779221 GT:AC:AF:NC 0:2:0.0869565217391:C=2,T=21, 0:4:0.0740740740741:C=4,T=50, +phiX174 1107 . T G 0 . AC=2;AF=0.0277777777778 GT:AC:AF:NC 0:0:0:T=22,G=1, 0:2:0.04:T=48,G=2, +phiX174 1110 . T C 0 . AC=2;AF=0.0277777777778 GT:AC:AF:NC 0:2:0.125:C=2,T=14, 0:0:0:C=1,T=56,G=1, +phiX174 1121 . T G 0 . AC=3;AF=0.041095890411 GT:AC:AF:NC 0:0:0:T=19, 0:3:0.0555555555556:T=51,G=3, +phiX174 1124 . T C 0 . AC=2;AF=0.0307692307692 GT:AC:AF:NC 0:0,2:0,0.117647058824:C=2,T=15,G=1, 0:0,0:0,0:T=48,G=1, +phiX174 1125 . T G 0 . AC=3;AF=0.0454545454545 GT:AC:AF:NC 0:3:0.176470588235:T=14,G=3, 0:0:0:T=49,G=1, +phiX174 1131 . T G 0 . AC=2;AF=0.025641025641 GT:AC:AF:NC 0:2:0.0952380952381:T=19,G=2, 0:0:0:T=57, +phiX174 1141 . T G 0 . AC=2;AF=0.0253164556962 GT:AC:AF:NC 0:0:0:T=20,G=1, 0:2:0.0338983050847:T=57,G=2, +phiX174 1145 . G C 0 . AC=2;AF=0.027397260274 GT:AC:AF:NC 0:2:0.0952380952381:C=2,G=19, 0:0:0:G=52, +phiX174 1147 . T C 0 . AC=2;AF=0.0281690140845 GT:AC:AF:NC 0:2:0.0952380952381:C=2,T=19, 0:0:0:T=50, +phiX174 1149 . T G,C 0 . AC=2,2;AF=0.027027027027,0.027027027027 GT:AC:AF:NC 0:2,2:0.1,0.1:C=2,T=16,G=2, 0:0,0:0,0:T=54, +phiX174 1156 . T G 0 . AC=2;AF=0.031746031746 GT:AC:AF:NC 0:0:0:C=1,T=15, 0:2:0.0416666666667:T=46,G=2, +phiX174 1159 . T C,G 0 . AC=2,2;AF=0.0298507462687,0.0298507462687 GT:AC:AF:NC 0:2,0:0.111111111111,0:C=2,T=16, 0:0,2:0,0.0408163265306:C=1,T=47,G=2, +phiX174 1164 . T G 0 . AC=4;AF=0.0666666666667 GT:AC:AF:NC 0:0:0:C=1,T=16,G=1, 0:4:0.0909090909091:T=40,G=4, +phiX174 1168 . T G 0 . AC=4;AF=0.0666666666667 GT:AC:AF:NC 0:2:0.0952380952381:C=1,T=19,G=2, 0:2:0.0512820512821:T=37,G=2, +phiX174 1171 . T G 0 . AC=2;AF=0.0327868852459 GT:AC:AF:NC 0:2:0.0869565217391:T=21,G=2, 0:0:0:T=38,G=1, +phiX174 1174 . C A 0 . AC=2;AF=0.0327868852459 GT:AC:AF:NC 0:2:0.0869565217391:A=2,C=21, 0:0:0:A=1,C=38, +phiX174 1178 . G C 0 . AC=2;AF=0.0307692307692 GT:AC:AF:NC 0:2:0.0740740740741:C=2,T=1,G=25, 0:0:0:G=38, +phiX174 1187 . T G 0 . AC=3;AF=0.04 GT:AC:AF:NC 0:3:0.0882352941176:C=1,T=31,G=3, 0:0:0:T=41, +phiX174 1189 . T C 0 . AC=3;AF=0.0379746835443 GT:AC:AF:NC 0:3:0.0857142857143:C=3,T=32,G=1, 0:0:0:T=44, +phiX174 1200 . T G 0 . AC=3;AF=0.0365853658537 GT:AC:AF:NC 0:3:0.0882352941176:T=31,G=3, 0:0:0:T=48,G=1, +phiX174 1202 . T G 0 . AC=3;AF=0.0365853658537 GT:AC:AF:NC 0:0:0:T=32, 0:3:0.06:T=47,G=3, +phiX174 1205 . A C 0 . AC=2;AF=0.0246913580247 GT:AC:AF:NC 0:2:0.0588235294118:A=32,C=2, 0:0:0:A=47,C=1,G=1, +phiX174 1207 . T G 0 . AC=2;AF=0.0246913580247 GT:AC:AF:NC 0:0:0:A=1,T=34, 0:2:0.0425531914894:T=45,G=2, +phiX174 1214 . G A 0 . AC=2;AF=0.0246913580247 GT:AC:AF:NC 0:0:0:T=1,G=34, 0:2:0.0425531914894:A=2,G=45, +phiX174 1216 . C A 0 . AC=2;AF=0.0246913580247 GT:AC:AF:NC 0:0:0:C=34, 0:2:0.0425531914894:A=2,C=45, +phiX174 1217 . C A 0 . AC=2;AF=0.025 GT:AC:AF:NC 0:2:0.0588235294118:A=2,C=32, 0:0:0:C=46, +phiX174 1219 . T G 0 . AC=2;AF=0.025641025641 GT:AC:AF:NC 0:0:0:C=1,T=33,G=1, 0:2:0.0444444444444:T=43,G=2, +phiX174 1223 . C T 0 . AC=2;AF=0.025 GT:AC:AF:NC 0:2:0.0555555555556:C=34,T=2, 0:0:0:C=44, +phiX174 1227 . A C 0 . AC=2;AF=0.027027027027 GT:AC:AF:NC 0:2:0.0588235294118:A=32,C=2, 0:0:0:A=40, +phiX174 1232 . T G 0 . AC=2;AF=0.0298507462687 GT:AC:AF:NC 0:0:0:T=34, 0:2:0.0606060606061:C=1,T=31,G=2, +phiX174 1236 . G T 0 . AC=2;AF=0.0294117647059 GT:AC:AF:NC 0:2:0.0606060606061:A=1,T=2,G=31, 0:0:0:G=35, +phiX174 1237 . T G 0 . AC=2;AF=0.0294117647059 GT:AC:AF:NC 0:2:0.0571428571429:A=1,T=33,G=2, 0:0:0:T=33,G=1, +phiX174 1240 . A C 0 . AC=2;AF=0.0294117647059 GT:AC:AF:NC 0:0,0:0,0:A=32,T=1, 0:0,2:0,0.0555555555556:A=34,C=2,T=1, +phiX174 1245 . G T 0 . AC=4;AF=0.0625 GT:AC:AF:NC 0:2:0.0666666666667:C=1,T=2,G=28, 0:2:0.0588235294118:T=2,G=32, +phiX174 1247 . A T 0 . AC=2;AF=0.031746031746 GT:AC:AF:NC 0:2:0.0714285714286:A=26,C=1,T=2,G=1, 0:0:0:A=35, +phiX174 1250 . A C 0 . AC=2;AF=0.0307692307692 GT:AC:AF:NC 0:2:0.0769230769231:A=24,C=2, 0:0:0:A=39, +phiX174 1251 . A G 0 . AC=3;AF=0.044776119403 GT:AC:AF:NC 0:0:0:A=29, 0:3:0.0789473684211:A=35,C=1,G=3, +phiX174 1253 . T G 0 . AC=4;AF=0.0526315789474 GT:AC:AF:NC 0:2:0.0645161290323:T=29,G=2, 0:2:0.0444444444444:T=43,G=2, +phiX174 1257 . T G 0 . AC=2;AF=0.0281690140845 GT:AC:AF:NC 0:0:0:T=26, 0:2:0.0444444444444:T=43,G=2, +phiX174 1259 . A C 0 . AC=2;AF=0.027397260274 GT:AC:AF:NC 0:0:0:A=28,T=1, 0:2:0.0444444444444:A=43,C=2,G=1, +phiX174 1274 . G C 0 . AC=2;AF=0.0240963855422 GT:AC:AF:NC 0:2:0.0512820512821:C=2,G=37, 0:0:0:G=44, +phiX174 1277 . A C 0 . AC=2;AF=0.025 GT:AC:AF:NC 0:2:0.05:A=38,C=2, 0:0:0:A=40,G=1, +phiX174 1279 . T G 0 . AC=2;AF=0.0253164556962 GT:AC:AF:NC 0:0:0:T=39,G=1, 0:2:0.05:C=1,T=38,G=2, +phiX174 1284 . T G,C 0 . AC=3,2;AF=0.038961038961,0.025974025974 GT:AC:AF:NC 0:3,2:0.075,0.05:C=2,T=35,G=3, 0:0,0:0,0:A=1,T=37, +phiX174 1289 . A C 0 . AC=3;AF=0.0454545454545 GT:AC:AF:NC 0:3:0.0789473684211:A=35,C=3, 0:0:0:A=28,G=1, +phiX174 1296 . A C 0 . AC=2;AF=0.03125 GT:AC:AF:NC 0:0:0:A=38, 0:2:0.0769230769231:A=24,C=2, +phiX174 1298 . A C 0 . AC=3;AF=0.0491803278689 GT:AC:AF:NC 0:3,0:0.0810810810811,0:A=34,C=3,G=1, 0:0,0:0,0:A=24,C=1,T=1,G=1, +phiX174 1306 . T G 0 . AC=5;AF=0.0943396226415 GT:AC:AF:NC 0:5:0.178571428571:T=23,G=5, 0:0:0:C=1,T=25, +phiX174 1307 . T G 0 . AC=3;AF=0.0545454545455 GT:AC:AF:NC 0:3:0.107142857143:T=25,G=3, 0:0:0:C=1,T=27, +phiX174 1314 . A G 0 . AC=5;AF=0.0877192982456 GT:AC:AF:NC 0:2,0:0.0740740740741,0:A=25,C=1,G=2, 0:3,0:0.1,0:A=27,C=1,G=3, +phiX174 1317 . A C 0 . AC=2;AF=0.0322580645161 GT:AC:AF:NC 0:2:0.0740740740741:A=25,C=2, 0:0:0:A=35, +phiX174 1320 . C G 0 . AC=2;AF=0.0322580645161 GT:AC:AF:NC 0:2:0.0769230769231:C=24,G=2, 0:0:0:C=36, +phiX174 1322 . G C 0 . AC=3;AF=0.0441176470588 GT:AC:AF:NC 0:3:0.107142857143:C=3,G=25, 0:0:0:G=40, +phiX174 1323 . C G 0 . AC=2;AF=0.0298507462687 GT:AC:AF:NC 0:0:0:C=28, 0:2:0.0512820512821:C=37,T=1,G=2, +phiX174 1329 . T C 0 . AC=2;AF=0.0277777777778 GT:AC:AF:NC 0:2:0.0666666666667:C=2,T=28, 0:0:0:T=42,G=1, +phiX174 1335 . C A 0 . AC=2;AF=0.025974025974 GT:AC:AF:NC 0:2:0.0625:A=2,C=30, 0:0:0:C=45, +phiX174 1338 . T A 0 . AC=2;AF=0.027027027027 GT:AC:AF:NC 0:0:0:T=32, 0:2:0.047619047619:A=2,T=40, +phiX174 1345 . T C 0 . AC=2;AF=0.0253164556962 GT:AC:AF:NC 0:0,0:0,0:T=36,G=1, 0:0,2:0,0.046511627907:C=2,T=41,G=1, +phiX174 1347 . A C 0 . AC=2;AF=0.0243902439024 GT:AC:AF:NC 0:2:0.0512820512821:A=37,C=2,G=1, 0:0:0:A=43, +phiX174 1348 . T G 0 . AC=2;AF=0.0240963855422 GT:AC:AF:NC 0:2:0.05:T=38,G=2, 0:0:0:T=43, +phiX174 1352 . A T 0 . AC=2;AF=0.027397260274 GT:AC:AF:NC 0:0:0:A=34,G=1, 0:2:0.0512820512821:A=37,T=2, +phiX174 1354 . T G,A 0 . AC=2,2;AF=0.027027027027,0.027027027027 GT:AC:AF:NC 0:2,0:0.0526315789474,0:T=36,G=2, 0:0,2:0,0.0555555555556:A=2,T=34,G=1, +phiX174 1359 . T G 0 . AC=3;AF=0.0405405405405 GT:AC:AF:NC 0:3:0.0714285714286:T=39,G=3, 0:0:0:T=32,G=1, +phiX174 1366 . G C 0 . AC=3;AF=0.0428571428571 GT:AC:AF:NC 0:3:0.0714285714286:C=3,T=1,G=39, 0:0:0:G=28, +phiX174 1370 . T G 0 . AC=2;AF=0.0289855072464 GT:AC:AF:NC 0:0:0:C=1,T=45,G=1, 0:2:0.0833333333333:A=1,T=22,G=2, +phiX174 1373 . T G 0 . AC=2;AF=0.027027027027 GT:AC:AF:NC 0:2:0.0408163265306:T=47,G=2, 0:0:0:T=25, +phiX174 1380 . G C 0 . AC=2;AF=0.0238095238095 GT:AC:AF:NC 0:0,2:0,0.0384615384615:C=2,T=1,G=50, 0:0,0:0,0:T=1,G=32, +phiX174 1381 . T G 0 . AC=3;AF=0.0361445783133 GT:AC:AF:NC 0:3:0.0576923076923:A=1,T=49,G=3, 0:0:0:T=31, +phiX174 1384 . T G 0 . AC=4;AF=0.0512820512821 GT:AC:AF:NC 0:4:0.0816326530612:T=45,G=4, 0:0:0:T=29,G=1, +phiX174 1385 . T G 0 . AC=4;AF=0.0512820512821 GT:AC:AF:NC 0:4:0.0816326530612:T=45,G=4, 0:0:0:C=1,T=29, +phiX174 1391 . A C 0 . AC=2;AF=0.0253164556962 GT:AC:AF:NC 0:2:0.0408163265306:A=47,C=2, 0:0:0:A=30,C=1, +phiX174 1393 . C A 0 . AC=2;AF=0.0253164556962 GT:AC:AF:NC 0:2:0.0425531914894:A=2,C=45, 0:0:0:C=32, +phiX174 1396 . T G 0 . AC=3;AF=0.04 GT:AC:AF:NC 0:3:0.0697674418605:T=40,G=3, 0:0:0:T=32, +phiX174 1397 . A C 0 . AC=2;AF=0.027027027027 GT:AC:AF:NC 0:0:0:A=42,T=1, 0:2:0.0625:A=30,C=2, +phiX174 1400 . A C 0 . AC=2;AF=0.027397260274 GT:AC:AF:NC 0:2:0.047619047619:A=40,C=2,G=1, 0:0:0:A=31, +phiX174 1401 . A C 0 . AC=2;AF=0.027397260274 GT:AC:AF:NC 0:2:0.047619047619:A=40,C=2, 0:0:0:A=31, +phiX174 1434 . C T 0 . AC=2;AF=0.0571428571429 GT:AC:AF:NC 0:2:0.142857142857:C=12,T=2, 0:0:0:C=21, +phiX174 1438 . G C 0 . AC=2;AF=0.0526315789474 GT:AC:AF:NC 0:0:0:G=15, 0:2:0.0869565217391:C=2,G=21, +phiX174 1450 . T G 0 . AC=3;AF=0.0789473684211 GT:AC:AF:NC 0:3:0.176470588235:A=1,C=1,T=14,G=3, 0:0:0:T=21,G=1, +phiX174 1459 . T G 0 . AC=3;AF=0.075 GT:AC:AF:NC 0:3:0.142857142857:C=1,T=18,G=3, 0:0:0:T=19, +phiX174 1470 . C A 0 . AC=2;AF=0.0384615384615 GT:AC:AF:NC 0:0:0:C=29, 0:2:0.0869565217391:A=2,C=21, +phiX174 1484 . C G 0 . AC=2;AF=0.0384615384615 GT:AC:AF:NC 0:2:0.0740740740741:A=1,C=25,G=2, 0:0:0:C=25, +phiX174 1487 . T G 0 . AC=2;AF=0.0363636363636 GT:AC:AF:NC 0:2:0.0689655172414:T=27,G=2, 0:0:0:T=26, +phiX174 1491 . G T 0 . AC=2;AF=0.037037037037 GT:AC:AF:NC 0:0:0:T=1,G=26, 0:2:0.0714285714286:T=2,G=26, +phiX174 1503 . A C 0 . AC=2;AF=0.0454545454545 GT:AC:AF:NC 0:0:0:A=20,C=1, 0:2:0.0833333333333:A=22,C=2, +phiX174 1536 . A C 0 . AC=3;AF=0.0545454545455 GT:AC:AF:NC 0:3:0.111111111111:A=24,C=3, 0:0:0:A=28,C=1, +phiX174 1541 . T G 0 . AC=2;AF=0.0363636363636 GT:AC:AF:NC 0:2:0.0740740740741:T=25,G=2, 0:0:0:T=28, +phiX174 1553 . A G 0 . AC=2;AF=0.0408163265306 GT:AC:AF:NC 0:2:0.0952380952381:A=19,C=1,G=2, 0:0:0:A=28, +phiX174 1566 . C G 0 . AC=2;AF=0.0307692307692 GT:AC:AF:NC 0:0:0:C=30,T=1, 0:2:0.0571428571429:C=33,G=2, +phiX174 1572 . T G 0 . AC=2;AF=0.03125 GT:AC:AF:NC 0:2:0.0625:T=30,G=2, 0:0:0:T=32, +phiX174 1574 . G C 0 . AC=3;AF=0.0441176470588 GT:AC:AF:NC 0:3:0.0882352941176:C=3,G=31, 0:0:0:G=34, +phiX174 1575 . A C 0 . AC=2;AF=0.0294117647059 GT:AC:AF:NC 0:2:0.0571428571429:A=33,C=2, 0:0:0:A=33,T=1, +phiX174 1577 . A C 0 . AC=3;AF=0.0441176470588 GT:AC:AF:NC 0:0:0:A=35, 0:3:0.0909090909091:A=30,C=3, +phiX174 1581 . T G 0 . AC=2;AF=0.0285714285714 GT:AC:AF:NC 0:2:0.0555555555556:C=1,T=34,G=2, 0:0:0:T=34,G=1, +phiX174 1591 . A C 0 . AC=3;AF=0.037037037037 GT:AC:AF:NC 0:3:0.0789473684211:A=35,C=3, 0:0:0:A=43,C=1, +phiX174 1597 . T G 0 . AC=2;AF=0.025641025641 GT:AC:AF:NC 0:2:0.0606060606061:T=31,G=2, 0:0:0:C=1,T=45, +phiX174 1603 . T G 0 . AC=2;AF=0.0243902439024 GT:AC:AF:NC 0:2:0.0625:T=30,G=2, 0:0:0:T=50, +phiX174 1606 . T G 0 . AC=2;AF=0.0243902439024 GT:AC:AF:NC 0:2:0.0666666666667:A=1,C=1,T=28,G=2, 0:0:0:T=52, +phiX174 1611 . A C 0 . AC=2;AF=0.0243902439024 GT:AC:AF:NC 0:2:0.0769230769231:A=24,C=2, 0:0:0:A=56,C=1, +phiX174 1613 . A G 0 . AC=2;AF=0.0240963855422 GT:AC:AF:NC 0:0:0:A=25,C=1,G=1, 0:2:0.0344827586207:A=56,G=2, +phiX174 1618 . C A 0 . AC=2;AF=0.0229885057471 GT:AC:AF:NC 0:2:0.0666666666667:A=2,C=28,G=1, 0:0:0:A=1,C=57, +phiX174 1619 . C A 0 . AC=2;AF=0.0229885057471 GT:AC:AF:NC 0:2:0.0666666666667:A=2,C=28, 0:0:0:C=57, +phiX174 1621 . A C 0 . AC=3;AF=0.0352941176471 GT:AC:AF:NC 0:0:0:A=31,T=1,G=1, 0:3:0.0555555555556:A=51,C=3, +phiX174 1622 . G T 0 . AC=2;AF=0.0224719101124 GT:AC:AF:NC 0:0:0:A=1,G=35, 0:2:0.037037037037:T=2,G=52, +phiX174 1624 . A G 0 . AC=2;AF=0.0235294117647 GT:AC:AF:NC 0:2:0.0555555555556:A=34,T=1,G=2, 0:0:0:A=49, +phiX174 1627 . T G,C 0 . AC=2,2;AF=0.0246913580247,0.0246913580247 GT:AC:AF:NC 0:0,0:0,0:T=37,G=1, 0:2,2:0.0454545454545,0.0454545454545:C=2,T=40,G=2, +phiX174 1630 . T G 0 . AC=2;AF=0.025641025641 GT:AC:AF:NC 0:2:0.0526315789474:C=1,T=36,G=2, 0:0:0:T=40, +phiX174 1632 . A C 0 . AC=2;AF=0.0246913580247 GT:AC:AF:NC 0:2:0.05:A=38,C=2, 0:0:0:A=41, +phiX174 1634 . T G 0 . AC=4;AF=0.0493827160494 GT:AC:AF:NC 0:2:0.047619047619:C=1,T=40,G=2, 0:2:0.0512820512821:T=37,G=2, +phiX174 1635 . T G 0 . AC=3;AF=0.0365853658537 GT:AC:AF:NC 0:3:0.0714285714286:T=39,G=3, 0:0:0:T=40, +phiX174 1639 . G A 0 . AC=3;AF=0.0405405405405 GT:AC:AF:NC 0:3:0.075:A=3,T=1,G=37, 0:0:0:G=34, +phiX174 1645 . T G 0 . AC=2;AF=0.027027027027 GT:AC:AF:NC 0:2:0.0425531914894:T=45,G=2, 0:0:0:T=27,G=1, +phiX174 1647 . A T,G,C 0 . AC=2,2,2;AF=0.0253164556962,0.0253164556962,0.0253164556962 GT:AC:AF:NC 0:0,2,2:0,0.0384615384615,0.0384615384615:A=48,C=2,G=2, 0:2,0,0:0.0740740740741,0,0:A=25,T=2, +phiX174 1652 . G A 0 . AC=2;AF=0.0240963855422 GT:AC:AF:NC 0:2:0.037037037037:A=2,G=52, 0:0:0:G=29, +phiX174 1656 . T G 0 . AC=3;AF=0.0348837209302 GT:AC:AF:NC 0:3:0.0555555555556:T=51,G=3, 0:0:0:T=32,G=1, +phiX174 1658 . A C 0 . AC=2;AF=0.0240963855422 GT:AC:AF:NC 0:2:0.04:A=48,C=2, 0:0:0:A=33, +phiX174 1659 . T G 0 . AC=3;AF=0.0365853658537 GT:AC:AF:NC 0:3:0.0625:C=1,T=45,G=3, 0:0:0:T=34, +phiX174 1663 . T G 0 . AC=4;AF=0.0454545454545 GT:AC:AF:NC 0:0:0:A=1,T=53, 0:4:0.114285714286:T=31,G=4, +phiX174 1665 . C A 0 . AC=2;AF=0.0212765957447 GT:AC:AF:NC 0:2:0.0384615384615:A=2,C=50,G=1, 0:0:0:C=42, +phiX174 1670 . G C 0 . AC=2;AF=0.0224719101124 GT:AC:AF:NC 0:0,2:0,0.0408163265306:A=1,C=2,T=1,G=47, 0:0,0:0,0:T=1,G=40, +phiX174 1671 . G T 0 . AC=2;AF=0.0212765957447 GT:AC:AF:NC 0:2:0.0384615384615:C=1,T=2,G=50, 0:0:0:G=42, +phiX174 1672 . A G 0 . AC=3;AF=0.0315789473684 GT:AC:AF:NC 0:3:0.0545454545455:A=52,G=3, 0:0:0:A=40,T=1, +phiX174 1673 . G T 0 . AC=2;AF=0.0212765957447 GT:AC:AF:NC 0:2:0.0377358490566:C=1,T=2,G=51, 0:0:0:G=41, +phiX174 1674 . G T 0 . AC=4;AF=0.04 GT:AC:AF:NC 0:4:0.0714285714286:T=4,G=52, 0:0:0:G=44, +phiX174 1677 . A C 0 . AC=3;AF=0.031914893617 GT:AC:AF:NC 0:3:0.0566037735849:A=50,C=3,N=1, 0:0:0:A=41, +phiX174 1678 . A C 0 . AC=2;AF=0.0212765957447 GT:AC:AF:NC 0:0:0:A=54,C=1,G=1, 0:2:0.05:A=38,C=2, +phiX174 1682 . T C 0 . AC=2;AF=0.0229885057471 GT:AC:AF:NC 0:2:0.0408163265306:C=2,T=47,G=1, 0:0:0:C=1,T=38, +phiX174 1685 . T G 0 . AC=2;AF=0.0238095238095 GT:AC:AF:NC 0:2:0.046511627907:A=1,T=41,G=2, 0:0:0:T=41,G=1, +phiX174 1686 . A C 0 . AC=2;AF=0.025641025641 GT:AC:AF:NC 0:0:0:A=39,C=1,T=1,G=1, 0:2:0.0512820512821:A=37,C=2, +phiX174 1687 . T C 0 . AC=2;AF=0.0243902439024 GT:AC:AF:NC 0:2:0.047619047619:C=2,T=40, 0:0:0:A=1,T=40, +phiX174 1688 . G C 0 . AC=2;AF=0.0246913580247 GT:AC:AF:NC 0:2:0.05:C=2,T=1,G=38, 0:0:0:G=41, +phiX174 1689 . A C 0 . AC=2;AF=0.0243902439024 GT:AC:AF:NC 0:2:0.0487804878049:A=39,C=2,T=1,G=1, 0:0:0:A=41,C=1, +phiX174 1691 . G C 0 . AC=2;AF=0.025974025974 GT:AC:AF:NC 0:2:0.05:C=2,T=1,G=38, 0:0:0:G=37, +phiX174 1693 . T G 0 . AC=2;AF=0.0238095238095 GT:AC:AF:NC 0:0:0:T=42,G=1, 0:2:0.047619047619:T=40,G=2, +phiX174 1695 . A C 0 . AC=2;AF=0.0232558139535 GT:AC:AF:NC 0:2:0.046511627907:A=41,C=2, 0:0:0:A=43,C=1, +phiX174 1702 . T C 0 . AC=2;AF=0.0246913580247 GT:AC:AF:NC 0:2:0.0512820512821:A=1,C=2,T=37, 0:0:0:C=1,T=42, +phiX174 1707 . T C 0 . AC=2;AF=0.0232558139535 GT:AC:AF:NC 0:2:0.05:C=2,T=38, 0:0:0:T=46,G=1, +phiX174 1708 . A C 0 . AC=2;AF=0.0238095238095 GT:AC:AF:NC 0:2:0.0526315789474:A=36,C=2, 0:0:0:A=46,C=1,G=1, +phiX174 1710 . T C 0 . AC=2;AF=0.025 GT:AC:AF:NC 0:0:0:C=1,T=34, 0:2:0.0434782608696:C=2,T=44, +phiX174 1723 . T C 0 . AC=2;AF=0.0243902439024 GT:AC:AF:NC 0:2:0.0588235294118:C=2,T=32, 0:0:0:T=48, +phiX174 1724 . A C 0 . AC=4;AF=0.0481927710843 GT:AC:AF:NC 0:2:0.0571428571429:A=33,C=2, 0:2:0.0416666666667:A=46,C=2, +phiX174 1725 . A C 0 . AC=2;AF=0.0243902439024 GT:AC:AF:NC 0:2:0.0606060606061:A=31,C=2, 0:0:0:A=49, +phiX174 1726 . T G 0 . AC=3;AF=0.037037037037 GT:AC:AF:NC 0:3:0.0882352941176:T=31,G=3, 0:0:0:C=1,T=47, +phiX174 1730 . T G 0 . AC=2;AF=0.027397260274 GT:AC:AF:NC 0:0:0:C=1,T=28, 0:2:0.0444444444444:T=43,G=2, +phiX174 1739 . G T 0 . AC=3;AF=0.0434782608696 GT:AC:AF:NC 0:3:0.09375:T=3,G=29, 0:0:0:G=37, +phiX174 1746 . A C 0 . AC=2;AF=0.0338983050847 GT:AC:AF:NC 0:2:0.0689655172414:A=27,C=2, 0:0:0:A=30,C=1,G=1, +phiX174 1749 . T G,C 0 . AC=3,2;AF=0.044776119403,0.0298507462687 GT:AC:AF:NC 0:3,2:0.0857142857143,0.0571428571429:C=2,T=30,G=3, 0:0,0:0,0:T=32,G=1, +phiX174 1755 . G C 0 . AC=2;AF=0.027397260274 GT:AC:AF:NC 0:2:0.05:C=2,G=38, 0:0:0:C=1,G=33, +phiX174 1756 . A G,C 0 . AC=2,2;AF=0.0281690140845,0.0281690140845 GT:AC:AF:NC 0:0,2:0,0.0526315789474:A=36,C=2,T=1,G=1, 0:2,0:0.0606060606061,0:A=31,C=1,G=2, +phiX174 1761 . A G 0 . AC=2;AF=0.0285714285714 GT:AC:AF:NC 0:2:0.0512820512821:A=37,C=1,G=2, 0:0:0:A=31, +phiX174 1764 . A C 0 . AC=2;AF=0.0294117647059 GT:AC:AF:NC 0:2:0.05:A=38,C=2, 0:0:0:A=28, +phiX174 1766 . A C 0 . AC=2;AF=0.0294117647059 GT:AC:AF:NC 0:2:0.0487804878049:A=39,C=2, 0:0:0:A=27, +phiX174 1769 . T G 0 . AC=2;AF=0.03125 GT:AC:AF:NC 0:2:0.0526315789474:A=1,C=1,T=36,G=2, 0:0:0:T=26,G=1, +phiX174 1771 . G C 0 . AC=2;AF=0.0298507462687 GT:AC:AF:NC 0:2:0.0526315789474:C=2,G=36, 0:0:0:G=29, +phiX174 1772 . T G 0 . AC=4;AF=0.0597014925373 GT:AC:AF:NC 0:4:0.108108108108:C=1,T=33,G=4, 0:0:0:T=30, +phiX174 1774 . A G 0 . AC=3;AF=0.0441176470588 GT:AC:AF:NC 0:3:0.0810810810811:A=34,G=3, 0:0:0:A=31, +phiX174 1775 . G A 0 . AC=2;AF=0.030303030303 GT:AC:AF:NC 0:2:0.0571428571429:A=2,G=33, 0:0:0:G=31, +phiX174 1776 . G T 0 . AC=2;AF=0.0289855072464 GT:AC:AF:NC 0:2:0.0540540540541:T=2,G=35, 0:0:0:G=32, +phiX174 1778 . C G 0 . AC=2;AF=0.0289855072464 GT:AC:AF:NC 0:2:0.0526315789474:C=36,G=2, 0:0:0:C=31, +phiX174 1779 . A C 0 . AC=3;AF=0.0441176470588 GT:AC:AF:NC 0:3:0.075:A=37,C=3, 0:0:0:A=28, +phiX174 1781 . T G 0 . AC=5;AF=0.0714285714286 GT:AC:AF:NC 0:2:0.0487804878049:C=1,T=39,G=2, 0:3:0.103448275862:T=26,G=3, +phiX174 1782 . T G 0 . AC=3;AF=0.0441176470588 GT:AC:AF:NC 0:3:0.075:T=37,G=3, 0:0:0:A=1,T=28, +phiX174 1783 . T C 0 . AC=2;AF=0.0298507462687 GT:AC:AF:NC 0:2:0.0512820512821:C=2,T=37, 0:0:0:T=28,G=1, +phiX174 1788 . G T 0 . AC=2;AF=0.0344827586207 GT:AC:AF:NC 0:0:0:A=1,T=1,G=29, 0:2:0.0689655172414:T=2,G=27, +phiX174 1792 . T G 0 . AC=2;AF=0.0344827586207 GT:AC:AF:NC 0:2:0.0606060606061:T=31,G=2, 0:0:0:T=25,G=1, +phiX174 1802 . A C 0 . AC=2;AF=0.030303030303 GT:AC:AF:NC 0:2:0.0571428571429:A=33,C=2, 0:0:0:A=31,C=1, +phiX174 1804 . C T 0 . AC=2;AF=0.0294117647059 GT:AC:AF:NC 0:2:0.0588235294118:A=1,C=32,T=2, 0:0:0:C=34, +phiX174 1805 . T G 0 . AC=3;AF=0.0441176470588 GT:AC:AF:NC 0:3:0.0857142857143:T=32,G=3, 0:0:0:T=33,G=1, +phiX174 1806 . A C 0 . AC=2;AF=0.0289855072464 GT:AC:AF:NC 0:2:0.0571428571429:A=33,C=2, 0:0:0:A=34,G=1, +phiX174 1824 . G C 0 . AC=2;AF=0.0281690140845 GT:AC:AF:NC 0:2:0.0714285714286:C=2,T=1,G=26, 0:0:0:G=43, +phiX174 1825 . T G 0 . AC=2;AF=0.0281690140845 GT:AC:AF:NC 0:2:0.0689655172414:T=27,G=2, 0:0:0:T=42,G=1, +phiX174 1833 . T G 0 . AC=2;AF=0.0289855072464 GT:AC:AF:NC 0:0:0:A=1,T=24, 0:2:0.0444444444444:T=43,G=2, +phiX174 1843 . T G 0 . AC=2;AF=0.031746031746 GT:AC:AF:NC 0:2:0.08:A=1,T=23,G=2, 0:0:0:T=38, +phiX174 1847 . A C 0 . AC=2;AF=0.0333333333333 GT:AC:AF:NC 0:2:0.0833333333333:A=22,C=2,T=1,G=1, 0:0:0:A=36, +phiX174 1854 . T G 0 . AC=2;AF=0.0298507462687 GT:AC:AF:NC 0:0:0:T=26, 0:2:0.0487804878049:T=39,G=2, +phiX174 1860 . T G 0 . AC=2;AF=0.0377358490566 GT:AC:AF:NC 0:0:0:A=1,T=20, 0:2:0.0606060606061:T=31,G=2, +phiX174 1863 . C G 0 . AC=2;AF=0.0384615384615 GT:AC:AF:NC 0:2:0.105263157895:A=1,C=17,G=2, 0:0:0:C=33, +phiX174 1866 . T G 0 . AC=2;AF=0.04 GT:AC:AF:NC 0:0:0:T=19,G=1, 0:2:0.0645161290323:T=29,G=2, +phiX174 1875 . T A 0 . AC=2;AF=0.0377358490566 GT:AC:AF:NC 0:2:0.0952380952381:A=2,T=19, 0:0:0:T=32,G=1, +phiX174 1876 . T G 0 . AC=3;AF=0.0588235294118 GT:AC:AF:NC 0:3:0.157894736842:C=1,T=16,G=3, 0:0:0:T=32,G=1, +phiX174 1882 . T G 0 . AC=2;AF=0.0377358490566 GT:AC:AF:NC 0:2:0.105263157895:T=17,G=2, 0:0:0:T=34, +phiX174 1889 . A C 0 . AC=2;AF=0.04 GT:AC:AF:NC 0:2:0.0869565217391:A=21,C=2, 0:0:0:A=27, +phiX174 1913 . G C 0 . AC=2;AF=0.04 GT:AC:AF:NC 0:0:0:G=27, 0:2:0.0869565217391:C=2,G=21, +phiX174 1915 . T G 0 . AC=2;AF=0.0425531914894 GT:AC:AF:NC 0:0:0:A=1,T=24, 0:2:0.0869565217391:T=21,G=2, +phiX174 1917 . A T 0 . AC=2;AF=0.0444444444444 GT:AC:AF:NC 0:2,0:0.0833333333333,0:A=22,T=2,G=1, 0:0,0:0,0:A=21,G=1, +phiX174 1930 . T A 0 . AC=2;AF=0.0526315789474 GT:AC:AF:NC 0:2:0.0909090909091:A=2,T=20, 0:0:0:T=16, +phiX174 1934 . A C 0 . AC=2;AF=0.0512820512821 GT:AC:AF:NC 0:2:0.0952380952381:A=19,C=2,T=1, 0:0:0:A=18,C=1, +phiX174 1940 . A C 0 . AC=2;AF=0.046511627907 GT:AC:AF:NC 0:2:0.0952380952381:A=19,C=2, 0:0:0:A=22,C=1, +phiX174 1942 . T G 0 . AC=2;AF=0.047619047619 GT:AC:AF:NC 0:0:0:T=19, 0:2:0.0869565217391:T=21,G=2, +phiX174 1945 . T G 0 . AC=3;AF=0.0769230769231 GT:AC:AF:NC 0:0:0:T=15,G=1, 0:3:0.125:T=21,G=3, +phiX174 1954 . T C 0 . AC=2;AF=0.0625 GT:AC:AF:NC 0:0,0:0,0:T=12,G=1, 0:0,2:0,0.1:C=2,T=18,G=1, +phiX174 1959 . T G 0 . AC=2;AF=0.0555555555556 GT:AC:AF:NC 0:2:0.153846153846:T=11,G=2, 0:0:0:T=23, +phiX174 1994 . A C 0 . AC=2;AF=0.0512820512821 GT:AC:AF:NC 0:2:0.153846153846:A=11,C=2, 0:0:0:A=26, +phiX174 1998 . A G 0 . AC=2;AF=0.05 GT:AC:AF:NC 0:2:0.142857142857:A=12,G=2, 0:0:0:A=26,
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/vcflib-test-genome-phix.fa Thu Mar 19 14:32:31 2015 -0400 @@ -0,0 +1,2 @@ +>phiX174 +GAGTTTTATCGCTTCCATGACGCAGAAGTTAACACTTTCGGATATTTCTGATGAGTCGAAAAATTATCTTGATAAAGCAGGAATTACTACTGCTTGTTTACGAATTAAATCGAAGTGGACTGCTGGCGGAAAATGAGAAAATTCGACCTATCCTTGCGCAGCTCGAGAAGCTCTTACTTTGCGACCTTTCGCCATCAACTAACGATTCTGTCAAAAACTGACGCGTTGGATGAGGAGAAGTGGCTTAATATGCTTGGCACGTTCGTCAAGGACTGGTTTAGATATGAGTCACATTTTGTTCATGGTAGAGATTCTCTTGTTGACATTTTAAAAGAGCGTGGATTACTATCTGAGTCCGATGCTGTTCAACCACTAATAGGTAAGAAATCATGAGTCAAGTTACTGAACAATCCGTACGTTTCCAGACCGCTTTGGCCTCTATTAAGCTCATTCAGGCTTCTGCCGTTTTGGATTTAACCGAAGATGATTTCGATTTTCTGACGAGTAACAAAGTTTGGATTGCTACTGACCGCTCTCGTGCTCGTCGCTGCGTTGAGGCTTGCGTTTATGGTACGCTGGACTTTGTGGGATACCCTCGCTTTCCTGCTCCTGTTGAGTTTATTGCTGCCGTCATTGCTTATTATGTTCATCCCGTCAACATTCAAACGGCCTGTCTCATCATGGAAGGCGCTGAATTTACGGAAAACATTATTAATGGCGTCGAGCGTCCGGTTAAAGCCGCTGAATTGTTCGCGTTTACCTTGCGTGTACGCGCAGGAAACACTGACGTTCTTACTGACGCAGAAGAAAACGTGCGTCAAAAATTACGTGCaGAAGGAGTGATGTAATGTCTAAAGGTAAAAAACGTTCTGGCGCTCGCCCTGGTCGTCCGCAGCCGTTGCGAGGTACTAAAGGCAAGCGTAAAGGCGCTCGTCTTTGGTATGTAGGTGGTCAACAATTTTAATTGCAGGGGCTTCGGCCCCTTACTTGAGGATAAATTATGTCTAATATTCAAACTGGCGCCGAGCGTATGCCGCATGACCTTTCCCATCTTGGCTTCCTTGCTGGTCAGATTGGTCGTCTTATTACCATTTCAACTACTCCGGTTATCGCTGGCGACTCCTTCGAGATGGACGCCGTTGGCGCTCTCCGTCTTTCTCCATTGCGTCGTGGCCTTGCTATTGACTCTACTGTAGACATTTTTACTTTTTATGTCCCTCATCGTCACGTTTATGGTGAACAGTGGATTAAGTTCATGAAGGATGGTGTTAATGCCACTCCTCTCCCGACTGTTAACACTACTGGTTATATTGACCATGCCGCTTTTCTTGGCACGATTAACCCTGATACCAATAAAATCCCTAAGCATTTGTTTCAGGGTTATTTGAATATCTATAACAACTATTTTAAAGCGCCGTGGATGCCTGACCGTACCGAGGCTAACCCTAATGAGCTTAATCAAGATGATGCTCGTTATGGTTTCCGTTGCTGCCATCTCAAAAACATTTGGACTGCTCCGCTTCCTCCTGAGACTGAGCTTTCTCGCCAAATGACGACTTCTACCACATCTATTGACATTATGGGTCTGCAAGCTGCTTATGCTAATTTGCATACTGACCAAGAACGTGATTACTTCATGCAGCGTTACCgTGATGTTATTTCTTCATTTGGAGGTAAAACCTCTTATGACGCTGACAACCGTCCTTTACTTGTCATGCGCTCTAATCTCTGGGCATCTGGCTATGATGTTGATGGAACTGACCAAACGTCGTTAGGCCAGTTTTCTGGTCGTGTTCAACAGACCTATAAACATTCTGTGCCGCGTTTCTTTGTTCCTGAGCATGGCACTATGTTTACTCTTGCGCTTGTTCGTTTTCCGCCTACTGCGACTAAAGAGATTCAGTACCTTAACGCTAAAGGTGCTTTGACTTATACCGATATTGCTGGCGACCCTGTTTTGTATGGCAACTTGCCGCCGCGTGAAATTTCTATGAAGGATGTTTTCCGTTCTGGTGATTCGTCTAAGAAGTTTAAGATTGCTGAGGGTCAGTGGTATCGTTATGCGCCTTCGTATGTTTCTCCTGCTTATCACCTTCTTGAAGGCTTCCCATTCATTCAGGAACCGCCTTCTGGTGATTTGCAAGAACGCGTACTTATTCGCCACCATGATTATGACCAGTGTTTCCAGTCCGTTCAGTTGTTGCAGTGGAATAGTCAGGTTAAATTTAATGTGACCGTTTATCGCAATCTGCCGACCACTCGCGATTCAATCATGACTTCGTGATAAAAGATTGAGTGTGAGGTTATAACGCCGAAGCGGTAAAAATTTTAATTTTTGCCGCTGAGGGGTTGACCAAGCGAAGCGCGGTAGGTTTTCTGCTTAGGAGTTTAATCATGTTTCAGACTTTTATTTCTCGCCATAATTCAAACTTTTTTTCTGATAAGCTGGTTCTCACTTCTGTTACTCCAGCTTCTTCGGCACCTGTTTTACAGACACCTAAAGCTACATCGTCAACGTTATATTTTGATAGTTTGACGGTTAATGCTGGTAATGGTGGTTTTCTTCATTGCATTCAGATGGATACATCTGTCAACGCCGCTAATCAGGTTGTTTCTGTTGGTGCTGATATTGCTTTTGATGCCGACCCTAAATTTTTTGCCTGTTTGGTTCGCTTTGAGTCTTCTTCGGTTCCGACTACCCTCCCGACTGCCTATGATGTTTATCCTTTGAATGGTCGCCATGATGGTGGTTATTATACCGTCAAGGACTGTGTGACTATTGACGTCCTTCCCCGTACGCCGGGCAATAAtGTTTATGTTGGTTTCATGGTTTGGTCTAACTTTACCGCTACTAAATGCCGCGGATTGGTTTCGCTGAATCAGGTTATTAAAGAGATTATTTGTCTCCAGCCACTTAAGTGAGGTGATTTATGTTTGGTGCTATTGCTGGCGGTATTGCTTCTGCTCTTGCTGGTGGCGCCATGTCTAAATTGTTTGGAGGCGGTCAAAAAGCCGCCTCCGGTGGCATTCAAGGTGATGTGCTTGCTACCGATAACAATACTGTAGGCATGGGTGATGCTGGTATTAAATCTGCCATTCAAGGCTCTAATGTTCCTAACCCTGATGAGGCCGCCCCTAGTTTTGTTTCTGGTGCTATGGCTAAAGCTGGTAAAGGACTTCTTGAAGGTACGTTGCAGGCTGGCACTTCTGCCGTTTCTGATAAGTTGCTTGATTTGGTTGGACTTGGTGGCAAGTCTGCCGCTGATAAAGGAAAGGATACTCGTGATTATCTTGCTGCTGCATTTCCTGAGCTTAATGCTTGGGAGCGTGCTGGTGCTGATGCTTCCTCTGCTGGTATGGTTGACGCCGGATTTGAGAATCAAAAAGAGCTTACTAAAATGCAACTGGACAATCAGAAAGAGATTGCCGAGATGCAAAATGAGACTCAAAAAGAGATTGCTGGCATTCAGTCGGCGACTTCACGCCAGAATACGAAAGACCAGGTATATGCACAAAATGAGATGCTTGCTTATCAACAGAAGGAGTCTACTGCTCGCGTTGCGTCTATTATGGAAAACACCAATCTTTCCAAGCAACAGCAGGTTTCCGAGATTATGCGCCAAATGCTTACTCAAGCTCAAACGGCTGGTCAGTATTTTACCAATGACCAAATCAAAGAAATGACTCGCAAGGTTAGTGCTGAGGTTGACTTAGTTCATCAGCAAACGCAGAATCAGCGGTATGGCTCTTCTCATATTGGCGCTACTGCAAAGGATATTTCTAATGTCGTCACTGATGCTGCTTCTGGTGTGGTTGATATTTTTCATGGTATTGATAAAGCTGTTGCCGATACTTGGAACAATTTCTGGAAAGACGGTAAAGCTGATGGTATTGGCTCTAATTTGTCTAGGAAATAACCGTCAGGATTGACACCCTCCCAATTGTATGTTTTCATGCCTCCAAATCTTGGAGGCTTTTTTATGGTTCGTTCTTATTACCCTTCTGAATGTCACGCTGATTATTTTGACTTTGAGCGTATCGAGGCTCTTAAACCTGCTATTGAGGCTTGTGGCATTTCTACTCTTTCTCAATCCCCAATGCTTGGCTTCCATAAGCAGATGGATAACCGCATCAAGCTCTTGGAAGAGATTCTGTCTTTTCGTATGCAGGGCGTTGAGTTCGATAATGGTGATATGTATGTTGACGGCCATAAGGCTGCTTCTGACGTTCGTGATGAGTTTGTATCTGTTACTGAGAAGTTAATGGATGAATTGGCACAATGCTACAATGTGCTCCCCCAACTTGATATTAATAACACTATAGACCACCGCCCCGAAGGGGACGAAAAATGGTTTTTAGAGAACGAGAAGACGGTTACGCAGTTTTGCCGCAAGCTGGCTGCTGAACGCCCTCTTAAGGATATTCGCGATGAGTATAATTACCCCAAAAAGAAAGGTATTAAGGATGAGTGTTCAAGATTGCTGGAGGCCTCCACTATGAAATCGCGTAGAGGCTTTaCTATTCAGCGTTTGATGAATGCAATGCGACAGGCTCATGCTGATGGTTGGTTTATCGTTTTTGACACTCTCACGTTGGCTGACGACCGATTAGAGGCGTTTTATGATAATCCCAATGCTTTGCGTGACTATTTTCGTGATATTGGTCGTATGGTTCTTGCTGCCGAGGGTCGCAAGGCTAATGATTCACACGCCGACTGCTATCAGTATTTTTGTGTGCCTGAGTATGGTACAGCTAATGGCCGTCTTCATTTCCATGCGGTGCAtTTTATGCGGACACTTCCTACAGGTAGCGTTGACCCTAATTTTGGTCGTCGGGTACGCAATCGCCGCCAGTTAAATAGCTTGCAAAATACGTGGCCTTATGGTTACAGTATGCCCATCGCAGTTCGCTACACGCAGGACGCTTTTTCACGTTCTGGTTGGTTGTGGCCTGTTGATGCTAAAGGTGAGCCGCTTAAAGCTACCAGTTATATGGCTGTTGGTTTCTATGTGGCTAAATACGTTAACAAAAAGTCAGATATGGACCTTGCTGCTAAAGGTCTAGGAGCTAAAGAATGGAACAACTCACTAAAAACCAAGCTGTCGCTACTTCCCAAGAAGCTGTTCAGAATCAGAATGAGCCGCAACTTCGGGATGAAAATGCTCACAATGACAAATCTGTCCACGGAGTGCTTAATCCAACTTACCAAGCTGGGTTACGACGCGACGCCGTTCAACCAGATATTGAAGCAGAACGCAAAAAGAGAGATGAGATTGAGGCTGGGAAAAGTTACTGTAGCCGACGTTTTGGCGGCGCAACCTGTGACGACAAATCTGCTCAAATTTATGCGCGCTTCGATAAAAATGATTGGCGTATCCAACCTGCA
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/tool-data/fasta_indexes.loc.sample Thu Mar 19 14:32:31 2015 -0400 @@ -0,0 +1,29 @@ +#This is a sample file distributed with Galaxy that enables tools +#to use a directory of Samtools indexed sequences data files. You will need +#to create these data files and then create a fasta_indexes.loc file +#similar to this one (store it in this directory) that points to +#the directories in which those files are stored. The fasta_indexes.loc +#file has this format (white space characters are TAB characters): +# +# <unique_build_id> <dbkey> <display_name> <file_base_path> +# +#So, for example, if you had hg19 Canonical indexed stored in +# +# /depot/data2/galaxy/hg19/sam/, +# +#then the fasta_indexes.loc entry would look like this: +# +#hg19canon hg19 Human (Homo sapiens): hg19 Canonical /depot/data2/galaxy/hg19/sam/hg19canon.fa +# +#and your /depot/data2/galaxy/hg19/sam/ directory +#would contain hg19canon.fa and hg19canon.fa.fai files. +# +#Your fasta_indexes.loc file should include an entry per line for +#each index set you have stored. The file in the path does actually +#exist, but it should never be directly used. Instead, the name serves +#as a prefix for the index file. For example: +# +#hg18canon hg18 Human (Homo sapiens): hg18 Canonical /depot/data2/galaxy/hg18/sam/hg18canon.fa +#hg18full hg18 Human (Homo sapiens): hg18 Full /depot/data2/galaxy/hg18/sam/hg18full.fa +#hg19canon hg19 Human (Homo sapiens): hg19 Canonical /depot/data2/galaxy/hg19/sam/hg19canon.fa +#hg19full hg19 Human (Homo sapiens): hg19 Full /depot/data2/galaxy/hg19/sam/hg19full.fa
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/tool-data/tool_data_table_conf.xml.sample Thu Mar 19 14:32:31 2015 -0400 @@ -0,0 +1,8 @@ +<!-- Use the file tool_data_table_conf.xml.oldlocstyle if you don't want to update your loc files as changed in revision 4550:535d276c92bc--> +<tables> + <!-- Location of SAMTools indexes and other files --> + <table name="fasta_indexes" comment_char="#"> + <columns>line_type, value, path</columns> + <file path="tool-data/fasta_indexes.loc" /> + </table> +</tables>
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/tool_dependencies.xml Thu Mar 19 14:32:31 2015 -0400 @@ -0,0 +1,6 @@ +<?xml version="1.0"?> +<tool_dependency> + <package name="vcflib" version="8a5602bf07"> + <repository changeset_revision="7e67466b033e" name="package_vcflib_8a5602bf07" owner="iuc" toolshed="https://toolshed.g2.bx.psu.edu" /> + </package> +</tool_dependency>
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/vcfgeno2haplo.xml Thu Mar 19 14:32:31 2015 -0400 @@ -0,0 +1,75 @@ +<tool id="vcfgeno2haplo" name="VCFgenotype-to-haplotype:" version="0.0.3"> + <description>Convert genotype-based phased alleles into haplotype alleles</description> + <macros> + <import>macros.xml</import> + </macros> + <expand macro="requirements"></expand> + <expand macro="stdio"></expand> + <command> + #set $reference_fasta_filename = "localref.fa" + #if str( $reference_source.reference_source_selector ) == "history": + ln -s "${reference_source.ref_file}" "${reference_fasta_filename}" && + #else: + #set $reference_fasta_filename = str( $reference_source.ref_file.fields.path ) + #end if + vcfgeno2haplo ${output_option} -w ${window_size} -r "${reference_fasta_filename}" "${input_vcf}" > "${out_file1}"</command> + <inputs> + <conditional name="reference_source"> + <param name="reference_source_selector" type="select" label="Choose the source for the reference genome"> + <option value="cached">Locally cached</option> + <option value="history">History</option> + </param> + <when value="cached"> + <param name="ref_file" type="select" label="Select reference genome"> + <options from_data_table="fasta_indexes"> + <!--<filter type="data_meta" key="dbkey" ref="input_bam" column="value"/>--> + </options> + <validator type="no_options" message="A built-in reference genome is not available for the build associated with the selected input file"/> + </param> + <param name="input_vcf" type="data" format="vcf" label="Select VCF dataset"> + <!-- Validators are commented to allow users apply too to any build. May need to be revised in the future + <validator type="unspecified_build" /> + <validator type="dataset_metadata_in_data_table" table_name="fasta_indexes" metadata_name="dbkey" metadata_column="1" message="Sequences are not currently available for the specified build." /> + --> + </param> + </when> + <when value="history"> <!-- FIX ME!!!! --> + <param name="ref_file" type="data" format="fasta" label="Using reference file" /> + <param name="input_vcf" type="data" format="vcf" label="Select VCF dataset" /> + </when> + </conditional> + <param name="window_size" type="text" size="4" value="30" label="compare records up to this many bp away" help="--window-size option (default = 30)" /> + <param name="output_option" type="boolean" truevalue="-o" falsevalue="" label="What to report?" help="--only-variants option" /> +<!-- <option value=" ">Output entire haplotype</option> + <option value="-o">Don't output the entire haplotype, just concatenate REF/ALT strings (delimited by ":")</option> + </param> --> + </inputs> + <outputs> + <data format="vcf" name="out_file1" /> + </outputs> + <tests> + <test> + <param name="reference_source_selector" value="history" /> + <param name="output_option" value="true" /> + <param name="window_size" value="5000" /> + <param name="input_vcf" value="vcflib-phix.vcf"/> + <param name="ref_file" value="vcflib-test-genome-phix.fa" /> + <output name="out_file1" file="vcfgeno2haplo-test1.vcf"/> + </test> + </tests> + <help> + +Convert genotype-based phased alleles within a window size specified by -w option into haplotype alleles. Will break haplotype construction when encountering non-phased genotypes on input. + +The options are:: + + -w, --window-size N Merge variants at most this many bp apart (default 30) + -o, --only-variants Don't output the entire haplotype, just concatenate + REF/ALT strings (delimited by ":") + +---- + +Vcfgeno2haplo @IS_PART_OF_VCFLIB@ +</help> + <expand macro="citations" /> +</tool>