BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,643,679	A→T	L209Q (CTG→CAG)	ydfU ←	protein YdfU

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,643,679	0	A	T	98.9%	286.9 / NA	93	L209Q (CTG→CAG)	ydfU	protein YdfU
Reads supporting (aligned to +/- strand): ref base A (1/0); new base T (56/36); total (57/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

ATGCGGATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGTGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATA‑CAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGATTTTGTTGAGCAAGCAGCTCCAGTGATTCCGATGTCTGGT  >  NC_000913/1643444‑1643925

atgCGGATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGTGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACg                                                                                                                                                                                                                                           >  1:363834/1‑250 (MQ=255)

        cgtcgGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGTGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCtt                                                                                                                                                                                                                                  >  1:328149/1‑251 (MQ=255)

        cgtcgGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGTGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATAAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCtt                                                                                                                                                                                                                                  <  2:489359/251‑1 (MQ=255)

                ggCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGTGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGCTGTCGCTTTCGCGGt                                                                                                                                                                                                                           >  2:815980/1‑250 (MQ=255)

                ggCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGTGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGt                                                                                                                                                                                                                           >  1:185889/1‑250 (MQ=255)

                             cagcaCACGGATGTGACTTAACCCAACGTGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGTGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTTGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGAGGCGGGCTGCTCACCCGGTACGATGTCGCTTTCGCGGTACCCCGAGCAggtt                                                                                                                                                                                                             >  2:639683/2‑248 (MQ=255)

                               gcaCACGGCTGTGTCTTAACCCAACGTGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGCtttt                                                                                                                                                                                                           >  1:114297/1‑251 (MQ=255)

                                 acacGGCTGTGTCTTAACCCAACGTGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGCTGCTTACCACCGTCTTTTCCTGTGGTGCTTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCCAGATTTTTTGTGCGTTGTTTCACTATGCTGGTGGCGGTCTGCTCTCCCGGTCCGCTGTCGCTTTCGCGGTACCCCGGGCAGCtttttt                                                                                                                                                                                                         >  2:65472/1‑251 (MQ=255)

                                         gtgtCTTAACCCAACGTGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTgg                                                                                                                                                                                                  <  1:415870/250‑1 (MQ=255)

                                          tgtCTTAACCCAACGTGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGt                                                                                                                                                                                                 >  1:470683/1‑250 (MQ=255)

                                           gtCTTAACCCAACGTGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTa                                                                                                                                                                                                >  1:59722/1‑250 (MQ=255)

                                               tAACCCAACGTGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATcc                                                                                                                                                                                            <  1:525984/250‑1 (MQ=255)

                                                aaCCCAACGTGTGTATTTCTCCTTAACCCAGCGGCGACGGTTAGGCAGCTTCATGAAAGATTCCGGAGACTTTGGATCAACGGTGATGCTTTCCACCGTCTTTTCCTGTGGTGATTTTTGGTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATccc                                                                                                                                                                                           <  2:747779/250‑1 (MQ=255)

                                                 aCCCAACTTGTGTATTTCGTCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAg                                                                                                                                                                                         <  2:617414/251‑1 (MQ=255)

                                                    cAACGTGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACCCCGAGCAGCTTTTTTCCGCTGGTAATCCCAGCg                                                                                                                                                                                       >  1:466139/1‑250 (MQ=255)

                                                     aaCGTGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGa                                                                                                                                                                                      >  1:755887/1‑250 (MQ=255)

                                                     aaCGTGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGa                                                                                                                                                                                      >  2:459748/1‑250 (MQ=255)

                                                     aaCGTGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGAGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGCTTTTTTGTGCGTTGTTTCACTATGCTGGTGGCGGTCTGCTCTCCCGGTACGCTGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGa                                                                                                                                                                                      >  2:620668/1‑250 (MQ=255)

                                                          gtgtATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAAGCCCAGCGAacgac                                                                                                                                                                                 <  2:865889/250‑1 (MQ=255)

                                                            gtATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGAcgcg                                                                                                                                                                              >  1:341888/1‑251 (MQ=255)

                                                                   ccctaatcccAGCGGCGTCGTTTAGGCAGCTTCATGAAAGATTCCAGAGACTCTGGATCAACGGTGATGCTTACCACCTTCTTTTCCTTTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTTCGTTGTTTCAATATGCTGGGGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACTCCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACTCGATA‑CAGc                                                                                                                                                                      <  2:100493/244‑1 (MQ=255)

                                                                    ccTTAACCCAGCGTCTACGTTTAGGCAGCTTCATGAAAGATTCCGGAGTCTCTGGATCAACGGTGATGCTTATCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGGCGCTTTCGCGGTACGCCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATA‑Ca                                                                                                                                                                        <  1:65472/248‑1 (MQ=255)

                                                                      ttAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATA‑CAGCCTc                                                                                                                                                                   >  1:868203/1‑251 (MQ=255)

                                                                       tcaCCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGGTTTTTTCCGCTGGTAATACCAGCGAACGACGCGATA‑CAGcatc                                                                                                                                                                   >  1:177058/3‑247 (MQ=255)

                                                                       tAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATA‑CAGCCt                                                                                                                                                                    <  2:520684/249‑1 (MQ=255)

                                                                       tAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTATGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATA‑CAGCCTc                                                                                                                                                                   <  1:593704/250‑1 (MQ=255)

                                                                        aaCCCAGCGGCGACGTTTAGGCAGCTTCATCAAAGATTCCGGAGACTCTGGCTCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGAGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGATGGTGGCGGGCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGGTTTTTTCAGCTGGTACTCCCCGCGACCGACGCGATaacgcctc                                                                                                                                                                    >  2:556095/1‑242 (MQ=255)

                                                                        aaCCCAGCGGCGACGTGGAGGCCGCTTCATGAAAGATTCCGGAGACTCTGGAGCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTGTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGGACACCGAGCAGATTTTTTCCGCTTGTAATCCCAGCGAACGACGCGATA‑CAGCCTCAg                                                                                                                                                                 <  2:301683/251‑1 (MQ=255)

                                                                          cccAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATA‑CAGCCTCAGGt                                                                                                                                                               >  1:792598/1‑251 (MQ=255)

                                                                            cAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATA‑CAGCCTCAGGt                                                                                                                                                               >  1:176692/1‑249 (MQ=255)

                                                                             aGCGGCGACGTTTAGGCAGCTTCATGACAGATTCCGGAGACTCTGGATCAACGGGGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATA‑CAGCCTCAGGTAGt                                                                                                                                                            <  2:588045/251‑1 (MQ=255)

                                                                              gcggcgACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCACTATGCTGGTGGCGGTGTGCTATCACGGTACGCTGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCAGCAGGTAATCCCAGCGAACGACGCTATA‑AAGACTCAGGTAGTg                                                                                                                                                           >  2:476295/1‑251 (MQ=255)

                                                                                  ggACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTTCCCCCGTCTTTACCTGTGGTTATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGAGTTTTTGTGCGTTGTTTCAAGAGGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATA‑CAGCCTCAGGTAGTGCATc                                                                                                                                                       <  2:862067/250‑1 (MQ=255)

                                                                                  cgACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATA‑CAGCCTCAGGTAGTGCATc                                                                                                                                                       >  2:647703/1‑251 (MQ=255)

                                                                                          aGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATA‑CAGCCTCAGGTAGTGCATCCACCACCt                                                                                                                                               >  1:501355/1‑251 (MQ=255)

                                                                                          aGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATA‑CAGCATCAGGAAGTGCATCCACCCCCt                                                                                                                                               >  2:134884/1‑251 (MQ=255)

                                                                                             cAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATA‑CAGCCTCAGGTAGTGCATCCACCACCTGAt                                                                                                                                            >  1:475914/1‑251 (MQ=255)

                                                                                             cAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATA‑CAGCCTCAGGTAGTGCATCCACCACCTGAt                                                                                                                                            >  1:346096/1‑251 (MQ=255)

                                                                                              agattCATGAAAGATTCCGGAAACTCTGGATCAACGGTGATGCTTACAACAGTCTTTTACTGGGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGAGCCAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGTCGGTCTGCCCTACCGATACGAGGTCAGTTTAGCGGAACACCGAGCAGCTTGTTTCCGCGGGTACTCCCAGCGAACGACGCGATC‑AAGCACCAGGAAGTGCAACCACCCCCTGAtt                                                                                                                                           >  2:733840/4‑251 (MQ=255)

                                                                                               gCTTCATGAAAGATTCCTGAGACTCTGGATCAACGGTGAGGCTTACCCCCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGGACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATA‑CAGCCTCAGGTAGTGCATCCACCACCTGATTg                                                                                                                                          <  2:340187/251‑1 (MQ=255)

                                                                                                cTTCATGAAAGATTCCGGAGACTCTGGATCAACGGGGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGAGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGCCGGTCTGCTCTCCCGGTGCGGAGTCGCTTTCGCGGTACACCGTGCAGATTTTTTCTGCTGTTAATTCCACCGACAGACGCGCTA‑CAGACTCAgggatagccgccaaccaccgattgc                                                                                                                                         >  2:203014/1‑228 (MQ=255)

                                                                                                 ttCATGCAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCACCGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTACCGGTACGATGTCGCTTTCGCGGTACACCGAGAAGATTTTTTCCGCTGGTAATCCCAGCGAACGAAGCGATA‑CAGCCTCAGGTAGTGCATCCACCACCGGATTGc                                                                                                                                         >  2:266076/1‑250 (MQ=255)

                                                                                                       aaaGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATA‑CAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGc                                                                                                                                  >  1:689699/1‑251 (MQ=255)

                                                                                                       aaaGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATA‑CAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGc                                                                                                                                  >  2:293555/1‑251 (MQ=255)

                                                                                                         aGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGGCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATA‑CAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGcca                                                                                                                                <  2:640471/251‑2 (MQ=255)

                                                                                                           aTTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATA‑CAGCCTCAGGTAGTGCATCCACCTCCTGATTGCAGACCGCccac                                                                                                                              <  2:328149/251‑1 (MQ=255)

                                                                                                           aTTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATA‑CAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCccac                                                                                                                              >  1:367154/1‑251 (MQ=255)

                                                                                                               cGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGGGCAGATTTTTTCCGCTGGTAATCCCAGCGACAGACCCGCAA‑CAGCCTCAGGTAGAGCATCCACCACCTGATTGCAGACAGCCCACCAgg                                                                                                                          >  2:92974/1‑251 (MQ=255)

                                                                                                               cGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATA‑CAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAgg                                                                                                                          <  1:620492/251‑1 (MQ=255)

                                                                                                                   gattctGGAACAACGGTGATGCTTACTACCGTCTTTTCCTGTGGTTATTTTTTTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGGGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATA‑CAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATa                                                                                                                       <  2:597561/247‑1 (MQ=255)

                                                                                                                     ctctGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATA‑CAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAAtt                                                                                                                    >  1:255817/1‑251 (MQ=255)

                                                                                                                       ctGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATA‑CAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTc                                                                                                                   >  1:815749/1‑250 (MQ=255)

                                                                                                                         ggATCAACGGTGATTCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATA‑CAGCCGCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAg                                                                                                                 <  2:524312/250‑1 (MQ=255)

                                                                                                                         ggATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATA‑CAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGc                                                                                                                <  1:459748/251‑1 (MQ=255)

                                                                                                                          gATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTAGCTTTCGCGGGACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCAATA‑CAGACTCAGGTCGTGCATCCACCACCGGATTGCAGACCACCCCCCAGGATAATTCAGcc                                                                                                               >  2:658493/1‑251 (MQ=255)

                                                                                                                           atcaccggaggtGCTTACCAACGTCTTTTCCTGGGGTGATTTTTTACGCTGGTGGGCGTAAGGCAAATGTGCCACATTTTTTGTGCGTTGTTTCAATATGCTGGTGGAGGTCTGCCCTCCAGGTACGATGTAGCTTTCGCGGTACAACGGGCAGCTTTTTTCCGCTGGCAATCCCAGCGAAAGACGAGATAACA‑CCTCAGGTCGGGCATCCACACCCTAATAGCAGCCCACCCACCAGGAAAAtgacacc                                                                                                               >  2:529428/12‑244 (MQ=255)

                                                                                                                           aTCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATGTTTTCCGCTGGTAATCCCAGCGAACGACGCGATA‑CAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCa                                                                                                              <  2:501355/251‑1 (MQ=255)

                                                                                                                                  gTGATGCTTTCCTCCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTTGGCGTAAGGCAACGGTGCAAGATTTTTTTTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATA‑CAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAAGTCAGCCAAAGATaa                                                                                                       <  2:523765/251‑1 (MQ=255)

                                                                                                                                         ttACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACCCCGAGCAGATTTTTTCCGCTGGTACTCCCAGCGATCGACGCGATA‑CAGCCTCAGGTAGTGCATCCACCACCTTATTGCAGACCGCCACCCAGGATAAATCAGCCAAAGATAATTCCCGc                                                                                                >  2:531076/1‑251 (MQ=255)

                                                                                                                                          tACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGGGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATA‑CAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGAAAATTCAGCCAAAGATAATTCCCGcc                                                                                               >  2:797145/1‑250 (MQ=255)

                                                                                                                                           accgcctgcTTTTCCGGTGTTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTTCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCGCCCGGTACGATGTCGCTTTCGCGTTACACCGAGCAGAATTTTTTCGCTGTTAATCCCAGCGAACGACGCGATA‑CAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCTCACCAGGATAATTCAGCCAAAGATAATTCCCGCTc                                                                                              <  1:787018/243‑1 (MQ=255)

                                                                                                                                              cccGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATA‑CAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTg                                                                                           <  1:143303/250‑1 (MQ=255)

                                                                                                                                              cccGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATA‑CAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTg                                                                                           <  1:297964/250‑1 (MQ=255)

                                                                                                                                               ccGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGCTTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATA‑CAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGCCCGCCCACCAGGATAATTAAGCCAAAGAAAATTCCCGCTCCTGc                                                                                          >  1:352882/1‑251 (MQ=255)

                                                                                                                                                 gTCTTTTCCTGTGGTGTTTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGTTGTCGCTTTCGCGGTACACCGAGCCGATTTTTTCCGCTTGTAATCCCCGCGAACGACGCGATA‑CAGCCTCAGGTCGTGCATCCCCTACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTACCGCTCCTGCGt                                                                                        >  2:262505/1‑251 (MQ=255)

                                                                                                                                                     tttCCTGTAGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGGGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATA‑CAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCttcgacccc                                                                                    >  2:59956/1‑243 (MQ=255)

                                                                                                                                                        ccTGTGGTTTTTTTTGGTGCTGGTGGGCGTAAGGCACCGGTGCAAGATTTTTTGTGCGTTGTTTAAATATGCTGGTGGCGGGCAGCACTCCCGGTACGATGTCGCTTTCGCGGTACACAGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATA‑CAGCCTCAGGTGATGAATCCACCACCTGATTGCAGAACGCCCACCAGGGTAATTCAACCAAAGATAATTCCCGCTCCTGAGTGCCGCtt                                                                                 >  2:298142/1‑251 (MQ=255)

                                                                                                                                                        ccTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATA‑CAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTAAGCCAAAGATAATTCCCGCTCCTGCGAACCGCtt                                                                                 >  2:336327/1‑251 (MQ=255)

                                                                                                                                                              gagATTTTTGTTGCTGGTTGGCGTAAGGCAACGGTGCAATATTTTTTGTGCGTTGTTTCAATATGCTGGTTGCGGTCTTCTCTCCCGGTACGATGGCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATA‑CAGCCTCAGGTAGTGCATCCACCACCTGATTGCTGGCCGCCCGCCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTTCGTACCGCTTATTGCg                                                                           <  1:447671/249‑1 (MQ=255)

                                                                                                                                                                          ggtggtggGCGGAAGGCAACGGTGCAAGATTTTTTGGGCGTTGTTTCAATATGCTGGTGGCGGTCTGTTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATA‑CAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATGCCCGCGCCTGCGTACCGCTTATTGCGTGACGGATGACg                                                               <  2:59722/249‑1 (MQ=255)

                                                                                                                                                                           cTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACAGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATA‑CAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGAACCGCTTATTGCGTGACGGCTGACGt                                                              >  1:428528/1‑251 (MQ=255)

                                                                                                                                                                           cTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGGACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGACCGACGCGATA‑CAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCCAAGATAATTCCCGCTCCTGCGTACCGCATATTGCGTGCCGGacgaagg                                                              >  2:469943/1‑245 (MQ=255)

                                                                                                                                                                                gggCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATA‑CAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAAt                                                          <  1:645732/250‑1 (MQ=255)

                                                                                                                                                                                  gCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTGTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGATTTCGCGGTACACCGAGCAGATTTTGTCCGCTGGTAATCCCAGCGAACGACGCGATA‑CAGCCTCAGGTAGTGCAGCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCAc                                                       <  1:414185/251‑1 (MQ=255)

                                                                                                                                                                                         gCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATA‑CAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAAGTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCt                                                <  1:522689/251‑1 (MQ=255)

                                                                                                                                                                                                  cAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATA‑CAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGAttt                                       <  1:138783/251‑1 (MQ=255)

                                                                                                                                                                                                   aaGATTTTTTTGGCGTTGTTTCAATATGCTGGTGGCGGGTTGCTCTCCCGGTACGAGGCCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATA‑CAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGAtttt                                      <  2:458232/251‑1 (MQ=255)

                                                                                                                                                                                                     gATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATA‑CAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGATTttgt                                    <  1:293555/251‑1 (MQ=255)

                                                                                                                                                                                                      tttttttGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCACAGCGAACGACGCGATA‑CAGCCTCAGGTAGTGCATCCACCACCTGCTTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGGCAATCACCCATGCTGTCAGATTttgtt                                   >  2:571489/2‑251 (MQ=255)

                                                                                                                                                                                                         ttttGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATA‑CAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGATTttgttg                                  >  1:25031/1‑249 (MQ=255)

                                                                                                                                                                                                               gCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATA‑CAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGATTTTGTTGagcaagc                           >  2:828370/1‑250 (MQ=255)

                                                                                                                                                                                                               gCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATA‑CAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGATTTTGTTGagcaagc                           >  1:194990/1‑250 (MQ=255)

                                                                                                                                                                                                                cGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATA‑CAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGATTTTGTTGagcaagc                           >  1:26862/1‑249 (MQ=255)

                                                                                                                                                                                                                 gttgttTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGAGTTACACCGAGCAGTTTTTTTCCGCTGGTAATCCCAGCGAACGCCGCGATA‑CAGCCTCAGGTAGTGCATCCCCCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTAAATACCCCATGCTGCCAGATTTTTTTGggcaagc                           >  2:590193/1‑248 (MQ=255)

                                                                                                                                                                                                                 ggtgttTCAATATGCTTGTGGCTGTCTGCCCTCCTGTTTCAAGTTCGCTTTCGCGTTACTCCGAGCAGATTTTTTCCTCTGGTAATCCCAGCGAACTCCGCGATA‑CAGCCTCAGGTAGTGCATCCACCACCTTCTTGCAGACCGCCCACCATGATAAGTCAGCCAAAGATAATTCCCGCACCTGCGTTCCGCTTATTGCGTTACGGATGACGTCTATCACCCATGCTGTCAGATTTTTTTGAGCAagcagc                        <  2:174584/249‑1 (MQ=255)

                                                                                                                                                                                                                         aaTATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCCGCGAACGACGCGATA‑CCGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGCCGGATGACGTCAATCACCCATGCTGTCAGATTTTGTTGAGCAAGCAGCTCCAGTgg                >  1:802753/1‑250 (MQ=255)

                                                                                                                                                                                                                          atatGCTGGTGGAGGGTTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATA‑CAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGATTTTGTTGAGCAAGCAGCTCCAGTGAt               <  2:470683/251‑1 (MQ=255)

                                                                                                                                                                                                                           tatGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATA‑CAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGATTTTGTTGAGCAAGCAGCTCCAGTGAtt              >  1:258512/1‑251 (MQ=255)

                                                                                                                                                                                                                             tGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATA‑CAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGATTTTGTTGAGCAAGCAGCTCCAGTGATTcc            >  1:325753/1‑251 (MQ=255)

                                                                                                                                                                                                                              gCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATA‑CAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGATTTTGTTGAGCAAGCAGCTCCAGTGATTCCg           >  1:233666/1‑251 (MQ=255)

                                                                                                                                                                                                                                gggtggCGGGCTGCTCTCCCGGTACGATGGCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATA‑CAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGATTTTGTTGAGCAAGCAGCTCCAGTGATTCCGAt         <  2:253610/250‑1 (MQ=255)

                                                                                                                                                                                                                                 ggtggCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATA‑CAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCTCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGATTTTGTTGAGCAAGCAGCTCCAGTGATTCCGATg        >  1:516898/1‑251 (MQ=255)

                                                                                                                                                                                                                                       ggTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGACCGACGCGATA‑CAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGACAGATTTTGTTGTGCAAGCAGCTCCCGTGATTCCGATGTCTGGt  >  2:301084/1‑251 (MQ=255)

ATGCGGATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGTGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATA‑CAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGATTTTGTTGAGCAAGCAGCTCCAGTGATTCCGATGTCTGGT  >  NC_000913/1643444‑1643925

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