BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	375,202	T→A	3.8%	intergenic (+321/‑257)	mhpE → / → mhpT	4‑hyroxy‑2‑oxovalerate/4‑hydroxy‑2‑oxopentanoic acid aldolase, class I/3‑hydroxyphenylpropionic transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	375,202	0	T	A	3.8%	405.4 / 4.4	123	intergenic (+321/‑257)	mhpE/mhpT	4‑hyroxy‑2‑oxovalerate/4‑hydroxy‑2‑oxopentanoic acid aldolase, class I/3‑hydroxyphenylpropionic transporter
Reads supporting (aligned to +/- strand): ref base T (53/60); new base A (3/2); total (59/64)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 124 total aligned reads displayed.

GGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATC‑TGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCAC  >  NC_000913/374968‑375439

ggCGTGAGAATTTGTAGGTCTGATAAGTCGCGTTAGCGTCGCATCAGGCATCTGCGCGCGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGc                                                                                                                                                                                                                                             <  1:1079896/239‑1 (MQ=255)

ggCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTCGCGTTGCATCAGGCATC‑TGCGCACGACCGCCGGATGCGGCGTGCACGCCTTATCCGGCCTCCGGGCGGCGGCGGAGATTGTAGCGCAGATAAGGAGCG‑TTAg                                                                                                                                                                                                                                              >  2:353089/1‑238 (MQ=11)

ggCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGGTGCGGCGTGAACGCCCTATCCGCCCTACGGTCTGCGCCGGCATTTGTAGGCCTGATAAGACTCG‑TTAGc                                                                                                                                                                                                                                             >  2:691107/1‑239 (MQ=11)

ggCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGc                                                                                                                                                                                                                                             >  1:736501/1‑239 (MQ=255)

   gTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGt                                                                                                                                                                                                                                           <  2:555064/238‑1 (MQ=40)

           ttGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCGCGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAgg                                                                                                                                                                                                                                  <  1:1070451/239‑1 (MQ=38)

           ttGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAgg                                                                                                                                                                                                                                  <  2:184658/239‑1 (MQ=38)

            tGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGc                                                                                                                                                                                                                                 >  1:396269/1‑239 (MQ=38)

               aGGACTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCAt                                                                                                                                                                                                                               <  1:73587/238‑1 (MQ=37)

                  tCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGGCGCATCAGGCATC‑Tgc                                                                                                                                                                                                                          >  1:536228/1‑239 (MQ=37)

                    tGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGGGAATTTGTAGGCCTGATAATACGCG‑TTACCGTCGCCTCAGGCATC‑Tgcg                                                                                                                                                                                                                         >  2:633452/1‑238 (MQ=19)

                      agaaGACGCGGTGGCGTGGCAGCAGGCAGCTGCGCACGGCTGCCGGATGCGGCGTGGACGCCTTGTGCGGCCTACGGAGGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCa                                                                                                                                                                                                                       <  2:497786/236‑1 (MQ=11)

                                ttAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTGACCGTCTTATCCGGCCCACGGGTGGCGCGGGCATTTGTAGGGCTGGTAAGACCCC‑TTAGCGGCCCATCAAGCCTC‑CGGGCAAGACTGCCggg                                                                                                                                                                                                            >  2:579494/1‑238 (MQ=11)

                                ttAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGa                                                                                                                                                                                                            >  1:1330350/1‑239 (MQ=38)

                                aaaGCGTCGCATCAGGCAT‑GGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACACTTCGC‑CCGTAAGCCGTAGGCCTGATAAGACGCGCAAAGCGTCGCATCAGGCAT‑‑GGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCCGTCTACACTTCGC‑CCATAAACCGTAGGCCTGATAAACCGCGCATAGCCTCGCATCAGGCCGC‑TCCGCACAACttttctttg                                                                                                                                                                                                          >  2:484127/3‑231 (MQ=21)

                                  aGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACACTTCGCCCG‑TAAGCCGTAGGCCTGATAAGACGCGCAAAGCGTCGCATCAGGCAT‑‑GGCGCACGACTGCCGGATgc                                                                                                                                                                                                         <  1:199757/239‑1 (MQ=11)

                                     gTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCC‑GATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATgcg                                                                                                                                                                                                        <  1:975449/237‑1 (MQ=38)

                                      tCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACGGCCCGACgcggc                                                                                                                                                                                                      >  2:317990/1‑239 (MQ=255)

                                      tCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGGGAATTTGTAGGCCTGATAAAACGCG‑GTAGCGCCCCATCCGGCATC‑TGCGCACGCCTGCCGGCTgcggc                                                                                                                                                                                                      >  2:862884/1‑239 (MQ=11)

                                      tCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCCGGCATC‑TGCGCACGACTGCCGGATgcggc                                                                                                                                                                                                      >  1:985744/1‑239 (MQ=39)

                                       cGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGCGTGGCGCGAGAATTTGTAGGCCTGTTAAGACGCG‑GTAGCGTCGCATGAGGCATA‑TGCGCACGACCGCCGGGTgcggcg                                                                                                                                                                                                     >  2:218342/1‑239 (MQ=11)

                                         cATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTa                                                                                                                                                                                                   >  2:214131/1‑239 (MQ=40)

                                          aTCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTaa                                                                                                                                                                                                  <  2:736501/239‑1 (MQ=255)

                                          aTCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTa                                                                                                                                                                                                   <  2:1108718/238‑1 (MQ=40)

                                                   cTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCCCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCCTC‑TGCGCCCGACCGCCGGATGCGGGGTAAACCCCTTa                                                                                                                                                                                          >  2:1096597/1‑238 (MQ=255)

                                                    tGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATc                                                                                                                                                                                        <  2:568472/239‑1 (MQ=40)

                                                         aCGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                                                                                                                                                   <  1:999503/239‑1 (MQ=38)

                                                              tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCCCGGCTGCCGGATGCCGCGTGAAGGCCTTACCCGGCCCACGGGGGGCGCGAGACTTTTTCATCCAGATAAGACCCG‑CTAGCGTCGCATCAGGCATC‑TGCGCCAGACTGGCGGGTGCGCGGTAAGCGGCCTCTCCGGcacacgt                                                                                                                                                                              >  2:609101/1‑233 (MQ=11)

                                                               gtcGGATGCGGAGTAGACGCCTCATCTGGCTCGCGGGTCGCGCGGGAAATTGTAGTATTGAAAAGAGGTG‑TTAGGGTGGCTACAGGCGTC‑TGAGCACGACTGCCGGATGCGGCGTGAACGCCTTAGCCGGCCTACGGGTGGGGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGa                                                                                                                                                                             <  2:869576/237‑1 (MQ=11)

                                                               gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGa                                                                                                                                                                             <  1:841517/239‑1 (MQ=39)

                                                                ccGGGTGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGa                                                                                                                                                                             <  2:741017/238‑1 (MQ=39)

                                                                ccGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGTGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCCCGCCTGCCGGATGCGGCGAAAACGCCTTATCCGGCCTACGGAt                                                                                                                                                                            >  2:832798/1‑239 (MQ=255)

                                                                ccGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTCGGCCTGATCAGACGCG‑TTAGCGTCGCATCCGGCATC‑CTCGCACCACTGCCGCGTGCGGAGTAAACGCCCTATCCGGCCTCCgtgc                                                                                                                                                                            >  2:18694/1‑236 (MQ=11)

                                                                          cgTGAGCGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTATGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTGAACGCCTTAGCCGGCCTACGGGTGGCGGGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGGACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGa                                                                                                                                                                   <  1:239800/238‑1 (MQ=255)

                                                                          cgTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAgg                                                                                                                                                                                                                                  <  2:13639/176‑1 (MQ=39)

                                                                          cgTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAgg                                                                                                                                                                                                                                  >  1:13639/1‑176 (MQ=39)

                                                                                      aTCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTg                                                                                                                                                      >  1:1011743/1‑239 (MQ=255)

                                                                                            ccTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAAGCCCTTATCCGGCCTACGGCTGGCGCGGGGATTTGTAGGCCTGCTAAg                                                                                                                                                 >  2:506981/1‑238 (MQ=255)

                                                                                            ccTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCCACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCCTCAGGCATC‑TGCGCCCGCCTGCCGGGTGCTGCGTAAACGCCTTTTTCGGCCTACGGCTGGCGACGGACTTTGCAGGCCTGCTAAGa                                                                                                                                                >  2:523490/1‑239 (MQ=14)

                                                                                             cTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGa                                                                                                                                                <  1:279317/238‑1 (MQ=255)

                                                                                             cTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGa                                                                                                                                                <  2:1330350/238‑1 (MQ=255)

                                                                                                     ggCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGc                                                                                                                                       >  1:164176/1‑239 (MQ=255)

                                                                                                      gcgcGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCg                                                                                                                                      >  1:297927/1‑239 (MQ=255)

                                                                                                      gcgcGGGAATTTGTAGGCCTGATAAGAAGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACCCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACCGCCGGATGCGGCGCAAACGCCTCATCCGGCCTACGGATTGCGCGGGCATTTGTATGCTCGATaaagacaggtagcg                                                                                                                                      >  2:687712/1‑227 (MQ=14)

                                                                                                        gcgGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGt                                                                                                                                     <  2:951779/238‑1 (MQ=40)

                                                                                                             aaTTTGTAGGCCTGAGAAGAGGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCAt                                                                                                                                <  2:446217/238‑1 (MQ=255)

                                                                                                              aTTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCa                                                                                                                              <  1:89345/239‑1 (MQ=38)

                                                                                                                    aGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCAt                                                                                                                         <  1:87472/238‑1 (MQ=38)

                                                                                                                    aGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCAt                                                                                                                         <  1:75407/238‑1 (MQ=38)

                                                                                                                       ccTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCAt                                                                                                                                                                                                                               <  1:576852/134‑1 (MQ=35)

                                                                                                                        cTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATC‑TGca                                                                                                                   <  2:600173/239‑1 (MQ=37)

                                                                                                                           aTAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACCCCTTATCCGGCCTACGGATGGCGCGGGGATTTGTAGGCCTGATAAGAAGCGGTGGCGGCGCATCCGGCATC‑TGccca                                                                                                                 >  2:221798/1‑238 (MQ=11)

                                                                                                                                aCGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGGCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGAAAACGCCTTATCCGGCCTAAGGATGGCGCGGGCATTTGTTCGCCTGATAAGACGCGTTCGCGTTCCACCAGGCATC‑TGCACACGACt                                                                                                            >  2:1302018/1‑238 (MQ=11)

                                                                                                                                      ttAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGGATTTGTGGGCCTGATAAGACGCGTTAGCGTCGCCTCAGGCGTC‑TGCCCAAGAGTGCCGGa                                                                                                      >  2:200620/1‑239 (MQ=255)

                                                                                                                                      ttAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATC‑TGCACACGACTGCCGGa                                                                                                      >  1:1210915/1‑239 (MQ=38)

                                                                                                                                      ccaGCGTCGCATCAGGCATCATGCAC‑CAATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGTT‑‑CG‑GACTATTTGTAGGCACGATAAGACGCG‑ACAGCGTCGCATCAGGCATCATGCAC‑CAATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGTT‑‑CG‑GACTATTTGTAGGCACGATAAGACGCGACAGCGTCGCATCAGGCATCATGCAC‑CAATTGCCGGATGCggc                                                                                                >  2:902169/3‑237 (MQ=17)

                                                                                                                                       tAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATC‑TGCACACGACTGCCGGa                                                                                                      <  2:1011743/238‑1 (MQ=38)

                                                                                                                                        aGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATC‑TGCACACGACTGCCGGAt                                                                                                     >  1:405362/1‑238 (MQ=38)

                                                                                                                                               cATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATC‑TGCACACGACTGCCGGATGCGATaa                                                                                              >  1:353616/1‑238 (MQ=40)

                                                                                                                                                       aTC‑TGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATC‑TGCACACGACTGCCGGATGCGATAAACGTCTTGt                                                                                     <  1:806881/239‑1 (MQ=255)

                                                                                                                                                        tC‑TGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATC‑TGCACACGACTGCCGGATGCGATAAACGTCTTGt                                                                                     >  1:817300/1‑238 (MQ=255)

                                                                                                                                                                  gACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATC‑TGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACa                                                                           >  1:1145363/1‑239 (MQ=255)

                                                                                                                                                                    cTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGt                                                                                                                                                                                                                                           <  1:629165/79‑1 (MQ=39)

                                                                                                                                                                    cTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGt                                                                                                                                                                                                                                           >  2:629165/1‑79 (MQ=39)

                                                                                                                                                                    cTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGGAAACGCCTTATCCGGCCTACGGATGGCGCGGGAAATTGTAAGCCTGATAAGACGCGTTAGCGTCGCGTCAGGCCAC‑GGCACACGACCGGCGGATGCGCTAAACGTCtttgtctggctacagt                                                                         >  2:500113/1‑225 (MQ=14)

                                                                                                                                                                     tGCCGGATGCGGTGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATC‑TGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACAttt                                                                        <  1:832798/239‑1 (MQ=255)

                                                                                                                                                                     tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATC‑TGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACAttt                                                                        <  2:405362/239‑1 (MQ=255)

                                                                                                                                                                              cggcgTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATC‑TGCACACGACTGCCGGATGCGCTAAACGTCTTGTCCGGCCTACATTTCGCCCGTAg                                                               >  2:950035/1‑239 (MQ=255)

                                                                                                                                                                                 cgTGACCGCCTTATCCGGCGTACGGGGGGCGGGAGAATTTGTAGGCCTGATAAGACGCG‑GTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATC‑TGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGc                                                             <  2:1132948/238‑1 (MQ=255)

                                                                                                                                                                                 cgTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATC‑TGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGc                                                             <  1:648440/238‑1 (MQ=255)

                                                                                                                                                                                 cgTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATC‑TGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGc                                                             <  1:200620/238‑1 (MQ=255)

                                                                                                                                                                                 cgTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATC‑TGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGc                                                             <  2:1210915/238‑1 (MQ=255)

                                                                                                                                                                                   tGAACGCCGTGTCCGGGCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACGGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATC‑TGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAg                                                           <  1:725918/238‑1 (MQ=255)

                                                                                                                                                                                     aaCGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATC‑TGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCa                                                        >  2:1388302/1‑239 (MQ=255)

                                                                                                                                                                                      aCGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATC‑TGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCa                                                        <  1:633452/238‑1 (MQ=255)

                                                                                                                                                                                          cTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATC‑TGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTaaa                                                   <  2:164176/239‑1 (MQ=255)

                                                                                                                                                                                          cTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATC‑TGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTaa                                                    <  1:886922/238‑1 (MQ=255)

                                                                                                                                                                                             aTCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGgcg                                                                                                                                                                        >  2:1114747/1‑120 (MQ=255)

                                                                                                                                                                                             aTCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGgcg                                                                                                                                                                        <  1:1114747/120‑1 (MQ=255)

                                                                                                                                                                                              tCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCa                                                                                                                              <  1:194486/161‑1 (MQ=38)

                                                                                                                                                                                                  gCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATC‑TGCACACGACTGCCGGATGCGATAAACGTCTTGCCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTg                                           >  1:1349864/1‑239 (MQ=255)

                                                                                                                                                                                                       cGGCTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATC‑TGCACACGACTGCCGGATGCGATAAACGTCTTGTCAGGCCTCCATTTCGCCCGTCGGCACTGATTAACTCGTTCTGATAAc                                      >  2:843225/1‑239 (MQ=17)

                                                                                                                                                                                                        gggTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATC‑TGca                                                                                                                   >  2:194486/1‑161 (MQ=39)

                                                                                                                                                                                                          gTGGCGCGAGAATTTGTAGGACTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCGTAGCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATC‑TGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTAcgc                                    <  2:517876/238‑1 (MQ=255)

                                                                                                                                                                                                            gtcgcgAGAATTTGTAGGCCTGATAAGACGCG‑TTAGAGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATC‑TGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGc                                 <  1:632977/237‑1 (MQ=255)

                                                                                                                                                                                                             gcgtGAGAATTTGTAGGTCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATC‑TGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTAcgcg                                   >  1:314879/1‑236 (MQ=255)

                                                                                                                                                                                                             gcgcGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATC‑TGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTAcgcg                                   >  1:523986/1‑236 (MQ=255)

                                                                                                                                                                                                             gcgcGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATC‑TGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTAcgcg                                   >  1:993165/1‑236 (MQ=255)

                                                                                                                                                                                                             gcgcGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATC‑TGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGc                                 <  2:353616/238‑1 (MQ=255)

                                                                                                                                                                                                             gcgcGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATC‑TGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGc                                 <  2:914907/238‑1 (MQ=255)

                                                                                                                                                                                                             gcgcGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATC‑TGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCCCCCGTAGGCAGTCCTTAACTAGTTCTGGTTAcccg                                   >  2:546547/1‑236 (MQ=255)

                                                                                                                                                                                                              cgcgAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATC‑TGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCt                                <  1:1096597/238‑1 (MQ=255)

                                                                                                                                                                                                                 gagaATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATC‑TGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGc                             <  1:354331/238‑1 (MQ=255)

                                                                                                                                                                                                                 gagaATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATC‑TGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGc                             <  1:506981/238‑1 (MQ=255)

                                                                                                                                                                                                                  agaATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATC‑TGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCAGCTGGCaa                           >  1:1059881/1‑239 (MQ=255)

                                                                                                                                                                                                                               cTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATC‑TGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTAcgc                                    >  2:488953/1‑217 (MQ=255)

                                                                                                                                                                                                                               cTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATC‑TGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTAcgc                                    <  1:488953/217‑1 (MQ=255)

                                                                                                                                                                                                                               cTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATC‑TGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTcgcg              >  1:1310021/1‑239 (MQ=255)

                                                                                                                                                                                                                                  agaaGACGCG‑GTAGCGTCGCAGCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATC‑GGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGAt            <  1:1323081/236‑1 (MQ=255)

                                                                                                                                                                                                                                             ttAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATC‑TGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCAc  >  2:298972/1‑238 (MQ=255)

GGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATC‑TGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATC‑TGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCAC  >  NC_000913/374968‑375439

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 12 ≤` `ATCG/ATCG` `< 14 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 39 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`