Mercurial > repos > iuc > vsnp_determine_ref_from_data
changeset 7:57bd5b859e86 draft default tip
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/vsnp commit c38fd63f7980c70390d104a73ba4c72b266444c3
| author | iuc |
|---|---|
| date | Fri, 10 Jun 2022 06:10:23 +0000 |
| parents | 532a11cdd818 |
| children | |
| files | macros.xml test-data/add_zc_metrics.tabular test-data/add_zc_metrics1.tabular test-data/add_zc_metrics2.tabular test-data/add_zc_metrics3.tabular test-data/add_zc_metrics4.tabular test-data/add_zc_metrics5.tabular test-data/r1_seqkit_fx2tab1.tabular test-data/r1_seqkit_fx2tab2.tabular test-data/r1_seqkit_stats1.tabular test-data/r1_seqkit_stats2.tabular test-data/r2_seqkit_fx2tab2.tabular test-data/r2_seqkit_stats2.tabular test-data/samtools_idxstats.tabular test-data/samtools_idxstats1.tabular test-data/samtools_idxstats2.tabular test-data/samtools_idxstats3.tabular test-data/samtools_idxstats4.tabular test-data/samtools_idxstats5.tabular test-data/statistics_output1.tabular test-data/statistics_output2.tabular test-data/statistics_output3.tabular vsnp_add_zero_coverage.py vsnp_build_tables.py vsnp_get_snps.py vsnp_statistics.py |
| diffstat | 26 files changed, 482 insertions(+), 283 deletions(-) [+] |
line wrap: on
line diff
--- a/macros.xml Mon Dec 06 18:29:21 2021 +0000 +++ b/macros.xml Fri Jun 10 06:10:23 2022 +0000 @@ -1,7 +1,6 @@ -<?xml version='1.0' encoding='UTF-8'?> <macros> - <token name="@TOOL_VERSION@">1.0</token> - <token name="@VERSION_SUFFIX@">2</token> + <token name="@TOOL_VERSION@">3.0.6</token> + <token name="@VERSION_SUFFIX@">0</token> <token name="@PROFILE@">21.01</token> <xml name="biopython_requirement"> <requirement type="package" version="1.79">biopython</requirement> @@ -13,7 +12,7 @@ <requirement type="package" version="3.0.9">openpyxl</requirement> </xml> <xml name="pandas_requirement"> - <requirement type="package" version="1.3.0">pandas</requirement> + <requirement type="package" version="1.4.2">pandas</requirement> </xml> <xml name="pysam_requirement"> <requirement type="package" version="0.15.4">pysam</requirement> @@ -24,6 +23,9 @@ <xml name="pyyaml_requirement"> <requirement type="package" version="5.4.1">pyyaml</requirement> </xml> + <xml name="seqkit_requirement"> + <requirement type="package" version="2.2.0">seqkit</requirement> + </xml> <xml name="xlrd_requirement"> <requirement type="package" version="2.0.1">xlrd</requirement> </xml>
--- a/test-data/add_zc_metrics.tabular Mon Dec 06 18:29:21 2021 +0000 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,3 +0,0 @@ -# File Number of Good SNPs Average Coverage Genome Coverage -MarkDuplicates on data 4_ MarkDuplicates BAM output 10.338671 98.74% -VCFfilter_ on data 7 611
--- a/test-data/add_zc_metrics1.tabular Mon Dec 06 18:29:21 2021 +0000 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,3 +0,0 @@ -# File Number of Good SNPs Average Coverage Genome Coverage -Mcap_Deer_DE_SRR650221_fastq_gz 0.439436 8.27% -Mcap_Deer_DE_SRR650221_fastq_gz 36
--- a/test-data/add_zc_metrics2.tabular Mon Dec 06 18:29:21 2021 +0000 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,3 +0,0 @@ -# File Number of Good SNPs Average Coverage Genome Coverage -MarkDuplicates on data 4_ MarkDuplicates BAM output 10.338671 98.74% -VCFfilter_ on data 7 611
--- a/test-data/add_zc_metrics3.tabular Mon Dec 06 18:29:21 2021 +0000 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,3 +0,0 @@ -# File Number of Good SNPs Average Coverage Genome Coverage -13-1941-6_S4_L001_R1_600000_fastq_gz 0.001252 0.13% -13-1941-6_S4_L001_R1_600000_fastq_gz 0
--- a/test-data/add_zc_metrics4.tabular Mon Dec 06 18:29:21 2021 +0000 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,3 +0,0 @@ -# File Number of Good SNPs Average Coverage Genome Coverage -Mcap_Deer_DE_SRR650221_fastq_gz 0.439436 8.27% -Mcap_Deer_DE_SRR650221_fastq_gz 36
--- a/test-data/add_zc_metrics5.tabular Mon Dec 06 18:29:21 2021 +0000 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,3 +0,0 @@ -# File Number of Good SNPs Average Coverage Genome Coverage -13-1941-6_S4_L001_600000_fastq 0.002146 0.16% -13-1941-6_S4_L001_600000_fastq 0
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/r1_seqkit_fx2tab1.tabular Fri Jun 10 06:10:23 2022 +0000 @@ -0,0 +1,101 @@ +#id seq qual length avg.qual +XXXXXXXXXX:7:1101:1582:1835#/1 NGCTTCAGTTGGTACAGGTGTGTCACCCGTAGTTTGTAGGCTTGGAAGCCTGTAAGCTGTAGGGAACGTGGCAAGAACCAGATCTGGCCTTCCATATCTC BP\ccccc`eeccdeef_`fgabeghehdfhb^caefhee^efcehaeehhWabehddeM_eehhdeebd\^^a_]`___aaa`c]`Y^aaaaa_]]]aY 100 52.19 +XXXXXXXXXX:7:1101:1610:1859#/1 NTTATGTGTTACGTGAATCTCCTGAAGACAGCAGAAACTAGGTTGGTGAATTCTTATCCATTCTGCCATTCTGTATCTTTTAAGTGGAGCATTTAGGCTA BP\cceceggggggiihiiiiiiiiiiihihhiiihihiiiihiiihhiiihiiihhiiiihhiiiifhiiihgigggggggeeeeeecccddcdccbcc 100 52.73 +XXXXXXXXXX:7:1101:1743:1871#/1 NCGGGATGGAGGCAGGGCAAGCGGATGCAGGCCCAGCCCTGCCCTTACAGGTAGCTGGAAGTGACGAAGGAGTGGGGCCCCCGGACTCTGCTCAGGTGGA BP\ccceeggggghiiidfffhghfhhfffhhhidgggffbbghhheggeg`ccbccb_bdbdbbccaV^aaX^a`]X[_cca]]X]bc]Y]b`X`Y`_` 100 52.42 +XXXXXXXXXX:7:1101:1536:1878#/1 NGGCCTTTTTGGACAGGTGGTGTGTGGTGGCCTTGGTATGTGCTTTCTCGTGTTACATCGCGCCATCATTGGTATATGGTTAGTGTGTTGGTTAGTAGGC BP\cccedggggbffhieghdbgdeghffffhih[aaefdhefhhihhiiidffdffhgbbagaeebbdcddb_bcccbR]YXT]_`ab``[^b^]bcbb 100 52.41 +XXXXXXXXXX:7:1101:2990:100153#/1 GCCAGGAGTGACTGTGTTGTTCTTTGCTTTATATACATAAGCGCATCTCTTGCCCAAATAGAATTCTGTTTCATCTCGGGCGTAAACACCTTCAATTTTA ___cccecaegceg[dgfhgfbgdfgfhhfhhdfedhffhhhhhh[fhbbd_S^aaegfg`dfhgghhbcgddd``geadcc_a_a_][_c]bccccbba 100 64.07 +XXXXXXXXXX:7:1101:1624:1906#/1 NCCAGGTTTAAGGCATTACATCGAGGGGCTGAACTTGAGACAGTCCCTGAAAATGTAAGAACAGGAGCTAGAACTTGCTTATGGGAAGGGCCTAGGGAGA BP\cceeegggggiiiiihiiiiiiiiiiiiiiiiihiihihhiiiiiiifhiihhihdbfhggggeeeeddbddbdccccccbcbccccccccccccc_ 100 52.70 +XXXXXXXXXX:7:1101:1666:1926#/1 NGCAAACTGGTCAGGACCTGGAGCTGCTTGCCTGGCCTTGGCACAAGAGGCCGACATATGCCCAGCCCCTAGGAGTGACAGATCCGTCCTCTGTGGGAAC BPaceeeegggggiiiiiiiiiiiiiiiiiiihihiiiiiihiiiiiiiiiiiiiiiggggggeeeeccccccbbcaccccccccccccccccccccabc 100 52.71 +XXXXXXXXXX:7:1101:2921:100163#/1 CTCCGGATGCCAAAGGATCTGACAGTGTATCTAGCTTTGGAAAACACAGGGGTCTGCCCTGTGAGCTGCTCCAACACCTTGGCTGCTCGCGTCAGTCTGT bbaeeeeegggggiiihihhhiiiigheegfhhgdhihhhhhihhfhhhhbgggghihhdggihggeggeeeeeeacccccbccccccccc_a__bbccc 100 67.82 +XXXXXXXXXX:7:1101:1513:1929#/1 NCACTATCAAGGTTAGTGTCTATATCTACTTGAGANAAAATACTGAAAATTCAGTGTATGAAGCCAAATCCTGATTTAACAAGTTATTGGTAGTATAAGT BPY`acccceeeY`bfefeggehe`gddgceh_eWBPYbedeehfdefhhh]fhehbffhh_eeehded`cfhhdheeeeeeeddccccac`b]]ZZa`T 100 49.64 +XXXXXXXXXX:7:1101:2759:100170#/1 CTGCTCATCACATTCCAGATGCCATCACAGGCAATGACCATGAATTCATGGTCGTCAGTGAGAGTCAGCACCTTGATGTCAGGAAGGGCTGAAATCATCT bbbeeeeegggggiiiiihiiiiiiiiiiiiiiiiiiiiiiiihiiihiiihiiiiiighhiiihiiiiiiigggggggeeeeeddcccccabccccccc 100 69.31 +XXXXXXXXXX:7:1101:1708:1937#/1 NGGTCGGTGTCCTTTTTTGGCGTCAAACTTATAGAGGGAAGCAGGGAATCTACCCAACCAAAGCCAATGAACAGTCAGTCAGTCCATCTAGTTAAGGTTC BP\accccecggghiihihifhihffhihf`fegfghf`\ebcefgiigdgeggdeeec`abcbcbbbcbccccc`bccbccb```bbcccccc`bb`bc 100 52.64 +XXXXXXXXXX:7:1101:2981:100211#/1 GGGGTCAGTGGGAAAGGTGTCAATGGGACTTCGGTCTCCATGCCGAAACACCAAAGTCACAAACTTCAACTCCTTGGCTAGTACACTTCGGTCTAGCCAG __beccceeeeggdhffghhhhiihhcghffhiefaghhhfffffffhhiigifhf\cdbdddeeeebcbbd_bb_^bcbaccbcbbbbaa^a_]b`b`a 100 66.08 +XXXXXXXXXX:7:1101:1688:1946#/1 NCCGCCTTCTTTCTCTACTTCGTGACCTTTTGTGGTCCCGAGACCTGCTGCATCTTCTGTCTGATGTTCGGCTTGAGTGTCTACTTCCTGAACGACTTCT BPYac``c[cccchdd]caddeddehhhhhd`cdb^b]cc[^c_^_acccddchdhdcd\bdd]bd^bbbU\^BBBBBBBBBBBBBBBBBBBBBBBBBBB 100 38.51 +XXXXXXXXXX:7:1101:2767:100225#/1 CGCCTGGTGTCAATGCGCACATCTGGAGTTCCCATCTCCTTCATGGCAAATTTCCGAATCTCTTTGAGTGCCCGAGGTGCACGCTTCTTGAAGCCCACTC ___eceeefgggeiihihhhhihihihfgac_cefffihehiiiihhhhhfghhfhiehiihihf_ddbbgddcb`c]`bbbccccc`bb`b_bcccccc 100 66.85 +XXXXXXXXXX:7:1101:1536:1959#/1 NGATACTGTTGGGTCAGGACCAGGGCCTGCACCACCGAAATTCTGTCTTTTACAAACTGGACCAAGCGGTTGAGAATACAGGGGCCGAAGAGTAGGATTA BSYaa^ccggcg`d^dba_\ffh\[Zbf_aZ[c_cfY_ee`aeaf[\eefeghhffg_]\\`db^cdb`[ZWZR[`b]bbbbWELT]__TX^SYY]]]bb 100 50.36 +XXXXXXXXXX:7:1101:2797:100234#/1 GTGCTGGTCTCTCTCCACGATCTCTTACGTGCTAGCAGGGTGAACAAGGACGCCTCGCTGTCATTGCACTCACGCCACAGGAGCCTGGCACATTTCCGCT JY^accWceeeaWdcZ_feW`R^eg^a[HPP^_ef[d[Z_eUQX^ae_^Q\`ZSbUadVLKV_Z_Z_`]Z_baa^\WW^TGXTXWR^GW[][``BBBBBB 100 44.12 +XXXXXXXXXX:7:1101:1552:1976#/1 NTGATGAGTTTGCTAATACAATGCCAGTCAGGCCACCTACGGTGAAAAGAAAGATGAATCCTAGGGCTCAGAGCACTGCAGCAGATCATTTCATATTGCT BP\cceeeegegfiiiihiiiiiiiiiiihhiefgfhhhihifghdhhiiighfhhifihhiiiaghhigggggcdebdbddddcccccdedddddddcc 100 52.72 +XXXXXXXXXX:7:1101:1748:1978#/1 NCGCTGAGCAGGGCAGGCAGTGGTGGAGGCAGCCCCCCAAGACCCCATCTCTATCACAGTGAGAGCATCTGAACTCCAAAGCCACGACAGAACAGCTCCC BP\ceeeeggggghiiiihihiiegihhihiifiihhi`gfffggggeeeeeeddcdbd_bccbcccbccc`bbbbccccbccccccccaccbbb_`bcc 100 52.66 +XXXXXXXXXX:7:1101:2779:100239#/1 TGGGTGCTCTCCTTGGTCAGCTGCTGCTCGAACAGAGCACTTCTAGATACAGAGTTGGACCCAGACACAGCAGGGGTGGCCATGGCGGCTGGGAAAGACC ___^cc^ceggeeghddghafhXbegefbaYYbfhfgf_ae[egegc_cccaacaefg[bZbZ\\ZbZZ^gZZ^[]WW\_^bBBBBBBBBBBBBBBBBBB 100 40.42 +XXXXXXXXXX:7:1101:1593:1980#/1 NGGCACCAGAAGCATGAGCATAGGCAGAGCCATTTCCACCATGCTGACGGCAGGCAAGAAGAAGACAAAAATGAACCCCTTTCATGCTTGTTTAGAAATG BP\ceeeeggggfiiiiiiiihiihiiihfihhiiiifhhhiiiiihhiiiiiihhigggggeeeeeddcccccccccccccccdcccdccccccccccc 100 52.71 +XXXXXXXXXX:7:1101:2946:100242#/1 TGGGGATGTTGTTAAAAAAAATACAGGCTCCCCCACAACTGGGGTGCCTGGGGGGAACTTGGTCTGCTTCAGCCCAAGAGGAATCAAAAGATCAAAAGCA ___eecccgggggihhfhihfdfhhiihhiifiiddbcd`\eghL\^a^bdccccTTWRY`bRYYXSYb``bbca`acW^a`^bbccccc_bb`bbcbac 100 58.29 +XXXXXXXXXX:7:1101:1987:1781#/1 NTTAGCAGGTCCTCGTAGCGTTTCCGCAACTCTCGGAATCTGGAGTCTTCGGAGTGCAACTCTGAGGGTCCCGGGAAACTTGCGCGGCTCGCCTCGGCCA BP\cceeegggggiiiiihiiiihiiiiiiiiiiiiiiifiiiihhghiiihiggfgeeeeeebddcccccccccaccccccccccccccaacccXa^a[ 100 52.65 +XXXXXXXXXX:7:1101:3046:100006#/1 CTGTACAGCATGGATGTGGACAGGGAAAACAGCAAGTGCACACAGTGGGACAGGCCACAGGACAGGCTGGGAGTGAATAAAGAGTTCACACTGCTTCCCT bbbeeeeegggggiiihiiiiiiiiiiiiiiihiiieghiihhihffghiiiiiihiihhiiiggggeeec\bZ^bcbccc_]^]_b`bc`bbccbcccc 100 66.70 +XXXXXXXXXX:7:1101:1900:1788#/1 NCCCGCTGCCCACGCCACTGTCCCGGCGCTGCCCTGGGCGCGGAATGCGCGTTAGGCAATCCCCGGTCTGCGCGTCCCACACGCAGATGTGGCCTGCCAG BP\ceeeefggfgiiihiihiiiifhhiiihhiihiiihgeecccccccccaccccccccccccac_accca]aa[aacccccccccccdbccaccccb_ 100 52.64 +XXXXXXXXXX:7:1101:3214:100027#/1 GGCTGAGGTACTGAATGAGGAACAGCAGGAGGTAACGTCTGAGGGTGCGGCTGGGGAACGGTTTGCACAGGGGTGGCTGTAGGAGCAGGAGGTGAAACAC ___eecee^cegghhhhhhfhhfgfhhfhhhh^eggfcefhffhbU^eghhehgggTXZ\aTTX\RRYTT[GWROTWW[GG]YGXW[baOQEQW]b`]XW 100 49.49 +XXXXXXXXXX:7:1101:1848:1879#/1 NCTTCTTCTTCCTCATGTTTTTTGATCAACTGCTCCACTTCTTTTTGTTTGAGAACTCTGATTACTGTCTTTCCATTCTCTCTTGTTAGTGTTGCAATTT BP\cceeegggggiiiiiiiiiihhihhiehfhihhhiiihiiiihfghihdggfhhiiihihh_ddcggeefdb`beeeeeddcddbdbbcc`bbcccc 100 52.70 +XXXXXXXXXX:7:1101:3237:100032#/1 CTCCCGAAGTGCTGGAATTACAGGCATGAGCCCCCACGCCTGGCCTGATCAAGCTGTTTCTCTTGCTTGGAATGCCCTTTGCTGCCTGGCCAGCTCCTAG _b_eeeeegegggiiiiiihiiffegfhhihiiiiiighiffhiiiiiiiiihhiibdeggggeeeeeebcbdbcbcccccccccccccccccccccccb 100 67.87 +XXXXXXXXXX:7:1101:3027:100049#/1 CACAGAAACACGGGTTCATGATGGTGAATGATATAGCGGGTGACAGATGGCTACAAACCAGTCATAGGCCATCACACTCAGGAGCATGTCTCTCTTCCAT bbbeeeeegggggidghiihhihidgfhiihihfiiiiii_eghiiiihihhiiiggggeebdddcddddcccccbbcccccccccccbbdccccccbbc 100 67.75 +XXXXXXXXXX:7:1101:1756:1891#/1 NCGATGTTGAAGCCTGAGACTAGTTCGGACTCCCCTTCGGCAAGGTCGAAGGGGGTTCGGTTGGTCTCTGCTAGTGTGGAGATAAATCATATTATGGCCA BP\cceeegggggiiiiihiiiihiiiiiiiiiiiiiiiiiiiiigiiegghige]_ccaaccc_bbbbcccccc_bbcbcbccccdccccddddccccc 100 52.67 +XXXXXXXXXX:7:1101:3238:100065#/1 GGATGATAGAGTGTGGGAGAGAAAAGGGAGAGAGACTTTGGAGCTGAGTGGGGACCGCCGGGGAGGGGAGGCCTCAGAGTGGCTCAGCTGTAGATGCGGC bb_eeeeegggggdfgifhghhiiiiiihghhfgdhihiiifhhiiiigghiiiiigeecccc_accc_ccccbcbbbc_bccccbcbcc_bccdcbcca 100 67.15 +XXXXXXXXXX:7:1101:1915:1901#/1 NTAGAGGCTTACTAGAAGTGTGAAAACGTAGGCTTGGATTAAGGCGACAGCGATTTCTAGGATAGTCAGTAGAATTAGAATTGTGAAGATGATAAGTATA BPY`cceegfggfihgghhhegfghfgihdfhhhhhhfeffdgighiiiiifhiifdgdgecddc_cbbbdbbdcb__`bcc__`_`bbcbbccccac`c 100 52.64 +XXXXXXXXXX:7:1101:3198:100082#/1 CCCCCGCTCAACTACAAATGTCAATATTGAATGAAGCATTAAAAGACAAACATAAAGTAACTTCAGCTTTTATCTAGCAATGCAGAATGAATACTAAAAT b__eeeeeggggfiiiiiiieghhffhfhiiiiiifhiiiiiihgfhhiiiiiiiiigghhihiiiifegfggeeeeededcdcddcccccddcccbbcd 100 68.63 +XXXXXXXXXX:7:1101:1964:1931#/1 NGGGGTCTTAGCTTTGGTTCTCCTTGCAAAGTTATTTCTAGTTAATTCATTATGCAGAAGGTATAGGGGTTAGTCCTTGCTATATTATGCTTGGTTATAA BP\`Z`cdggggfhhhhghfghhhhhhhhhhffhhhhhhhhdggfhhhhhhhhghhhfgfceghfhhhhdgdgdggbggfgfebdeeeedcddbbccccc 100 52.70 +XXXXXXXXXX:7:1101:3088:100091#/1 TCCCACGTTTCCCACTTTAAAACAGAAGACTGCAGAGGATGTGGTAGAGACTGCAACTGCCAATGTAGCATCCGTTTTCCCAATTTCTCAGTAACAGGAT bbbeeeecgggggiiiiiiiiihffhihiiiiifgdghbghggighhhihhiiihiiiihiiiiihiiiiggggebcecdcbddcccdcdbccccccc`b 100 68.56 +XXXXXXXXXX:7:1101:1840:1948#/1 NCAGGGTGAGGGGCTTCGTCTGCTGGATGTACTTGTTGCCGAACTGAGTGACCCCCCGGGGCCAGCCGGTCTGCGAGCGATTGGGCGGTACCACAGACAT BP\ceecegggggiiihhhhiiiiiighhfhhiifhiiiiiiiiiiiihiiifhhgccccccacccccccccccccccccccccccccacaccccccccc 100 52.69 +XXXXXXXXXX:7:1101:3105:100094#/1 GGGAACTCCAGGGATGAACATGAGAAGGCCCTTCTGCTGTAGTACTAACTGCAATCAGTGTCTCGGGCTTCATGGGCTTGGCTTTACGTAACAAGCCAGG ^[[^ccccc^eccgf\R`bggfffbffhhfhbXc]Xcfcg]ece_cegfadafhgfheNWW\_]eZFZW`^HVZ]UZ`ZKZGT\b_b[KYXGY_XQWbW^ 100 50.33 +XXXXXXXXXX:7:1101:1958:1952#/1 NTCATGGTAGGGGTAAAAGGAGGGCAATTTCTAGATCAAATAATAAGAAGGTAATAGCTACTAAGAAGAATTTTATGGAGAAAGGGACGCGGGCGGGGGA BPY`cccc^aaecf`e[dbeeedaaeVefhhhdeXc_fehdb_eedhhhed\acbefdhdddee`bbeeeeeacca]_aa__]]_aaaaaa]aaa_]XaT 100 52.40 +XXXXXXXXXX:7:1101:3190:100106#/1 CACGGTTAGGGTACCGCGGCCGTTAAACATGTGTCACTGGGCAGGCGGTGCCTCTAATACTGGTGATGCTAGAGGTGATGTTTTTGGTAAACAGGCGGGG bbbeeceegggggiihhhihifhiihiiiiifgeghihhhhhfhgfecR]_bcccccdcccdc]b]bbccccbccc^`bc_bcccaccccccccacccac 100 64.50 +XXXXXXXXXX:7:1101:1993:1999#/1 GTTTAATTGTGGGAGGCTGAGAGCAGCATTTTTACCAAATTGGTGTAAAAATGAAACGGGGTCCAGACGTGTACACTGACAGTTACAGTAAATTTCAAAA bbbeeeeegecgfdhhiiihihiiiiiiiihiighihhfhihhcf_ehhhiiiiihihhggadgggeeeccddddbccccccccccccedccecccdccc 100 68.03 +XXXXXXXXXX:7:1101:3117:100110#/1 GTTTTATGTTGTTAATGTGGTGGGTGAGTGAGCCCCATTGTGTTGTGGTAAATATGTAGAGGGAGTATAGGGCTGTGACTAGTATGTTGAGTCCTGTAAG ___eeeeeggggfiiiigihghiibfcgefdgfhii`ghiefefhghhaedghdhibfghfhiff^`gbggeeeedbcbdbbbccdedcccbbbcbbbcb 100 67.19 +XXXXXXXXXX:7:1101:2159:1798#/1 NCTTCCAGTCCCCTCCTCGATTTCACACCGGGACACACACCTGGGTCAGGTGAGCATGACACACAGCCCCCAGGATGGGCACCCCACCCTGCAGAGCAGA BP\ceeeegggggiiiiiiiiiihiiiiiiiihhiiiifiiiiiifgiiiggggggeeeeeddccccaccccccccccccccccccccccccccbccccb 100 52.69 +XXXXXXXXXX:7:1101:3147:100111#/1 AGAAATTATTTTTAAATGTCAGAGAGAATGCTATTATTGGCCAACTTGAAATTTCCTAATTAGCTATCTGAGCTATTTAAGTCTGTACTTCAGTGTCTTC bbbeeeeegggggiiiiifiiiiiiiiiiiiihhiiiihhiiiiiiiiiiiiiiiihigiiiiiiiiiiiihihfcghhihbggggfgeeeeea_bdddd 100 69.60 +XXXXXXXXXX:7:1101:2152:1838#/1 NGGGAAACGCCATATCGGGGGCACCGATTATTAGGGGAACTAGTCAGTTGCCAAAGCCTCCGATTATGATGGGTATTACTATGAAGAAGATTATTACAAA BP\cceeegggggiiiiiiiiiiiiighigggggeeecccdccccdcceccccccccccccccccccccdccc^bbccccddccdcccccbcccddcccc 100 52.66 +XXXXXXXXXX:7:1101:3065:100152#/1 CACATAGTCAGCCTGCTTGATCACGTCGATCTTGTCCAGGAAGTACTGGGCACAGTCAATCAGCTGGTACTCGTTGGAGCGTTCGTAGCAGGCACGCACT ___ecacccggcgfhhhhhhh_gfdeghhhhhhhfhhdcdhXfX^f_eghhfUegbfh_be__dhbb`g`bga_aZ^]bdaX\^a^`aa`_[_aW]][]a 100 62.63 +XXXXXXXXXX:7:1101:2180:1843#/1 NTAGGGTGTAGCCTGAGAATAGGGGAAATCAGTGAATGAAGCCTCCTATGATGGCAAATACAGCTCCTATTGATAGGACATAGTGGAAGTGGGCTACAAC BP\cceecgggfehiehihhihhhhdffhihifdhhifghhhiiihihihiiihiiiiihiiihfhiffefgfdegeecac`d]bddbb_bcaacccccc 100 52.70 +XXXXXXXXXX:7:1101:3154:100159#/1 GTTTCATCATGTTGGCCAGGATGGTCTCAATCTCTTGACCTCGTGATCAGCCTGCCTTGGCCTCCCAAAGTACTGGGATTATAGGCGTGAGCCACCGCGT bbbeeeeegggggiiiiiiigiiiggiiiiififhiggiiiiifhiiiighiiiiiiiiihiiiiihhffffgggdgeeeeeedbdcaccaccccccaca 100 68.83 +XXXXXXXXXX:7:1101:2125:1861#/1 NGTGGGAGCAGCTGGCTGTGTGGCACCTGCTGTGGTGGCTGCTGGCGTGGCCGGAGTGAAGGGCAACGTGGCAGGTGCCGGGGGCTGGGCTGGATTCCCT BS\cceeegggggiiiiighghihfhifhhihfgfeghfbffgffaf_efffegX]Z]]_bccPZ^`aX`a]`acBBBBBBBBBBBBBBBBBBBBBBBBB 100 38.83 +XXXXXXXXXX:7:1101:3198:100173#/1 TGGCAGGTGGTGGTGGTGCTCAAAGAGCAGGTCCAGGTGGAAAGTGAGCACCGACAATGGCTGCAGCAGGAGCAGCAGCTCTGTGGACAGGGAGGGACAA ___cccccece`gbZaffgadfgg_ddffghbafffdaffffhhedghhhhhffhgggadc_`bbbba^a\\XXT^[^[Y_]`Y]b]`b__aaBBBBBBB 100 44.48 +XXXXXXXXXX:7:1101:2076:1911#/1 NGGAGTAGGAGTACCCAGTCTCTAGTGCCAGCTGCAAGCGTTAATTTGGTTAAAGTCTCCTTGATTGTCCTGTTCATTCTTTCTACCTGACCTGAACTCT BS\ceeeegggggiihiiiiihiihgiiihiiiiiiiihhhihghiiiihhiiiihihiihiihiiihhiiihfggggggeeeeabcdbdddcccccccc 100 52.81 +XXXXXXXXXX:7:1101:3166:100190#/1 CCGGCCGAAATTGTTTTTTAAATTTACGTGTCACCCTTCATGCAGGGCCCTGCTACTCTTCGTAGCCCTCCTCCATGTTTAGTGTGGGTGCGGCTGCAGC bbbeeeeeffgegdgigfiihhiiihiifggghffc`gfffgheifffagdbddgdbeced`bdcaccbaccccb`b`bccbY`bccaX]_accaa^R^b 100 63.87 +XXXXXXXXXX:7:1101:2196:1920#/1 NTGGGGTGGGTATGGCGGGGCCAGCCTTGCTCAGTGCTGTGGTGATGGCCACAGTAGCCTCGGGGGGCACAGCCACGGCTGGGGGAGCAGGCCCTTCCTT BP\ccc`ceeceehhhhhhh[aZZZ^U]LZ]]ZTZ_]]TTY_KTT]]`[a]_[[R]]R__aR[aaaaTEW^XP[GRLTTOOXaaBBBBBBBBBBBBBBBB 100 40.29 +XXXXXXXXXX:7:1101:3225:100207#/1 GCTATTTTCTGCTAGGGGGTGGAAGCGGATGAGTAAGAAGATTCCTGCTACAACTATAGTGCTTGAGTGGAGTAGGGCTGAGACTGGGGTGGGGCCTTCT bbbeeeeeggfgghiiihiYcghiihiihiififhfghfhhhhhhggggggeeeeeddddddcdcac]`ababbbbacccccccbcccc[aaacaccccc 100 66.30 +XXXXXXXXXX:7:1101:2115:1927#/1 NCTCCTTCTTGGCCTGGAGGCGCTCTTCACGGCGCTTGCGTGCTTCCTTGGTCTTAGACCTGCGGGCCTCAGCCTGGTCAGCCAGGAGCTTCTTGCGGGC BS\ccecegeeggiiiffggfhhiiiihhdhihfdghddfggeeeeecdddZ]bcdcbcccbccacaaccccccc^`YYYX`bW`_`a[`_b``_baccT 100 52.51 +XXXXXXXXXX:7:1101:3019:100219#/1 GGCGGGTTTTAGGGGCTCTTTGGTGAAGAGTTTTATGGCGTCAGCGAAGGGTTGTAGTAGCCCGTAGGGGCCTACAACGTTGGGGCCTTTGCGTAGTTGT _[^ccc_cggccgghdgcbgdf_WcXae[cbefffgbefh_d`dgddb]aaZTZW__cc`bb]ZZ_^aaZ]a^`b_Y__ab[`_XT_^bbbbaBBBBBBB 100 44.46 +XXXXXXXXXX:7:1101:2179:1937#/1 NAGGAATACATATAGTCATCATTGCCAGACTTAATATGAGAGGTGAAATGTTCGATCCAATTTTCCTTCCTGGATAAGTTTTTCTTTCCTATCCCTCTCA BP\cceeeggggghigfghiiiiiiihiiiiiiiiiiifhhhheghhiihhhihiiiiiiiiifhghhiiiiihhgffdgfgegeeeeeeeddcddcccb 100 52.73 +XXXXXXXXXX:7:1101:3202:100230#/1 CTGAAAAACAAAACAAAACTTAGCTAAGTTTTATGGCAAGTGGGGGCAGGGTCAGGTGAAGAGGTAGGAATTGCCTGAGGTTGGGGGTCCCTAAAACATA __aeeeeegggggihhhiiiiihihhhhfghhhihhhiihfffgdfhhihiddggeceeeeedcbbdcbbbcbcc]bbcc^abccccT[acbbcccb^`b 100 65.28 +XXXXXXXXXX:7:1101:2149:1945#/1 NGTGTGGTGGCCTTGGTATGTGCTTTCTCGTGTTACATCGCGCCATCATTGGTATATGGTTAGTGTGTTGGTTAGTAGGCCTAGTATGAGGAGCGTTATG BPYaccccegggghhhceghghhhhhhhhhfhhhhhhhhhhhhhffhghhhdffghghhbfghgfddgddcdddbebcbbbbbbbdcdbbbbaaaaaaab 100 52.67 +XXXXXXXXXX:7:1101:3211:100242#/1 CTGTTCCCCAATACCTGCAATTCACAGATATAATTTAGCATTTTTGTTTTTATATACTAATTTAGGCAAAACCCCCCCGCCACTGAATCGTCCAACAAAA ___`acacgggcghfdddc_gbggghddgfgfhhccdab[e[_Xbe^^cacdfgbaa[caSWWbfbf`efc_Zbg[XZaccaacccbcc[[WW^bbW[[a 100 61.18 +XXXXXXXXXX:7:1101:2169:1964#/1 NGGGAGCAAAAAGCGTGGGTACAGAAGCGAGAAGCGAGCTGAGTGGTTAGTTCAAATCAAGGCTTGGGGTATTTCCTGGTCATTTGGGGAACCTGAAACT BP\ccceeggggfehfhhh_fhifbgf[effhh`cfV__`gefccg`ddg`d`beeeadddbcccaaacNWY_bZ]`bb_YY__aGSTXGTW^R^R_b_` 100 50.12 +XXXXXXXXXX:7:1101:3168:100244#/1 CCCTTCCTGTGTCCATGTGATCTCATTGTTCAATTCCCACCTATGAGTGAGAATATGCGGTGTTTGGTTTTTTGTCCTTGCGATAGTTTACTGAGAATGA ___eeeeegegeghiihfbgghhhhfhaggiiiidbfhfhfbdegdfcffgdgfhhfhhfceefghicfffeeceddbddba^acc_bccccc[bbbcc_ 100 66.77 +XXXXXXXXXX:7:1101:3005:100246#/1 CAAAAGTTAAATGAAAAACCACAAGTTTAGTTTGACTCCGTCTCCTAGGGTGGATTTCATTCAGATATTTGTTCCATATTATAGGAGGGTGGATCCTAGC bbbeeeeegggggihhhihihiiiiiiiiihhihhhihiighiiiiihii^efaegfffhhfhhhiihhiidgggggdeeeeeeedccaZ\_accccccc 100 67.69 +XXXXXXXXXX:7:1101:2313:1789#/1 NTGCGCCTGGAAGCCGCCAGGCCGCAGCGCGGGGAACAGGGCGGCCACCCGGCTCAGCAGCCGCTCGCAGGGTCCCGGACCCGCCGCAGCCTCGGCTTCC BP\cceeegggggiihiiiiiihiihiiifhheccccccaccccccaccccccccccbcccccc_a]Zaaccbb`aaacaccaccccaaacccX]^^abc 100 52.58 +XXXXXXXXXX:7:1101:3253:100014#/1 GCGGATGATTCAGCCATAATTTACGTCTCGAGTGATGTGGGCGATTGATGAAAAGGCGGTTGAGGCGTCTGGTGAGTAGTGCATGGCTAGGAATAGTCCT bbbeceeegggggiiihihiiiiiiigfhiafbeccb_fdgaffhhhhhihihihhiieT]abcccccccccbb_WY]_bbbbccccccccbcbb`Y`bb 100 64.18 +XXXXXXXXXX:7:1101:2361:1794#/1 NGTTAGGCTCTCCACGACTGAGCTGGTCTCGGCACCAAGGGGGTCAGATTGCTTAAAGGCAAGGAGTTCAGGACCAGCCTGGGCAATGTGGGGAGACCTC BPYacceeegggfhhhhfhhfhhb_eaefdgghhfgfffgbfcfdeggfgeeeeeddbd\``aaabbbbdbbbb^[[^aba_aababb`_Y[aaaR[a`b 100 52.45 +XXXXXXXXXX:7:1101:2337:1794#/1 NGGGAAGGTGTCTGTCATCGTCCCCAGGAGGGAGGTGGCGTTAGCTGTTTGTCCTGTCATTAGGTTGGTAACTCTCCAAGTAACATTGAAGATCTGGTGA BP\c`ceceggggiihiiihiffffghi`ghhggi^ce`f_faghihbghigegdgR_ceeeeeU]`c^`_bcbbccbbb]`bcbcb`b`cbbbc_bbcc 100 52.47 +XXXXXXXXXX:7:1101:3284:100039#/1 TTGGGCGACAAGGACATGGAGCTGGAGGAGAGCGTGGAGGAGCGGCTGTTGAGTTCAGAGCCCGAGTCCACTGCCACCGAGGCCGGGCTGGCGGCCAGGG bbbeeeeegggggiihiiiiiiiiihhihiffhifhiihhifhigeec`db`bbccdcccbcccccacccccccc`bcc_acacccccccca_a[_BBBB 100 46.93 +XXXXXXXXXX:7:1101:2477:1795#/1 NGGCGGTGCACCGAGCTGGAGGCTGCCGCAGAGGAGGCGGCCGCCGCTGCGGCCATGTACGTGTAGGGGTAAGGGAACAGGCTTCCGAAAGGGGACATGG BP\ccccefeeggffghffhhihhhhhfifhhiidghfddcca_aWaWac_accabccdcaX`R^bba_X_bcccW`^acc_^[X`_[[[a_acacc]][ 100 52.30 +XXXXXXXXXX:7:1101:3310:100056#/1 GGGCGATTGATGAAAAGGCGGTTGAGGCGTCTGGTGAGTAGTGCATGGCTAGGAATAGTCCTGTGGTGATTTGGAGGATCAGGCAGGCGCCAAGGAGTGA bbbeeeeegffgfiihiihiighiiihiihhiii\edgfhiddggeggeeeedbddcdbdcccbbc^a`bcdccb`aabbaccb^acccacccca`cX`b 100 65.89 +XXXXXXXXXX:7:1101:2355:1821#/1 NTGTATTGAAGCCAAATTTCGTGCCAGGGTGGTGTTAGAATGTGAGCTATTTTCATTATATTTGGACTAGAATTTTTTCTTTTAGGAGAATGGTGGGGAA BS\cceeegggggiiiiiiihfgiihhii^egbffhfhfhhiehhhiihiiiiifgh`ffghiiighfhhgggfgggcdddddcbbbc`abcc_bcccac 100 52.78 +XXXXXXXXXX:7:1101:3420:100060#/1 GTGAAAAATGTTTAATTTACTTCAAAAATATCTTGCATAAGTATCTCATCCCTAAGCTCAAACAAACAAGAAAATTAAAGGCGGAAAGAATACATAAAAT ab_eeeeeggffgihhhihgfdhhhhihf]ghiihefhhig^fgbghbfdhhiihbffhhiihiiihhhhhcfhigfhdg`dFXZaccbcbcccddcbcc 100 56.56 +XXXXXXXXXX:7:1101:2418:1834#/1 NCGGTTTCAGTCGTCACTGCTGCCCCAGCCTTTCGCTGTCCTTCTTGCAAAAAGCTGCTCCCATCGGTATACCAGGTGTGGTCGGCGTCTGGGATGGGCT BPYcaceeggggfiiiiiiiihifiiihiiiiihdfhhfhfhiiiiihiihhhhhihihhiiegggeeecbeddcabbbcccccccc[aaccc`bccacc 100 52.67 +XXXXXXXXXX:7:1101:3267:100061#/1 CTGGGTAGTTCTGGCTCCGGGTCTCTCATGAGGTTGCAGTCACGATATTGGTTGGAGTTGCAGCCATCAGAAGGCTAGACTGCGTGTGAAAGACTCACTT _bbeeccegggggiiiiiihigghihifhiiiiffhiiihhihidghiihheghhhighgggfeeeeeeccddcccccccccccaacccdcccccccccc 100 68.01 +XXXXXXXXXX:7:1101:2378:1838#/1 NGCCCAACCAACAGCTTCAACAGCTATTGCAAGGTGCGCCAAAGGCATCTTGTCATCCCTCACACGAATCTCACTTCCTGTGAATTTGCAGGGAGGCAGA BP\ceeeegggggiiiihiiiiiiiiiiiiiiihfhiihiiiiiiiiiiiiiiiiiiiiiiggggeeebcddccccccccccddddcccdcccccccccc 100 52.71 +XXXXXXXXXX:7:1101:3416:100083#/1 GCCCAGTCCTGGGGGCAGGGGAGAAAGGGAGGGGTGAAATGAGAAGAGCCCCGTTCTTCCTCAGTGCTGGGGGGTTGGAGACACCACCCCCTTCTGCTGC ab_eeeeeggggfiiiiiiiighghhiiihhhiiL\bgeeeeeeddbddcccc^acacccccccbccccccccc]X_abcbcccaccacccccccccccc 100 61.17 +XXXXXXXXXX:7:1101:2481:1853#/1 NGCAGGTAGAGGCTTACTAGAAGTGTGAAAACGTAGGCTTGGATTAAGGCGACAGCGATTTCTAGGATAGTCAGTAGAATTAGAATTGTGAAGATGATAA BS\cceccgggggiiiiiiihiihhfhhihiiifhgiiiiiihiiiiiiiiiiiiiieggeeeeeebdddbddccbbccccccbbbccbcbcccccccdb 100 52.78 +XXXXXXXXXX:7:1101:3411:100111#/1 GTGGTGGCCTTGGTATGTGCTTTCTCGTGTTACATCGCGCCATCATTGGTATATGGTTAGTGTGTTGGTTAGTAGGCCTAGTATGAGGAGCGTTATGGAG babeeeeegggggghiihiiiiiiiihdfgiiiiiifiiiiiiihiiiiefhhiiiegfhggghgggfbddeceeedbdbb_bedcccccccccccccba 100 68.20 +XXXXXXXXXX:7:1101:3258:100128#/1 GCTGTGGCTTTTTTAAGGATAACTACCTTGGGGGCCTTTTCATTGTTTTCCAACTCCGGGATCTGGTCACGCAGGGCAAAAAAGCTCCGTTTTAGCTCGT _b_ecceeggggfihghifffhhfghidghhiihighiihighiighhiifdffggggeecccccbccbccccca[_aaccccccbbbcacbcbbcccc[ 100 66.63 +XXXXXXXXXX:7:1101:2252:1856#/1 NCAGAAATATGGCAGTGCCACAGAGAGTGTGTCCTCATTGGAGTACTTGCGCTTAAATTCATCCAACACAAAGGTACTCTTGGGCAGGTGAGCGAAGGGG BPYcaccce^ecc^e`begg_a_fXQ`[Q_eb_b`d]ee^[ac^e^e[cbe^edeedVWW_bb\_\dd`bdGV^aaab``UUG\\\]_T_]aaaZO[T]a 100 50.04 +XXXXXXXXXX:7:1101:3428:100129#/1 GTTGGAGTGTAGGATAAATCATGCTAAGGCGAGGATGAAACCGATATCGCCGATACGGTTGTATAGGATTGCTTGAATGGCTGCTGTGTTGGCATCTGCT ___eeac\cecggfhhihihfhhiiiihiiihfhffihhhhihhidg`fgfgeeccccZ^^\]bbbc_bb`bbccabcbb`ba[]`YbS]Q[RSRX`bb_ 100 60.12 +XXXXXXXXXX:7:1101:2394:1871#/1 NTTTATTTCTAGGCCTACTCAGGTAAAAAATCAGTGCGAGCTTAGCGCTGTGATGAGTGTGCCTGCAAAGATGGTAGAGTAGATGACGGGTTGGGCCAGG BP\ceeeegggggiiiiiiiiiigghiiiiiiiigihihhiiiiiiiihighiiiiihgeggggeeeeedddddb_bccbcccdcccccc_ccccccccc 100 52.70 +XXXXXXXXXX:7:1101:3387:100138#/1 TTGAGGGAGGGCAGACACAGGCTGAGGAGAGTGGATAGGAGCAGGAGAACCTGCCACACCCCACACCCCACTGAGAGCTCATCAGCTGCAGCCCTGGGCC bbbeeeeegggggiiiiiiiiPdgfhihhfgXafffhfhhhiiiicfgifhhhegegceeccccccccc_aacc``__cccccbbcccc`bcccc`accc 100 63.66 +XXXXXXXXXX:7:1101:2269:1904#/1 NGCACCTGTAATCCCAGCTACTCAGGAGGCTGAGGCAGGAGAATCACTTGAACCCGGGAGGTGGAGGTTGCAGTGAGCCAGGGTCAGGCCATTGCACTCC BPY`cceceaga^dadghhhgfhhbfghhhhhhhhfhgad`ec[ecedbbbecbdghgaaXKU^GZ\GTX`bRY]__bbT^aaRRWbb^a^]bbbbbb^` 100 49.98 +XXXXXXXXXX:7:1101:3444:100163#/1 TGTGGCTGTTCTATGAACCACCTTCTTCTTTCTGCGAGCAGTTCCTTTCCCACCAATGCGCACTTGTGCCTGCAGTTTGGCGAGGTTTTCCTGGGTCATG __\c\cc`^`c^be``bad`ddegaX`dgegfaIXW^^OW^O^ce]a^ae^egffgd_VZ_FT\ZV_]c`bb_bX]Z]R^baQREHW`WY_bYXGX^^BB 100 46.94 +XXXXXXXXXX:7:1101:2259:1943#/1 NTCTTCCTCATCTTCTTGTTCTTTTCTTTTTCCTTCTCCTTCTCCTTTTTCTTTTGTGACATGATTTCACTCTTGTTACCTACGCTGGAGAGCAGTGGTG BP\accc`b[QbKb`cd_``dedddcddYXecSbdI^P^XcSXbX^Y^^c^a^cc_HHXacIXWM\WWWa\HW_]NV\``_BBBBBBBBBBBBBBBBBBB 100 39.63 +XXXXXXXXXX:7:1101:3371:100179#/1 GGTTCTTCGAATGTGTGGTAGGGTGGGGGGCATCCATATAGTCACTCCAGGTTTATGGAGGGTTCTTCTACTATTAGGACTTTTCGCTTCGAAGCGAAGG abaeeeeegggggfgegheghii[eghiicaacccccccdccddccccccc`bbcccccccc[aacccbccdcddbccbccccccccccccccccc^aaa 100 65.77 +XXXXXXXXXX:7:1101:2371:1957#/1 NCCGTCGGAAATGGTGAAGGGAGACTCGAAGTACTCTGAGGCTTGGAGGAGGGGAAAATAGAGACCCAGTAAAATTGTAATAAGCAGTGCTTGAATTATT BP\cccccc^ggghbfdcghf^ecfghhhhhaffbgfd[egfcd_GWbfZcagU^N^`__ZZ_b^aaaa]]`bbc_RT]]b_b_SYXGJS_b]^bb`]_a 100 49.81 +XXXXXXXXXX:7:1101:3311:100186#/1 GAAGTGTCTGACTAAAGTCCGTTGCCCATTTTTTAATCAGGTTGTTCATTTTCTTACTGTTCAGTTTTAAGAGTTCTTTGTGTATTTTGTATGACAGTCC bbbeceeegggggiiihiiiihiiiiiiihiiiihhhhhhiefhhiihiiiiiiiihiiiiihhhghihhhffdddggggeedeeeedbdbddccccbcb 100 69.12 +XXXXXXXXXX:7:1101:2394:1961#/1 NCTGGAAGTGATGGTCCGGAACCGTTCCTGGCCCGCTGTGTCCCAGATCTGAAGTTTGATAGTTTTGCCATCCAGCTCGATGGTTCGGATCTTGAAGTCC BP\cceeeggeggifhiiiiihifhhihhffhiifgfffgfghihhhfhiihiigggefgeeceeeebbddbdcc`bcccccaY``aaacabcb[`bbcc 100 52.65 +XXXXXXXXXX:7:1101:3438:100192#/1 GGCGGATGAAGCAGATAGTGAGGAAAGTTGAGCCAATAATGACGTGAAGTCCGTGGAAGCCTGTGGCTACAAAAAATGTTGAGCCGTAGATGCCGTCGGA bbbeeeeegggggihiiifghhiihiighihiiiiihiiihihihiiiiiiiihiiifhhiefgggeeeedcdcc_acabbcccccacccbbbbcXa]a] 100 66.94 +XXXXXXXXXX:7:1101:2333:1962#/1 NGGTTGCTTGCGTGAGGAAATACTTGATGGCAGCTTCTGTGGAACGAGGGTTTATTTTTTTGGTTAGAACTGGAATAAAAGCTAGCATGTTTATTTCTAG BPYacccca]ee_c\egafbgfgghf]fhfdbg__a^efghaffcf`eghefhdghhdg^`[]a`bab]`bYb]b__bbbbabb__bcbbddbbbdbd`] 100 52.52 +XXXXXXXXXX:7:1101:3479:100209#/1 TAGAGGTGGGAACAAGAAAACATTGATATAGTTTTGAATTAACATTAAAGAGGCAAAATTTCACATTTAAAATGGCCCAAAATATCCCACTATTCATTTT abbeeeacgggggiiihiiiiiiiiiihiiiihihiiihiiiiiiiiiiiiiiighifgiiihiiiiiihiggggdgeeedcddeddccccccddddded 100 69.20 +XXXXXXXXXX:7:1101:2459:1990#/1 CGGGGAGGTAGCGATGAGAGTAATAGATAGGGCTCAGGCGTTTGTGTATGATATGTTTGCGGTTTCGATGATGTGGTCTTTGGAGTAGAAACCTGTGAGG a_beeeee`cggghiiiihhdfgfdhhiihiiihhiiihiffhifgefhifhhihfhggegeabcdaaddcdc`_`bcccccc[bbccccccccabbcb^ 100 66.80 +XXXXXXXXXX:7:1101:3417:100210#/1 CCTGGATATTGTTTTTTCTAATAATGCCTATCACAAAAATGTATAAATGCCTTCGTATCATACTTTTTTAAAAACAACTGTCAATAATGAAAACTGAATG aa_ceeeegggggiiiiiiiihiiigiiiiiibgghiiiiheffhiiihehiihighhhiihiiiiiiggggggeaabddbdbcddccccccccbc_bcc 100 68.25 +XXXXXXXXXX:7:1101:2372:1994#/1 GCGGTGATGGTAGCCTTTCTGCCCAGCGCGTGCCACAGAGAAGGCTACACGAGCAGGATGCCATGCCCCAATACAGGCCACCTTGCGCAGGCCTCGGTGG ___ecccceceeggggfhghiiifggfgfigghf_fdgdgiffffffhhhhaggeeecdbbZ`]bbcbcc_cccbccac^^^`bcccac[_X[_^acW^a 100 64.27 +XXXXXXXXXX:7:1101:3452:100214#/1 GAAAAGGATGGAGGCAAATTTTAGCAGCAACCTCCACCTGACTGCTTGACTGGAGTGCTCTGAATTTTAACATTGGACTTCTCAGCACCACCAGCTGTTG a__eeecagegffhhihiiiiifffhiiihhhfhihhihbggegbbffghihihibgg`gfhhhecffgggdgegedeeedcdcddccccccccccc]Y] 100 66.98 +XXXXXXXXXX:7:1101:2677:1830#/1 NTCATTTTGTTTAAAAACTCGCTTGTCCCAGATTTGGCTATGGTACCAAGGTTTTTAACCAACTCTTCTCTGGTCATTCCTACACCGGTGTCTGTGACAT BP\cceeeggggehhiiiiiiiiihfhhiihiiifffiiifhfeeghgiiicfhihhhhhihhffhhhhggedggceeeeeeccddcccccbbcbccccc 100 52.71 +XXXXXXXXXX:7:1101:3354:100219#/1 TCTGGCTCGTCTGAGGGCACGGATGTCCTTCCCAATGCCTCGTGTGAGTTGGTAAGTATTGCCAGCTAAGGAAAGAAAGTTGGATATGTTCTCACCCAAA aabeeeeefegeehhiiihiiiihfgfhggdgighhfdhdfdfg\e_eeghhhhhideggeggeeeeedcbZ`b^^`bb_`bcbbbbcZ_`cbdbcc[`_ 100 65.63 +XXXXXXXXXX:7:1101:2603:1846#/1 NTCCACTCCTAGGAGGTCCCCAGAGTCCTCAGATTCAGAGACAGAAAGGATGGGGGACGGGCACAGTGGCTCACGCCTGTACTCCCAGCACTTTGGGAGG BP\cceeeeggfghfhcegihihffeffdhihhiihhhhhghfhhfhhib`fgffh`geecccccbbc^_accb``_^``bb`]^`^`acc`bbY^]`aa 100 52.59 +XXXXXXXXXX:7:1101:3600:100019#/1 GGAGAGGTTAAAGGAGCCACTTATTAGTAATGTTGATAGTAGAATGATGGCTAGGGTGACTTCATATGAGATTGTTTGGGCTACTGCTCGCAGTGCGCCG aa_eeeeegggggiiihiiiiiiiihighihihhiiihifhhiiiiiiiiiiihiigiiiiiihhhhighfghggcggddgeeeedddcccccccccccc 100 69.05 +XXXXXXXXXX:7:1101:2535:1848#/1 NGTAGCGTACCAGCCCAGGATGAGCTCATTTGGAGAAACTTTTTTATGCAGTTCATACATATTCTTAGCAAATTCCATGTCAACAGCCACTTCATCTTCT BSYcceeegggggiiiiihghiifhhhiiiiiifhhihghihiiihihfhiiihiiihiiihhigggfgggeeeeeeecdddddcaccccccccdcbcdc 100 52.79
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/r1_seqkit_fx2tab2.tabular Fri Jun 10 06:10:23 2022 +0000 @@ -0,0 +1,101 @@ +#id seq qual length avg.qual +XXXXXXXXXX:7:1101:1582:1835#/1 NGCTTCAGTTGGTACAGGTGTGTCACCCGTAGTTTGTAGGCTTGGAAGCCTGTAAGCTGTAGGGAACGTGGCAAGAACCAGATCTGGCCTTCCATATCTC BP\ccccc`eeccdeef_`fgabeghehdfhb^caefhee^efcehaeehhWabehddeM_eehhdeebd\^^a_]`___aaa`c]`Y^aaaaa_]]]aY 100 52.19 +XXXXXXXXXX:7:1101:1610:1859#/1 NTTATGTGTTACGTGAATCTCCTGAAGACAGCAGAAACTAGGTTGGTGAATTCTTATCCATTCTGCCATTCTGTATCTTTTAAGTGGAGCATTTAGGCTA BP\cceceggggggiihiiiiiiiiiiihihhiiihihiiiihiiihhiiihiiihhiiiihhiiiifhiiihgigggggggeeeeeecccddcdccbcc 100 52.73 +XXXXXXXXXX:7:1101:1743:1871#/1 NCGGGATGGAGGCAGGGCAAGCGGATGCAGGCCCAGCCCTGCCCTTACAGGTAGCTGGAAGTGACGAAGGAGTGGGGCCCCCGGACTCTGCTCAGGTGGA BP\ccceeggggghiiidfffhghfhhfffhhhidgggffbbghhheggeg`ccbccb_bdbdbbccaV^aaX^a`]X[_cca]]X]bc]Y]b`X`Y`_` 100 52.42 +XXXXXXXXXX:7:1101:1536:1878#/1 NGGCCTTTTTGGACAGGTGGTGTGTGGTGGCCTTGGTATGTGCTTTCTCGTGTTACATCGCGCCATCATTGGTATATGGTTAGTGTGTTGGTTAGTAGGC BP\cccedggggbffhieghdbgdeghffffhih[aaefdhefhhihhiiidffdffhgbbagaeebbdcddb_bcccbR]YXT]_`ab``[^b^]bcbb 100 52.41 +XXXXXXXXXX:7:1101:2990:100153#/1 GCCAGGAGTGACTGTGTTGTTCTTTGCTTTATATACATAAGCGCATCTCTTGCCCAAATAGAATTCTGTTTCATCTCGGGCGTAAACACCTTCAATTTTA ___cccecaegceg[dgfhgfbgdfgfhhfhhdfedhffhhhhhh[fhbbd_S^aaegfg`dfhgghhbcgddd``geadcc_a_a_][_c]bccccbba 100 64.07 +XXXXXXXXXX:7:1101:1624:1906#/1 NCCAGGTTTAAGGCATTACATCGAGGGGCTGAACTTGAGACAGTCCCTGAAAATGTAAGAACAGGAGCTAGAACTTGCTTATGGGAAGGGCCTAGGGAGA BP\cceeegggggiiiiihiiiiiiiiiiiiiiiiihiihihhiiiiiiifhiihhihdbfhggggeeeeddbddbdccccccbcbccccccccccccc_ 100 52.70 +XXXXXXXXXX:7:1101:1666:1926#/1 NGCAAACTGGTCAGGACCTGGAGCTGCTTGCCTGGCCTTGGCACAAGAGGCCGACATATGCCCAGCCCCTAGGAGTGACAGATCCGTCCTCTGTGGGAAC BPaceeeegggggiiiiiiiiiiiiiiiiiiihihiiiiiihiiiiiiiiiiiiiiiggggggeeeeccccccbbcaccccccccccccccccccccabc 100 52.71 +XXXXXXXXXX:7:1101:2921:100163#/1 CTCCGGATGCCAAAGGATCTGACAGTGTATCTAGCTTTGGAAAACACAGGGGTCTGCCCTGTGAGCTGCTCCAACACCTTGGCTGCTCGCGTCAGTCTGT bbaeeeeegggggiiihihhhiiiigheegfhhgdhihhhhhihhfhhhhbgggghihhdggihggeggeeeeeeacccccbccccccccc_a__bbccc 100 67.82 +XXXXXXXXXX:7:1101:1513:1929#/1 NCACTATCAAGGTTAGTGTCTATATCTACTTGAGANAAAATACTGAAAATTCAGTGTATGAAGCCAAATCCTGATTTAACAAGTTATTGGTAGTATAAGT BPY`acccceeeY`bfefeggehe`gddgceh_eWBPYbedeehfdefhhh]fhehbffhh_eeehded`cfhhdheeeeeeeddccccac`b]]ZZa`T 100 49.64 +XXXXXXXXXX:7:1101:2759:100170#/1 CTGCTCATCACATTCCAGATGCCATCACAGGCAATGACCATGAATTCATGGTCGTCAGTGAGAGTCAGCACCTTGATGTCAGGAAGGGCTGAAATCATCT bbbeeeeegggggiiiiihiiiiiiiiiiiiiiiiiiiiiiiihiiihiiihiiiiiighhiiihiiiiiiigggggggeeeeeddcccccabccccccc 100 69.31 +XXXXXXXXXX:7:1101:1708:1937#/1 NGGTCGGTGTCCTTTTTTGGCGTCAAACTTATAGAGGGAAGCAGGGAATCTACCCAACCAAAGCCAATGAACAGTCAGTCAGTCCATCTAGTTAAGGTTC BP\accccecggghiihihifhihffhihf`fegfghf`\ebcefgiigdgeggdeeec`abcbcbbbcbccccc`bccbccb```bbcccccc`bb`bc 100 52.64 +XXXXXXXXXX:7:1101:2981:100211#/1 GGGGTCAGTGGGAAAGGTGTCAATGGGACTTCGGTCTCCATGCCGAAACACCAAAGTCACAAACTTCAACTCCTTGGCTAGTACACTTCGGTCTAGCCAG __beccceeeeggdhffghhhhiihhcghffhiefaghhhfffffffhhiigifhf\cdbdddeeeebcbbd_bb_^bcbaccbcbbbbaa^a_]b`b`a 100 66.08 +XXXXXXXXXX:7:1101:1688:1946#/1 NCCGCCTTCTTTCTCTACTTCGTGACCTTTTGTGGTCCCGAGACCTGCTGCATCTTCTGTCTGATGTTCGGCTTGAGTGTCTACTTCCTGAACGACTTCT BPYac``c[cccchdd]caddeddehhhhhd`cdb^b]cc[^c_^_acccddchdhdcd\bdd]bd^bbbU\^BBBBBBBBBBBBBBBBBBBBBBBBBBB 100 38.51 +XXXXXXXXXX:7:1101:2767:100225#/1 CGCCTGGTGTCAATGCGCACATCTGGAGTTCCCATCTCCTTCATGGCAAATTTCCGAATCTCTTTGAGTGCCCGAGGTGCACGCTTCTTGAAGCCCACTC ___eceeefgggeiihihhhhihihihfgac_cefffihehiiiihhhhhfghhfhiehiihihf_ddbbgddcb`c]`bbbccccc`bb`b_bcccccc 100 66.85 +XXXXXXXXXX:7:1101:1536:1959#/1 NGATACTGTTGGGTCAGGACCAGGGCCTGCACCACCGAAATTCTGTCTTTTACAAACTGGACCAAGCGGTTGAGAATACAGGGGCCGAAGAGTAGGATTA BSYaa^ccggcg`d^dba_\ffh\[Zbf_aZ[c_cfY_ee`aeaf[\eefeghhffg_]\\`db^cdb`[ZWZR[`b]bbbbWELT]__TX^SYY]]]bb 100 50.36 +XXXXXXXXXX:7:1101:2797:100234#/1 GTGCTGGTCTCTCTCCACGATCTCTTACGTGCTAGCAGGGTGAACAAGGACGCCTCGCTGTCATTGCACTCACGCCACAGGAGCCTGGCACATTTCCGCT JY^accWceeeaWdcZ_feW`R^eg^a[HPP^_ef[d[Z_eUQX^ae_^Q\`ZSbUadVLKV_Z_Z_`]Z_baa^\WW^TGXTXWR^GW[][``BBBBBB 100 44.12 +XXXXXXXXXX:7:1101:1552:1976#/1 NTGATGAGTTTGCTAATACAATGCCAGTCAGGCCACCTACGGTGAAAAGAAAGATGAATCCTAGGGCTCAGAGCACTGCAGCAGATCATTTCATATTGCT BP\cceeeegegfiiiihiiiiiiiiiiihhiefgfhhhihifghdhhiiighfhhifihhiiiaghhigggggcdebdbddddcccccdedddddddcc 100 52.72 +XXXXXXXXXX:7:1101:1748:1978#/1 NCGCTGAGCAGGGCAGGCAGTGGTGGAGGCAGCCCCCCAAGACCCCATCTCTATCACAGTGAGAGCATCTGAACTCCAAAGCCACGACAGAACAGCTCCC BP\ceeeeggggghiiiihihiiegihhihiifiihhi`gfffggggeeeeeeddcdbd_bccbcccbccc`bbbbccccbccccccccaccbbb_`bcc 100 52.66 +XXXXXXXXXX:7:1101:2779:100239#/1 TGGGTGCTCTCCTTGGTCAGCTGCTGCTCGAACAGAGCACTTCTAGATACAGAGTTGGACCCAGACACAGCAGGGGTGGCCATGGCGGCTGGGAAAGACC ___^cc^ceggeeghddghafhXbegefbaYYbfhfgf_ae[egegc_cccaacaefg[bZbZ\\ZbZZ^gZZ^[]WW\_^bBBBBBBBBBBBBBBBBBB 100 40.42 +XXXXXXXXXX:7:1101:1593:1980#/1 NGGCACCAGAAGCATGAGCATAGGCAGAGCCATTTCCACCATGCTGACGGCAGGCAAGAAGAAGACAAAAATGAACCCCTTTCATGCTTGTTTAGAAATG BP\ceeeeggggfiiiiiiiihiihiiihfihhiiiifhhhiiiiihhiiiiiihhigggggeeeeeddcccccccccccccccdcccdccccccccccc 100 52.71 +XXXXXXXXXX:7:1101:2946:100242#/1 TGGGGATGTTGTTAAAAAAAATACAGGCTCCCCCACAACTGGGGTGCCTGGGGGGAACTTGGTCTGCTTCAGCCCAAGAGGAATCAAAAGATCAAAAGCA ___eecccgggggihhfhihfdfhhiihhiifiiddbcd`\eghL\^a^bdccccTTWRY`bRYYXSYb``bbca`acW^a`^bbccccc_bb`bbcbac 100 58.29 +XXXXXXXXXX:7:1101:1987:1781#/1 NTTAGCAGGTCCTCGTAGCGTTTCCGCAACTCTCGGAATCTGGAGTCTTCGGAGTGCAACTCTGAGGGTCCCGGGAAACTTGCGCGGCTCGCCTCGGCCA BP\cceeegggggiiiiihiiiihiiiiiiiiiiiiiiifiiiihhghiiihiggfgeeeeeebddcccccccccaccccccccccccccaacccXa^a[ 100 52.65 +XXXXXXXXXX:7:1101:3046:100006#/1 CTGTACAGCATGGATGTGGACAGGGAAAACAGCAAGTGCACACAGTGGGACAGGCCACAGGACAGGCTGGGAGTGAATAAAGAGTTCACACTGCTTCCCT bbbeeeeegggggiiihiiiiiiiiiiiiiiihiiieghiihhihffghiiiiiihiihhiiiggggeeec\bZ^bcbccc_]^]_b`bc`bbccbcccc 100 66.70 +XXXXXXXXXX:7:1101:1900:1788#/1 NCCCGCTGCCCACGCCACTGTCCCGGCGCTGCCCTGGGCGCGGAATGCGCGTTAGGCAATCCCCGGTCTGCGCGTCCCACACGCAGATGTGGCCTGCCAG BP\ceeeefggfgiiihiihiiiifhhiiihhiihiiihgeecccccccccaccccccccccccac_accca]aa[aacccccccccccdbccaccccb_ 100 52.64 +XXXXXXXXXX:7:1101:3214:100027#/1 GGCTGAGGTACTGAATGAGGAACAGCAGGAGGTAACGTCTGAGGGTGCGGCTGGGGAACGGTTTGCACAGGGGTGGCTGTAGGAGCAGGAGGTGAAACAC ___eecee^cegghhhhhhfhhfgfhhfhhhh^eggfcefhffhbU^eghhehgggTXZ\aTTX\RRYTT[GWROTWW[GG]YGXW[baOQEQW]b`]XW 100 49.49 +XXXXXXXXXX:7:1101:1848:1879#/1 NCTTCTTCTTCCTCATGTTTTTTGATCAACTGCTCCACTTCTTTTTGTTTGAGAACTCTGATTACTGTCTTTCCATTCTCTCTTGTTAGTGTTGCAATTT BP\cceeegggggiiiiiiiiiihhihhiehfhihhhiiihiiiihfghihdggfhhiiihihh_ddcggeefdb`beeeeeddcddbdbbcc`bbcccc 100 52.70 +XXXXXXXXXX:7:1101:3237:100032#/1 CTCCCGAAGTGCTGGAATTACAGGCATGAGCCCCCACGCCTGGCCTGATCAAGCTGTTTCTCTTGCTTGGAATGCCCTTTGCTGCCTGGCCAGCTCCTAG _b_eeeeegegggiiiiiihiiffegfhhihiiiiiighiffhiiiiiiiiihhiibdeggggeeeeeebcbdbcbcccccccccccccccccccccccb 100 67.87 +XXXXXXXXXX:7:1101:3027:100049#/1 CACAGAAACACGGGTTCATGATGGTGAATGATATAGCGGGTGACAGATGGCTACAAACCAGTCATAGGCCATCACACTCAGGAGCATGTCTCTCTTCCAT bbbeeeeegggggidghiihhihidgfhiihihfiiiiii_eghiiiihihhiiiggggeebdddcddddcccccbbcccccccccccbbdccccccbbc 100 67.75 +XXXXXXXXXX:7:1101:1756:1891#/1 NCGATGTTGAAGCCTGAGACTAGTTCGGACTCCCCTTCGGCAAGGTCGAAGGGGGTTCGGTTGGTCTCTGCTAGTGTGGAGATAAATCATATTATGGCCA BP\cceeegggggiiiiihiiiihiiiiiiiiiiiiiiiiiiiiigiiegghige]_ccaaccc_bbbbcccccc_bbcbcbccccdccccddddccccc 100 52.67 +XXXXXXXXXX:7:1101:3238:100065#/1 GGATGATAGAGTGTGGGAGAGAAAAGGGAGAGAGACTTTGGAGCTGAGTGGGGACCGCCGGGGAGGGGAGGCCTCAGAGTGGCTCAGCTGTAGATGCGGC bb_eeeeegggggdfgifhghhiiiiiihghhfgdhihiiifhhiiiigghiiiiigeecccc_accc_ccccbcbbbc_bccccbcbcc_bccdcbcca 100 67.15 +XXXXXXXXXX:7:1101:1915:1901#/1 NTAGAGGCTTACTAGAAGTGTGAAAACGTAGGCTTGGATTAAGGCGACAGCGATTTCTAGGATAGTCAGTAGAATTAGAATTGTGAAGATGATAAGTATA BPY`cceegfggfihgghhhegfghfgihdfhhhhhhfeffdgighiiiiifhiifdgdgecddc_cbbbdbbdcb__`bcc__`_`bbcbbccccac`c 100 52.64 +XXXXXXXXXX:7:1101:3198:100082#/1 CCCCCGCTCAACTACAAATGTCAATATTGAATGAAGCATTAAAAGACAAACATAAAGTAACTTCAGCTTTTATCTAGCAATGCAGAATGAATACTAAAAT b__eeeeeggggfiiiiiiieghhffhfhiiiiiifhiiiiiihgfhhiiiiiiiiigghhihiiiifegfggeeeeededcdcddcccccddcccbbcd 100 68.63 +XXXXXXXXXX:7:1101:1964:1931#/1 NGGGGTCTTAGCTTTGGTTCTCCTTGCAAAGTTATTTCTAGTTAATTCATTATGCAGAAGGTATAGGGGTTAGTCCTTGCTATATTATGCTTGGTTATAA BP\`Z`cdggggfhhhhghfghhhhhhhhhhffhhhhhhhhdggfhhhhhhhhghhhfgfceghfhhhhdgdgdggbggfgfebdeeeedcddbbccccc 100 52.70 +XXXXXXXXXX:7:1101:3088:100091#/1 TCCCACGTTTCCCACTTTAAAACAGAAGACTGCAGAGGATGTGGTAGAGACTGCAACTGCCAATGTAGCATCCGTTTTCCCAATTTCTCAGTAACAGGAT bbbeeeecgggggiiiiiiiiihffhihiiiiifgdghbghggighhhihhiiihiiiihiiiiihiiiiggggebcecdcbddcccdcdbccccccc`b 100 68.56 +XXXXXXXXXX:7:1101:1840:1948#/1 NCAGGGTGAGGGGCTTCGTCTGCTGGATGTACTTGTTGCCGAACTGAGTGACCCCCCGGGGCCAGCCGGTCTGCGAGCGATTGGGCGGTACCACAGACAT BP\ceecegggggiiihhhhiiiiiighhfhhiifhiiiiiiiiiiiihiiifhhgccccccacccccccccccccccccccccccccacaccccccccc 100 52.69 +XXXXXXXXXX:7:1101:3105:100094#/1 GGGAACTCCAGGGATGAACATGAGAAGGCCCTTCTGCTGTAGTACTAACTGCAATCAGTGTCTCGGGCTTCATGGGCTTGGCTTTACGTAACAAGCCAGG ^[[^ccccc^eccgf\R`bggfffbffhhfhbXc]Xcfcg]ece_cegfadafhgfheNWW\_]eZFZW`^HVZ]UZ`ZKZGT\b_b[KYXGY_XQWbW^ 100 50.33 +XXXXXXXXXX:7:1101:1958:1952#/1 NTCATGGTAGGGGTAAAAGGAGGGCAATTTCTAGATCAAATAATAAGAAGGTAATAGCTACTAAGAAGAATTTTATGGAGAAAGGGACGCGGGCGGGGGA BPY`cccc^aaecf`e[dbeeedaaeVefhhhdeXc_fehdb_eedhhhed\acbefdhdddee`bbeeeeeacca]_aa__]]_aaaaaa]aaa_]XaT 100 52.40 +XXXXXXXXXX:7:1101:3190:100106#/1 CACGGTTAGGGTACCGCGGCCGTTAAACATGTGTCACTGGGCAGGCGGTGCCTCTAATACTGGTGATGCTAGAGGTGATGTTTTTGGTAAACAGGCGGGG bbbeeceegggggiihhhihifhiihiiiiifgeghihhhhhfhgfecR]_bcccccdcccdc]b]bbccccbccc^`bc_bcccaccccccccacccac 100 64.50 +XXXXXXXXXX:7:1101:1993:1999#/1 GTTTAATTGTGGGAGGCTGAGAGCAGCATTTTTACCAAATTGGTGTAAAAATGAAACGGGGTCCAGACGTGTACACTGACAGTTACAGTAAATTTCAAAA bbbeeeeegecgfdhhiiihihiiiiiiiihiighihhfhihhcf_ehhhiiiiihihhggadgggeeeccddddbccccccccccccedccecccdccc 100 68.03 +XXXXXXXXXX:7:1101:3117:100110#/1 GTTTTATGTTGTTAATGTGGTGGGTGAGTGAGCCCCATTGTGTTGTGGTAAATATGTAGAGGGAGTATAGGGCTGTGACTAGTATGTTGAGTCCTGTAAG ___eeeeeggggfiiiigihghiibfcgefdgfhii`ghiefefhghhaedghdhibfghfhiff^`gbggeeeedbcbdbbbccdedcccbbbcbbbcb 100 67.19 +XXXXXXXXXX:7:1101:2159:1798#/1 NCTTCCAGTCCCCTCCTCGATTTCACACCGGGACACACACCTGGGTCAGGTGAGCATGACACACAGCCCCCAGGATGGGCACCCCACCCTGCAGAGCAGA BP\ceeeegggggiiiiiiiiiihiiiiiiiihhiiiifiiiiiifgiiiggggggeeeeeddccccaccccccccccccccccccccccccccbccccb 100 52.69 +XXXXXXXXXX:7:1101:3147:100111#/1 AGAAATTATTTTTAAATGTCAGAGAGAATGCTATTATTGGCCAACTTGAAATTTCCTAATTAGCTATCTGAGCTATTTAAGTCTGTACTTCAGTGTCTTC bbbeeeeegggggiiiiifiiiiiiiiiiiiihhiiiihhiiiiiiiiiiiiiiiihigiiiiiiiiiiiihihfcghhihbggggfgeeeeea_bdddd 100 69.60 +XXXXXXXXXX:7:1101:2152:1838#/1 NGGGAAACGCCATATCGGGGGCACCGATTATTAGGGGAACTAGTCAGTTGCCAAAGCCTCCGATTATGATGGGTATTACTATGAAGAAGATTATTACAAA BP\cceeegggggiiiiiiiiiiiiighigggggeeecccdccccdcceccccccccccccccccccccdccc^bbccccddccdcccccbcccddcccc 100 52.66 +XXXXXXXXXX:7:1101:3065:100152#/1 CACATAGTCAGCCTGCTTGATCACGTCGATCTTGTCCAGGAAGTACTGGGCACAGTCAATCAGCTGGTACTCGTTGGAGCGTTCGTAGCAGGCACGCACT ___ecacccggcgfhhhhhhh_gfdeghhhhhhhfhhdcdhXfX^f_eghhfUegbfh_be__dhbb`g`bga_aZ^]bdaX\^a^`aa`_[_aW]][]a 100 62.63 +XXXXXXXXXX:7:1101:2180:1843#/1 NTAGGGTGTAGCCTGAGAATAGGGGAAATCAGTGAATGAAGCCTCCTATGATGGCAAATACAGCTCCTATTGATAGGACATAGTGGAAGTGGGCTACAAC BP\cceecgggfehiehihhihhhhdffhihifdhhifghhhiiihihihiiihiiiiihiiihfhiffefgfdegeecac`d]bddbb_bcaacccccc 100 52.70 +XXXXXXXXXX:7:1101:3154:100159#/1 GTTTCATCATGTTGGCCAGGATGGTCTCAATCTCTTGACCTCGTGATCAGCCTGCCTTGGCCTCCCAAAGTACTGGGATTATAGGCGTGAGCCACCGCGT bbbeeeeegggggiiiiiiigiiiggiiiiififhiggiiiiifhiiiighiiiiiiiiihiiiiihhffffgggdgeeeeeedbdcaccaccccccaca 100 68.83 +XXXXXXXXXX:7:1101:2125:1861#/1 NGTGGGAGCAGCTGGCTGTGTGGCACCTGCTGTGGTGGCTGCTGGCGTGGCCGGAGTGAAGGGCAACGTGGCAGGTGCCGGGGGCTGGGCTGGATTCCCT BS\cceeegggggiiiiighghihfhifhhihfgfeghfbffgffaf_efffegX]Z]]_bccPZ^`aX`a]`acBBBBBBBBBBBBBBBBBBBBBBBBB 100 38.83 +XXXXXXXXXX:7:1101:3198:100173#/1 TGGCAGGTGGTGGTGGTGCTCAAAGAGCAGGTCCAGGTGGAAAGTGAGCACCGACAATGGCTGCAGCAGGAGCAGCAGCTCTGTGGACAGGGAGGGACAA ___cccccece`gbZaffgadfgg_ddffghbafffdaffffhhedghhhhhffhgggadc_`bbbba^a\\XXT^[^[Y_]`Y]b]`b__aaBBBBBBB 100 44.48 +XXXXXXXXXX:7:1101:2076:1911#/1 NGGAGTAGGAGTACCCAGTCTCTAGTGCCAGCTGCAAGCGTTAATTTGGTTAAAGTCTCCTTGATTGTCCTGTTCATTCTTTCTACCTGACCTGAACTCT BS\ceeeegggggiihiiiiihiihgiiihiiiiiiiihhhihghiiiihhiiiihihiihiihiiihhiiihfggggggeeeeabcdbdddcccccccc 100 52.81 +XXXXXXXXXX:7:1101:3166:100190#/1 CCGGCCGAAATTGTTTTTTAAATTTACGTGTCACCCTTCATGCAGGGCCCTGCTACTCTTCGTAGCCCTCCTCCATGTTTAGTGTGGGTGCGGCTGCAGC bbbeeeeeffgegdgigfiihhiiihiifggghffc`gfffgheifffagdbddgdbeced`bdcaccbaccccb`b`bccbY`bccaX]_accaa^R^b 100 63.87 +XXXXXXXXXX:7:1101:2196:1920#/1 NTGGGGTGGGTATGGCGGGGCCAGCCTTGCTCAGTGCTGTGGTGATGGCCACAGTAGCCTCGGGGGGCACAGCCACGGCTGGGGGAGCAGGCCCTTCCTT BP\ccc`ceeceehhhhhhh[aZZZ^U]LZ]]ZTZ_]]TTY_KTT]]`[a]_[[R]]R__aR[aaaaTEW^XP[GRLTTOOXaaBBBBBBBBBBBBBBBB 100 40.29 +XXXXXXXXXX:7:1101:3225:100207#/1 GCTATTTTCTGCTAGGGGGTGGAAGCGGATGAGTAAGAAGATTCCTGCTACAACTATAGTGCTTGAGTGGAGTAGGGCTGAGACTGGGGTGGGGCCTTCT bbbeeeeeggfgghiiihiYcghiihiihiififhfghfhhhhhhggggggeeeeeddddddcdcac]`ababbbbacccccccbcccc[aaacaccccc 100 66.30 +XXXXXXXXXX:7:1101:2115:1927#/1 NCTCCTTCTTGGCCTGGAGGCGCTCTTCACGGCGCTTGCGTGCTTCCTTGGTCTTAGACCTGCGGGCCTCAGCCTGGTCAGCCAGGAGCTTCTTGCGGGC BS\ccecegeeggiiiffggfhhiiiihhdhihfdghddfggeeeeecdddZ]bcdcbcccbccacaaccccccc^`YYYX`bW`_`a[`_b``_baccT 100 52.51 +XXXXXXXXXX:7:1101:3019:100219#/1 GGCGGGTTTTAGGGGCTCTTTGGTGAAGAGTTTTATGGCGTCAGCGAAGGGTTGTAGTAGCCCGTAGGGGCCTACAACGTTGGGGCCTTTGCGTAGTTGT _[^ccc_cggccgghdgcbgdf_WcXae[cbefffgbefh_d`dgddb]aaZTZW__cc`bb]ZZ_^aaZ]a^`b_Y__ab[`_XT_^bbbbaBBBBBBB 100 44.46 +XXXXXXXXXX:7:1101:2179:1937#/1 NAGGAATACATATAGTCATCATTGCCAGACTTAATATGAGAGGTGAAATGTTCGATCCAATTTTCCTTCCTGGATAAGTTTTTCTTTCCTATCCCTCTCA BP\cceeeggggghigfghiiiiiiihiiiiiiiiiiifhhhheghhiihhhihiiiiiiiiifhghhiiiiihhgffdgfgegeeeeeeeddcddcccb 100 52.73 +XXXXXXXXXX:7:1101:3202:100230#/1 CTGAAAAACAAAACAAAACTTAGCTAAGTTTTATGGCAAGTGGGGGCAGGGTCAGGTGAAGAGGTAGGAATTGCCTGAGGTTGGGGGTCCCTAAAACATA __aeeeeegggggihhhiiiiihihhhhfghhhihhhiihfffgdfhhihiddggeceeeeedcbbdcbbbcbcc]bbcc^abccccT[acbbcccb^`b 100 65.28 +XXXXXXXXXX:7:1101:2149:1945#/1 NGTGTGGTGGCCTTGGTATGTGCTTTCTCGTGTTACATCGCGCCATCATTGGTATATGGTTAGTGTGTTGGTTAGTAGGCCTAGTATGAGGAGCGTTATG BPYaccccegggghhhceghghhhhhhhhhfhhhhhhhhhhhhhffhghhhdffghghhbfghgfddgddcdddbebcbbbbbbbdcdbbbbaaaaaaab 100 52.67 +XXXXXXXXXX:7:1101:3211:100242#/1 CTGTTCCCCAATACCTGCAATTCACAGATATAATTTAGCATTTTTGTTTTTATATACTAATTTAGGCAAAACCCCCCCGCCACTGAATCGTCCAACAAAA ___`acacgggcghfdddc_gbggghddgfgfhhccdab[e[_Xbe^^cacdfgbaa[caSWWbfbf`efc_Zbg[XZaccaacccbcc[[WW^bbW[[a 100 61.18 +XXXXXXXXXX:7:1101:2169:1964#/1 NGGGAGCAAAAAGCGTGGGTACAGAAGCGAGAAGCGAGCTGAGTGGTTAGTTCAAATCAAGGCTTGGGGTATTTCCTGGTCATTTGGGGAACCTGAAACT BP\ccceeggggfehfhhh_fhifbgf[effhh`cfV__`gefccg`ddg`d`beeeadddbcccaaacNWY_bZ]`bb_YY__aGSTXGTW^R^R_b_` 100 50.12 +XXXXXXXXXX:7:1101:3168:100244#/1 CCCTTCCTGTGTCCATGTGATCTCATTGTTCAATTCCCACCTATGAGTGAGAATATGCGGTGTTTGGTTTTTTGTCCTTGCGATAGTTTACTGAGAATGA ___eeeeegegeghiihfbgghhhhfhaggiiiidbfhfhfbdegdfcffgdgfhhfhhfceefghicfffeeceddbddba^acc_bccccc[bbbcc_ 100 66.77 +XXXXXXXXXX:7:1101:3005:100246#/1 CAAAAGTTAAATGAAAAACCACAAGTTTAGTTTGACTCCGTCTCCTAGGGTGGATTTCATTCAGATATTTGTTCCATATTATAGGAGGGTGGATCCTAGC bbbeeeeegggggihhhihihiiiiiiiiihhihhhihiighiiiiihii^efaegfffhhfhhhiihhiidgggggdeeeeeeedccaZ\_accccccc 100 67.69 +XXXXXXXXXX:7:1101:2313:1789#/1 NTGCGCCTGGAAGCCGCCAGGCCGCAGCGCGGGGAACAGGGCGGCCACCCGGCTCAGCAGCCGCTCGCAGGGTCCCGGACCCGCCGCAGCCTCGGCTTCC BP\cceeegggggiihiiiiiihiihiiifhheccccccaccccccaccccccccccbcccccc_a]Zaaccbb`aaacaccaccccaaacccX]^^abc 100 52.58 +XXXXXXXXXX:7:1101:3253:100014#/1 GCGGATGATTCAGCCATAATTTACGTCTCGAGTGATGTGGGCGATTGATGAAAAGGCGGTTGAGGCGTCTGGTGAGTAGTGCATGGCTAGGAATAGTCCT bbbeceeegggggiiihihiiiiiiigfhiafbeccb_fdgaffhhhhhihihihhiieT]abcccccccccbb_WY]_bbbbccccccccbcbb`Y`bb 100 64.18 +XXXXXXXXXX:7:1101:2361:1794#/1 NGTTAGGCTCTCCACGACTGAGCTGGTCTCGGCACCAAGGGGGTCAGATTGCTTAAAGGCAAGGAGTTCAGGACCAGCCTGGGCAATGTGGGGAGACCTC BPYacceeegggfhhhhfhhfhhb_eaefdgghhfgfffgbfcfdeggfgeeeeeddbd\``aaabbbbdbbbb^[[^aba_aababb`_Y[aaaR[a`b 100 52.45 +XXXXXXXXXX:7:1101:2337:1794#/1 NGGGAAGGTGTCTGTCATCGTCCCCAGGAGGGAGGTGGCGTTAGCTGTTTGTCCTGTCATTAGGTTGGTAACTCTCCAAGTAACATTGAAGATCTGGTGA BP\c`ceceggggiihiiihiffffghi`ghhggi^ce`f_faghihbghigegdgR_ceeeeeU]`c^`_bcbbccbbb]`bcbcb`b`cbbbc_bbcc 100 52.47 +XXXXXXXXXX:7:1101:3284:100039#/1 TTGGGCGACAAGGACATGGAGCTGGAGGAGAGCGTGGAGGAGCGGCTGTTGAGTTCAGAGCCCGAGTCCACTGCCACCGAGGCCGGGCTGGCGGCCAGGG bbbeeeeegggggiihiiiiiiiiihhihiffhifhiihhifhigeec`db`bbccdcccbcccccacccccccc`bcc_acacccccccca_a[_BBBB 100 46.93 +XXXXXXXXXX:7:1101:2477:1795#/1 NGGCGGTGCACCGAGCTGGAGGCTGCCGCAGAGGAGGCGGCCGCCGCTGCGGCCATGTACGTGTAGGGGTAAGGGAACAGGCTTCCGAAAGGGGACATGG BP\ccccefeeggffghffhhihhhhhfifhhiidghfddcca_aWaWac_accabccdcaX`R^bba_X_bcccW`^acc_^[X`_[[[a_acacc]][ 100 52.30 +XXXXXXXXXX:7:1101:3310:100056#/1 GGGCGATTGATGAAAAGGCGGTTGAGGCGTCTGGTGAGTAGTGCATGGCTAGGAATAGTCCTGTGGTGATTTGGAGGATCAGGCAGGCGCCAAGGAGTGA bbbeeeeegffgfiihiihiighiiihiihhiii\edgfhiddggeggeeeedbddcdbdcccbbc^a`bcdccb`aabbaccb^acccacccca`cX`b 100 65.89 +XXXXXXXXXX:7:1101:2355:1821#/1 NTGTATTGAAGCCAAATTTCGTGCCAGGGTGGTGTTAGAATGTGAGCTATTTTCATTATATTTGGACTAGAATTTTTTCTTTTAGGAGAATGGTGGGGAA BS\cceeegggggiiiiiiihfgiihhii^egbffhfhfhhiehhhiihiiiiifgh`ffghiiighfhhgggfgggcdddddcbbbc`abcc_bcccac 100 52.78 +XXXXXXXXXX:7:1101:3420:100060#/1 GTGAAAAATGTTTAATTTACTTCAAAAATATCTTGCATAAGTATCTCATCCCTAAGCTCAAACAAACAAGAAAATTAAAGGCGGAAAGAATACATAAAAT ab_eeeeeggffgihhhihgfdhhhhihf]ghiihefhhig^fgbghbfdhhiihbffhhiihiiihhhhhcfhigfhdg`dFXZaccbcbcccddcbcc 100 56.56 +XXXXXXXXXX:7:1101:2418:1834#/1 NCGGTTTCAGTCGTCACTGCTGCCCCAGCCTTTCGCTGTCCTTCTTGCAAAAAGCTGCTCCCATCGGTATACCAGGTGTGGTCGGCGTCTGGGATGGGCT BPYcaceeggggfiiiiiiiihifiiihiiiiihdfhhfhfhiiiiihiihhhhhihihhiiegggeeecbeddcabbbcccccccc[aaccc`bccacc 100 52.67 +XXXXXXXXXX:7:1101:3267:100061#/1 CTGGGTAGTTCTGGCTCCGGGTCTCTCATGAGGTTGCAGTCACGATATTGGTTGGAGTTGCAGCCATCAGAAGGCTAGACTGCGTGTGAAAGACTCACTT _bbeeccegggggiiiiiihigghihifhiiiiffhiiihhihidghiihheghhhighgggfeeeeeeccddcccccccccccaacccdcccccccccc 100 68.01 +XXXXXXXXXX:7:1101:2378:1838#/1 NGCCCAACCAACAGCTTCAACAGCTATTGCAAGGTGCGCCAAAGGCATCTTGTCATCCCTCACACGAATCTCACTTCCTGTGAATTTGCAGGGAGGCAGA BP\ceeeegggggiiiihiiiiiiiiiiiiiiihfhiihiiiiiiiiiiiiiiiiiiiiiiggggeeebcddccccccccccddddcccdcccccccccc 100 52.71 +XXXXXXXXXX:7:1101:3416:100083#/1 GCCCAGTCCTGGGGGCAGGGGAGAAAGGGAGGGGTGAAATGAGAAGAGCCCCGTTCTTCCTCAGTGCTGGGGGGTTGGAGACACCACCCCCTTCTGCTGC ab_eeeeeggggfiiiiiiiighghhiiihhhiiL\bgeeeeeeddbddcccc^acacccccccbccccccccc]X_abcbcccaccacccccccccccc 100 61.17 +XXXXXXXXXX:7:1101:2481:1853#/1 NGCAGGTAGAGGCTTACTAGAAGTGTGAAAACGTAGGCTTGGATTAAGGCGACAGCGATTTCTAGGATAGTCAGTAGAATTAGAATTGTGAAGATGATAA BS\cceccgggggiiiiiiihiihhfhhihiiifhgiiiiiihiiiiiiiiiiiiiieggeeeeeebdddbddccbbccccccbbbccbcbcccccccdb 100 52.78 +XXXXXXXXXX:7:1101:3411:100111#/1 GTGGTGGCCTTGGTATGTGCTTTCTCGTGTTACATCGCGCCATCATTGGTATATGGTTAGTGTGTTGGTTAGTAGGCCTAGTATGAGGAGCGTTATGGAG babeeeeegggggghiihiiiiiiiihdfgiiiiiifiiiiiiihiiiiefhhiiiegfhggghgggfbddeceeedbdbb_bedcccccccccccccba 100 68.20 +XXXXXXXXXX:7:1101:3258:100128#/1 GCTGTGGCTTTTTTAAGGATAACTACCTTGGGGGCCTTTTCATTGTTTTCCAACTCCGGGATCTGGTCACGCAGGGCAAAAAAGCTCCGTTTTAGCTCGT _b_ecceeggggfihghifffhhfghidghhiihighiihighiighhiifdffggggeecccccbccbccccca[_aaccccccbbbcacbcbbcccc[ 100 66.63 +XXXXXXXXXX:7:1101:2252:1856#/1 NCAGAAATATGGCAGTGCCACAGAGAGTGTGTCCTCATTGGAGTACTTGCGCTTAAATTCATCCAACACAAAGGTACTCTTGGGCAGGTGAGCGAAGGGG BPYcaccce^ecc^e`begg_a_fXQ`[Q_eb_b`d]ee^[ac^e^e[cbe^edeedVWW_bb\_\dd`bdGV^aaab``UUG\\\]_T_]aaaZO[T]a 100 50.04 +XXXXXXXXXX:7:1101:3428:100129#/1 GTTGGAGTGTAGGATAAATCATGCTAAGGCGAGGATGAAACCGATATCGCCGATACGGTTGTATAGGATTGCTTGAATGGCTGCTGTGTTGGCATCTGCT ___eeac\cecggfhhihihfhhiiiihiiihfhffihhhhihhidg`fgfgeeccccZ^^\]bbbc_bb`bbccabcbb`ba[]`YbS]Q[RSRX`bb_ 100 60.12 +XXXXXXXXXX:7:1101:2394:1871#/1 NTTTATTTCTAGGCCTACTCAGGTAAAAAATCAGTGCGAGCTTAGCGCTGTGATGAGTGTGCCTGCAAAGATGGTAGAGTAGATGACGGGTTGGGCCAGG BP\ceeeegggggiiiiiiiiiigghiiiiiiiigihihhiiiiiiiihighiiiiihgeggggeeeeedddddb_bccbcccdcccccc_ccccccccc 100 52.70 +XXXXXXXXXX:7:1101:3387:100138#/1 TTGAGGGAGGGCAGACACAGGCTGAGGAGAGTGGATAGGAGCAGGAGAACCTGCCACACCCCACACCCCACTGAGAGCTCATCAGCTGCAGCCCTGGGCC bbbeeeeegggggiiiiiiiiPdgfhihhfgXafffhfhhhiiiicfgifhhhegegceeccccccccc_aacc``__cccccbbcccc`bcccc`accc 100 63.66 +XXXXXXXXXX:7:1101:2269:1904#/1 NGCACCTGTAATCCCAGCTACTCAGGAGGCTGAGGCAGGAGAATCACTTGAACCCGGGAGGTGGAGGTTGCAGTGAGCCAGGGTCAGGCCATTGCACTCC BPY`cceceaga^dadghhhgfhhbfghhhhhhhhfhgad`ec[ecedbbbecbdghgaaXKU^GZ\GTX`bRY]__bbT^aaRRWbb^a^]bbbbbb^` 100 49.98 +XXXXXXXXXX:7:1101:3444:100163#/1 TGTGGCTGTTCTATGAACCACCTTCTTCTTTCTGCGAGCAGTTCCTTTCCCACCAATGCGCACTTGTGCCTGCAGTTTGGCGAGGTTTTCCTGGGTCATG __\c\cc`^`c^be``bad`ddegaX`dgegfaIXW^^OW^O^ce]a^ae^egffgd_VZ_FT\ZV_]c`bb_bX]Z]R^baQREHW`WY_bYXGX^^BB 100 46.94 +XXXXXXXXXX:7:1101:2259:1943#/1 NTCTTCCTCATCTTCTTGTTCTTTTCTTTTTCCTTCTCCTTCTCCTTTTTCTTTTGTGACATGATTTCACTCTTGTTACCTACGCTGGAGAGCAGTGGTG BP\accc`b[QbKb`cd_``dedddcddYXecSbdI^P^XcSXbX^Y^^c^a^cc_HHXacIXWM\WWWa\HW_]NV\``_BBBBBBBBBBBBBBBBBBB 100 39.63 +XXXXXXXXXX:7:1101:3371:100179#/1 GGTTCTTCGAATGTGTGGTAGGGTGGGGGGCATCCATATAGTCACTCCAGGTTTATGGAGGGTTCTTCTACTATTAGGACTTTTCGCTTCGAAGCGAAGG abaeeeeegggggfgegheghii[eghiicaacccccccdccddccccccc`bbcccccccc[aacccbccdcddbccbccccccccccccccccc^aaa 100 65.77 +XXXXXXXXXX:7:1101:2371:1957#/1 NCCGTCGGAAATGGTGAAGGGAGACTCGAAGTACTCTGAGGCTTGGAGGAGGGGAAAATAGAGACCCAGTAAAATTGTAATAAGCAGTGCTTGAATTATT BP\cccccc^ggghbfdcghf^ecfghhhhhaffbgfd[egfcd_GWbfZcagU^N^`__ZZ_b^aaaa]]`bbc_RT]]b_b_SYXGJS_b]^bb`]_a 100 49.81 +XXXXXXXXXX:7:1101:3311:100186#/1 GAAGTGTCTGACTAAAGTCCGTTGCCCATTTTTTAATCAGGTTGTTCATTTTCTTACTGTTCAGTTTTAAGAGTTCTTTGTGTATTTTGTATGACAGTCC bbbeceeegggggiiihiiiihiiiiiiihiiiihhhhhhiefhhiihiiiiiiiihiiiiihhhghihhhffdddggggeedeeeedbdbddccccbcb 100 69.12 +XXXXXXXXXX:7:1101:2394:1961#/1 NCTGGAAGTGATGGTCCGGAACCGTTCCTGGCCCGCTGTGTCCCAGATCTGAAGTTTGATAGTTTTGCCATCCAGCTCGATGGTTCGGATCTTGAAGTCC BP\cceeeggeggifhiiiiihifhhihhffhiifgfffgfghihhhfhiihiigggefgeeceeeebbddbdcc`bcccccaY``aaacabcb[`bbcc 100 52.65 +XXXXXXXXXX:7:1101:3438:100192#/1 GGCGGATGAAGCAGATAGTGAGGAAAGTTGAGCCAATAATGACGTGAAGTCCGTGGAAGCCTGTGGCTACAAAAAATGTTGAGCCGTAGATGCCGTCGGA bbbeeeeegggggihiiifghhiihiighihiiiiihiiihihihiiiiiiiihiiifhhiefgggeeeedcdcc_acabbcccccacccbbbbcXa]a] 100 66.94 +XXXXXXXXXX:7:1101:2333:1962#/1 NGGTTGCTTGCGTGAGGAAATACTTGATGGCAGCTTCTGTGGAACGAGGGTTTATTTTTTTGGTTAGAACTGGAATAAAAGCTAGCATGTTTATTTCTAG BPYacccca]ee_c\egafbgfgghf]fhfdbg__a^efghaffcf`eghefhdghhdg^`[]a`bab]`bYb]b__bbbbabb__bcbbddbbbdbd`] 100 52.52 +XXXXXXXXXX:7:1101:3479:100209#/1 TAGAGGTGGGAACAAGAAAACATTGATATAGTTTTGAATTAACATTAAAGAGGCAAAATTTCACATTTAAAATGGCCCAAAATATCCCACTATTCATTTT abbeeeacgggggiiihiiiiiiiiiihiiiihihiiihiiiiiiiiiiiiiiighifgiiihiiiiiihiggggdgeeedcddeddccccccddddded 100 69.20 +XXXXXXXXXX:7:1101:2459:1990#/1 CGGGGAGGTAGCGATGAGAGTAATAGATAGGGCTCAGGCGTTTGTGTATGATATGTTTGCGGTTTCGATGATGTGGTCTTTGGAGTAGAAACCTGTGAGG a_beeeee`cggghiiiihhdfgfdhhiihiiihhiiihiffhifgefhifhhihfhggegeabcdaaddcdc`_`bcccccc[bbccccccccabbcb^ 100 66.80 +XXXXXXXXXX:7:1101:3417:100210#/1 CCTGGATATTGTTTTTTCTAATAATGCCTATCACAAAAATGTATAAATGCCTTCGTATCATACTTTTTTAAAAACAACTGTCAATAATGAAAACTGAATG aa_ceeeegggggiiiiiiiihiiigiiiiiibgghiiiiheffhiiihehiihighhhiihiiiiiiggggggeaabddbdbcddccccccccbc_bcc 100 68.25 +XXXXXXXXXX:7:1101:2372:1994#/1 GCGGTGATGGTAGCCTTTCTGCCCAGCGCGTGCCACAGAGAAGGCTACACGAGCAGGATGCCATGCCCCAATACAGGCCACCTTGCGCAGGCCTCGGTGG ___ecccceceeggggfhghiiifggfgfigghf_fdgdgiffffffhhhhaggeeecdbbZ`]bbcbcc_cccbccac^^^`bcccac[_X[_^acW^a 100 64.27 +XXXXXXXXXX:7:1101:3452:100214#/1 GAAAAGGATGGAGGCAAATTTTAGCAGCAACCTCCACCTGACTGCTTGACTGGAGTGCTCTGAATTTTAACATTGGACTTCTCAGCACCACCAGCTGTTG a__eeecagegffhhihiiiiifffhiiihhhfhihhihbggegbbffghihihibgg`gfhhhecffgggdgegedeeedcdcddccccccccccc]Y] 100 66.98 +XXXXXXXXXX:7:1101:2677:1830#/1 NTCATTTTGTTTAAAAACTCGCTTGTCCCAGATTTGGCTATGGTACCAAGGTTTTTAACCAACTCTTCTCTGGTCATTCCTACACCGGTGTCTGTGACAT BP\cceeeggggehhiiiiiiiiihfhhiihiiifffiiifhfeeghgiiicfhihhhhhihhffhhhhggedggceeeeeeccddcccccbbcbccccc 100 52.71 +XXXXXXXXXX:7:1101:3354:100219#/1 TCTGGCTCGTCTGAGGGCACGGATGTCCTTCCCAATGCCTCGTGTGAGTTGGTAAGTATTGCCAGCTAAGGAAAGAAAGTTGGATATGTTCTCACCCAAA aabeeeeefegeehhiiihiiiihfgfhggdgighhfdhdfdfg\e_eeghhhhhideggeggeeeeedcbZ`b^^`bb_`bcbbbbcZ_`cbdbcc[`_ 100 65.63 +XXXXXXXXXX:7:1101:2603:1846#/1 NTCCACTCCTAGGAGGTCCCCAGAGTCCTCAGATTCAGAGACAGAAAGGATGGGGGACGGGCACAGTGGCTCACGCCTGTACTCCCAGCACTTTGGGAGG BP\cceeeeggfghfhcegihihffeffdhihhiihhhhhghfhhfhhib`fgffh`geecccccbbc^_accb``_^``bb`]^`^`acc`bbY^]`aa 100 52.59 +XXXXXXXXXX:7:1101:3600:100019#/1 GGAGAGGTTAAAGGAGCCACTTATTAGTAATGTTGATAGTAGAATGATGGCTAGGGTGACTTCATATGAGATTGTTTGGGCTACTGCTCGCAGTGCGCCG aa_eeeeegggggiiihiiiiiiiihighihihhiiihifhhiiiiiiiiiiihiigiiiiiihhhhighfghggcggddgeeeedddcccccccccccc 100 69.05 +XXXXXXXXXX:7:1101:2535:1848#/1 NGTAGCGTACCAGCCCAGGATGAGCTCATTTGGAGAAACTTTTTTATGCAGTTCATACATATTCTTAGCAAATTCCATGTCAACAGCCACTTCATCTTCT BSYcceeegggggiiiiihghiifhhhiiiiiifhhihghihiiihihfhiiihiiihiiihhigggfgggeeeeeeecdddddcaccccccccdcbcdc 100 52.79
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/r1_seqkit_stats1.tabular Fri Jun 10 06:10:23 2022 +0000 @@ -0,0 +1,2 @@ +file format type num_seqs sum_len min_len avg_len max_len Q1 Q2 Q3 sum_gap N50 Q20(%) Q30(%) +13-1941-6_S4_L001_R1_600000_fastq FASTQ DNA 100 10000 100 100.0 100 50.0 100.0 50.0 0 100 100.00 100.00
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/r1_seqkit_stats2.tabular Fri Jun 10 06:10:23 2022 +0000 @@ -0,0 +1,2 @@ +file format type num_seqs sum_len min_len avg_len max_len Q1 Q2 Q3 sum_gap N50 Q20(%) Q30(%) +13-1941-6_S4_L001_R1_600000_fastq FASTQ DNA 100 10000 100 100.0 100 50.0 100.0 50.0 0 100 100.00 100.00
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/r2_seqkit_fx2tab2.tabular Fri Jun 10 06:10:23 2022 +0000 @@ -0,0 +1,101 @@ +#id seq qual length avg.qual +XXXXXXXXXX:7:1101:1582:1835#/2 GCAATGATGGGGCCGCTCCCGGGGCCACACCTGCTTCTGGATCCGCTGCTTCCAAGTTCCTTTGCCTCCCTCCTGGCATCACTGTCTGCGACTCAGGCCG ___ccc`cgegfedfhfhhhfdfec_eg`_fhfhiihiid`ggdadc`bd`b_bcccb__bbbbcbc_`cX^abbccbbccccccccbbXXaaccb`baa 100 64.29 +XXXXXXXXXX:7:1101:1610:1859#/2 CAATAACACAGTCCAAAATATATGTATGTATGTATATATATATGCAGGCAACAAATAGCACAATGAATAAAATAGTACTTCATCTCAGTACTAACATTGA b_aeeeeeggfgchhiiiiiihhighiighihiihfhihiiiiiehhiihhhhhiiihiiiiiiiihiiihiiiidfhhiihigggggbceeeeeedddd 100 69.64 +XXXXXXXXXX:7:1101:1743:1871#/2 GCTTCCAAGGGACAAAGTCCTGCCCTTGGGTGTGGAAGACACCGTGGACAAGCTCAAGATGCTGGAAGGCCGCAAGACCAGCATCCGCAAGTCAGTGCAG bbbeeeeegggggiiiiihhigdgihiiihdcfhhhihfffhihfdhhhhihhhhhgfggggeeeeeeeca]a_accbcccccccba^_acbbccccccc 100 67.25 +XXXXXXXXXX:7:1101:1536:1878#/2 TGCCTATCATATAGTAAAACCCAGCCCATGACCCCTAACAGGGGCCCTCTCAGCCCTCCTAATGACCTCCGGCCTAGCCATGTGATTTCACTTCCACTCC ___ceeeegfegfhiiiiiii`ecghbcegfhihhiddhifhifgihiihe```efaffeeegfgececeabacccc[`bcccb`cdccccbcY]bccbb 100 65.99 +XXXXXXXXXX:7:1101:2990:100153#/2 TGGAGGCCTGCTGGGAACGGGACTTCTAAAAGGAACTATGTCTGGAAGGCTGTGGTCCAAGGCCATTTTTGCTGGCTATAAGCGGGGTCTCCGGAACCAA ___cccceggecgiihffghighiiiiifffgidfhhiiihiiihiihh^\eeghgffhihggggeeb_aabbdcccccdbbcacacOXbccccaaa^^P 100 61.60 +XXXXXXXXXX:7:1101:1624:1906#/2 GGCTGGTCTGTCCCCTCTGCATCTGGGGCTCCTCCTGGCCCTTCTCCCGCTGCCTGTCCTGCCTCCTGGCCCTGCAGCAGCACCTTGCCCAGTGATCTTT bbbeeeeegfgggiiiiiiihiiihiiihihiihiihihiiiiiiffhiiefhifgghegggggeececddcccc^bcccccccccbbcbccbbbcbbbc 100 67.90 +XXXXXXXXXX:7:1101:1666:1926#/2 GGCTGTGGGAACTCTAGCTGCTGTCAGCCCCTCCTGTAGGAGCTGGCCCTGCCCAGGCTCCTGACTTCCCTCAGGAAGTCTTTCTGTCTTTCTCCATCAG bbbeeeeeggeggiihiiiiiihiiiiiiiiiiiiiiiihhiiiihiiiihiihiihhiiihgbggfgeeeeadd]_a]bcccbccccccccc`bccbbc 100 67.77 +XXXXXXXXXX:7:1101:2921:100163#/2 GAAAAGGAGAACCCCATGCGGGAACTTCGCATCCGCAAACTCTGTCTCAACATCTGTGTTGGGGAGAGTGGAGACAGACTGACGCGAGCAGCCAAGGTGT bbbeeeeeggfggiihiiiiiihhiiihhhfcghfhiiihii]ggfdggfdgeeeebcddcbacccaa]b_ccbcccccccccca_ac_acccacccY__ 100 66.10 +XXXXXXXXXX:7:1101:1513:1929#/2 AAAATGAGGGCCCGTAACAGAACCAGTGTGTGTATAACGAAAACCATGTATAAAATGGGCCTATCACCCTTGTCAGAGATATAAATTACCACATTTGCCT ^__ZZc]cegYZ^decYe[ed^ZbaeYcbf^ecS^ehfghfhiifefhef_effcghghfggfeeeeedbdc`cbc`b`]b``c`bccbcbb[`bccccc 100 62.67 +XXXXXXXXXX:7:1101:2759:100170#/2 GGTCACCATGGATGGGCGAGTCAACGGGGGCCTCAACCTCTCCAGAGCCATTGGGGACCACTTCTATAAGAGAAACAAGAACCTGCCACCTGAGGAACAG bb_eeeeegggggiiiiihifhhiiiiiihiiffggbdeeeeddddccccccccccccccccccccddccccbccccccccccccccccccccbcb^`ac 100 66.80 +XXXXXXXXXX:7:1101:1708:1937#/2 GAATACTTTGTGTTTACCCACTGCCCCCTGACCACTGCTGAAGGCAGAAAGGGACGATCACCTACAGTACCTGGTTTGGGTCTTTATTCTCTCATTCTAG a_aeeeeeggbeghihgfhfhhiiffhiifhhhfhffhiidhfhggiiiihghhiihigggfeadeebbddcbcZ_`^aa^`bbbbbbbb`bc_c`c`_b 100 66.38 +XXXXXXXXXX:7:1101:2981:100211#/2 TTCCTCCTAACTCCTGCCAGAAACAGCTCTCCTCAACATGAGAGCTGCACCCCTCCTCCTGGCCAGGGCAGCAAGCCTTAGCCTTGGCTTCTTGTTTCTG __aeeeeegacgghghhhhfff`cffhfhdghhhhiihhhd[ag`ggghfefhiiiihhghacffd`deeeZ]bb`abcbccb`bcbb``^_bcSY[_`b 100 64.03 +XXXXXXXXXX:7:1101:1688:1946#/2 GAGAAAGGGAAGAAAGAGAAAGGAAAAGACGAAGGGAAGAGGAAGAAAGAGAAAAAGAAAGGGCTCGTGACAGAGAAAGAAGAAAGAGAAGTCGTTCACG aa_`ceeegccg`ffgfgffgdgfhiihhihhiiii[egZefcdfdfgfhegdfhieeggeceedccaccccbbb`bbc`acccccbbbbcbbc_a^bcc 100 66.02 +XXXXXXXXXX:7:1101:2767:100225#/2 TGCAGAATGGCTCCCGCAAAGAAGGGTGGCGAGAAGAAAAAGGGCCGTTCTGCCATCAACGAAGTGGTAACCCGAGAATACACCATCAACATTCACAAGC ___eeeee[eeggffhihhfffhhh]HY^eeYefffhihihgbgggeabdZ`bcccccb`_[_caccbccbccccccccc_bbbc``_]`bb]]`__bP[ 100 57.32 +XXXXXXXXXX:7:1101:1536:1959#/2 GAGACAAAAATTGTTTGAATCAGGACAAGGGAGGTTTGAGGGACTGTTTAACAGGTCCCCATGGTTCACGACCCTGATATCCACCATTATGGGCCCTCTG _W_PY^Y\cc]c[bdJ`YJb^J`eJJJJQWPH^cHY^gc\efUGG^WW\beebgfV\_gg_X^RMHVVVQ^RKXZ`GGGT]ZT]aa``Y]bbc^W_acX` 100 46.90 +XXXXXXXXXX:7:1101:2797:100234#/2 GGAAACGCTGGAGAAGACAGTGTTGCTAGTCTGTTCCTGTTCGTGAGAATGCTTTGAATGGAGGCTCACTGTGCTGCGGCGGCTGCGGCCGTAGCGGGCG aP\JY`JO`YY`QJQ`Q`R`[JQ[^X[IPPXYXYPbbYccee[HWOONNXN]Wd\[HNW___UbGU\\VRV\\]]_BBBBBBBBBBBBBBBBBBBBBBBB 100 38.76 +XXXXXXXXXX:7:1101:1552:1976#/2 TGTGAGCACACCATATATTTACAGTAGGAATAGACGTAGACACACGAGCATATTTCACCTCCGCTACCATAATCATCGCTATCCCCACCGGCGTCAAAGT bbbeeeeeggggghhfhffhgihfggihdbfffgihiiiiiiiiifghiiffgfhhdghhichiegeeeee_bbbddbcccccbccca_accaaaccbc] 100 67.40 +XXXXXXXXXX:7:1101:1748:1978#/2 CAAGGAGCTGTTCTGTCGTGGCTTTGGAATTCAGATGCTCTCACTGTGATAGAGATGGGGTCTTGGGGGGCTGGGAGCTGTTCTGTCGTGGCTTTGGAGT a__eeeeegggggiihhighihiiiiihhiihiiiiihiiiihiiifgeffgghhhihii_fhhiihgecccccc`aacbbccdccca^accccbbcbc` 100 67.54 +XXXXXXXXXX:7:1101:2779:100239#/2 GCTGTGGTTACCTGTCCCCTGCCCCTGGCCTCCAGCCAGCTGCCAGTGGCCCCTCAGCTGCCTGTGGTGAGGCAATACTTAGCTAGGTTCTTAGACGGCC ^[^aaccde`c`^dbbeefdgchcaZZeaf]fcgefff`edgd^e_W^]eZbebd`cddf_ddg`deacdddcV`]b`bcdbc]b_`GT]]b_]`BBBBB 100 45.65 +XXXXXXXXXX:7:1101:1593:1980#/2 GGCTGTTCTCGAACTCCTGACCTCAGGTGATCTGCCTGCCTTGGCCTCCCAAAGTGCTGGGATTACAGGCGTGAGCCACTGCGCCTAGCCGTCACATTTC bbbeeeeegggggiiiiiiiiiiiiiieghhiiiiihhiiiiiiiiiiiiiiihghiiiiighiiiiiiigdddeddddcdbcccccccccacccbccdd 100 68.79 +XXXXXXXXXX:7:1101:2946:100242#/2 AATCCCTTCACCCCAGTTGCCAAACAGACCCCCCACCCCCTGGATTTTCCTTCTCCCTCCATCCCTTGACGGTTCTGGCCTTCCCAAACTGCTTTTGATC b_beeeeegggeegagggdhiihhhihhibffhhfhgghihf_cZcbd]]^`ddee__caGMZ]^Zbb_]\FKT[bb]`cbcbbc`]a^bcc]`bbbO]b 100 52.75 +XXXXXXXXXX:7:1101:1987:1781#/2 GCCGCAACCTGCACAGCCATGCCCGGGCAAGAACTCAGGACGCTGAATGGCTCTCAGATGCTCCTGGTGTTGCTGGTGCTCTCGTGGCTGCCGCATGGGG bbbeeeeegggggiiiiiihiiiiiiiiiiiiiiiiiihhihiiiiigggggeeeeeeddddddcccbbacccccccccccccccccccccc[_aacaac 100 67.28 +XXXXXXXXXX:7:1101:3046:100006#/2 GACGAGATGGAGTTCACCGAGGCCGAGAGCAACATGAATGACCTGGTGTCCGAGTACCAGCAGTACCAGGATGCCACAGCCGAGGAGGAGGGCGAGTTCG bbbeeeeeggggehiiiiihhihiihghhhiiihhiiifhiiihiihigfdgeeedddddcccccccccccbccccccccccccccccaccacc[aX_ba 100 66.53 +XXXXXXXXXX:7:1101:1900:1788#/2 CGGCTATGCGCCACCCGAGACGGAGATCGTGCCGCTTGTGCTGCGCGGCCACCTCATGGACATCGAGTGCCTGGCCAGCGACGGCATGCTGCTGGTGAGC bbbeeeeefggggiiiiiiggiiheghiigghhhiiihhhiiiiigeccc__ccccbcccccccccccccccccccccccccccccccccccbcc[^b_b 100 66.67 +XXXXXXXXXX:7:1101:3214:100027#/2 CACAACAACTGCTACTTCAGTTGCGCAGACAGTATCAACACCCACAACACAAGATCAGACCCCAAGTTCTGCTGTTTCAGTTGCCACGCCTACAGTTAGT _bbeeeeefeggghhiiihhghhhhhi[fhhf^cghiifhhiagiiiiihfbffggiiggggecdeZ`bdddcccccbdccccbcbcccccccbbcccc] 100 66.59 +XXXXXXXXXX:7:1101:1848:1879#/2 GGAGAAATGACCTTGAAGTCAGCACTTGCTTTAGCTATCAAAGTACTAAATAAGACCATGGATGTTAGTAAACTCTCTGCTGAAAAAGTGGAAATTGCAA bbbeeeeegggfgiiihighiiihiiihiiiihhiiiiiiiiighhiiihegghiihhiiiiiiiiiiiiiiiiiiiiiiigfgfgee^beddcdbdbcc 100 69.46 +XXXXXXXXXX:7:1101:3237:100032#/2 CGCTTATGCAGACAGAGCACAGCTGACTGCACGAACCAAAGGAGACGCATGAGCCAGAGAGACAAGGAAGGCTTTCAAAGGAAGGGATGGTCAAGCTGGG bbbeeeeefggggihiiihiiiiihiiiiiiiiiihihiiiiiiiihihifhdfgfgggfeeeeddcdcccccccccccccccaccccccbbbbcbbccc 100 67.97 +XXXXXXXXXX:7:1101:3027:100049#/2 CTCCACCACGGTCCCCAAGATGATTGTGGACATCCAGTCTCACAGCAGAGTCATCTCCTATGCAGGCTGCCTGACTCAGATGTCTCTCTTTGCCATTTTT bbbeeeeegggggiiiiihhihhhihfhhhiiihhiicfgiiiiidffdhighihhiiifihiiifgggfeeeeecdbdddcdcdccccccccc_bcddc 100 68.43 +XXXXXXXXXX:7:1101:1756:1891#/2 CTCATATGAAGTCACCCTAGCCATCATTCTACTATCAACATTACTAATAAGTGGCTCCTTTAACCTCTCCACCCTTATCACAACACAAGAACACCTCTGA bbbeeeeeggggfiiiiiiiiiiiihiiiiiiiiihhihiiiiiighhiihhhiiiiiiiihhghhiiiiiiiiiifgggggeeeeccdddccccccbbc 100 69.48 +XXXXXXXXXX:7:1101:3238:100065#/2 GGGGAAGGAAAACAGGGTATGGGGGAAAGGACCCTGGGGAGCGAAGTGGAGGATACAACCTTGGGCCTGCAGGCCAGGCTACCTACCCACTTGGAAACCC bbbeeeeegggggiiiidebhihihfihiihiiiiiiihhiihigeceeebcbdcccccccccccccccccccccccaccbccbcbccccacccbbb`ac 100 67.08 +XXXXXXXXXX:7:1101:1915:1901#/2 TACTCATTCAACCAATAGCCCTGGCCGTACGCCTAACCGCTAACATTACTGCAGGCCACCTACTCATGCACCTAATTGGAAGCGCCACCCTAGCAATATC _bbeedeegeeggihhicdhhhidghheegihhhhiiiiihihifiihhigiihiggdeeeedcdbdbdcccbcccdc`bbbccccccccccbbcbcccd 100 67.56 +XXXXXXXXXX:7:1101:3198:100082#/2 GTCAGCTTACATCCCTGAGCAGGAAAGTTTACCCATGAAGATTGGTGGGATTTTTTGTTTGTTTGTTTTGTTTTGTTTGTTGTTTGTTGTTTGTTTTTTT babeeeeeeegggihihiiiiiiiiiifhihhhihhhgffgihhhX^ege]eghiihihffdggfgededdbddccccbccab`b^^bccccccbcbaca 100 66.54 +XXXXXXXXXX:7:1101:1964:1931#/2 CCCAAACCCACTCCACCTTACTACCAGACAACCTTAGCCAAACCATTTACCCAAATAAAGTATAGGCGATAGAAATTGAAACCTGGCGCAATAGATATAG a__ecece_agfefgfghhhhhehhfeeghfffffhfgffcffddfghihfgffghhihhefddffhgd]cdddcdc`bcccab^^^_acaccccc_bcc 100 66.53 +XXXXXXXXXX:7:1101:3088:100091#/2 AACATCTTAAGTCTTACCAGCCAGAATGTAGATACATAACTACAGTTAGTTGTGGAAGCTTTGGGATTTTAGTTGAGTGCTATCATCCCCAAATGTCATT __aeeeeegfgeghfihiiiiihiighigfggifhiiiiihihhhfhiihgghghhihbdgiihiacehhghhhhhfZbdeggeeedacbdadbbbcdcd 100 67.89 +XXXXXXXXXX:7:1101:1840:1948#/2 GGCGGGCAGTAGCCGCTGAGGGGATTGCAGATAACCGCTTCCCGCACGGGGAAAGTCTACCCTGCCTGCCACTTTCTGCTCGCCGTCAGCGCCGGAGCTC bbbeeeeegggggiiiiifghiighiiiihhfhigggeeecddccccccccccccccdcccccccccccccccccccccccccccccccccacccaca^b 100 66.86 +XXXXXXXXXX:7:1101:3105:100094#/2 GCTTTGTCAAGCCCGCTCTGCTGCATCAGTAGTGAATGGAGAGAATTATTTTGGCAGGGCTGGGATGGAGTTCAGGTGCCCTTGCTTGTTACGTGAAGTC ^V^ccce`egeegh[b[dbafbcagacce^SXe`_caee^afbf]ac[I^dfaeefZaBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB 100 36.74 +XXXXXXXXXX:7:1101:1958:1952#/2 CTTACTACTAATAATTATTACATTTTGACTACCACAACTCAACGGCTACATAGAAAAATCCACCCCTTACGAGTGCGGCTTCGACCCTATATCCCCCGCC bb_eeeeefgfgehhfhhihg]gfgfffgdghicegdhhhffgefhihiifhfgiiiif`fhffhh\dbdb^a^^b`_^aacN[_aaccbbd[]`_aX_a 100 62.07 +XXXXXXXXXX:7:1101:3190:100106#/2 AAAACACTGAACTGACAATTAACAGCCCAATATCTACAATCAACCAACAAGTCATTATTACCCTCACTGTCAACCCAACACAGGCATGCTCATAAGGAAA bbbeeeeegggggiiihiiiiiiiiiiiiihhiiihhiiifhiiihiiiieghhhihiiihfhhihiihhiiggggeeeeccdbcccccbc`bcbcbbbb 100 68.77 +XXXXXXXXXX:7:1101:1993:1999#/2 GCCCCCCGCACAGCCACCCGTGGGGTGCTGCCCTTATGCCCGGCGCCGGGCACCAACTCCATGTTTGGTGTTTGTCTGTGTTTGTTTTTCAAGAAATGAT bbbeeeeeggfggiiiiihhcffhh_e_dgfhgffgdcebecccccccccccccccccbcbccccdcc^^`acc`ccb_b]bbbP[^bc^_bbbcbcccb 100 63.20 +XXXXXXXXXX:7:1101:3117:100110#/2 AGGACTTCAAACTCTACTCCCACTAATAGCTTTTTGATGACTTCTAGCAAGCCTCGCTAACCTCGCCTTACCCCCCACTATTAACCTACTGGGAGAACTC bbbceeeefgfggiiiiihiihiihhiiiiiiihihefhiiihihihihifdhiihhiiiiiiiiiggeeedecccccccbcddccccccc\accccccc 100 68.17 +XXXXXXXXXX:7:1101:2159:1798#/2 CGGGCAGGGACCTGCAGGAGGCTTTCCCAGCAGGTGCACAGAAGGAAGCGCTGAGCCCACTGTGTGCAGAGGGTGTGTGCCCAGGGCACTGCCACGGCCC bb_eeeeegggggiiiiigiihiiiiiihiiiiiafghhehiiiihiiiifgigfggededdbdddcdccccc[_aacabbbcccccccbccccccccc] 100 67.04 +XXXXXXXXXX:7:1101:3147:100111#/2 CAAAATTTTGACAGGGAAGGAAGTAACTGCAAAACACTAGGCTTTAGTAGGTACTTATATAAAATCTAGTCCAGTTCTCTCATTTAAAAAAATGAAGACA bbbeeeeeggegghiighhihiibgfhiihhiiiiiiihihihiihhhhieeghiiihhiiihihihhhhiiiihfgcgggeeeeeeeddccccccbc`b 100 69.07 +XXXXXXXXXX:7:1101:2152:1838#/2 ATTCGGCGCATGAGCTGGAGTCCTAGGCACAGCTCTAAGCCTCCTTATTCGAGCCGAGCTGGGCCAGCCAGGCAACCTTCTAGGTAACGACCACATCTAC bbbeeeeegggggiiiiiiighiiiiiiiiiiiiiiiihiiiiiiiiiifhffhiggeccddccccccccccccaccccccccccbcccccacccccbbc 100 67.83 +XXXXXXXXXX:7:1101:3065:100152#/2 GGAGTGCGTGCCTGCTACGAACGCTCCAACGAGTACCAGCTGATTGACTGTGCCCAGTACTTCCTGGACAAGATCGACGTGATCAAGCAGGCTGACTATG __^cacccbeegghhiihhihihhhhchifdhfbgfhddffiibedffhhgffhiiiegfdggeeeca`bcddccaccaaaX^bbbbcc[^TW`W^``_b 100 63.71 +XXXXXXXXXX:7:1101:2180:1843#/2 ATCTGCTGCAGTGCTCTGAGCCCTAGGATTCATCTTTCTTTTCACCGTAGGTGGCCTGACTGGCATTGTATTAGCAAACTCATCACTAGACATCGTACTA bbbeeeeefgggghhhghihifiihgf_gfcdhhhiifhiihihiiighhhfffgdgfhiiihhhdffbfgffgeg^dcbeeeecdbbbdccccc`^`a_ 100 67.67 +XXXXXXXXXX:7:1101:3154:100159#/2 GTTTAGATAACCTCTAAGATCCAGCTTTAAAATTCTTGACTGAAGGCCTAACCTTTAGAAATATCTCGTTAGAGGCCAGACGCGGTGGCTCACGCCTATA babeeeeegggggiiiiihiiihiihiiiiiiiiiiiiiiiiiiiihiiihihiiiiiiiiiiiihihhiiihihggegeeeccc_cccccbccccaccb 100 69.02 +XXXXXXXXXX:7:1101:2125:1861#/2 CCCTCCACCACCTCTGTGGCTCCAGCTACCACATCTGGAGGCTTCTCATTCACTGGTGGAAGCACGGCCCAACCCTCCGGTTTCAACATTGGCTCAGCAG bbbeeeceggggfihihihhiiihhfihiihiihhhiihfhgfhhiififfhhfhhfghbghhcgggadcab^_^`ccaaW_]bbcbbbbbccbb]bcbc 100 66.16 +XXXXXXXXXX:7:1101:3198:100173#/2 CGCTTCCATGCCTTTATCCTGGGCCTCCTCAACACCAAGCAGTTGGAGCTGTGGTTTTCCAGTCTCCAGGAAGATGCAGGCCTGCTCTCCCTCCTGTACC _b_eccec`ggeehhihdfdgdf`dfhihfhhhhiiiihihhdgdhdfg_eecdebegfbdeadfdfhfgddaaadc_`accc\`c]bbcb^^abbYb`b 100 65.68 +XXXXXXXXXX:7:1101:2076:1911#/2 CCCAGGTATTGGGGACAGATAATGGGCCTGCCTTCATCTCCCAGGTAAGTCAGTCGGTGGCCAAGTTGCTGGGGATTGATTGGAAACTGCATCGTGCTTA abbeeeacggfggihihihiiiiihiiiiiiiiihhiiifiiiiiefhgeffgaehh`fhihgdgbgeeeeeccbacccd`bcbbbcbcbbbcaa_cccc 100 67.53 +XXXXXXXXXX:7:1101:3166:100190#/2 TGTGGTGGTTTCTCACAGCAGCCATGGGGAACGGAAGCGTCCCCAGGCCACAGCCTCCCCGCAGTGTGCGGACTCGGCCTGCTGAGCCTGCATGTTGGCG _a_ccacceeecghgfhhhcgfhigghhhiOb`g`ccbfcghhfbfcdgdcgbdegbcceccccW`Z]`^acccc^^a_accccccb^][[b_]SSGW__ 100 56.13 +XXXXXXXXXX:7:1101:2196:1920#/2 ACGGGGCTGTCTGCGTAAAACTTGCCACTTAGAGCCCTGCCTGCACCGTCGTCACGGCCCAGCCTGCAGGGTCGGGTGTGCTGACAGCCGGCTGTTAAGC abbeeeeegggggiihiifhiiiiiiiiiiihiiiiiiiiiifghiiiiifhhggfeeccacccccacccc`bacc^a_ccccccccccccccccccdbb 100 67.38 +XXXXXXXXXX:7:1101:3225:100207#/2 CTCCAACTCATGAGACCCACAACAAATAGCCCTTCTAAACGCTAATCCAAGCCTCACCCCACTACTAGGCCTCCTCCTAGCAGCAGCAGGCAAATCAGCC bbbeeeeegfgggiiihiiiiiiiihihfhhifghihihhiiiiiiiiiihghhiiihihhgggeeeeeedcdcdcbcccccccaacccbcaa^bcc__b 100 67.82 +XXXXXXXXXX:7:1101:2115:1927#/2 CGTGAATCTAAGAAGATCGATCGCCACATGTATCACAGCCTGTACCTGAAGGTGAAGGGGAATGTGTTCAAAAACAAGCGGATTCTCATGGAACACATCC bbaeeeadgggfghhiiiiiiihghiiiihehfgfhiffhifhhfiiigfegWbghfdghfhiigggfggeeeeecccccaacccccbccbcbc`bcccb 100 67.01 +XXXXXXXXXX:7:1101:3019:100219#/2 CCAACCTCCTACTCCTCATTGTACCCATTCTAATCGCAATGGCATTCCTAATGCTTACCGAACGAAAAATTCTAGGCTATATACAACTACGCAAAGGCCC a__accacagfaehihiffgh`efdcfff]gb`fgcfh_fffdfgg[abcdeghfhhfhcffffcedec_b_`dY_dcccdccccccbccaca[Xb`WW^ 100 63.93 +XXXXXXXXXX:7:1101:2179:1937#/2 CAAGAACAATAAATGTCTATCCTAATTTCAGACCCACTCCTAAAAACTCACTCATGGGAGCTCTGTGTGGATTTGGGCCCCTCATCTTCATTTATTATAT bbbeeeeegggggiiiiiiiihiiiiiiiihdhhhiiiiiiiiiiiihihgghiiiiiiiiiiiigiggidhggggfcedccccccccd_ceececdede 100 69.03 +XXXXXXXXXX:7:1101:3202:100230#/2 TGTACAAAGAGTAGCTGCCCCTTCTCAGGGCAGATGATGCTTTGAGAACATACTTTGGCCATTACCCCCAGCTTTGAGGAAAATGGGCTTTGGATGATTA _a_ceccea`aceggfhfihihiidhdh]e`eg^efgdhdhhihfhacgdefcgfhacbbegdffecga`ggeedb_c[`]_bbccaccccbccccbcbc 100 65.88 +XXXXXXXXXX:7:1101:2149:1945#/2 GCCTATCATATAGTAAAACCCAGCCCATGACCCCTAACAGGGGCCCTCTCAGCCCTCCTAATGACCTCCGGCCTAGCCATGTGATTTCACTTCCACTCCA bbbeeeeeggggfegghiiiihhiihififfhiiihihihiiihhihiihhhhiihhiigeffggeeecbbccccccccccdcccccbdcccccccb``b 100 68.03 +XXXXXXXXXX:7:1101:3211:100242#/2 TGAGGGCGTTGGGATGGCTTGTAGGATTTTTAAGAGAGGGTGTCTACAGACAAGCATTTTCTCTGTAGAGCAGCCACACGTTGTATAAACATAAACTTTA _a^ccccc_eeaeZ^bgg_cac[^^X^ccffg_a_fffGVaPZ_b]]e`c`cd``]db^]ZbZaZ]_b]]bbbBBBBBBBBBBBBBBBBBBBBBBBBBBB 100 38.61 +XXXXXXXXXX:7:1101:2169:1964#/2 CTAGGACAAGGGCCAAGAACCGATGGTACCCACCTGGGCCCCGGCTCAGGGCCAAGAACAGATGATACCCAGATAAGGCGGAACCAGCAACAGTTTCTGA bbbeeeeegfggfiiiiiiiihhiihhhiiefhihhiiiihiigiggegedceccccbccccccbcdcbbccacccc]]acccaccaccccccbcdcccc 100 66.99 +XXXXXXXXXX:7:1101:3168:100244#/2 ACTGGATTAAGAAAATGTGGCACATACACACCATGGAATACTATGCAGCCATAAAAAATGATGAGTTCATGTCCTTTGTAGGGACATGGATGAAATTGGA _bbee`eeeegfgiihigfghiefghgifhihdhiihiidghhffhcgehfd_^eafgiiihifhcgghhggeggfgeacdecbbcbccbbbbccccccc 100 67.61 +XXXXXXXXXX:7:1101:3005:100246#/2 CAGGGAGGGAGACACCTAGATTAGCAGCTCAATTTGTACTACTTCAGCCAATCTGTGAATTTAAAAACTACACTGTTGCCTTGCTAGGATCCACCCTCCT aabeeeee`gggghihiiihiiiiiiiiihiihiiihiifhiiiihiiiighiiiiihihiiiiihiigggggeeeeeddcdddddbaccbcbc`_cccc 100 68.54 +XXXXXXXXXX:7:1101:2313:1789#/2 GCTCGCCGCTCGCTATGGCGTCGCTCACCGTGAAGGCCTACCTTCTGGGCAAGGAGGACGCGGCGCGCGAGATTCGCCGCTTCAGCTTCTGCTGCAGCCC bbbeeeeeggggfiiihhihfhiiihiiihifdhghggggeeeeeddcdcdccccc\[`ccccccV_aaXZ]_bbccaa___ab_bbcccbb_Y_bbb`a 100 64.11 +XXXXXXXXXX:7:1101:3253:100014#/2 CAAAACTAACCCCCTAATAAAATTAATTAACCACTCGTTCATCGACCTCCCCACCCCATCCAACATCTCCGCATGATGAAACTTCGGCTCACTCCTTGGC bbbeeeeegcgggiiiiiiiihiiiiiiiiiiiiihVafhiiidhh[eggghdghhigfgggeeeeedddccccc_bcdcccccccacaccccccc_`bc 100 66.29 +XXXXXXXXXX:7:1101:2361:1794#/2 CGTGGTTCACTGCAGCCTCAACCTCCTAGGCTCAAGCAATTCTCCTACTTTCCTGAGTAGCTGGGACTACAGGCATGCGCCACCATACCCACCTAATTTT __^c^`c`cgeegghhhi`gghagihddgfegfhiihihfhihiiifghihghihhifhfbcggcdecghagggggebccc`a`^b_bb`bcaaRYY]bB 100 52.68 +XXXXXXXXXX:7:1101:2337:1794#/2 CGACAACTCCTCCGGCCGGAACAGCATGGAAGGTCCAGCGTTCTCAAAATCCCTTAAAGATAAGATTAACCCGTGGGGCCCCCTAATAGTTATGGGGATC ab_eececggegehhhhhfghfffgbgffihefWbafbg`_cghhdbdgg]_bccdcdbb`_bb]b`_bcc^\WNZLXTT[acaabb]`GYY`]`W^GEW 100 51.43 +XXXXXXXXXX:7:1101:3284:100039#/2 TTCCTCAATCTGAACACCATGCGCCCGCGGCTGCGCTACAGCCCCTACTCCATCCCGGTGCCGGTCCCGGACGGCAGCAGTCTGCTCACCACCGCCCTGC bbbeeeeegggggiiiiiiihhghghhfiihiidgfhfgfegeeeecdddbccccccccaacccZaccccacccccccc_bccccccccccacccccccc 100 66.60 +XXXXXXXXXX:7:1101:2477:1795#/2 GGCCTGGATTCCACGGCCATGGCCTCTGCCGCTGCGGCGCAGGGACTGTCCGGGGCGTCCGCGGCCACCCTGCCCTTCCACCTCCAGCAGCACGTCCTGG ___ceee^ceegfifhfhcfhhiihhihhffgghiaffZd]a```cccccbaaacacaZ_acVZ]aac^ZWX`aabb[X]Y^aaW``]b[`a_TW`^WSX 100 60.46 +XXXXXXXXXX:7:1101:3310:100056#/2 CGCAAAACTAACCCCCTAATAAAATTAATTAACCACTCATTCATCGACCTCCCCACCCCATCCAACATCTCCGCATGATGAAACTTCGGCTCACTCCTTG bbbeeeeegggggiiiihiihiiihihfhhiiihiiiiihiiiiihhiiihhfhhiiiigggdeeeedddddccccccccdcccbcbcccccbbb]bcbc 100 68.01 +XXXXXXXXXX:7:1101:2355:1821#/2 GTCTTCGTCATGCTCAATGAAAATGTACTAAGTTCTCCCTCTGTGCTAGTAACAGTGATTCCCCACCATTCTCCTAAAAGAAAAAATTCTAGTCCAAATA aaaeeeecgggggiiiiiiiiihiighiihiighhiiiiiiiihiihiffeghihghiiiiiiiiiiiiiihiiiifggggggeecddcdccdccccccb 100 69.29 +XXXXXXXXXX:7:1101:3420:100060#/2 CCCAGAAAGGCTATTTTCAGATGAAATCGATATTAGAAGCTATATTAGCTGAAACAACTCCTTTTACTGCGTAGAACCTATATCGAGAGTGTGTGTATAT bbbeeeeecfgggehifhiiihhfhihhiffhhhhhdfgfihhihiiihgfhhihiiiihiihhiiihhi]dgeg`gceeeedccca^aZ]_bQ^[_acb 100 64.38 +XXXXXXXXXX:7:1101:2418:1834#/2 GGTGGCCCTCAATCCTGCCACCTTGCTCCCTCTACCGGAAAAAGGAGCCCCCCATGATTGCCTCGAGATCTTGGCTGAAACGCATGGAACCAGACCGGAT ab_ceeeeggggghiihihiiiiiifhiiihffgiiiifhiihfiiiiiighgaeeeedddddcccccaccccccccacccccaac^]_bcccacca[aX 100 66.17 +XXXXXXXXXX:7:1101:3267:100061#/2 GCCAGCTGCATGTTGTGAGGACACTCAAGCAGCTTTTTGGAGAGGTCCCATGGTGAGGTATTGAGATTTCCTGCCAATGTCCACATGAGTGAGCTCCCTT bbbeeeeeggggegigfihidfdiiiiiihfhhhhiihighhhhighiiifghfghghbeghihihegfbgggeeeeedddcdbbccccbbcccbccccc 100 68.35 +XXXXXXXXXX:7:1101:2378:1838#/2 GCTTGCTGCTGCTGGAGGTGTTTCCCATGATGAATTGCTTGACTTAGCAAAGTTTCATTTCGGTGACTCTTTATGCACACACAAAGGAGAAATACCAGCT bbbeeeeeggfggiighi`deghiiiiiiiiiiiiihihiihiiiffihiigfghihiifhhihihhhiiiihgggfggeeeccddcccccccccccccb 100 68.88 +XXXXXXXXXX:7:1101:3416:100083#/2 TGATCAACATCTTTTCTTGCCTCTGTCCCCTTCTCTCATCTCTTAGCTCCCCTCCAACCTGGGGGGCAGTGGTGTGGAGAAGCCACAGGCCTGAGATTTC b_beeeeegggggiiiiihiiihfhiiihiiiiiiiiiiiiiiiihhhiffghiiiiiiiihia`]accbccbaaa^bacccccbbaac_^ccccbbbcb 100 67.08 +XXXXXXXXXX:7:1101:2481:1853#/2 GCCTACTCATTCAACCAATAGCCCTGGCCGTACGCCTAACCGCTAACATTACTGCAGGCCACCTACTCATGCACCTAATTGGAAGCGCCACCCTAGCAAT bbbeeeeegggggiiiiiiiiiiiiiiiiihiiiiiihiiiihiiiiiiihiiiiggggeeeedddcdcccccccccccdcccccacccccccccccccc 100 68.19 +XXXXXXXXXX:7:1101:3411:100111#/2 CACATGCCTATCATATAGTAAAACCCAGCCCATGACCCCTAACAGGGGCCCTCTCAGCCCTCCTAATGACCTCCGGCCTAGCCATGTGATTTCACTTCCA bbbeeeeegggggihihiggiiiiihiiiiiiihhihihfhiihhfiihdaghiiiiiiihggggeeeeeeddcccccccccccccbbc_bc__bbbbcc 100 67.87 +XXXXXXXXXX:7:1101:3258:100128#/2 CCAGCCCCAGGTCCTCGGACACCGAGGAGAATGTCAAGAGGCGAACACACAACGTCTTGGAGCGCCAGAGGAGGAACGAGCTAAAACGGAGCTTTTTTGC _bbeceeegfggghfhhiihhhiidgf`gfhdh^effcebeghgdcedcc`^_aaa[^bbbbca]_Vaacaccc^W^_a_X_b`b``[^_aT[]`bcaX` 100 62.43 +XXXXXXXXXX:7:1101:2252:1856#/2 GAAGGGTTCACGGGAAGAGAAGCAGAAGCCCCAGGCTGAGCGGAAGGAGGAGAAAAAGGCGGCTGCCCCTGCTCCTGAGGAGGAGATGGATGAATGTGAG __Zc``cccagaefhihicggffffffdf]bcdP^aXaQaeghhghcddgPTZZ__GX\ZX_Zaa^`^[W[abb]]]GRWY[XGGGRY]RY_Y_Y_bc[b 100 50.14 +XXXXXXXXXX:7:1101:3428:100129#/2 CAGTTCTTCAAATATCTACTCATCTTCCTAATTACCATACTAATCTTAGTTACCGCTAACAACCTATTCCAACTGTTCATCGGCTGAGAGGGCGTAGGAA a_bceceeegegghhhhiiffffhfiiifihhihhhhffghiiibghidgccgeghhhhidghffhiifihhhfhgegcfggeeeccccca_aa_\]^b` 100 67.40 +XXXXXXXXXX:7:1101:2394:1871#/2 GCTAAATAAGCTATCGGGCCCATACCCCGAAAATGTTAGTTATACCCTTCCCGTACTAATTAATCCCCTGGCCCAACCCGTCATCTACTCTACCATCTTT bbbeeeeegggggiiiiiiiiiiiiiiiiiiiiiifhiighiiiiiiiiighiggggggeeeeeeddcdccccccccccccccccbbcdcccccccccbc 100 68.16 +XXXXXXXXXX:7:1101:3387:100138#/2 GCCAGGCCTGCCAGGGAGGCCAGGGCTGACTGCTCTTCCCCTCCCATGGCTGGTGCCCAGCCCCCCCCCCCCCCCCCCCACAAACAAGGGGCCCGGCCCC bbbeeeeegggggihifhihihiiihiiihefhfhhhhiiiiiihfdW\ebdgZc`ggbdgeeecccccccBBBBBBBBBBBBBBBBBBBBBBBBBBBBB 100 38.37 +XXXXXXXXXX:7:1101:2269:1904#/2 AAAGGAGGCAGATGGGTCCTAGAGGCAGCCCATGGCAAGGGAGAGGTCTGCTTTTTTCTTTTTTTTTGAGACGGAGTTTCACTTTTGTTGCCCAGGCTGG _[_c`Za]ecceghfeeae`Xecee^bghffccgfdf]egfggOacM\Zeefbd]_[bg]bcc[]aaZ^ccY[OETJR]Y]_YS]`YJRWSSGJQ[[`aB 100 48.61 +XXXXXXXXXX:7:1101:3444:100163#/2 CGTGTGTGCGCCTAATCTCAGGTGGTCCACCCGAGACCCCTTGAGCACCAACCCTAGTCCCCCGCGCGGCCCCTTATTCGCTCCGACAAGATGAAAGAAA aaYY\YQ\\_Y[[behifh`afQXebdefd^^U__U_cac^\aeg_edgf\Z_\^`]\bRZZ\[_^WWZTTXX`ab]bY^T_`Xa[ET[Q^_`bb]]b`[ 100 53.54 +XXXXXXXXXX:7:1101:2259:1943#/2 AAAAATTAGCTGGGCATGCTGGCGCACACCTGTGGCTCTAGCTACTCGAGGAGATGAGGTAGAAGGATCACTTGAGCCCAGGAGTTGGAGCCACAGTGAG ^__aaca`beacg`^`bb_dg]bcgdgHPY^e^cgcaefffhdc`c[Ze``TbUV\HV`M\a^aZaUZ`]]bbcR]]^^^`_[WKY]]b^BBBBBBBBBB 100 42.61 +XXXXXXXXXX:7:1101:3371:100179#/2 CCGGAATGCCCCGACGTTACTCGGACTACCCCGATGCATACACCACATGAAACATCCTATCATCTGTAGGCTCATTCATTTCTCTAACAGCAGTAATATT bbbeeeeegggggiiihiiiiiiihiiiiifhiehhiiiiiiiigggfeeeeeeddddddccdcccccbccccccdddddeddccdbccccccbcc`bcb 100 67.68 +XXXXXXXXXX:7:1101:2371:1957#/2 CGCTAAATCCCCTAGAAGTCCCACTCCTAAACACATCCGTATTACTCGCATCAGGAGTATCAATCACCTGAGCTCACCATAGTCTAATAGAAAACAACCG _b_eeeeecefgghfhiihhfhiiidfhhhdfhhii_ghhhiihieihhhdfghihheghiihhfgggfb]aceab`a]]bb`ccccdcb`bbccccccc 100 66.91 +XXXXXXXXXX:7:1101:3311:100186#/2 AAGTAATTGATAAGCTGGACTGTGCTAAAATTAGAAATATCTGCTCTGTGAAAGACAATGTCAAGAAAATGTAAAGACAAGCTACAGAATGGGAGAGAAT _b_cceeegggggiiiiifhiighihiiiiiihiiiiiiiiiiihiiigifghiihiiiiiiiiiiihihigghhihiiiigfgggebeeeeddcbcacc 100 69.35 +XXXXXXXXXX:7:1101:2394:1961#/2 CCCGAATATGACTACCTGTTTAAGCTGCTTTTGATTGGCGACTCAGGCGTGGGCAAGTCATGCCTGCTCCTGCGGTTTGCTGATGACACGTACACAGAGA b__eeeecgggggiiiiihhifhghhhhhiiiifiiffhihiihfhaffff]`gbgggggfgecbeeddddbbac^\_abb_bbbbbbcc`^accccb^[ 100 66.13 +XXXXXXXXXX:7:1101:3438:100192#/2 CCTGAGCTCACCATAGTCTAATAGAAAACAACCGAAACCAAATAATTCAAGCACTGCTTATTACAATTTTACTGGGTCTCTATTTTACCCTCCTACAAGC bbbeeeeegfggfiiifgiiiiiiiiiihiiiihgfhiiiihiiiiiiiihhfhhiiiiiiiihifheggggfeeeccbdddcddccccccccbcbcbcc 100 68.74 +XXXXXXXXXX:7:1101:2333:1962#/2 CAACCCGTCATCTACTCTACCATCTTTGCAGGCACACTCATCACAGCGCTAAGCTCGCACTGATTTTTTACCTGAGTAGGCCTAGAAATAAACATGCTAG ^a_ccc^O[beeedfh`e_gdffgdc]cedcehdR^e^RcfhhdcX^e`gafbfffa_FTS^RMMM]_aZ]Z]``XY]RTGTGY[bbb`bcb`__b`Y`_ 100 51.00 +XXXXXXXXXX:7:1101:3479:100209#/2 CCAGAAACTGCTAAAGAACATGGCCTGTTTGACATGTTCATGAGTCACCTGACCCACAGCATATATGCTTATGACTAAACCCTCCACTCCTGATTCTCAA bbbeeeeegggggiiiiihihiihiifeghfgiihhchghiiiheghiiihiiiiiiiiiiihiiiiiihiihiicgggggeeecddddcccbcccccdc 100 69.18 +XXXXXXXXXX:7:1101:2459:1990#/2 CCACAACCTTAACAATGAACAAGATATTCGAAAAATAGGAGGACTACTCAAAACCATACCTCTCACTTCAACCTCCCTCACCATTGGCAGCCTAGCATTA bbbeeeeegggggiihhiihiiiihfghhiiiigegfhfghiiiiiiiihihhhhiihihhhiiiiiggfdggeeeedddcccccccccccccccccccc 100 68.75 +XXXXXXXXXX:7:1101:3417:100210#/2 GCAGACTACTACTATGGATAACTACAGTGCGATTCGTTCATGACTGAGCGTGGAGCAGAGTCCACACAGAAGTTGATCGTGGTACCCAAGGCTGCTTCTC bbbeeeeegggggiiiiiiiiiiiiiighiiiihiihiiiiiiiihifiiighhiiihiifgiiiggggeee_cedcdcdcc``ccccccccbccccccc 100 68.39 +XXXXXXXXXX:7:1101:2372:1994#/2 GTTTGGGCAGGATGAGATGATCGACGTCATCGGGGTGACCAAGGGCAAAGGCTACAAAGGGGTCACCAGTCGTTGGCACACCAAGAAGCTGCCCCGCAAG _^^c`ccabee`ca]f[dd[]J`dPbbedX_gghiOVUbf]_bfg_d\db^Z^_`_`b_b__NFXT_]`TYWWXX`YY]X^^^`W^`bccbbcccX]cBB 100 48.35 +XXXXXXXXXX:7:1101:3452:100214#/2 CGTTTTTGGAAACCAGTGCTAAGAATGCAACGAATGTAGAACAGTCTTTCATGACGATGGCAGCTGAGATTAAAAAGCGAATGGGTCCCGGAGCAACAGC bbbcecccgeggghiidcggfhhdfhhbffehhhhi]ccfghiicfffffbcfbegba_efgfhhgc^dggccee[R`[_^NGYQYR[a^a]_[`X^`^a 100 55.72 +XXXXXXXXXX:7:1101:2677:1830#/2 TGAAAATGCTCTTTCTGGAAATGAGGAACTAACAGTCAAAATTAAGTGTGATAAGGAGAAGAACCTGCTGCATGTCACAGACACCGGTGTAGGAATGACC bbbeeeeegggggiiiiheghhhihhiihhhhfhighiiiiihifhefhghhiiiiiiiiihiihihhhiiiiiiigfggfeeeecc^a_ccccbccccc 100 68.81 +XXXXXXXXXX:7:1101:3354:100219#/2 ACAGCCGCTTTGACCGGGATTACCAGCAGTACCATGGACTACGGAAAGAAGTGGTGGACACAAGCCCAAGCCCACGACCTGGTCATCAAAAGCCTTGACA bbbeeeeegggggiiiiihiiiiiiiihighhiiiihfiihiihfghibgggggdgeceeedcdbaaaaaccc`a__accccbbccccccccccccbcbb 100 67.23 +XXXXXXXXXX:7:1101:2603:1846#/2 GCTAATTTTTGTATTTTTCATAGAGACAGGGTTTCACCATGTTGGCCAGGCTGGTCTCGAACTCCTGATCTTAGGTGAGCCGCCCGCCTCAGCCTCCCAA bbbeeeeeggfegiihhhhifihiifiiiffbfhihdcghfehheghhhhhfhh`efhhaegihihfcfhiigffdgdceecccc\^cccccbbbcccbc 100 67.68 +XXXXXXXXXX:7:1101:3600:100019#/2 CCCCTGGTCAACCTCAACCTAGGCCTCCTATTTATTCTAGCCACCTCTAGCCTAGCCGTTTACTCAATCCTCTGATCAGGGTGAGCATCAAACTCAAACT bbbeeeeeggfggiiiihiiiihhiiiiihiiiiiiiihihihhghihiigghhhhhhfgfhihiiiiiggggggeeeecc]^`dbbccccccccccccc 100 68.39 +XXXXXXXXXX:7:1101:2535:1848#/2 CCCGAGTTATCGGGACCCTGTTGGGAACTGTCGACAAACACTCAGTGGAGGTCACCAATTGCTTTTCAGTGCCGCACAATGAGTCAGAAGATGAAGTGGC _bbeeeceggggfhihiihihiihiihiiihhhiiiiiiiihiiifhhhigggggggeeedeedddccddcccccZZ\X`bbcbccbccccccc_Y_bcc 100 65.68
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/r2_seqkit_stats2.tabular Fri Jun 10 06:10:23 2022 +0000 @@ -0,0 +1,2 @@ +file format type num_seqs sum_len min_len avg_len max_len Q1 Q2 Q3 sum_gap N50 Q20(%) Q30(%) +13-1941-6_S4_L001_R2_600000_fastq FASTQ DNA 100 10000 100 100.0 100 50.0 100.0 50.0 0 100 100.00 100.00
--- a/test-data/samtools_idxstats.tabular Mon Dec 06 18:29:21 2021 +0000 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,2 +0,0 @@ -NC_002945.4 4349904 45 4047 -* 0 0 5
--- a/test-data/samtools_idxstats1.tabular Mon Dec 06 18:29:21 2021 +0000 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,2 +0,0 @@ -NC_002945.4 4349904 17063 0 -* 0 0 223
--- a/test-data/samtools_idxstats2.tabular Mon Dec 06 18:29:21 2021 +0000 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,2 +0,0 @@ -NC_002945.4 4349904 45 4047 -* 0 0 5
--- a/test-data/samtools_idxstats3.tabular Mon Dec 06 18:29:21 2021 +0000 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,2 +0,0 @@ -NC_002945.4 4349904 24 0 -* 0 0 2
--- a/test-data/samtools_idxstats4.tabular Mon Dec 06 18:29:21 2021 +0000 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,2 +0,0 @@ -NC_002945.4 4349904 17063 0 -* 0 0 223
--- a/test-data/samtools_idxstats5.tabular Mon Dec 06 18:29:21 2021 +0000 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,2 +0,0 @@ -NC_002945.4 4349904 46 2 -* 0 0 4
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/statistics_output1.tabular Fri Jun 10 06:10:23 2022 +0000 @@ -0,0 +1,2 @@ +FASTQ File Size Read Count Length Sum Min Length Ave Length Max Length Q1 Q2 Q3 Sum Gap N50 Passing Q20 Passing Q30 Read Quality Ave +13-1941-6_S4_L001_R1_600000_fastq_gz 8.7 KB 100 10000 100 100.0 100 50.0 100.0 50.0 0 100 100.00 100.00 101.01
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/statistics_output2.tabular Fri Jun 10 06:10:23 2022 +0000 @@ -0,0 +1,2 @@ +R1 FASTQ R1 File Size R1 Read Count R1 Length Sum R1 Min Length R1 Ave Length R1 Max Length R1 Q1 R1 Q2 R1 Q3 R1 Sum Gap R1 N50 R1 Passing Q20 R1 Passing Q30 R1 Read Quality Ave R2 FASTQ R2 File Size R2 Read Count R2 Length Sum R2 Min Length R2 Ave Length R2 Max Length R2 Q1 R2 Q2 R2 Q3 R2 Sum Gap R2 N50 R2 Passing Q20 R2 Passing Q30 R2 Read Quality Ave +13-1941-6_S4_L001_R1_600000_fastq_gz 8.7 KB 100 10000 100 100.0 100 50.0 100.0 50.0 0 100 100.00 100.00 101.01 13-1941-6_S4_L001_R2_600000_fastq_gz 8.5 KB 100 10000 100 100.0 100 50.0 100.0 50.0 0 100 100.00 100.00 101.01
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/statistics_output3.tabular Fri Jun 10 06:10:23 2022 +0000 @@ -0,0 +1,2 @@ +R1 FASTQ R1 File Size R1 Read Count R1 Length Sum R1 Min Length R1 Ave Length R1 Max Length R1 Q1 R1 Q2 R1 Q3 R1 Sum Gap R1 N50 R1 Passing Q20 R1 Passing Q30 R1 Read Quality Ave R2 FASTQ R2 File Size R2 Read Count R2 Length Sum R2 Min Length R2 Ave Length R2 Max Length R2 Q1 R2 Q2 R2 Q3 R2 Sum Gap R2 N50 R2 Passing Q20 R2 Passing Q30 R2 Read Quality Ave +Unnamed Collection_R1 8.7 KB 100 10000 100 100.0 100 50.0 100.0 50.0 0 100 100.00 100.00 101.01 Unnamed Collection_R2 8.5 KB 100 10000 100 100.0 100 50.0 100.0 50.0 0 100 100.00 100.00 101.01
--- a/vsnp_add_zero_coverage.py Mon Dec 06 18:29:21 2021 +0000 +++ b/vsnp_add_zero_coverage.py Fri Jun 10 06:10:23 2022 +0000 @@ -33,15 +33,17 @@ def get_zero_df(reference): # Create a zero coverage dictionary. zero_dict = {} + reference_length = 0 for record in SeqIO.parse(reference, "fasta"): chrom = record.id total_len = len(record.seq) + reference_length = reference_length + len(record.seq) for pos in list(range(1, total_len + 1)): zero_dict["%s-%s" % (str(chrom), str(pos))] = 0 # Convert it to a data frame with depth_x # and depth_y columns - index is NaN. zero_df = pandas.DataFrame.from_dict(zero_dict, orient='index', columns=["depth"]) - return zero_df + return zero_df, reference_length def output_zc_vcf_file(base_file_name, vcf_file, zero_df, total_zero_coverage, output_vcf): @@ -82,25 +84,28 @@ return good_snp_count -def output_metrics_file(base_file_name, average_coverage, genome_coverage, good_snp_count, output_metrics): - bam_metrics = [base_file_name, "", "%4f" % average_coverage, genome_coverage] - vcf_metrics = [base_file_name, str(good_snp_count), "", ""] - metrics_columns = ["File", "Number of Good SNPs", "Average Coverage", "Genome Coverage"] +def output_metrics_file(base_file_name, average_coverage, genome_coverage, good_snp_count, dbkey, reference, + reference_length, total_zero_coverage, percent_ref_with_zero_coverage, output_metrics): + columns = ["BAM File", "Reference", "reference Length", "Genome with Coverage", "Average Depth", + "No Coverage Bases", "Percent Ref with Zero Coverage", "Quality SNPs"] + values = [base_file_name, dbkey, str(reference_length), genome_coverage, average_coverage, + str(total_zero_coverage), percent_ref_with_zero_coverage, str(good_snp_count)] with open(output_metrics, "w") as fh: - fh.write("# %s\n" % "\t".join(metrics_columns)) - fh.write("%s\n" % "\t".join(bam_metrics)) - fh.write("%s\n" % "\t".join(vcf_metrics)) + fh.write("# %s\n" % "\t".join(columns)) + fh.write("%s\n" % "\t".join(values)) -def output_files(vcf_file, total_zero_coverage, zero_df, output_vcf, average_coverage, genome_coverage, output_metrics): +def output_files(vcf_file, total_zero_coverage, percent_ref_with_zero_coverage, zero_df, output_vcf, + average_coverage, genome_coverage, output_metrics, reference, reference_length, dbkey): base_file_name = get_sample_name(vcf_file) good_snp_count = output_zc_vcf_file(base_file_name, vcf_file, zero_df, total_zero_coverage, output_vcf) - output_metrics_file(base_file_name, average_coverage, genome_coverage, good_snp_count, output_metrics) + output_metrics_file(base_file_name, average_coverage, genome_coverage, good_snp_count, dbkey, reference, + reference_length, total_zero_coverage, percent_ref_with_zero_coverage, output_metrics) -def get_coverage_and_snp_count(bam_file, vcf_file, reference, output_metrics, output_vcf): +def get_coverage_and_snp_count(bam_file, vcf_file, dbkey, reference, output_metrics, output_vcf): coverage_df = get_coverage_df(bam_file) - zero_df = get_zero_df(reference) + zero_df, reference_length = get_zero_df(reference) coverage_df = zero_df.merge(coverage_df, left_index=True, right_index=True, how='outer') # depth_x "0" column no longer needed. coverage_df = coverage_df.drop(columns=['depth_x']) @@ -109,19 +114,22 @@ coverage_df = coverage_df.fillna(0) coverage_df['depth'] = coverage_df['depth'].apply(int) total_length = len(coverage_df) - average_coverage = coverage_df['depth'].mean() + average_coverage = "{:.2f}".format(coverage_df['depth'].mean()) zero_df = coverage_df[coverage_df['depth'] == 0] total_zero_coverage = len(zero_df) + percent_ref_with_zero_coverage = "{:.6%}".format(total_zero_coverage / reference_length * 100) total_coverage = total_length - total_zero_coverage genome_coverage = "{:.2%}".format(total_coverage / total_length) # Output a zero-coverage vcf fil and the metrics file. - output_files(vcf_file, total_zero_coverage, zero_df, output_vcf, average_coverage, genome_coverage, output_metrics) + output_files(vcf_file, total_zero_coverage, percent_ref_with_zero_coverage, zero_df, output_vcf, + average_coverage, genome_coverage, output_metrics, reference, reference_length, dbkey) if __name__ == '__main__': parser = argparse.ArgumentParser() parser.add_argument('--bam_input', action='store', dest='bam_input', help='bam input file') + parser.add_argument('--dbkey', action='store', dest='dbkey', help='bam input dbkey') parser.add_argument('--output_metrics', action='store', dest='output_metrics', required=False, default=None, help='Output metrics text file') parser.add_argument('--output_vcf', action='store', dest='output_vcf', required=False, default=None, help='Output VCF file') parser.add_argument('--reference', action='store', dest='reference', help='Reference dataset') @@ -129,4 +137,4 @@ args = parser.parse_args() - get_coverage_and_snp_count(args.bam_input, args.vcf_input, args.reference, args.output_metrics, args.output_vcf) + get_coverage_and_snp_count(args.bam_input, args.vcf_input, args.dbkey, args.reference, args.output_metrics, args.output_vcf)
--- a/vsnp_build_tables.py Mon Dec 06 18:29:21 2021 +0000 +++ b/vsnp_build_tables.py Fri Jun 10 06:10:23 2022 +0000 @@ -1,6 +1,7 @@ #!/usr/bin/env python import argparse +import itertools import multiprocessing import os import queue @@ -153,11 +154,6 @@ return base_file_name -def output_cascade_table(cascade_order, mqdf, group, annotation_dict): - cascade_order_mq = pandas.concat([cascade_order, mqdf], join='inner') - output_table(cascade_order_mq, "cascade", group, annotation_dict) - - def output_excel(df, type_str, group, annotation_dict, count=None): # Output the temporary json file that # is used by the excel_formatter. @@ -183,19 +179,6 @@ excel_formatter(json_file_name, excel_file_name, group, annotation_dict) -def output_sort_table(cascade_order, mqdf, group, annotation_dict): - sort_df = cascade_order.T - sort_df['abs_value'] = sort_df.index - sort_df[['chrom', 'pos']] = sort_df['abs_value'].str.split(':', expand=True) - sort_df = sort_df.drop(['abs_value', 'chrom'], axis=1) - sort_df.pos = sort_df.pos.astype(int) - sort_df = sort_df.sort_values(by=['pos']) - sort_df = sort_df.drop(['pos'], axis=1) - sort_df = sort_df.T - sort_order_mq = pandas.concat([sort_df, mqdf], join='inner') - output_table(sort_order_mq, "sort", group, annotation_dict) - - def output_table(df, type_str, group, annotation_dict): if isinstance(group, str) and group.startswith("dataset"): # Inputs are single files, not collections, @@ -233,10 +216,6 @@ except queue.Empty: break newick_file, json_file, json_avg_mq_file = tup - avg_mq_series = pandas.read_json(json_avg_mq_file, typ='series', orient='split') - # Map quality to dataframe. - mqdf = avg_mq_series.to_frame(name='MQ') - mqdf = mqdf.T # Get the group. group = get_sample_name(newick_file) snps_df = pandas.read_json(json_file, orient='split') @@ -245,67 +224,53 @@ line = re.sub('[:,]', '\n', line) line = re.sub('[)(]', '', line) line = re.sub(r'[0-9].*\.[0-9].*\n', '', line) + line = re.sub("'", '', line) line = re.sub('root\n', '', line) sample_order = line.split('\n') - sample_order = list([_f for _f in sample_order if _f]) + sample_order = list(filter(None, sample_order)) sample_order.insert(0, 'root') - tree_order = snps_df.loc[sample_order] + tree_order_df = snps_df.loc[sample_order] + # Output the sorted table. + output_table(tree_order_df, "sort", group, annotation_dict) # Count number of SNPs in each column. snp_per_column = [] - for column_header in tree_order: + for column_header in tree_order_df: count = 0 - column = tree_order[column_header] + column = tree_order_df[column_header] for element in column: - if element != column[0]: + # column[0] is top row/root/reference, + # element is everything below it. + if element != column[0] and element != '-': count = count + 1 snp_per_column.append(count) - row1 = pandas.Series(snp_per_column, tree_order.columns, name="snp_per_column") + row1 = pandas.Series(snp_per_column, tree_order_df.columns, name="snp_per_column") # Count number of SNPS from the # top of each column in the table. snp_from_top = [] - for column_header in tree_order: + for column_header in tree_order_df: count = 0 - column = tree_order[column_header] - # for each element in the column - # skip the first element - for element in column[1:]: - if element == column[0]: - count = count + 1 - else: - break - snp_from_top.append(count) - row2 = pandas.Series(snp_from_top, tree_order.columns, name="snp_from_top") - tree_order = tree_order.append([row1]) - tree_order = tree_order.append([row2]) - # In pandas=0.18.1 even this does not work: - # abc = row1.to_frame() - # abc = abc.T --> tree_order.shape (5, 18), abc.shape (1, 18) - # tree_order.append(abc) - # Continue to get error: "*** ValueError: all the input arrays must have same number of dimensions" - tree_order = tree_order.T - tree_order = tree_order.sort_values(['snp_from_top', 'snp_per_column'], ascending=[True, False]) - tree_order = tree_order.T + column = tree_order_df[column_header] + index_list_of_ref_differences = [] + for ind, list_item in enumerate(column[1:].to_list()): + if list_item not in [column[0], '-']: + index_list_of_ref_differences.append(ind) + c = itertools.count() + val = max((list(g) for _, g in itertools.groupby(index_list_of_ref_differences, lambda x: x - next(c))), key=len) + # Starting row number with longest continous SNPs in column + snp_from_top.append(val[0]) + row2 = pandas.Series(snp_from_top, tree_order_df.columns, name="snp_from_top") + tree_order_df = tree_order_df.append([row1]) + tree_order_df = tree_order_df.append([row2]) + tree_order_df = tree_order_df.T + tree_order_df = tree_order_df.sort_values(['snp_from_top', 'snp_per_column'], ascending=[True, False]) + tree_order_df = tree_order_df.T # Remove snp_per_column and snp_from_top rows. - cascade_order = tree_order[:-2] + cascade_order_df = tree_order_df[:-2] # Output the cascade table. - output_cascade_table(cascade_order, mqdf, group, annotation_dict) - # Output the sorted table. - output_sort_table(cascade_order, mqdf, group, annotation_dict) + output_table(cascade_order_df, "cascade", group, annotation_dict) task_queue.task_done() -def set_num_cpus(num_files, processes): - num_cpus = len(os.sched_getaffinity(0)) - if num_files < num_cpus and num_files < processes: - return num_files - if num_cpus < processes: - half_cpus = int(num_cpus / 2) - if num_files < half_cpus: - return num_files - return half_cpus - return processes - - if __name__ == '__main__': parser = argparse.ArgumentParser() @@ -357,7 +322,6 @@ queue1 = multiprocessing.JoinableQueue() queue2 = multiprocessing.JoinableQueue() num_files = len(newick_files) - cpus = set_num_cpus(num_files, args.processes) # Set a timeout for get()s in the queue. timeout = 0.05 @@ -367,7 +331,7 @@ queue1.put((newick_file, json_file, json_avg_mq_file)) # Complete the preprocess_tables task. - processes = [multiprocessing.Process(target=preprocess_tables, args=(queue1, annotation_dict, timeout, )) for _ in range(cpus)] + processes = [multiprocessing.Process(target=preprocess_tables, args=(queue1, annotation_dict, timeout, )) for _ in range(args.processes)] for p in processes: p.start() for p in processes:
--- a/vsnp_get_snps.py Mon Dec 06 18:29:21 2021 +0000 +++ b/vsnp_get_snps.py Fri Jun 10 06:10:23 2022 +0000 @@ -21,18 +21,6 @@ return datetime.fromtimestamp(time.time()).strftime('%Y-%m-%d %H-%M-%S') -def set_num_cpus(num_files, processes): - num_cpus = len(os.sched_getaffinity(0)) - if num_files < num_cpus and num_files < processes: - return num_files - if num_cpus < processes: - half_cpus = int(num_cpus / 2) - if num_files < half_cpus: - return num_files - return half_cpus - return processes - - def setup_all_vcfs(vcf_files, vcf_dirs): # Create the all_vcfs directory and link # all input vcf files into it for processing. @@ -465,7 +453,6 @@ multiprocessing.set_start_method('spawn') queue1 = multiprocessing.JoinableQueue() num_files = len(vcf_files) - cpus = set_num_cpus(num_files, args.processes) # Set a timeout for get()s in the queue. timeout = 0.05 @@ -495,7 +482,7 @@ queue1.put(vcf_dir) # Complete the get_snps task. - processes = [multiprocessing.Process(target=snp_finder.get_snps, args=(queue1, timeout, )) for _ in range(cpus)] + processes = [multiprocessing.Process(target=snp_finder.get_snps, args=(queue1, timeout, )) for _ in range(args.processes)] for p in processes: p.start() for p in processes:
--- a/vsnp_statistics.py Mon Dec 06 18:29:21 2021 +0000 +++ b/vsnp_statistics.py Fri Jun 10 06:10:23 2022 +0000 @@ -1,25 +1,29 @@ #!/usr/bin/env python import argparse -import gzip import os -from functools import partial - -import numpy -import pandas -from Bio import SeqIO class Statistics: - def __init__(self, reference, fastq_file, file_size, total_reads, mean_read_length, mean_read_quality, reads_passing_q30): - self.reference = reference - self.fastq_file = fastq_file + def __init__(self, file_name, file_size, seq_type, num_seqs, sum_len, min_len, avg_len, + max_len, q1, q2, q3, sum_gap, n50, pass_q20, pass_q30, read_quality_average): + self.file_name = file_name self.file_size = file_size - self.total_reads = total_reads - self.mean_read_length = mean_read_length - self.mean_read_quality = mean_read_quality - self.reads_passing_q30 = reads_passing_q30 + self.seq_type = seq_type + self.num_seqs = num_seqs + self.sum_len = sum_len + self.min_len = min_len + self.avg_len = avg_len + self.max_len = max_len + self.q1 = q1 + self.q2 = q2 + self.q3 = q3 + self.sum_gap = sum_gap + self.n50 = n50 + self.pass_q20 = pass_q20 + self.pass_q30 = pass_q30 + self.read_quality_average = read_quality_average def nice_size(size): @@ -43,75 +47,19 @@ return '??? bytes' -def get_statistics(dbkey, fastq_file, gzipped): - sampling_size = 10000 - # Read fastq_file into a data fram to - # get the phred quality scores. - _open = partial(gzip.open, mode='rt') if gzipped else open - with _open(fastq_file) as fh: - identifiers = [] - seqs = [] - letter_annotations = [] - for seq_record in SeqIO.parse(fh, "fastq"): - identifiers.append(seq_record.id) - seqs.append(seq_record.seq) - letter_annotations.append(seq_record.letter_annotations["phred_quality"]) - # Convert lists to Pandas series. - s1 = pandas.Series(identifiers, name='id') - s2 = pandas.Series(seqs, name='seq') - # Gather Series into a data frame. - fastq_df = pandas.DataFrame(dict(id=s1, seq=s2)).set_index(['id']) - # Starting at row 3, keep every 4 row - # random sample specified number of rows. - file_size = nice_size(os.path.getsize(fastq_file)) - total_reads = len(seqs) - # Mean Read Length - if sampling_size > total_reads: - sampling_size = total_reads - try: - fastq_df = fastq_df.iloc[3::4].sample(sampling_size) - except ValueError: - fastq_df = fastq_df.iloc[3::4].sample(sampling_size, replace=True) - dict_mean = {} - list_length = [] - i = 0 - for id, seq, in fastq_df.iterrows(): - dict_mean[id] = numpy.mean(letter_annotations[i]) - list_length.append(len(seq.array[0])) - i += 1 - mean_read_length = '%.1f' % numpy.mean(list_length) - # Mean Read Quality - df_mean = pandas.DataFrame.from_dict(dict_mean, orient='index', columns=['ave']) - mean_read_quality = '%.1f' % df_mean['ave'].mean() - # Reads Passing Q30 - reads_gt_q30 = len(df_mean[df_mean['ave'] >= 30]) - reads_passing_q30 = '{:10.2f}'.format(reads_gt_q30 / sampling_size) - stats = Statistics(dbkey, os.path.basename(fastq_file), file_size, total_reads, mean_read_length, - mean_read_quality, reads_passing_q30) - return stats - - -def accrue_statistics(dbkey, read1, read2, gzipped): - read1_stats = get_statistics(dbkey, read1, gzipped) - if read2 is None: - read2_stats = None - else: - read2_stats = get_statistics(dbkey, read2, gzipped) - return read1_stats, read2_stats - - -def output_statistics(read1_stats, read2_stats, idxstats_file, metrics_file, output_file): +def output_statistics(read1_stats, read2_stats, output_file): paired_reads = read2_stats is not None if paired_reads: - columns = ['Read1 FASTQ', 'File Size', 'Reads', 'Mean Read Length', 'Mean Read Quality', - 'Reads Passing Q30', 'Read2 FASTQ', 'File Size', 'Reads', 'Mean Read Length', 'Mean Read Quality', - 'Reads Passing Q30', 'Total Reads', 'All Mapped Reads', 'Unmapped Reads', - 'Unmapped Reads Percentage of Total', 'Reference with Coverage', 'Average Depth of Coverage', - 'Good SNP Count', 'Reference'] + columns = ['R1 FASTQ', 'R1 File Size', 'R1 Read Count', 'R1 Length Sum', 'R1 Min Length', + 'R1 Ave Length', 'R1 Max Length', 'R1 Q1', 'R1 Q2', 'R1 Q3', 'R1 Sum Gap', + 'R1 N50', 'R1 Passing Q20', 'R1 Passing Q30', 'R1 Read Quality Ave', 'R2 FASTQ', + 'R2 File Size', 'R2 Read Count', 'R2 Length Sum', 'R2 Min Length', 'R2 Ave Length', + 'R2 Max Length', 'R2 Q1', 'R2 Q2', 'R2 Q3', 'R2 Sum Gap', 'R2 N50', 'R2 Passing Q20', + 'R2 Passing Q30', 'R2 Read Quality Ave'] else: - columns = ['FASTQ', 'File Size', 'Mean Read Length', 'Mean Read Quality', 'Reads Passing Q30', - 'Total Reads', 'All Mapped Reads', 'Unmapped Reads', 'Unmapped Reads Percentage of Total', - 'Reference with Coverage', 'Average Depth of Coverage', 'Good SNP Count', 'Reference'] + columns = ['FASTQ', 'File Size', 'Read Count', 'Length Sum', 'Min Length', 'Ave Length', + 'Max Length', 'Q1', 'Q2', 'Q3', 'Sum Gap', 'N50', 'Passing Q20', 'Passing Q30', + 'Read Quality Ave'] with open(output_file, "w") as outfh: # Make sure the header starts with a # so # MultiQC can properly handle the output. @@ -119,97 +67,102 @@ line_items = [] # Get the current stats and associated files. # Get and output the statistics. - line_items.append(read1_stats.fastq_file) + line_items.append(read1_stats.file_name) line_items.append(read1_stats.file_size) - if paired_reads: - line_items.append(read1_stats.total_reads) - line_items.append(read1_stats.mean_read_length) - line_items.append(read1_stats.mean_read_quality) - line_items.append(read1_stats.reads_passing_q30) - if paired_reads: - line_items.append(read2_stats.fastq_file) - line_items.append(read2_stats.file_size) - line_items.append(read2_stats.total_reads) - line_items.append(read2_stats.mean_read_length) - line_items.append(read2_stats.mean_read_quality) - line_items.append(read2_stats.reads_passing_q30) - # Total Reads + line_items.append(read1_stats.num_seqs) + line_items.append(read1_stats.sum_len) + line_items.append(read1_stats.min_len) + line_items.append(read1_stats.avg_len) + line_items.append(read1_stats.max_len) + line_items.append(read1_stats.q1) + line_items.append(read1_stats.q2) + line_items.append(read1_stats.q3) + line_items.append(read1_stats.sum_gap) + line_items.append(read1_stats.n50) + line_items.append(read1_stats.pass_q20) + line_items.append(read1_stats.pass_q30) + line_items.append(read1_stats.read_quality_average) if paired_reads: - total_reads = read1_stats.total_reads + read2_stats.total_reads - else: - total_reads = read1_stats.total_reads - line_items.append(total_reads) - # All Mapped Reads - all_mapped_reads, unmapped_reads = process_idxstats_file(idxstats_file) - line_items.append(all_mapped_reads) - line_items.append(unmapped_reads) - # Unmapped Reads Percentage of Total - if unmapped_reads > 0: - unmapped_reads_percentage = '{:10.2f}'.format(unmapped_reads / total_reads) - else: - unmapped_reads_percentage = 0 - line_items.append(unmapped_reads_percentage) - # Reference with Coverage - ref_with_coverage, avg_depth_of_coverage, good_snp_count = process_metrics_file(metrics_file) - line_items.append(ref_with_coverage) - line_items.append(avg_depth_of_coverage) - line_items.append(good_snp_count) - line_items.append(read1_stats.reference) + line_items.append(read2_stats.file_name) + line_items.append(read2_stats.file_size) + line_items.append(read2_stats.num_seqs) + line_items.append(read2_stats.sum_len) + line_items.append(read2_stats.min_len) + line_items.append(read2_stats.avg_len) + line_items.append(read2_stats.max_len) + line_items.append(read2_stats.q1) + line_items.append(read2_stats.q2) + line_items.append(read2_stats.q3) + line_items.append(read2_stats.sum_gap) + line_items.append(read2_stats.n50) + line_items.append(read2_stats.pass_q20) + line_items.append(read2_stats.pass_q30) + line_items.append(read2_stats.read_quality_average) outfh.write('%s\n' % '\t'.join(str(x) for x in line_items)) -def process_idxstats_file(idxstats_file): - all_mapped_reads = 0 - unmapped_reads = 0 - with open(idxstats_file, "r") as fh: +def get_statistics(fastq_file, seqkit_stats_file, seqkit_fx2tab_file): + file_size = nice_size(os.path.getsize(fastq_file)) + # SeqKit statistics. + with open(seqkit_stats_file, "r") as fh: + # This is a 2-line file for i, line in enumerate(fh): - line = line.rstrip('\r\n') - items = line.split("\t") if i == 0: - # NC_002945.4 4349904 213570 4047 - all_mapped_reads = int(items[2]) - elif i == 1: - # * 0 0 82774 - unmapped_reads = int(items[3]) - return all_mapped_reads, unmapped_reads - - -def process_metrics_file(metrics_file): - ref_with_coverage = '0%' - avg_depth_of_coverage = 0 - good_snp_count = 0 - with open(metrics_file, "r") as ifh: - for i, line in enumerate(ifh): - if i == 0: - # Skip comments. + # Skip header continue line = line.rstrip('\r\n') items = line.split("\t") - if i == 1: - # MarkDuplicates 10.338671 98.74% - ref_with_coverage = items[3] - avg_depth_of_coverage = items[2] - elif i == 2: - # VCFfilter 611 - good_snp_count = items[1] - return ref_with_coverage, avg_depth_of_coverage, good_snp_count + file_name = fastq_file + seq_type = items[2] + num_seqs = items[3] + sum_len = items[4] + min_len = items[5] + avg_len = items[6] + max_len = items[7] + q1 = items[8] + q2 = items[9] + q3 = items[10] + sum_gap = items[11] + n50 = items[12] + try: + pass_q20 = items[13] + except IndexError: + pass_q20 = 0 + try: + pass_q30 = items[14] + except IndexError: + pass_q30 = 0 + # Average read quality is not normalized on length. + avg_sum = 0 + with open(seqkit_fx2tab_file, "r") as fh: + for i, line in enumerate(fh): + if i == 0: + # Skip header + continue + line = line.rstrip('\r\n') + items = line.split("\t") + avg_sum += float(items[3]) + read_quality_average = "{:.2f}".format(avg_sum / float(i - 1)) + return Statistics(file_name, file_size, seq_type, num_seqs, sum_len, min_len, avg_len, + max_len, q1, q2, q3, sum_gap, n50, pass_q20, pass_q30, read_quality_average) parser = argparse.ArgumentParser() -parser.add_argument('--dbkey', action='store', dest='dbkey', help='Reference dbkey') -parser.add_argument('--gzipped', action='store_true', dest='gzipped', required=False, default=False, help='Input files are gzipped') parser.add_argument('--output', action='store', dest='output', help='Output Excel statistics file') parser.add_argument('--read1', action='store', dest='read1', help='Required: single read') parser.add_argument('--read2', action='store', dest='read2', required=False, default=None, help='Optional: paired read') -parser.add_argument('--samtools_idxstats', action='store', dest='samtools_idxstats', help='Output of samtools_idxstats') -parser.add_argument('--vsnp_azc_metrics', action='store', dest='vsnp_azc_metrics', help='Output of vsnp_add_zero_coverage') +parser.add_argument('--read1_seqkit_stats', action='store', dest='read1_seqkit_stats', help='Output of SeqKit statistics for forward read') +parser.add_argument('--read2_seqkit_stats', action='store', dest='read2_seqkit_stats', required=False, default=None, help='Output of SeqKit statistics for reverse read') +parser.add_argument('--read1_seqkit_fx2tab', action='store', dest='read1_seqkit_fx2tab', help='Output of SeqKit fx2tab for forward read') +parser.add_argument('--read2_seqkit_fx2tab', action='store', dest='read2_seqkit_fx2tab', required=False, default=None, help='Output of SeqKit fx2tab for reverse read') args = parser.parse_args() -stats_list = [] -idxstats_files = [] -metrics_files = [] -# Accumulate inputs. -read1_stats, read2_stats = accrue_statistics(args.dbkey, args.read1, args.read2, args.gzipped) -output_statistics(read1_stats, read2_stats, args.samtools_idxstats, args.vsnp_azc_metrics, args.output) +read1_stats = get_statistics(args.read1, args.read1_seqkit_stats, args.read1_seqkit_fx2tab) +if args.read2 is None: + read2_stats = None +else: + read2_stats = get_statistics(args.read2, args.read2_seqkit_stats, args.read2_seqkit_fx2tab) + +output_statistics(read1_stats, read2_stats, args.output)