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	97.8%	141.4 / NA	45	L209Q (CTG→CAG)	ydfU	protein YdfU
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (20/24);  minor base G (1/0);  total (21/24)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGATGATGCGGATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGTGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGATTTTGTTGAGCAAGCAGCTCCAGTGATTCCGATGTCT  >  NC_000913/1643439‑1643922                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    aGATGATGCGGATCGTCGGCGGGCATACCGCCACAAGCACACGGCTGTGTCTTAACCCAACGTGTGTATTTCTCCTTACCCCAGCGGCGACGTTTAGACAGCATCATGAAAGATTCCGGTGACTATGGATCAACGGTGATGCTTACCACAGTATTTTCCGGTGGTGATTTTTGTTGCTGGTGGGCAAAAGGCAACGGTGCAAAATTTTTTGTGCGATGTTTCCGTCTGCTGGTGGCGGTCGGCTCTCCCgg                                                                                                                                                                                                                                           >  2:383893/1‑251 (MQ=255)     gatgCGGATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAGCCCAACGTGTGTATGTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCATCGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTAc                                                                                                                                                                                                                                        <  2:425371/250‑1 (MQ=255)               gtcgGATGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGTGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGAGTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTTCTCTCCCGGTACGATGTCGCttt                                                                                                                                                                                                                             <  2:194890/251‑1 (MQ=255)                     cgtcattccGCATCAAGCACACGGCTGATTTTTAACCCAACTTGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGGCTCGTGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGGTGCTGGTGGGCGTAATGCAACGGTGCAAGATTTTTTGTGCTTTGTTTCAATATGCTGGTGGCGGTCTGCGCTCCCGGTACGATGTCGCTTTCGCGGt                                                                                                                                                                                                                       <  2:17244/244‑1 (MQ=255)                        cATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGTGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTaca                                                                                                                                                                                                                    <  2:89799/251‑1 (MQ=255)                        cATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGTGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTaca                                                                                                                                                                                                                    <  1:491133/251‑1 (MQ=255)                                  caagcaCACGGCTGTGTCTTAACCCAACGTGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGGGCAGAt                                                                                                                                                                                                          <  1:148327/251‑1 (MQ=255)                                  caagcaCACGGCTGTGTCTTAACCCAACGTGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGc                                                                                                                                                                                                              >  1:93700/1‑247 (MQ=255)                                      cacaCGGCTGTGTCTTAACCCAAAGTGTGAATTTCTCCTTACCCCAGCGGCGCCGTTTAGACAGCTTCATGAAAGATTCCGGAAACTCTGGATCAACGGTGATGCTTACCACCGTCATTTCCTGTGGCGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGCGCCAGATTTTTTGTGCGTTGTTTCAACATGCTGGTGGCGATCTGCTCTACCGGTACGATGTCGCTTTCGCGGTACACCTAGCCGAttttt                                                                                                                                                                                                      >  2:192708/1‑251 (MQ=255)                                               gtgtCTTAACCCAACGTGTGTATTTCTCCTTAACCAAGCGGCGCCGTTTAGGCCGCTTCATGAAAGCTTCCGGCGCCTCTGGCTCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGt                                                                                                                                                                                             >  2:482311/1‑251 (MQ=255)                                                        cccAACGTGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAATGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGc                                                                                                                                                                                    <  1:355685/251‑1 (MQ=255)                                                                gtgtATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAacgacg                                                                                                                                                                            <  1:384341/251‑1 (MQ=255)                                                                           cTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCACCGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCCGATTTTTTCCGCTGGAACTCCCAGCGAACGACGCGCTACAGCCt                                                                                                                                                                 >  2:8904/1‑251 (MQ=255)                                                                              aaCCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGAATCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTc                                                                                                                                                                >  1:340965/1‑249 (MQ=255)                                                                                   aGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAgg                                                                                                                                                         >  1:539981/1‑250 (MQ=255)                                                                                       gcgACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGAGCAACGGTGATGCTTACCTCCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGTTGCAAGAGTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTGTCGCTGTACACCGAGCAGATTGTTTCCGCTGTTAATCCCAGCGAACGACTCTATACAGCCTCAGTTATTGCAt                                                                                                                                                     <  2:308441/251‑1 (MQ=255)                                                                                         gACGTTTAGGAAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACAACCGTCTTTTCCTGTGGTGATTTTTGATGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGATTTCGAGGTACACCGAGCAGATTTTTTCCGCTGGTAATCACAGCGAACGACGCGATACAGCCTCAGGTAGTGCATcc                                                                                                                                                   >  2:212244/1‑251 (MQ=255)                                                                                              ttAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCTCAGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTTGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATccacca                                                                                                                                               <  1:482527/250‑1 (MQ=255)                                                                                                 ggCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGCTTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTg                                                                                                                                           >  2:152009/1‑251 (MQ=255)                                                                                                  gCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGCTGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCAACCACCTg                                                                                                                                           >  1:337292/1‑250 (MQ=255)                                                                                                   cAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGAt                                                                                                                                         <  2:107951/251‑1 (MQ=255)                                                                                                     gCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTg                                                                                                                                       <  2:160846/251‑1 (MQ=255)                                                                                                             aaaGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGc                                                                                                                               <  2:16663/251‑1 (MQ=255)                                                                                                              aaGAGTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCGGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGca                                                                                                                              <  2:210696/251‑2 (MQ=255)                                                                                                                gATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAgcccgcccc                                                                                                                            >  2:264845/1‑250 (MQ=255)                                                                                                                    ccGGAGACGCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCccacca                                                                                                                         <  2:185097/250‑1 (MQ=255)                                                                                                                       gagaCTCTTGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTGCGCGGTACTCCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGGCTCAGGTAGTGCATCCCCCACCTTATTGCAGACCGCCCTCCAGGAt                                                                                                                     <  2:340965/251‑1 (MQ=255)                                                                                                                                          gATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAAt                                                                                                   <  2:212913/250‑1 (MQ=255)                                                                                                                                                   caccGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCCCCAGGATAATTCAGCCAAAGATAATTCCCGCTCCt                                                                                         <  1:212244/251‑1 (MQ=255)                                                                                                                                                         cTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTAc                                                                                   >  2:49615/1‑251 (MQ=255)                                                                                                                                                           tttCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGCTAATTCAGCCAAAGATAATTCCCGCTCCTGCGTAccc                                                                                 >  2:154258/1‑250 (MQ=255)                                                                                                                                                                         ttttGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGAc                                                                    <  2:93700/250‑1 (MQ=255)                                                                                                                                                                           ttgttgCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTGCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGTATAATTCAGCCAAAGATAAGTCCCGCTCCTGCGTACCGCTTATTGCGTGACgg                                                                  <  1:482311/250‑1 (MQ=255)                                                                                                                                                                            tgttgCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGAt                                                                >  1:377549/1‑251 (MQ=255)                                                                                                                                                                                     tggGCGTAAGGCAACGGGGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATCATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCaa                                                        >  2:8532/1‑250 (MQ=255)                                                                                                                                                                                     ggggCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAAt                                                       <  1:174495/250‑1 (MQ=255)                                                                                                                                                                                        gCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTAAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGAAGTCAATCAc                                                    >  2:485886/1‑251 (MQ=255)                                                                                                                                                                                             aGGCAACGGTGCAAGATTTTTTGGGCGTTGTTTCAATATGCAGGTGGCGGTCTGCTCTCCAGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATATTTTCAGCTGGTAATCCCAGCGAACGAAGCGATACAGCCTCAGGTAGTGAATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGCTGACGTCAATCACCCATg                                               >  1:387983/1‑251 (MQ=255)                                                                                                                                                                                                  aCGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTc                                          <  1:31818/251‑1 (MQ=255)                                                                                                                                                                                                         aaGATTTTTTGTGCGTTGTTTCAATGTGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATGCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGAtttt                                   <  1:383893/251‑1 (MQ=255)                                                                                                                                                                                                                           tttCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGATTTTGTTGAGCAAGCAGCTcc                  >  1:206048/1‑250 (MQ=255)                                                                                                                                                                                                                              cAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGATTTTGTTGAGCAAGCAGCTCCAGt               <  2:377549/250‑1 (MQ=255)                                                                                                                                                                                                                              cAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGATTTTGTTGAGCAAGCAGCTCCAGTg              >  1:356703/1‑251 (MQ=255)                                                                                                                                                                                                                                       ggtggCGGTATGCTCTCCCGGAACGATGTCGCTTTCGCGGTACACCGAGCAAATTTTTTCCGCTGGTAATCCAAGCGAACGACGCGATACAGCCTAAGGTAGCGCATCCACCACCTGATTGCAGACCGCACACCCGGATACTTCAGCCACAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCCGCTTTTGTTGAGCAAACATCTCCAGtatttacgc       >  2:373691/1‑241 (MQ=255)                                                                                                                                                                                                                                          ggCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGATTTTGTTGAGCAAGCAGCTCCAGTGATTCCGATGTCt  >  2:245140/1‑251 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    AGATGATGCGGATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGTGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGATTTTGTTGAGCAAGCAGCTCCAGTGATTCCGATGTCT  >  NC_000913/1643439‑1643922

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