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	96.5%	190.2 / NA	58	L209Q (CTG→CAG)	ydfU	protein YdfU
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base T (24/32);  total (26/32)
Fisher's exact test for biased strand distribution p-value = 1.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

CAGATGATGCGGATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGTGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGAT‑CAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGATTTTGTTGAGCAAGCAGCTCCAGTGATTCCGATGT  >  NC_000913/1643438‑1643920                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                  cAGATGCTGCGGATCGTAGGCTGGCATACCGCAGCAAGAACACGGCTGTGTCTTAACCCAACGTGTGTATTTATCCTTAACCCAGCGGCGACGATTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGAT‑CAACGGTGATGCTTACCACAGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGGCTGCTCTCCCg                                                                                                                                                                                                                                          >  2:525497/1‑251 (MQ=255)        ggCGGATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGTGTGTATGTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTACATGAAAGATTCCGGAGACTCTGGAT‑CAACGGTGAGGCTTACCGCCGTCGTTTCCTGTGGGGATTTTTGTTGCGGGTGGGCGTAAGGCAACGGTGCAAGATGTTTTGTGCGTTGTGTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGAt                                                                                                                                                                                                                                   <  2:568816/250‑1 (MQ=255)               cgtcgGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGTGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGAT‑CAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCtt                                                                                                                                                                                                                            <  1:529572/251‑1 (MQ=255)                  cgGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGTGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGAT‑CAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTGGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTCCGATGTCGCTTTc                                                                                                                                                                                                                          <  1:201791/250‑1 (MQ=255)                      tGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGTGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGAT‑CAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGt                                                                                                                                                                                                                     >  1:502250/1‑251 (MQ=255)                        gCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGTGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGAT‑CAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTac                                                                                                                                                                                                                   >  1:97378/1‑251 (MQ=255)                           tACCGCAGCAAGCACACGACTGTGGCTTAACCCACCGTGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGTT‑CAACGGTGATGCTTACAACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTACGGCAACGGTGCAAGAGTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTacaac                                                                                                                                                                                                                >  2:49209/1‑249 (MQ=255)                                   caagcaCACGGCTGTGTCTTAACCCAACGTGTGTATTTCTCCTGAACCCAGCGGCGTCGTTTAGGCAGCTTCATGAAATATTCCGGAGACTCTGGAT‑CAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTTATTTTTGTTGCTGGTGGTCGTAAGGCAACGTTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCGCTCGGTTCGATGTCGCTTGCGCTGTACTCCGAGCAGa                                                                                                                                                                                                         <  2:225637/250‑1 (MQ=255)                                   caagcaCACGGCTGTGTCGTAACCCAACGTGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGAT‑CAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGAt                                                                                                                                                                                                        <  1:560979/251‑1 (MQ=255)                                       cacaCGACTGTGTCTCAACACAACGTATGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCCGCTTCATGAAAGATTCCGGAGACTCTGGAT‑CAACGGTGCTGCTTACCACCGTCTTTTCCTGTGGCGATTTTTGTTGCTGGTGGGCGTAAGCCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGCTGTCGCTTTCGCGGTACCCCGAGCAGAttttt                                                                                                                                                                                                    >  2:83229/1‑251 (MQ=255)                                           cGGCTGTGTCTTAACCCAACGTGTGTATTTCTCCTTAACCCAGCGTCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGAT‑CAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGGGCAAGATTTTTTGTGCGTTGTTTCAATCTGCTGGTGGCGGTCTGCTCTCCCGGTACGCTGTCGCTTTCGCGGTACACCGAGCCGATTTTTTccc                                                                                                                                                                                                >  1:382189/1‑250 (MQ=255)                                                gtgtCTTAACCCAACGTTTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGAT‑CAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCTTGCAGATTTTTTCCGCTGGt                                                                                                                                                                                           <  2:464184/251‑1 (MQ=255)                                                             aCGTGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGAT‑CAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAac                                                                                                                                                                              <  2:502250/251‑1 (MQ=255)                                                               gtgtgtATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGATACTCTGGAT‑CAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTGCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAacga                                                                                                                                                                            <  2:84305/251‑1 (MQ=255)                                                               gtgtgtATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGAT‑CAACGGTGATGCTTACCACCGTCTTTTACTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGAGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCAGGTACAATGGCGCTTTCGCTGAACAACGATCAGATTTTTTCCGCTGGTAAtacccgcgcgcgc                                                                                                                                                                            >  2:473744/1‑239 (MQ=255)                                                                 gtgtATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGAT‑CAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAacgac                                                                                                                                                                           <  1:449212/250‑1 (MQ=255)                                                                 gtgtATTTCTCCTTAACCCAGCGGCGACGGTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGAT‑CAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTGTTTCCGCTGGTAATCCCAGCGAacgacg                                                                                                                                                                          <  2:420954/251‑1 (MQ=255)                                                                   gtATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGAT‑CAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGAGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGAcgc                                                                                                                                                                         <  1:304481/250‑1 (MQ=255)                                                                               aaCCCAGCGGCGGCGTTTAGGCAGCTTCATGAAAGATTCCAGCGAGTCTGGATACAA‑GGTGCAGCTTACCACTGTCTTTTCCTGTGGGGATTTTTGTTGCTGGTGGGCGCACTGCAACGGTGCACGATATTTTGTGCGTTGTGTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTATTCCCAGCGAACGACGCGATACAGCCTc                                                                                                                                                              >  2:269752/1‑249 (MQ=255)                                                                               aaCCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGAT‑CAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAATGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGGCTGCTCGCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCGGGTAATCCCAGCGAACGACTCGATACAGCCTCa                                                                                                                                                             <  1:471287/250‑1 (MQ=255)                                                                                aCCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGAT‑CAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAgg                                                                                                                                                           >  2:520378/1‑251 (MQ=255)                                                                                 cccAGCGACGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGAT‑CAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGACCGACGCGTTACAGCCTAAGGt                                                                                                                                                          >  2:502245/1‑251 (MQ=255)                                                                                    aGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGAT‑CAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGt                                                                                                                                                       >  1:542231/1‑251 (MQ=255)                                                                                       ggcgACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGAT‑CAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGc                                                                                                                                                     >  2:96538/1‑250 (MQ=255)                                                                                          gACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGAT‑CAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACAcc                                                                                                                                                                                                                <  1:532013/188‑1 (MQ=255)                                                                                          gACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGAT‑CAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACAcc                                                                                                                                                                                                                >  2:532013/1‑188 (MQ=255)                                                                                          gACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGAT‑CAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATcc                                                                                                                                                 >  2:209245/1‑251 (MQ=255)                                                                                             gTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGAT‑CAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATccacc                                                                                                                                              >  1:39924/1‑251 (MQ=255)                                                                                                  ggCAGCTTCATGAAAGATTCCGGAGACTCTGGAT‑CAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTg                                                                                                                                         >  1:146826/1‑251 (MQ=255)                                                                                                   gCAGCTTCATGAAAGATTCCGGAGACTCTGGAT‑CAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTg                                                                                                                                         >  1:454370/1‑250 (MQ=255)                                                                                                          cATGAAAGATTCCGGAGACTCTGGAT‑CAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGc                                                                                                                                    >  1:532707/1‑248 (MQ=255)                                                                                                             gAAAGATTCCGGAGACTCTGGAT‑CAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGAcc                                                                                                                               <  2:146826/250‑1 (MQ=255)                                                                                                                         cgaCTCTGGAT‑CAACGGTGATGCTTCCCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGACAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGCCCGCCCACCAGGAt                                                                                                                   >  2:258227/2‑250 (MQ=255)                                                                                                                             tctGGAT‑CAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTc                                                                                                              >  1:448437/1‑251 (MQ=255)                                                                                                                                ggAT‑CAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCGCCTGATTGCAGACCGCCCACCAGGATAATTc                                                                                                              <  1:487779/248‑1 (MQ=255)                                                                                                                                          gTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATaa                                                                                                  <  2:532707/251‑1 (MQ=255)                                                                                                                                               gCTTACCATCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAAGCCCAGCGAACGACGCGAGACAGCCGCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTcc                                                                                              <  1:541947/250‑1 (MQ=255)                                                                                                                                                    ccaccGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTc                                                                                         <  1:278601/250‑1 (MQ=255)                                                                                                                                                     caccGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCt                                                                                       <  1:54476/251‑1 (MQ=255)                                                                                                                                                      accGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTg                                                                                      <  1:209245/251‑1 (MQ=255)                                                                                                                                                                ccTGTGGTGATTTTTGTTGCTGGTTGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGGACACCGAGCAGATTTTTTCCGCTGTTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCt                                                                             <  2:307592/250‑1 (MQ=255)                                                                                                                                                                 cTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACCCCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCCCCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTACCGCTCCTGCGTACCGCttt                                                                           >  2:51739/1‑250 (MQ=255)                                                                                                                                                                     ggtgATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGc                                                                       <  1:51739/251‑1 (MQ=255)                                                                                                                                                                        gatgtttgatgcaggtggtgggaaGGCGACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGGTGGTGGCGGTCTGCTCGCCCGGTACGAGGTCGCTTGCGCGGTACACCGGGCAGATTGTTGTCGCTGGTAATCCCAGCGAACGACGCGATATAGCCTCAGGTAGTGCAGTCACCACGGGATTGCAGACCGCCCACGTGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGt                                                                     <  1:240218/228‑1 (MQ=255)                                                                                                                                                                        gATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGTATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTg                                                                    <  1:525497/251‑1 (MQ=255)                                                                                                                                                                               gttgCTGGTGGGCGTACGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGAACACCGAGCAGATTTTTTCCGCTGGTAATCCCATCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAAACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCACTTATTGCGTGACGGa                                                               >  2:383818/1‑249 (MQ=255)                                                                                                                                                                                                    aCGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGACGACGTCAATCACCCATGCTGTc                                        >  2:437235/1‑251 (MQ=255)                                                                                                                                                                                                     tggTGCAAGATTGTTTGTGCGTTGTTTCAATATGCTGGTGTCGGTCTGCTCGCTCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCGGCTGGTAATCCCAGCGAACGTCGCGTTACAGCCGCAGGTAGTGCATCGACCACCTGATTGCAGACCGCCCGCCAGGATAATTTAGCCGAAGATAGTTCCCGGTCCTGCGTACCGCTTAGTGCGTGACGGGTGACGTCAATCACCCATGCTGTc                                        <  2:382189/249‑1 (MQ=255)                                                                                                                                                                                                        tGCAAGATTTTTTATGACTTGCTTAAACATGATGGCGGCGGCCACATCTCCCGGTACGACATCGCTTTCGCGGAACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGAGAAGGACGCGATACAGCCTCCGGTAGCGCATACAACACCTGATCGCCGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCACTCCTGAGTACCGCTTATTGAGTGACGGATGACGTCAATCACCCATGCTGTc                                        >  2:566517/1‑247 (MQ=255)                                                                                                                                                                                                         gCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGAt                                    <  1:520378/250‑1 (MQ=255)                                                                                                                                                                                                            aGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCTGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGATTttg                                <  1:383818/251‑1 (MQ=255)                                                                                                                                                                                                                  tttGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGATTTTGTTGag                           <  2:97378/250‑1 (MQ=255)                                                                                                                                                                                                                     gtgCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGATTTTGTTGagcaag                       <  1:473019/251‑1 (MQ=255)                                                                                                                                                                                                                     gtgCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGATTTTGTTGagcaa                        <  1:187379/250‑1 (MQ=255)                                                                                                                                                                                                                        cGTTGTTTCAATATGCTGGTGGCGGTCTGCTCGCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTTCAGACCGCCCACCAGGATAAGTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGTATGACGTCAATCACCCATGCTGTCAGATTTTGTTGAGCAagcag                    <  1:262360/251‑1 (MQ=255)                                                                                                                                                                                                                                    atGCTGGCGGCGGACAGCTCTACCAGTACAATGTCGCTTTCGCGACAAACCAAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACCCGCCACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCCGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCACTCACCCATTCTGTCAGATTTTGTTGAGCAAGCAGATCCAGTGATTc        >  2:326467/1‑251 (MQ=255)                                                                                                                                                                                                                                       cTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGATTTTGTTGAGCAAGCAGCTCCAGTGATTCCg      >  2:156199/1‑250 (MQ=255)                                                                                                                                                                                                                                          gtggCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGGACACCGAGCAGATTTGTTCCGCTGGTAATCCCAGAGAACGACTCGATACAGCCTCAGGTAGTGCATCCACCTCCTGATTGCATACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGCTTTTGTTGAGCAAGCAGCTCCAGTGATTCCGATGt  <  1:45332/251‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                  CAGATGATGCGGATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGTGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGAT‑CAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGATTTTGTTGAGCAAGCAGCTCCAGTGATTCCGATGT  >  NC_000913/1643438‑1643920

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