search

amnona / q2-dbbact

A Qiime2 plugin for dbBact
MIT License
1 stars 1 forks source link

Plugin error from dbbact #1

Open anilchauhanhp9 opened 2 years ago

anilchauhanhp9 commented 2 years ago

Hi,

I ran q2-dbbact using qiime2-2021.4 and got the following error

Plugin error from dbbact: Duplicate observation IDs

The content of the log file is following

trimming primers if needed trim-primers converting hashes to sequences using repseqs file testing kmer head for region on 100 sequences: No region exact match detected. No sequence matched any of the regions. Testing 3973 sequences for embedded primers (3) Trimming primers before test to length 10 Trimmed primers are: {'C[AC]TGG[CT]TCAG': 'v1', 'G[ACGT][CGT]GC[AT][CG]CAG': 'v3', 'GCCGCGGTAA': 'v4'} Found maximal matching primer: G[ACGT][CGT]GC[AT][CG]CAG with enough matches (0.593003) found embedded primer G[ACGT][CGT]GC[AT][CG]CAG for region v3. removing it. trimming 3973 sequences with primer G[ACGT][CGT]GC[AT][CG]CAG trimmed 2357 sequences, 1616 not found Traceback (most recent call last): File "/home/anil/miniconda3/envs/qiime2-2021.4/lib/python3.8/site-packages/q2cli/commands.py", line 329, in call results = action(arguments) File "", line 2, in enrich_pipeline File "/home/anil/miniconda3/envs/qiime2-2021.4/lib/python3.8/site-packages/qiime2/sdk/action.py", line 244, in bound_callable outputs = self._callableexecutor(scope, callable_args, File "/home/anil/miniconda3/envs/qiime2-2021.4/lib/python3.8/site-packages/qiime2/sdk/action.py", line 484, in _callableexecutor outputs = self._callable(scope.ctx, view_args) File "/home/anil/miniconda3/envs/qiime2-2021.4/lib/python3.8/site-packages/q2_dbbact/methods.py", line 288, in enrich_pipeline table, = trim_func(table=table, repseqs=repseqs) File "", line 2, in trim_primers File "/home/anil/miniconda3/envs/qiime2-2021.4/lib/python3.8/site-packages/qiime2/sdk/action.py", line 244, in bound_callable outputs = self._callableexecutor(scope, callable_args, File "/home/anil/miniconda3/envs/qiime2-2021.4/lib/python3.8/site-packages/qiime2/sdk/action.py", line 390, in _callableexecutor output_views = self._callable(**view_args) File "/home/anil/miniconda3/envs/qiime2-2021.4/lib/python3.8/site-packages/q2_dbbact/methods.py", line 405, in trim_primers table.update_ids(trimmed_seqs_map, axis='observation', strict=False, inplace=True) File "/home/anil/miniconda3/envs/qiime2-2021.4/lib/python3.8/site-packages/biom/table.py", line 1316, in update_ids errcheck(result) File "/home/anil/miniconda3/envs/qiime2-2021.4/lib/python3.8/site-packages/biom/err.py", line 474, in errcheck raise ret biom.exception.TableException: Duplicate observation IDs

I do not know what the error is please help

amnona commented 2 years ago

Hi, Could you provide the exact command you used when the error occurred? Also, could you elaborate a bit more about the dataset you are trying to analyze? Is is a 16S rRNA amplicon dataset? What region/primers were used for the amplification (i.e V4 515f etc.)? How were the reads processed (deblur/dada2/another method)? Did you demultiplex based on barcodes? Looking at the error message, it seems q2-dbbact identified the reads as coming from the V3 region. It detected the V3 primers as part of the read, and tried to remove them. However, some sequences had different primers but the same sequence, which caused the error.

Cheers, Amnon

On Sun, May 15, 2022 at 6:02 PM anilchauhanhp9 @.***> wrote:

Hi,

I ran q2-dbbact using qiime2-2021.4 and got the following error

Plugin error from dbbact: Duplicate observation IDs

The content of the log file is following

trimming primers if needed trim-primers converting hashes to sequences using repseqs file testing kmer head for region on 100 sequences: No region exact match detected. No sequence matched any of the regions. Testing 3973 sequences for embedded primers (3) Trimming primers before test to length 10 Trimmed primers are: {'C[AC]TGG[CT]TCAG': 'v1', 'G[ACGT][CGT]GC[AT][CG]CAG': 'v3', 'GCCGCGGTAA': 'v4'} Found maximal matching primer: G[ACGT][CGT]GC[AT][CG]CAG with enough matches (0.593003) found embedded primer G[ACGT][CGT]GC[AT][CG]CAG for region v3. removing it. trimming 3973 sequences with primer G[ACGT][CGT]GC[AT][CG]CAG trimmed 2357 sequences, 1616 not found Traceback (most recent call last): File "/home/anil/miniconda3/envs/qiime2-2021.4/lib/python3.8/site-packages/q2cli/commands.py", line 329, in call results = action(arguments) File "", line 2, in enrich_pipeline File "/home/anil/miniconda3/envs/qiime2-2021.4/lib/python3.8/site-packages/qiime2/sdk/action.py", line 244, in bound_callable outputs = self.callable_executor(scope, callable_args, File "/home/anil/miniconda3/envs/qiime2-2021.4/lib/python3.8/site-packages/qiime2/sdk/action.py", line 484, in callable_executor outputs = self._callable(scope.ctx, view_args) File "/home/anil/miniconda3/envs/qiime2-2021.4/lib/python3.8/site-packages/q2_dbbact/methods.py", line 288, in enrich_pipeline table, = trim_func(table=table, repseqs=repseqs) File "", line 2, in trim_primers File "/home/anil/miniconda3/envs/qiime2-2021.4/lib/python3.8/site-packages/qiime2/sdk/action.py", line 244, in bound_callable outputs = self.callable_executor(scope, callable_args, File "/home/anil/miniconda3/envs/qiime2-2021.4/lib/python3.8/site-packages/qiime2/sdk/action.py", line 390, in callable_executor output_views = self._callable(**view_args) File "/home/anil/miniconda3/envs/qiime2-2021.4/lib/python3.8/site-packages/q2_dbbact/methods.py", line 405, in trim_primers table.update_ids(trimmed_seqs_map, axis='observation', strict=False, inplace=True) File "/home/anil/miniconda3/envs/qiime2-2021.4/lib/python3.8/site-packages/biom/table.py", line 1316, in update_ids errcheck(result) File "/home/anil/miniconda3/envs/qiime2-2021.4/lib/python3.8/site-packages/biom/err.py", line 474, in errcheck raise ret biom.exception.TableException: Duplicate observation IDs

— Reply to this email directly, view it on GitHub https://github.com/amnona/q2-dbbact/issues/1, or unsubscribe https://github.com/notifications/unsubscribe-auth/ABMQB4Q7ZJUBA7EXCNGOUU3VKEGXLANCNFSM5V7GVVFA . You are receiving this because you are subscribed to this thread.Message ID: @.***>

anilchauhanhp9 commented 2 years ago

Hi, Thanks for the reply... The command i ran was as follows "qiime dbbact enrich-pipeline --i-table table.qza --i-repseqs rep-seq.qza --m-metadata-file metadata.tsv --p-field deglaciation-phase --output-dir cfs-pipeline"

The dataset is 16S rRNA amplicons of glacier soil samples consisting of V3-V4 regions. The reads were demultiplexed, and processed by using DADA2.

Before running Qiime2 I already removed the sequencing primers using cutadapet tool.

amnona commented 2 years ago

It seems the cutadapt program did not fully remove the adaptor from all the sequences. q2-dbbact tries to remove the adaptor (as it detects it in some sequences) but does not handle this situation well and hence the error. Could you attach the first 1000 lines of the fastq file of one sample prior to and following cutadapt? What were the cutadapt parameters you used? Also, I think q2-dbbact will work well if you apply it to the reads without the cutadapt step (as it has its own cutadapt routine which is run by default).

Thanks Amnon

On Tue, May 17, 2022 at 7:15 AM anilchauhanhp9 @.***> wrote:

Hi, The command i ran was as follows "qiime dbbact enrich-pipeline --i-table table.qza --i-repseqs rep-seq.qza --m-metadata-file metadata.tsv --p-field deglaciation-phase --output-dir cfs-pipeline"

The dataset is 16S rRNA amplicons of glacier soil samples consisting of V3-V4 regions. The reads were demultiplexed, and processed by using DADA2.

Before running Qiime2 I already removed the sequencing primers using cutadapet tool.

— Reply to this email directly, view it on GitHub https://github.com/amnona/q2-dbbact/issues/1#issuecomment-1128390703, or unsubscribe https://github.com/notifications/unsubscribe-auth/ABMQB4RI2E76AIIYGKCDAD3VKMMMPANCNFSM5V7GVVFA . You are receiving this because you commented.Message ID: @.***>

amnona commented 2 years ago

@anilchauhanhp9 Created a new release trying to address this issue. Can you try to update your q2-dbbact plugin from pypi and let me know if it works with the cutadapt files you used? Also, if you can attach the first 1000 lines of a pre- and post- cutadapt file you used?

Thanks

anilchauhanhp9 commented 2 years ago

Hi, Sorry for the late reply. I'll do the same and inform you Thanks

anilchauhanhp9 commented 2 years ago

It seems the cutadapt program did not fully remove the adaptor from all the sequences. q2-dbbact tries to remove the adaptor (as it detects it in some sequences) but does not handle this situation well and hence the error. Could you attach the first 1000 lines of the fastq file of one sample prior to and following cutadapt? What were the cutadapt parameters you used? Also, I think q2-dbbact will work well if you apply it to the reads without the cutadapt step (as it has its own cutadapt routine which is run by default). Thanks Amnon … On Tue, May 17, 2022 at 7:15 AM anilchauhanhp9 @.> wrote: Hi, The command i ran was as follows "qiime dbbact enrich-pipeline --i-table table.qza --i-repseqs rep-seq.qza --m-metadata-file metadata.tsv --p-field deglaciation-phase --output-dir cfs-pipeline" The dataset is 16S rRNA amplicons of glacier soil samples consisting of V3-V4 regions. The reads were demultiplexed, and processed by using DADA2. Before running Qiime2 I already removed the sequencing primers using cutadapet tool. — Reply to this email directly, view it on GitHub <#1 (comment)>, or unsubscribe https://github.com/notifications/unsubscribe-auth/ABMQB4RI2E76AIIYGKCDAD3VKMMMPANCNFSM5V7GVVFA . You are receiving this because you commented.Message ID: @.>

Hi, the fastq file lines before running cutadapt is as follows

@M05193:157:000000000-JHBCT:1:1101:19441:1844 1:N:0:GAATTCGT+AGGCGAAG TCAATAACTGTGTAGTAAGCTGCCTTCTCAATTGGTGTTCTATGTCATATCTAAGCATTTCACCGCTACATGACATATTCCGCTTACCTCCACAATATTCAAGACAAATAGTATCCATGGCAGTTTCTGAGTTAAGCTCAGAGATTTCACCACGGACTTACCTGCCCGCCTACGCACCCTTTAAACCCAGTGAATCCGGATAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCG + AAA1AB1FFBFFGCFGFB1FGCCFGGC1B3EFG11DFHFFHCGFHFBGFHHHFBBFBGHHHBHGCEGGGHHHFFFGEGHHHGGGGCGHGHHHHHFHEGHHFFHGHHEFHHHHHHHHHHHBHHFHHHHHFHHHHHHHHHHHGHHHHHHHHHHHHGGGGGHFHHHFHHHGGGGGHGGCGCGHHHHHHHHGGHHHHHHHHGGGGGGHGGGGGHHHHGHHGGEFFGEGGB@9@@AAFFF/9FF?-=9AFF/BBA @M05193:157:000000000-JHBCT:1:1101:16906:1883 1:N:0:GAATTCGT+AGGCGAAG GACTACGAGGGTATCTAATCCCATTTGCTCCCCTAGCTTTCGTCTCTCAGTGTCAGTAATGGCCCAGTAGAGTGCCTTCGCCATCGGTGGTCTTTCCAATATCTACGCATTTCACCGCTCCACTGGAAATTCCCTCTACCCCTACCATACTCTAGTTTGAAAGTTTCCACGGCATCTTTGAAGTTAAGCCTCAATTTTTAACAGTGGACTTTTCGAACCACCTACAGACGCTTTACGCCCAGTGATTCCG + AAAAAFBAABBBGGGGF5GGGGCGGHBFHHHGGHGHHHHHHGHGHHHGFFGFHHHHHHHHGHGHEEHHHHHHHHHHHHHHDEEGHEDGGGGHHHHHHHHHHHHHHHGGGGGHHHHHGGGGGHHHHHHHHHHHHHHHHHHHHHGGHHHHHHHHHHHHHHHHGHHHHHHHHHGGGGGHHHHHHHHHHHHHHHHHHHHHHHHHHGHHHHHHHHHHHHGGHF@FGGGGGGGGGGGGGGGF9DDCFEFFFFFFFB @M05193:157:000000000-JHBCT:1:1101:17994:1899 1:N:0:GAATTCGT+AGGCGAAG GACTACGGGGGTATCTAATCCTGTTTGCTCCCCACGCTTTCGCTCCTCAGCGTCAGTTACTGCCCAGAGACCCGCCTTCGCCACCGGTGTTCCTCCTGATATCTGCGCATTTCACCGCTACACCAGGAATTCCAGTCTCCCCTACAGCACTCTAGTCTGCCCGTATCGACTGCAGGCTCGGGGTTGAGCCCCGAGTTTTCACAGCCGACGCGACAAACCGCCTACGAGCTCTTTACGCCCAATAATTCCG + BBBBBFBBBBBBGGGGGGGGGGHHHHFHHHHGGHGGGGGGHEGGHGHHHFGEGGGGHHHHHHHHHHGFHHHGHGGGGGHEGGGGHGGGGGGHHHHHHHHHHHHHHHGGGGGHHHHHGGGGGHHHHHGHHHHHHHHHHHHHHGHGHHHHHHHHHHHHHHHHHHGGGHHGGGGGGGGGGGGFGGGFFFFFFFFFFFFEFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF?FFFF?FFFFFFFB @M05193:157:000000000-JHBCT:1:1101:15487:1918 1:N:0:GAATTCGT+AGGCGAAG AGTGGGGAATATTGCACAATGGGCGCAAGCCTGATGCAGCCATGCCGCGTGTATGAAGAAGGCCTTCGGGTTGTAAAGTACTTTCAGCGGGGAGGAAGGCGATAAGGTTAATAACCTTGTCGATTGACGTTACCCGCAGAAGAAGCACCGGCTAACTCCGTGCCAGCAGCCGCGGTAATACGGAGGGTGCAAGCGTTGTCCGGATTTATTGGGTTTAAAGGGTGCGTAGGCGGCTGATTAAGTCTGGGGT +

AADAAFFFDGGGGFGGGGHF2E2AEHGH2FGHHFGGHHHHHHAEEGGGHHHHHDGHHGHHHHHHGGHGGGHHGHHHHHHHHHHGGGGCGGGGHGGHGGGGGHHHFHHHHHHHHHHHHGHGHGFHGGGGGHHGGGGGFHHHHHHHHCGGGBGGGGGGGGGGGGGGGGGGGAADFFFFFFF.AFFAEF/;FFFFFFFFFFF.DFFFFFFFFFFFFFFFFFFFFFFFADFFF=;F.BF/B/FFFEFE-D @M05193:157:000000000-JHBCT:1:1101:17224:1943 1:N:0:GAATTCGT+AGGCGAAG CCTACGGGGGGCAGCAGTAAGGAATATTGGTCAATGGACGCAAGTCTGAACCAGCCATGCCGCGTGAAGGATTAAGGCCCTCTGGGTTGTAAACTTCTTTTATCTGGGACGAAAAACGGGTATTCATACTCACTTGACGGTACCAGAGGAATAAGCACCGGCTAACTCCGTGCCAGCAGCCGCGGTAATACGGAGGGTGCAAGCGTTATCCGGATTTACTGGGTTTAAAGGGTGCGTAGGTGGGTCTGTA + CCCCCCCCCCCCGGGGGGGGGGHHHHHHHGHHHHHHHHGGGGGGHHHHHGHHHHHGHHHHHGGGGGGGHHHHHHHGHHHGHHHHHHHFGGHHHHHHHHHHHHHHHHGGGGGGGGGGGGGGDGGGGFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFEFFFFFFFFEFFFFFFFFFFFFFFFF @M05193:157:000000000-JHBCT:1:1101:11649:1945 1:N:0:GAATTCGT+AGGCGAAG CCTACGGGTGGCAGCAGTGGGGAATATTGGACAATGGGGGCAACCCTGATCCAGCCATGCCGCGTGGGTGAAGAAGGCCTTCGGGTTGTAAAGCCCTTTTGTACGGGAAGAAAAGCTGCGGGTGAACAGCCTGTAGTCATGACGGTACCGTAAGAATAAGCACCGGCTAACTTCGTGCCAGCAGCCGCGGTAATACGAAGGGTGCAAGCGTTACTCGGAATTACTGGGCGTAAAGCGTGCGTAGGCGGTG + ABBBBBBBBD@FGGGGCGGBECGCGFGHDDEGHHHFFHA0AEFGHGHGFGHHHFHGHHFHHGGCGEGGGFHGFHHFFHHHHHEFCHG?EHHGHFHHHHHHFGHHH??DFDHBGHFGH1GD@CDGGGHFHHHHHHHHHHHHFFGG?EEFHGGHGHHFFFHGGGGGGGGGGFGGGGGGGGGGGGFFFFFFFFFFEFBFFEEFFFFFFEFFFFFFAEFF?DFEEFEFFFFEFFFFFABF?F?D-ED-.BD-9- @M05193:157:000000000-JHBCT:1:1101:12689:1962 1:N:0:GAATTCGT+AGGCGAAG GACTACAGGGGTATCTAATCCTGTTCGCTCCCCACGCTTTCGCTCCTCAGCGTCAGTTACTGCCCAGAGACCCGCCTTCGCCACCGGTGTTCCTCCTGATATCTGCGCATTTCACCGCTACACCAGGAATTCCAGTCTCCCCTGCAGTACTCAAGTCAGCCCGTATCGACTGCAAGCCCACCGTTGAGCGCTGGGTTTTCACAGCCGACGTGACAAACCGCCTACGAGCTCTTTACGCCCAGTAATTCCG + BBBBBFFFFBBBGFGGGBGGGGHHHHGGGHGGGHEGGGGGHGGGHGHHHHHGGGGGHGFFHHHHFHHAHHHFHGGGFGHGGGGGHGGGGGGFHHHHHHHHHHFGHHGGGGGGHGHHGGGCGFGHGFG/FHHHHFHGGHHHEGHGHHGHHGGHHHHHGHGHHGFGGGHGGGGDHHEH.EHHGGGFFGGGGGGGGGBFGGGGGFFFE?B;FFFFFFF/FFAD;BFFEDDDFDFFFFFFFFC;.EFBFBFFFB @M05193:157:000000000-JHBCT:1:1101:16633:1968 1:N:0:GAATTCGT+AGGCGAAG CCTACGGGTGGCAGCAGTGGGGAATATTGCGCAATGGGCGAAAGCCTGACGCAGCGACGCCGCGTGAGGGACGAAGGCCTTCGGGTTGTAAACCTCTTTCAGCAGGGACGAAGCGAAAGTGACGGTACCTGCAGAAGAAGCACCGGCCAACTACGTGCCAGCAGCCGCGGTAATACGTAGGGTGCAAGCGTTGTCCGGAATTATTGGGCGTAAAGGGTGCGTAGGCGGTGCGGTAAGTGTGCTGTGAAAT + CCCCCBCCDBCFGGGGGGGGGGGGHHHGHHGGGGGHHGHGGGGGHHHHHHGGGGGGGGGGGGGGGGGGHGGGGGGGGHHHHHHGGHGGGHHHHGHHHHHHHHHHHHGGGGGGGGGGGGGGHHHGGGHGHHGGGGGGGGGGGGGGGGGGGGGGFFFFFFFFFFFFFFFFFFDFFFFFFFFFFFFEFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFEFFFFFFBCFFF?FEFFFFFFFFFFFFFF @M05193:157:000000000-JHBCT:1:1101:15694:1989 1:N:0:GAATTCGT+AGGCGAAG CCTACGGGTGGCAGCAGTGGGGAATATTGCGCAATGGGCGAAAGCCTGACGCAGCGACGCCGCGTGAGGGATGAAGGCCTTCGGGTTGTAAACCTCTTTCAGCAGGGACGAAAACAGACGGTACCTGCAGAAGAAGCACCGGCCAACTACGTGCCAGCAGCCGCGGTAATACGTAGGGTGCAAGCGTTGTCCGGAATTATTGGGCGTAAAGAGCTCGTAGGCGGTCTGTCACGTCGGCTGTGAAAACTCG + BBBBBBB?ADBFGGEGGGGGGGCGHGHHHHGGGGGHHHHGGGGGHHGHHHGGGGGGGGGGGGGGGGGGHGGHGHGHHHHHGHHGGHGGGGHHHHHHHHHHHHHHHHGGGGGGGGHHHGHGGDHGGHHHGGBGHHHHHHHHGGGGGGGGGGGFFGGGGGGGGGFAFDAAEFFFF9BFFF/ADFFFFFFFFFFFFF?=EFFFFFFFFFFAF?F/B.FFF@DEFFFFF;DFBFFFBFEFFFEFFFFFFFFFFF @M05193:157:000000000-JHBCT:1:1101:12987:2015 1:N:0:GAATTCGT+AGGCGAAG GCGGGTATCTAATCCTGTTTGATCCCCACGCTTTCGTGCCTCAGCGTCAATATATGCGTAGTAAGCTGCCTTCGCAATAGGTGTTCTATGTCATATCTAAGCATTTCACCGCTACATGACATATTCCGCTTACCTCCACAATATTCAAGACAGATAGTATCAATGGCAGTTCTGGAGTTGAGCTCCAGGATTTCACCACTGACTTACCTGCCCGCCTACGCACCCTTTAAACCCAGTGAATCCGGACAAC + ABCCBDBFFFFFGGGGGGGGGGHHHHGHGGGGGGHGGGGHHHHHHGGGGGGHHHHHHGGGGGFFGHHHHHHHHHGGGGHHHHHHGHHHHHHHHHHHHHHHHHHHHHHHHHGGGGGHHHGHHHHHHHHHGGGGHHHHHHHHHHHHHHHHHHHHHGHHHHHHHHHHHHHHGHHHHHHFHHHHHHHHHHHGHHHHHHHHHHHHHHHHHHHHHHHHHFGGGGGGGGGGGGGGGFGGGGGGGGFGGGGFFFFFFF @M05193:157:000000000-JHBCT:1:1101:15270:2040 1:N:0:GAATTCGT+AGGCGAAG GACTACGGGGGTATCTAATCCTGTTTGCTCCCCACGCTTTCGTGCATGAGCGTCAGTGTTATCCCAGGGGGCTGCCTTCGCCATCGGTATTCCTCCACATATCTACGCATTTCACTGCTACACGTGGAATTCTACCCCCCTCTGACACACTCTAGCCATGCAGTCACAAATGCAGTTCCCAGGTTAAGCCCGGGCATTTCACATCTGTCTTACATAACCGCCTGCGCACGCTTTACGCCCAGTAATTCCG + AAAA1CAADDDDGGGDDBEGGB1BGGBFHFHGGHEGGECEDEEFF/FFF1FEE>EE22@F22GEF1@/EE?EECAC1GHF?E//F?EEHEDFDGGHGHFGHHH22FF?ACGFGH2FBGFFDGFCCGEGHGGGDGGH0G-@-CFCFHHFHGHHHFGHFHFGFFGHCHH0FFFBFFGGGGGF0FGF0C009FGGG-@@-BFB99B/BF//BFFFFFFBFF/-@@BFF@@@??@A-EBFEF@@FF/FFF/FFF @M05193:157:000000000-JHBCT:1:1101:19202:2062 1:N:0:GAATTCGT+AGGCGAAG CCTACGGGTGGCAGCAGTAAGGAATATTGGTCAATGGACGCAAGTCTGAACCAGCCATGCCGCGTGGAGGATTAAGGCCCTCTGGGTTGTAAACTCCTTTTATCTGGGAAGAAATGCATTTTTTCTAAAGTGCTTGACGGTACCAGAGGAATAAGCACCGGCTAACTCCGTGCCAGCAGCCGCGGTAATACGGAGGGTGCAAGCGTTATCCGGATTCACTGGGTTTAAAGGGTGCGTAGGCGGGCAGGTA + ABBBBBBBDBBFGGGGGGGGGGHHHGHHHFHHHGHHHHHGGGCGHHHHGFHHHHHGHHFHHGGGGGGGGEGGHHHHGHHGHHHHHEFF@EHHHHHHHHHHHHHHHHGGHHGHHHGHGHHHHHHGHHGFHHHHGHHGFHGGGHGHGHHEEHHFCFDHHHGCGGGGHHHHHGGEGHH0EHHHE?CAADDFGGG?EGGCFF.D9FFFF?DDFFFBF.AAFFFFFEEFFDFBB/FBFAAFFFADFF-=-=DAAB @M05193:157:000000000-JHBCT:1:1101:19663:2066 1:N:0:GAATTCGT+AGGCGAAG GACTACGCGGGTATCTAATCCTGTTTGCTCCCCACGCTTTCGCACCTCAGCGTCAGTATTGGGCCAGTGAGCCGCCTTCGCCACTGGTGTTCTTGCGAATATCTACGAATTTCACCTCTACACTCGCAGTTCCACTCACCTCTCCCAAACTCCAGATTCCCAGTATCAAAGGCAGTTCCAAGGTTGAGCCTTGGGATTTCACCTCTGACTTAAAAATCCGCCTACGTGCGCTTTACGCCCAGTAATTCCG + ABBBBFBDBBBBGGGGGGGGGGGHHHHHHGHGGHGGGGGGHGG1AFHGGHHEEG1EGHHFHGHHHGHFBGFFHGGGGGHEFGGGHHFHHHGHHH3F//>CGHHFHHAEHHHHFHHHHHHHHHG1?F@D?AGHHHHHGFGHHHHHHHE/GGHF0CGDF000DGHHEGHHHH0;GHDHHHHF0CEGG/CFGGGFGEFGGGGGGGGGGGGGGGGG000C9FFFFFDDF.;FF-9@FFFFFFFFF.FFFFFFFF @M05193:157:000000000-JHBCT:1:1101:16820:2070 1:N:0:GAATTCGT+AGGCGAAG GACTACCGGGGTATCTAATCCTGTTTGCTCCCCACGCTTTCGTGCATGAGCGTCAGTCTTGACCCAGGGGGCTGCCTTCGCCATCGGTATTCCTCCACATATCTACGCATTTCACTGCTACACGTGGAATTCTACCCCCCTCTGCCAAACTCTAGCCATGCAGTCACAAATGCAGTTCCCAGGTTGAGCCCGGGCATTTCACATCTGTCTTACATAACCGCCTGCGCACGCTTTACGCCCAGTAATTCCG + AAA?>CCDDDD>EGFGGBFAGFHHHHHHFHFFECFEEGGGHGFGGC2FGHHGEGGGHH1EBFGHHFGFCEGCGGFAHHHHCEFGFGGGHGGHFHHCGF<1GGBDDFFGCGGBFHFEDGHGBGB?<<FGHHGHFHHH0?CGCCHFHHHFH<GHHHHEHHGH/GGHHFHHHHFHHFHHFHFCGCEFGGGGGF?AGGA@GGFGGGGGGGGGFFFFFFFFFFF@@?FFF@@@@<@FEFFFF?@@BEFFFFFFFF @M05193:157:000000000-JHBCT:1:1101:11504:2121 1:N:0:GAATTCGT+AGGCGAAG CCTACGGGGGGCAGCAGTGGGGAATATTGGACAATGGGCGAAAGCCTGATCCAGCAATGCCGCGTGAGTGATGAAGGCCTTAGGGTTGTAAAGCTCTTTTACCCGGGAAGATAATGACTGTACCGGGAGAATAAGCCCCGGCTAACTCCGTGCCAGCAGCCGCGGTAATACGGAGGGGGCTAGCGTTGTTCGGAATTACTGGGCGTAAAGCGCGCGTAGGCGGCTTTGTAAGTTAGAGGTGAAAGCCCGG + CCCCCCBCCCCCGGGGGGGGGGGGHHHHHHHHHHGHHHHGGGGGHHHHHHHHHHHHHHHHHGGGGGCGHHHHHHHHHHHHHHHGGHHGGHHHHHHHHHGGGGGGGGGGGGGGGGGGGGGGGGGGFFFCFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFBBFFFFBBFFFFFFFBFFFFF @M05193:157:000000000-JHBCT:1:1101:18681:2125 1:N:0:GAATTCGT+AGGCGAAG CCTACGGGTGGCAGCAGTGGGGAATTTTGGACAATGGGGGCAACCCTGATCCAGCAATGCCGCGTGAGTGAAGAAGGCCTTCGGGTTGTAAAGCTCTTTTGTCAGGGAAGAAACGGTTTGGGTTAATACCTTGAACTAATGACGGTACCTGAAGAATAAGCACCGGCTAACTACGTGCCAGCAGCCGCGGTAATACGTAGGGTGCAAGCGTTAATCGGAATTACTGGGCGTAAAGCGTGCGCAGGCGGTT + BBBBBBBBDADFGGGGGGGGGGGDHHHHHGHHHHHHGH2EFGGHHGHHHHHHHHHFHHHHHGGGGGGGHHHHHHHHHHHHHHHGGHGGGGEHHGHHHHHH2BGHHHGGHGHHHFGFDEEDFFGGHHHHHHHHGFHHFB1FHHHGCDHGHHHCGGHHHHHHHHHHGGGGGHHHGGGGGGGGGGGGGGGGGGFBFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFBFFEFFFFFFFFFFFFFFFFFFF-@ @M05193:157:000000000-JHBCT:1:1101:19265:2137 1:N:0:GAATTCGT+AGGCGAAG GACTACCGGGGTATCTAATCCTGTTTGCTCCCCACGCTTTCGCGCCTCAGCGTCAGTAACGGGCCAGTGAGTCGCCTTCGCCACTGGTGTTCTTCCGAATATCTACGAATTTCACCTCTACACTCGGAATTCCACTCACCTCTCCCGTACTCAAGATGTCCAGTATCAAAGGCAGTTCCAAGGTTGAGCCCTGGGCTTTCACCTCTGACTTAAACATCCGCCTACGCGCCCTTTACGCCCAGTAATTCCG + BCCCCFFCCCCCGGGGGGGGGGHHHHHHHHHGGHGGGGGGHGGGGGGGHHHGGGGGHHHHGGCGGGHHHHHHHGGGGGHHGGGGHHHHHHGHHHHHGEGGGHHHHHGGHHHHHHHHHHHHHHHHHGGGGGHHHHHHHHHHHHHHHHGHGGHHHHHGHHGHHHHHGHHHBGGHHHGHHHHHHHGHHHGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGBFFFFFFFFFFFFFFFFEFFFFFFFFFFFFF @M05193:157:000000000-JHBCT:1:1101:14288:2162 1:N:0:GAATTCGT+AGGCGAAG CCTACGGGGTGCAGCAGTTGGGAATCTTGCACAATGGGGGAAACCCTGATGCAGCGACGCCGCGTGGACGATGAAGGCCTTAGGGTTGTAAAGTCCTTTTTTGGGGAAAGACTTAGGACGGTACCCCAAGAATAAGCTCCGGCTAACTACGTGCCAGCAGCCGCGGTAAGACGTAGGGAGCAAGCGTTGTCCGGATTTACTGGGCGTAAAGAGCTCGTAGGCGGTTCGGTAAGTCGGGTGTGAAAACTCT + CABCCCCCCDDCGGGGGGGGGGHHHHHHHHHHHHHHHHGGGGHHHGHHHHHHHHHGGGGGGGGGGGGGHGGGGHHHHHHHHHHHHHHGGHHHHHHHHHHHHGGGGGGGHHHHHGHHHHHGGGHGHHHGGHHHHHHHHHHHGGGGGHHHHHGGGGGGGGGGGGGGGGFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFDFFFFFFFFFFCFFFFFFBCFFFAFFFFFFFFEDDFFFFFFFFF @M05193:157:000000000-JHBCT:1:1101:20376:2173 1:N:0:GAATTCGT+AGGCGAAG GACTACGCGGGTATCTAATCCTGTTTGATCCCCACGCTTTCGCTCCTCAGCGTCAGTTACTGCCCAGAGACCCGCCTTCGCCACCGGTGTTCCTCCTGATATCTGCGCATTTCACCGCTACACCAGGAATTCCAGTCTCCCCTACAGCACTCTAGTCTGCCCGTATCGACTGCAGGCTCGGGGTTGAGCCCCGAGTTTTCACAGCCGACGCGACAAACCGCCTACGAGCTCTTTACGCCCAATAATTCCG + BCCDCFBCCCCBGGGGGGGGGGHHHHHHHHHHGHGGGGGGHGGGHGHHHHHGGGGGHHHHHHHHHHHHHHHGHGGGGGHHGGGGHGGGGGGHHHHHHHHHHHHHHHGGGGGHHHHHGGGGGHHHHHGHHHHHHHHHHHHHHGHGHHHHHHHHHHHHHHHHHHGGGHHGHGGGGGGGGGGGGGGGGFFFFFFFFFFEFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFHHHFFFFFFFFFFFFFF @M05193:157:000000000-JHBCT:1:1101:11103:2183 1:N:0:GAATTCGT+AGGCGAAG GACTACACGGGTATCTAATCCTGTTTGCTCCCCACGCTTTCGCGTCTCAGCGTCAGTTACCGTCCAGATGGCCGCCTTCGCCACCGGTGTTCCTCTGCATATCTACGCATTTCACTGCTACACGCGGAATTCCATCCCCCTCTACCGTACTCTAGCTATACAGTCACAGATGCAATTCCCAGGTTGAGCCCGGGGATTTCACAACTGTCTTATATAACCGCCTGCGCACGCTTTACGCCCAGTAATTCCG + ABABAFFFBDBDGGGGGGGGGGHHHHHHHHHGGHGGGGGGFGDEGGEHHHHGGGGDGHHHHGGGGFFHHHHFHGGGGGHHGGGGHGGGGGGHHHHHHGHHHHHHHHGGGGGHHHHHHHHHGFGGGGGGGGHHHHHHFHGDGGHHHHGGHGHHHHHHGHHHHHHHHFHHFHHHHGGGGGFGGEGCFFGGGGGGGGGFFFFFFFFFFFFFFBFFFFFFFFFFFFFEFD;BFFFFFEFFFFAFFFFFFEFFFF @M05193:157:000000000-JHBCT:1:1101:14294:2193 1:N:0:GAATTCGT+AGGCGAAG CCTACGGGGGGCAGCAGTGGGGAATATTGCACAATGGGCGGAAGCCTGATGCAGCGACGCCGCGTGAGGGATGAAGGCCTTCGGGTTGTAAACCTCTTTCAGCAGGGAAGAAGCGGTCGGGTTTTCTTGGTCGTGACGGTACCTGCAGAAGAAGCACCGGCTAACTACGTGCCAGCAGCCGCGGTAATACGTAGGGTGCGAGCGTTGTCCGGAATTATTGGGCGTAAAGAGCTTGTAGGCGGTTTGTCGC + BAACCCBCCCCCGGGGGGGGGGGGHHHHHHGHHHHHHHHGGGGGHHHHHHHHHHHGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGFFFFFFFFFFFFFFFFFFFFFFDFFFF9EFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFADFFFFFFFFFFBFFFFFACADFFFFFFABFFFFFFFFFFBFFFFFFFFFFFFFFFFF-@DFFFFF @M05193:157:000000000-JHBCT:1:1101:18285:2208 1:N:0:GAATTCGT+AGGCGAAG GACTACGCGGGTATCTAATCCCGTTTGCTCCCCTAGCTTTCGTCTCTCAGTGTCAGTCGTGGCCCAGCAGAGTGCCTTCGCCATTGGTGTTCTTCCCGATCTCTACGCATTTCACCGCTACACCGGGAATTCCCTCTGCCCCTACCAGACTCCAGCTTTTAAGTTTCTTCTGCCGACCCAGGGTTAAGCCCTGATCTTTGACAGAAGACTTTAAAAACCACCTACAGACGCTTTACGCCCAATCATTCCG + BBBBBFA3DBBBGGGGGGGGGFGHHHGHHHHGGFEHHHBHHGFEGHHHHHHHHHHHFAEF1EFHHGHHHHHE@FHHHHHHGGGEGHHH/GGHHHHHHGGGGHHHHHGGGGGFHHHHGGGGGHHHCGGGGHHHHGHHHHHHHGHGHFGGCEHHGHHHHHHHGHHHFCFGFG0GHFGGCGGEGGGFGEGGGGGGG/BFGGGFGGFFFFF9FFFFFBFFFFFFFFFFFFBFFAFFFFFFFFFFFEFFFFFFFF @M05193:157:000000000-JHBCT:1:1101:17873:2233 1:N:0:GAATTCGT+AGGCGAAG GACTACGAGGGTATCTAATCCTGTTTGCTCCCCACGCTGTCGCGCCTCAGTGTCAGTAACTGCCCAGCAGGCCGCCTTCGCCACCGGTCTTCCTCCGGATCTCTACGCATTCCACCGCTACACCCGGAATTCCGCCTGCCCCTACAGTACTCTAGCCTAACAGTCGGCACGGCCGACTCGGGGTTGAGCCCCGAGCTTTCACCGCACGCTTAAGCCACCACCTACGCGCCCTTTACGCCCAGTGATTCCG + BCCCCFCCCCCCGGGGGGGGGGHHHHHHHHHGGHGGGGGGHGHGGGGGHHHHHHHHHHHHHHHHHHHHHHHGGGGGGGHHGGGGHGGGGGHHHHHHGGGGGHHHHHGGGGGHHHHHGGGGGHHHGGGGGGHHHHGGGHHHHGHGHHHHHHHHHHHHHHHHHHHHHGCGGGGGGGGGGGGGFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF @M05193:157:000000000-JHBCT:1:1101:19778:2274 1:N:0:GAATTCGT+AGGCGAAG CAAAGGCAGTTCCGGGGTTGAGCCCCGGGCTTTCACCTCTGACTTGAATAGCCGCCTACGTGCGCTTTACGCCCAGTAATTCCGATTAACGCTCGCACCCTACGTATTACCGCGGCTGCTGGCAAGATCGGAAGAGCACACGTCTGAACTCCAGTCACGAATTCGTATCTCGTATGCCGTCTTCTGCTTGAAAAAAAAAAAAAAGCGAAAAACGCTCTTAGGCGGCCCTGTAACTCAAATTTGAAATCCC + BBBB@FBFBFFFGGGGGGGGGGHHGGG2EEGHHHGHHHHHHHHGHFHHHGHHHGGGGGH1GGHGGGGGGHGGGGGHHHGHHHHHGGGGHHEEFGEGFGGGHHHGHHHHHGGFGGGGGGGHHHGHHGHHHHGGGGHGHHHHHGHHHGHHHHHGGGFGGGGGDGGGEBFFFFFGADFFGGGDFAEFF/BFFF;/BFFF=B-9----....---..-;...9/9//-----.99/9///;//:////9//BB/ @M05193:157:000000000-JHBCT:1:1101:18722:2299 1:N:0:GAATTCGT+AGGCGAAG GACTACGCGGGTATCTAATCCTGTTTGCTCCCCACGCTTTCGTGCATGAGCGTCAGTACAGGTCCAGGGGATTGCCTTCGCCATCGGTGTTCCTCCGCATATCTACGCATTTCACTGCTACACGCGGAATTCCATCCCCCTCTACCGTACTCTAGCCATGCAGTCACAAATGCAGTTCCCAGGTTGAGCCCGGGGATTTCACATCTGTCTTGCAAAGCCGCCTGCGCACGCTTTACGCCCAGTAATTCCG + ABCCCFCCCCCCGGGGGGGGGGHHHHHHHHHGGHGGGGGGHGHGGHHHHHHGGGGGHHHHHHHHHHHHGGGGHGHHHHHHGGGGHGGGGGGHHHHHGGGGGHHHHHGGGGGHHHHHHHHHHHHGGGGGGGHHHHHHHHGGGGHHHHGGHGHHHHHHHHHHHHHHHHHGGGGGGGGGGGGGGGGGGGGGGGGGGGFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF @M05193:157:000000000-JHBCT:1:1101:15595:2353 1:N:0:GAATTCGT+AGGCGAAG GACTACTGGGGTATCTAATCCTGTTTGCTCCCCACGCTTTCGTGCATGAGCGTCAGTGTTATCCCAGGAGGCTGCCTTCGCCATCGATGTTCCTCCACATATCTACGCATTTCACTGCTACACGTGGAATTCCACCTCCCTCTGACACACTCGAGCCTTGCAGTCCCAAATGCAGTTCCCAGGTTAAGCCCGGGGATTTCACATCTGGCTTACAAAACCGCCTGCGCACGCTTTACGCCCAGTAATTCCG + B@BCCFFFFCCCGGGGGGGGGGHHHHHHGHHGGHGGGGGGHFHGGHHHHHGGGGGGHHHHHHHHHHFGHGGGHHGHGHHHGGGGHGHHHHHHHHHHHHHHHHHHHHGGGGGEHHHHHHHHHHHGHHGGHHGHHHHHHHHGHHHHHHHHHGHHGHGGGHHGHHFHHHHHHHHHGHHHHFHHHHGHHGHHHGHGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGFFFFFFFFFFFFFFFFFFFFFFFFFFF @M05193:157:000000000-JHBCT:1:1101:15964:2367 1:N:0:GAATTCGT+AGGCGAAG CCTACGGGGGGCAGCAGTGGGGAATTTTGCGCAATGGGCTAACGCCTGACGCAGCAACGCCGCGTGGAGGATGAAGGTCTTTGGATTGTAAACTCCTGTCAGCAGGGAAAAAGCAACTTGACCTCATACGTCAGGGAGTTGATTGTACCTGCAGAGGAAGCCCCGGCTAACTCTGTGCCAGCAGCCGCGGTAATACAGAGGGGGCAAGCGTTATTCGGAATTATTGGGCGTAAAGGGCGCGTAGGCGGCT + BBBBBBBBB2>DGGGGGGGGGGGGHCHHDHCGGGGHHGHBGH<C/CHHHFGADGGHHHGGGGGGFGG?GDGFHHHHH9CFG0FGFGBFGGGGGFBGGGGGGGBFF?DGGFGFDFBF.BFFFB9BFFFFFFFAAAFFEDDFFF/FFFFB9FFFF/FFAFFFFFFFFFFFFFFFBBFFFFFFFFFFFFFFFFFFFFFFEFFFFFFFFFFFFFAD?EFFFFFEFFFFFFFFFFFFFFFEFADFF;DFFFFFFF @M05193:157:000000000-JHBCT:1:1101:13577:2375 1:N:0:GAATTCGT+AGGCGAAG GACTACGAGGGTATCTAATCCTGTTTGCTCCCCACGCTTTCGTGCCTCAGCGTCAGTACTGGTCCAGAGAGCCGCCTTCGCAACTGATGTTCCTCCCGATCTCTACGCATTTCACTGCTACACCGGGAATTCCACTCTCCTCTACCGTACTCTAGCTCGCCAGTTTCCATCGCCATTCCCAGGTTGAGCCCGGGGCTTTCACGACAGACTTAACAAACCGCCTACGCACGCTTTACGCCCAGTAATTCCG + ABBBAFBBBBBAGGGGGGGGGGHHHHHHHHHGGHGGGGGGHGHGGHHGHHHGGGGGHHHHHHHHFHHHHGGAGCFGGEHHGGGGHHFHHHHFHHHHGGGGGHHHHHGGGGGHHHHHHHHHGHHHGGGGGHGHHHHHHHHGHHHHFHGGHGHHHGHHHHGGGFGHHHHHHHHGGGGGHHHHHGGHGGGG0FGGGGGGGGFBFFFFDF;FFFFFFFFFFFFFFDDFFFFFFF-DFEFFEFFFFFFFFEBFFF @M05193:157:000000000-JHBCT:1:1101:15104:2429 1:N:0:GAATTCGT+AGGCGAAG GACTACCGGGGTATCTAATCCCGTTTGCTCCCCTAGCCTTCGTGCCTCAGTGTCAGGAACTGTCCAGAGACCCGCCTTCGCCACTGGTCTTCCTTCCGATCTCTACGCATTTCACTGCTACACCGTAAATTCCAGTCTCCTCTCCAGTCCTCAAGCAAGGCAGTATCGAATGCAGTTTCGGAGTTAAGCTCCGAAATTTCACACCCGACTTACCCCGCCACCTACGCACCCTTTACGCCCAATGAATCCG + A3AABFCDDDBDFCGGGGGGFFEHHH?GEGGGGB2GG3GGHGGGGHHGGF5DFFGF2FBF3F5BGDGHHB?FF@EE1EHGEFGEBCBG3GFHHHB3B/?F/0FGHFGEG?AGFHHFGFHHHHHC@GEGHHHF2F2@@2FFHFDG0F111FDGHHHHH00FHHDGDGH.>FHHHFD0<0=<-.<:CGHFH00:ED.-;CFGGGFG.C-AFGGE0F99?BF@AF?/FA-AF-;EF/9B.A=BDEFBFFBFF. @M05193:157:000000000-JHBCT:1:1101:20437:2430 1:N:0:GAATTCGT+AGGCGAAG GACTACACGGGTATCTAATCCTGTTTGCTCCCCACGCTTTCGTGCATGAGCGTCAGTGTTATCCCAGGGGGCTGCCTTCGCCATCGGTATTCCTCCCGATATCTACGCATTTCACTGCTACACCGGGAATTCCACTACCCTCTACCGCACTCTAGCCTGCCAGTATCCAATGCAATTCCCAGGTTGAGCCCGGGGATTTCACGACAGACTTACAAAACCGCCTGCGCACGCTTTACGCCCAGTAATTCCG + ABBBBFFFBBBDGGGGGGGGGGHGHHBHHHHGGHGGGGGGHGHGGHHHHBGGEGGGHHHHHHHHHHHGGGGGGGGHHHHHGGGGGGGGGGHGHHFHGGGGGGHHHGGFGGGHHHHHHHHHHHHHGGCFGHHHHHHHHHHHHHGHHHGGGGGHFHHHGGHHHHHHHHHHHHHHHHGHHHHGGGGGGEGGGGGGGGGGFFFFFFFFFF=FFFFFFFBFFFFFFFFFFFFFFFFFFFFF?DBFFFFFFFFFFF @M05193:157:000000000-JHBCT:1:1101:19654:2504 1:N:0:GAATTCGT+AGGCGAAG CCTACGGGTGGCAGCAGTGGGGGATCTTGCGCAATGGGCGAAAGCCTGACGCAGCCATGCCGCGTGGATGATGAAGGTCTTAGGATTGTAAAATCCTTTCACCGGGGACGATAATGACGGTACCCGGAGAAGAAGCCCCGGCTAACTTCGTGCCAGCAGCCGCGGTAATACGAAGGGGGCTAGCGTTGCTCGGAATTACTGGGCGTAAAGGGAGCGTAGGCGGATCGTTAAGTTGGGGGTGAAAGCCCGG + CBCCCCCCCCCFGGGGGGGGGGGGGGHHHHGGGGGHHHHGGGGGHHHHHHGGGGGHHHHHHGGGFGGGHHHHHHHHHGHHHGHHHHHHHHHHHHHHHHHHHHHGGGGGGGGGGGHHHHGGGHGHHGGGGGHHHHHHHHGGGGGGGGGGGGFGGGGGGGGGGGGFFDFFFFFFFFFFFFFFFFFFFFFFFFFFFFEEFFFFFFFFFFFFFFFFFFFFDFFFFFF;FFFFFFFFFFFFFFFFFFFFFFFFFF @M05193:157:000000000-JHBCT:1:1101:13802:2509 1:N:0:GAATTCGT+AGGCGAAG CCTACGGGGGGCAGCAGTAGGGAATCTTCCACAATGGACGAAAGTCTGATGGAGCAACGCCGCGTGAGTGAAGAAGGCCTTCGGGTTGTAAAGCTCTTTCGGTGGGGAAGAAACGGTTACGGTGAATACCCGTGACTAATGACGGTACCCAAAGAAGAAGCACCGGCTAACTACGTGCCAGCAGCCGCGGTAATACGTAGGGTGCGAGCGTTAATCGGAATTACTGGGCGTAAAGCGTGCGCAGGCGGTT + A3ABCBCCCCCCGGGGGGGGGGHHHHHHHHHHHHHHGHHGGGGGHHHHHHHHHHHHHHGGGGGGGGGGHHHHHHHHHHHHHHHGGHGGGHHHHGGGGGGGGGAFGGGGGGGGGGGFFFFFFFFFFFFFFFFFDFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFBFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFDDDFFFFFFFF9@ @M05193:157:000000000-JHBCT:1:1101:10738:2509 1:N:0:GAATTCGT+AGGCGAAG GACTACGAGGGTATCTAATCCCATTCGCTCCCCTAGCTTTCGTCCCTGAGTGTCAGTTATGGTCCAGCAAAGCGCCTTCGCCACCGATGTTCTTCCTGATCTCTACGCATTTCACCGCTACACCAGGAATTCCCGTGACCCCTACCAAACTCTAGCCTTGTAGTTTCCATCGCCTGTATGGGGTTGAGCCCCATTCTTTAACGACAGACTTACAGAGCCACCTGCGGACGCTTTACGCCCAATAATTCCG + ABBCCFCBCCCCGGGGGGGGGGHHHHHGGHGGGHGHHHHHHGHGGGHHHHHHHHHHHHHHHHHHHHHGHHHHHGGGGGHHGGGGHG/EGHHHHHHHHHHHHHHHHHGGGGGHHHHHGGGGGHHHHHGHHHHHHHGHGGGHHGGGHHHHHHHHHHHHHHHHHHHHHHHHHHHGGGGHHHGHHHGGGGHHHHGGGHGGGGGGGGGEGGGFGGGGGGGGGGGGGGGFFFFFFFFFFFFFFFFFFFFFFFFFFF @M05193:157:000000000-JHBCT:1:1101:14394:2524 1:N:0:GAATTCGT+AGGCGAAG GACTACTAGGGTATCTAATCCTGTTTGCTCCCCACGCTTTCGTGCCTCAGCGTCAGTACTGGTCCAGGTGGCCGCCTTCGCCACAGATGTTCCTCCCGATCTCTACGCATTTCACCGCTACACCGGGAATTCCACCACCCTCTACCGTACTCTAGCGACCCAGTTTCCAATGCCATTCCCAGGTTGAGCCCGGGGATTTCACATCAGACTTAAATCACCGCCTACGCACGCTTTACGCCCAGTAATTCCG + AAABBFB5DFBCGBGFGFGGGGHHHCHHHHFF?GEAEEEGHEB2FEHHHHHGGEEGFH55AFBEAGHFG3FAGEEGEGHHGGGFHH3GFHDHHGBGG?GGCHDHHHEGGDFHFFH3/<CCCHHFDCGGGHGHH<GGHHE?GHCHFHGGHGHGHGHH--<CGHHFHH00;;GCHH0GFHHHHBE9BFF9FFFG@B@BAFGFFFFFF/;BFFFFF/9/BBF.BBFDFDB.@@B?EEFFFFBBFFBFFFFFFF @M05193:157:000000000-JHBCT:1:1101:16497:2533 1:N:0:GAATTCGT+AGGCGAAG GACTACGCGGGTATCTAATCCTGTTCGCTCCCCACGCTTTCGTGCCTCAGCGTCAGTACCAGCCTAGTCAGCTGCCTTCGCAATCGGGGTTCTGGACGGTATCTATGCATTTCACCGCTACACCGCCCATTCCGCCGACCTCGTCTGGACTCAAGCCCCGCAGTTTCCAAGGCAGTTCCGTTGTTGAGCAACGGGCTTTCACCCCGGACTTACGAAGCCGCCTACGCACCCTTTAAACCCAATAAATCCG + CCBCCFCBCCCBGGGGGGGGGGHHHHGGGHGGGHGGGGGGHGHGGHHHHHHGGGGGHHHHHGGGHHHHHHHHHHHHHHHHGGGGHGGGGGGHHHHHHGGFHGHHHHHHHHHHHHHHGGGGGHHHGGGGGHHHHHGGGGGGGHGHGHGHHHHHHHGGGGGGGGGGGGGGGGGGGGGGFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF @M05193:157:000000000-JHBCT:1:1101:20181:2540 1:N:0:GAATTCGT+AGGCGAAG CCTACGGGGGGCAGCAGTGGGGAATATTGCGCAATGGGCGAAAGCCTGACGCAGCGACGCCGCGTGAGGGATGAAGGCCTTCGGGTTGTAAACCTCTTTCAGCAGGGACGAAAACAGACGGTACCTGCAGAAGAAGCACCGGCCAACTACGTGCCAGCAGCCGCGGTAATACGTAGGGTGCAAGCGTTGTCCGGAATTATTGGGCGTAAAGAGCTCGTAGGCGGTTTGTCACGTCGGCTGTGAAAACTCG + BBBCBBBCBCCCGGGGGGGGGGGGHCHHHHGGGGGHGHHGCGGGHHHHHGADGGCGFDGGGGGGGC.DFGGFGGGGGGGGGGGGGGGGGGGGFFFFFFFFFFFFFFFFEFFFFFFFFFFFDAEDFFFFFBFFFFFFFFFEF;FFFFEFFFF?FAFFFEFFFE>FFF.;EB9/./:.A./:AFFFE?;DFA9FA-D=EFFFFFEEFF@;.?FFBBFFF>DEDFC-:BD?:/:FD?A9BD9-;9BFFFFFB. @M05193:157:000000000-JHBCT:1:1101:11032:2544 1:N:0:GAATTCGT+AGGCGAAG GATTTCACCGCTGACTTACCGACCCACCTACGCACCCTTTAAACCCAGTGAATCCGGATAACGCTTGCACCCTCCGTATTACCGCGGCTGCTGGCACGGAGTTAGCCGGTGCTTATTCCTCTGGTACCGTCAACCACATTAGAAAACGTGGAGTTCTTCCCAGATAAAAGGAGTTTACAACCCAGAGGGCCTTAATCCTCCACGCGGCATGGCTGGTTCAGACTTGCGTCCATTGACCAATATTCCTTAC + CCCCCFFFFBCBGGGGGGGGGGGGGHGGHHHGGGGGGHHHHHHHHGGHHHHHHHHGFGGGHHGGGGGHHHHGHHGGG1GH4GHGGGGGGGGHHHGHHGG?FGHHHHHGGFGGHHGHHHHHHHHHHHHHGGGGGHHGGHHHHGHHHHHGHHGGGHHHHHHHHHHHHHHHHHGGGGGGGGGGGGGGGGGGGGGGGGGGHGGFFFFFFFFFFFFFFFFFFFFFFBFFFFFFFFFFFFFFFFF?FFFFFFFF0B @M05193:157:000000000-JHBCT:1:1101:17408:2548 1:N:0:GAATTCGT+AGGCGAAG GACTACAGGGGTATCTAATCCTGTTTGCTCCCCACGCTTTCGCTCCTCAGCGTCAGTTACTGCCCAGAGACCCGCCTTCGCCACCGGTGTTCCTCCTGATATCTGCGCATTTCACCGCTACACCAGGAATTCCAGTCTCCCCTACCAGACTCAAGTCATGACAGTATCGAGTGGCGTCTCCAGGTTGAGCCTGGAGTTTTCACACCCGACTTACCCCGCCACCTACGCACCCTTTACGCCCAATGAATCC + CCCCCFFFFCCBEGGGGGGGGGHHHHHHHHHGGHGGGGGGHGGGHGHHHHHGGGGGHHHHHHHHHHHHHFHGHGGGGGHHGGGGHGGGGFGHHHHHHHHHHHHHHHGGGGGHHHHHGGGGGHHHHHGGHHHHGHHHHHHHHGHGHHHHHGHHHHHHHHHHHHHHHHHEFGGHHGGGGGCHHHHHGGHGGGGGGGGGFGEGGGGGFGGGGGGGGFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF @M05193:157:000000000-JHBCT:1:1101:15460:2557 1:N:0:GAATTCGT+AGGCGAAG GACTACGCGGGTATCTAATCCTGTTTGCTCCCCACGCTTTCGCTCCTCAGCGTCAGTTACGGCCCAGAGACCCGCCTTCGCCACCGGTGTTCCTCCTGATATCTGCGCATTTCACCGCTACACCAGGAATTCCAGTCTCCCCTACCGCACTCTAGTCTGCCCGTATCGCCTGCAGACCCGCAGTTGAGCCACGGGATTTCACAGCCGACGCGACAAACCGCCTACGAGCTCTTTACGCCCAATAATTCCG + A?AAAFADBBBBGGGGGFGFGGEHFGCGHHFF2FEEEEGGHECGGGHEGCAEEEE1DDGDFEFEE11?E?2EE>EEGEG?>E?GGE/E/?CFD?FGH22GFHHFF2@CCDGBD1<1//>ACGHGHHEGFGFHBHFGFDHHF</CFHCFGGGHHB0CGF:CGGB99CFA.CA:FG/BG.E?--9DF/9BF//B;B>?DBBF/9//9.@-9-9-;B@@FFF-;@A.:A9.-.;/B/;;.;-999FFFFFFB/ @M05193:157:000000000-JHBCT:1:1101:12922:2561 1:N:0:GAATTCGT+AGGCGAAG GACTACTGGGGTATCTAATCCTGTTTGCTCCCCACGCTTTCGTGCATGAGCGTCAGTACAGGTCCAGGGGATTGCCTTCGCCACCGGTGTTCCTCCTGATATCTGCGCATTTCACCGCTACACCAGGAATTCCAGTCTCCCCTACCGCACTCTAGTCTGCCCGTATCGACTGCAGGCTCGGGGTTGAGCCCCGAGTTTTCACAGCCGACGTGACAGACCGCCTACGAGCTCTTTACGCCCAATAATTCCG + CCBCBFFFFCCCGGGGGGGGGGHHHHHHHHHGGHGGGGGGHGHGGHHHHHHGGGGGHHHGHHHHHHHHGGGGHHHHHHHHGGGGHGGGGGGHHHHHHHHHHHHHHHGGGGGHHHHHGGGGGHHHHHGHHHHHHHHHHHHHHGHGHGGGGGGHHHHHHHHHHHGGGHHGHGGGHHHHGHHGGGGGGGGGGGGGGGGCFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF @M05193:157:000000000-JHBCT:1:1101:19466:2569 1:N:0:GAATTCGT+AGGCGAAG CCTACGGGTGGCAGCAGTGGGGAATTTTGGACAATGGGCGCAAGCCTGATCCAGCAATGCCGCGTGAGTGAAGAAGGCCTTCGGGTTGTAAAGCTCTTTTGTCAGGGAAGAAACGGTTCTGGCTAATACCTGGGGCTAATGACGGTACCTGAAGAATAAGCACCGGCTAACTACGTGCCAGCAGCCGCGGTAATACGTAGGGTGCAAGCGTTAATCGGAATTACTGGGCGTAAAGCGTGCGCAGGCGGTT + A?ABBABBDBBFGGGGGGGGGGGGHHHHHGHHHHHGGGHGEGGGHGHHHHHGHHHHHHHFHGGG>>EGGHHFHHHHHHGHHHHGGHGGGHGHHGHHHHHHHGHHHHHGHGHHGGEGEGGHHHHGHHFHHHHHHGGGGHHHHHHGGGHEHHHHFHHGHGHHGHHGGGGGGHFGHGF9FGGFGGGFGFACBGEFFEFFFBFFFFBEDFFFFFADFFFFFFDFEFFFFFFFECFFDFBBD.9AFF;=FFFF@D @M05193:157:000000000-JHBCT:1:1101:16624:2570 1:N:0:GAATTCGT+AGGCGAAG CCTACGGGTGGCAGCAGTGGGGAATCTTGCGCAATGGGCGGAAGCCTGACGCAGCAACGCCGCGTGGGGGACGAAGGCTCTCGGGTTGTAAACCCCTTTCAGTAGGGACGATTGTGACGGTACCTACAGAAGAAGCACCGGCTAACTACGTGCCAGCAGCCGCGGTAATACGTAGGGTGCGAGCGTTGTCCGGATTTATTGGGCGTAAAGAGCTCGTAGGCGGTTCAGTAAGTCGGTTGTGAAAACTCTA + CCCCDCDDBBCFGGGGGGGGGGGGHHHHHHGGGGGHHHHGGGGGHHHHHHGGGGGHHHGGGGGGGGGGGGGGGGGGGHHHHHGGGHGGGHHHHHGGGHHHHHHHHHGHGGGGGGFHHHGGFHGGGGGGGGGGGGGGGGGGGGGGGGGGFFFFFFFFFFFFFFFFFADFFFFFFFFFFBEFFFFBFDFFFFFFFFDFFFFFFFFFFFFFFHFFFFFFFFFFFFF@FFFFFFFFBFFFBFFDFFFFFFFFFF @M05193:157:000000000-JHBCT:1:1101:16163:2599 1:N:0:GAATTCGT+AGGCGAAG GACTACGAGGGTATCTAATCCTGTTCGCTCCCCATGCTTTCGCTTCTCAGTGTCAGTAATGGCCCAGAGACCTGCCTTCGCCACTGGTGTTCTTCCGAATATCTACGAATTTCACCTCTACACTCGGAATTCCACTCACCTCTCCGACGTTCTAGCGATGCAGTCTTAAAGGCAATTCCGGAGTTGAGCTCCGGGCTTTCACCTCTAACTTACAAAGCCGCCTACGTGCGCTTTACGCCCAGTAATTCCG + ABBAAFBABBBBGGGGGGGGGGHHHHGGGHGGGHHFHHHFHGGGHGHHHGGHHHHGHHHHHHH3GEGHHHHGHGHGHHFHGGGGHGHHHFGHHHHH?AEFGHHHHHEGHHHHHHFHHFHGHHGHHGGGFGHHHHHHHHHHHHHHHGGGGDEFGHHHGGHGGHHHHHHHHHHHHHHHGHHGGCFGHHGGHHHHGGGGGGGGGGGGGGGGGGGGGGGGGGGCGFFFFFFFFFFFFFFFEFFFFFFFFFFFFF @M05193:157:000000000-JHBCT:1:1101:19273:2663 1:N:0:GAATTCGT+AGGCGAAG GACTACGAGGGTATCTAATCCTGTTTGCTACCCACGCTTTCGTGCCTCAGCGTCAGTTATGGTCTAGAAAGCCGTTTACACCACCGGTGTTCCTCCTGATATCTACGCATTTCACCGCTACACCAGGAATTCCACTTTCCTCTACCGTACTCTAGCGACCCAGTTTCCAATGCCATTCCCAGGTTGAGCCCGGGGATTTCACATCAGACTTAAATCACCGCCTACGCACGCTTTACGCCCAGTAATTCCG + CCCCCF3ABCCBGGGGGGGGGGHHHHHHHFHHHHGGGGGGHGEFGHHHHHHGGGGGHBHHHHHHFHHHFHHHHG>EHGHFHHHHGGGGGGGHHHHHHHHHHHHHGHGGGGGHHHHHGGGGGHHHHHGHHHHHHHHHFHHHHHHHHHGFGGHHHHHHFGGGGHHHHHHGFHHHHFDGHHHHHHGHHHHHHHGGGGGGGGGGGGGGFGFFGGGGGGGFGGFFFFFFFFFDFF=DEFFFEFFFFFFFFFFFFF @M05193:157:000000000-JHBCT:1:1101:15738:2674 1:N:0:GAATTCGT+AGGCGAAG CAAGCCAAACAGTATCGAAGGCAATTCTGTGGTTGAGCCACAGGCTTTCACCCCCGACTTGAATGGCCGCCTACGCGCCCTTTACGCCCAGTAAATCCGAGCAACGCTAGCCCCCTTCGTATTACCGCGGCTGCTGGCACGAAGTTAGCCGGGGCTTATTCTTCCGGTACCGTCATTATCGTCCCGGAAAAAGAGCTTTACAACCCTAAGGCAGATCGGAAGAGCACACGTCTGAACTCCAGTCACGAAT + CCCBCCFFFBFFGGGGGGFGGGHHHHHHHFHHGG2FGFHFHHGHGHHHHHHHHGGGGGGGHGGHHHGHGGGGGHGGGGGGGHHGHGGCGGEGGHHHHHG/>C/BFDGGGGFHHGGGHHHHHHHHHHGGGGGGGGHHHGHHGGGGEHHFHHGGGFGFFGGGGGGFFGGGGFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFDDF @M05193:157:000000000-JHBCT:1:1101:19263:2679 1:N:0:GAATTCGT+AGGCGAAG GACTACACGGGTATCTAATCCTGTTCGCTCCCCATGCTTTCGCTTCTCAGTGTCAGTAATGGCCCAGAGACCTGCCTTCGCCATCGGTGTTCCTCCTGATATCTGCGCATTTCACCGCTACACCAGGAATTCCAGTCTCCCCTACCACACTCTAGCTATCTAGTTTCAAAGGCAGTTCCGGGGTTGAGCCCCGGGCTTAGATCGGAAGAGCACACGTCTGAACTCCAGTCACGAATTCGTATCTCGTATG + CCDDCFFFCCCCGGGGGGGGGGHHHHGGGHGGGHHHHHHHHGGGHGHHHHHHHHHHHHHHHHHHHGHHHHFGHHHHHHHHGGGGHGGGGGGHHHHHHHHHHHHHHHGGGGGHHHHHGGGGGHHHHHGHHHHHHHHHHHHHHGHGHHHHHGHHHHHHHHHHFHHHHHHHHHHHHGGHHHHGGGGGGGHHHHGGGGGGGGGGGGGGDGGGGGGFGGFFFFFFFFFEFFFFFFFFFFFFFHFFFBFFFFFFFF @M05193:157:000000000-JHBCT:1:1101:15665:2684 1:N:0:GAATTCGT+AGGCGAAG CCTACGGGTTGCAGCAGTGGGGAATATTGGACAATGGGCGCAAGCCTGATCCAGCCATGCCGCGTGTGTGAAGAAGGCCTTCGGGTTGTAAAGCACTTTTGTTGGGGAAGAAAAGCACTGGATTAATACCTTGGTGTTCTGACGGTACCCAAAGAATAAGCACCGGCTAACTTCGTGCCAGCAGCCGCGGTAATACGAAGGGTGCAAGCGTTACTCGGAATTACTGGGCGTAAAGCGTGCGTAGGCGGTT + AABABABBA2AFGGGGGGGGGGGGHHHHHBGHHHHHCHHGGGGGGGHHHHHHHHHGHHHGHGGCEEEHHHHHHHHHHHHHFHGGGGGGGHHHHHHHHHHFD2BG/CDG?GHHB2FG2G2@GBGHHFHHHCHHDHHGE/FGHHHG-CGGFHBCGHHHHHHHHHHHGGGGGHHGGHGGGGGGGGGGGGGGGCEGGFFFFAEDFFFEFFFFFFFFFEFFFFFFEFEBBFF.A=FFFFFFFD?DF.ABF/;DBD @M05193:157:000000000-JHBCT:1:1101:20024:2725 1:N:0:GAATTCGT+AGGCGAAG CCTACGGGAGGCAGCAGTGGGGAATATTGGACAATGGGCGCAAGCCTGATCCAGCAATGCCGCGTGGATGAAGAAGGTTTTCGGATTGTAAAGTCCTTTCGGCGGGGACGATGATGACGGTACCCGCAGAAGAAGCCCCGGCTAACTTCGTGCCAGCAGCCGCGGTAATACGAAGGGGGCTAGCGTTGCTCGGAATGACTGGGCGTAAAGGGCGCGTAGGCGGATTGATCAGTCGGGCGTGAAATTCCTG + CCBCCCCCCCCCGGGGGGGGGGGGHHHHHHHHHHHHHHHGGGGGHGHHHHHHHHHHHHHHHGGGGGGGHHHHHHHHHCGHGHGGGGGHHHHHHHHHHHHHGGGGGGGC@-@C<D<00CF:C.CGA@-CBFBCF//FFF?FA?CFBFFFABA;AFF?FAFF;B9C.;AFB/9:EAFFFF;9FABC.ABAFFACAE.BB/BBEFEF-:DFFFF.-B->=D/@C--ABBFBF//BDFF-@9B9AFBBBB9/ @M05193:157:000000000-JHBCT:1:1101:14065:2751 1:N:0:GAATTCGT+AGGCGAAG GACTACGGGGGTATCTAATCCTGTTTGCTCCCCACGCTTTCGTGTCTCAGCGTCAGTATCGGGCCAGTAAGCCGCCTTCGCCACTGGTGTTCTTCCGAATATCTACGAATTTCACCTCTACACTCGGAGTTCCGCTTACCTCTCCCGAACTCAAGATACCCAGTATTAAAGGCTATTCCAGGGTTGAGCCCTGGGCTTTCACCTCTAACTTAAGTATCCGCCTACACACCCTTTACGCCCAGTAATTCCG + ?ABBBFBBBBBBGGGGFGGGGGHHHHHHHHHGGHGGGFGGHGHGHHHHHHHGGGGGHHHHFGGGGGHHFHHHHGGGGGHHGGGGHHHHHHGHHHHHDGGGGHHHHHGGHHHHHHHHHHHHHHHHHGGGGGHHHGGGGGHHHHHHHHGGGGHHHHHHHHHHHHHHHHGGGGGGGGGGGGGGGGGGGGGGGGGFGGGGGFFFFFFFFFFFFFFBFFBBFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF @M05193:157:000000000-JHBCT:1:1101:14767:2768 1:N:0:GAATTCGT+AGGCGAAG CCTACGGGGGGCAGCAGTGGGGAATATTGGACAATGGGCGAAAGCCTGATCCAGCAATGCCGCGTGAGTGATGAAGGCCTTAGGGTTGTAAAGCTCTTTTACCCGGGAAGATAATGACTGTACCGGGAGAATAAGCCCCGGCTAACTCCGTGCCAGCAGCCGCGGTAATACGGAGGGGGCTAGCGTTGTTCGGAATTACTGGGCGTAAAGCGCACGTAGGCGGCTTTGTAAGTTAGAGGTGAAAGCCCGG + AAABBBAB@DBBG??GGGGGGCGG2GHHHHHHHHHGHGHGGGGGHHHHHHHGHHHHHHHHHGGGFGGGHGHHHHGGHAGHFHCFFCGGGHHHHHHHHBFGGGFFEACGGDFGGFGGGGGGGGGGFFFBDFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFEDCFF.AEFFFE.ADFFFFFCFFFFDDFEEFDBAA?FFFFFFFFFFFDFFFFFFFF.FFFFFFFB-FFBFFFB/FFFFFFFFFFFFFECF @M05193:157:000000000-JHBCT:1:1101:13059:2801 1:N:0:GAATTCGT+AGGCGAAG GGGAGACTGGAATTCCTGGTGTAGCGGTGAAATGCGCAGATATCAGGAGGAACACCGATGGCGAAGGCAGGGCTCTGGGCCATTACTGACACTGAGAAGCGAAAGCATGGGGAGCGAACAGGATTAGATACCCCCGTAGTCAGATCGGAAGAGCACACGTCTGAACTCCAGTCACGAATTCGTATCTCGTATGCCGTCATCTGCTTGAAAAAAAACGAAATGAATGGGCGTAACAGCTCCGCACTCGACT + 1>A111ADFAF11GA33111AAF1FEEE?EFDFH1EEGA/AG21DFAGE//0FGCFE/E//A>/>E//0BE/E0?>B10/?/BEGGHH2B1B1FFC11F<//EEC<1>B</<//?C//B/CG011?1GH1?<1<->C-.<=GDG0/.AC.EHH/:CB.ACG/CC0CF00::;FB0;.9.C;9///;9;F..9BFG@F@-9BBFFFFFFBFB-9---9----/;//;B--9--;-//-;//A-@--/9--A @M05193:157:000000000-JHBCT:1:1101:14254:2805 1:N:0:GAATTCGT+AGGCGAAG CCTACGGGGGGCAGCAGTAAGGAATATTGGTCAATGGACGCAAGTCTGAACCAGCCATGCCGCGTGGAGGATGAAGGTCCTCTGGATTGTAAACTTCTTTTATCTGGGAAGAAATGCATTTTTTCTAAAGTGTTTGACGGTACCAGAAGAATAAGCACCGGCTAACTCCGTGTCAGCAGCCGCGGTAATACGGAGGGTGCAAGCGTTATCCGGATTCACTGGGCTTAAAGGGTGCGTAGGTGGGCAGTTA + AAAAAAAD10ADGGGGGGFGGGHHHFHHHF1FHHHHHFHG@GGG>FHHHBFGHHHFH<GCGGGGGGCGECGHHHHHB<D0DHFHHFHHFCGHCCGC0CH:0:00;GHCGBA.FGGGB0BBGFB?/;0F;B0C9;CF0FC@-9-9BF-FFBBB/B/:99-9-;-:BF///:-9-;9FFFFF9=@9@BA9FFFFFFB;;-9-//BEFFABEF/A@-9ABB//BFFEFE//9BBBFEB@E?@F99-9?FBFFF @M05193:157:000000000-JHBCT:1:1101:17404:2813 1:N:0:GAATTCGT+AGGCGAAG CCTACGGGTGGCAGCAGTGGGGAATATTGCACAATGGGCGCAAGCCTGATGCAGCGACGCCGCGTGAGGGATGACGGCCTTTGGGTTGTAAACCTCTTTCAGCAGGGAAGAAGCGCAAGTGACGGTACCTGCAGAAGAAGCACCGGCTAACTACGTGCCAGCAGCCGCGGTAATACGTAGGGTGCAAGCGTTGTCCGGAATTATTGGGCGTAAAGAGCTTGTAGGTGGTTTGTCGCGTCTGTTGTGAAAT + ABCCCCCCBCCFGGGGGGGGGGGGHHHHHHHHHHHHHHHGGGGGHGHHHHGHHHHGGGGGGGGGGGDGHGGHGHHGGGFGHHHHGGGGGHHHHHHHHHHHHHHGGHGGHGHHGHGGFCFDGHHDGECEHHHHGHHHHFGCGHHHFCAFGFGGGGFGGGGGGGGGGGGGF=.;90BFE:9ADEBEFFFBADFFFFFFFFF;EFFFFFEFFFF>DFFFFBFFFFFFE0;E9AFFFFFFF@B?EFFFFDBBFF @M05193:157:000000000-JHBCT:1:1101:13265:2835 1:N:0:GAATTCGT+AGGCGAAG