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	100.0%	51.4 / NA	17	L209Q (CTG→CAG)	ydfU	protein YdfU
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (9/8);  total (9/8)

GCACACGGCTGTGTCTTAACCCAACGTGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGATTTTGTT  >  NC_000913/1643475‑1643892                                                                                                                                                                                                             |                                                                                                                                                                                                                      gcaCACGGCTGTGTCTTAACCCAACGTGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCCTGAAAGATTCCGGAGACTCTGGATCACCGGTGATGCTTACCCCCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGACGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTATGCTCTCCAGGTACGATGTCGCTTTCGCGGTCCACCGAGCAGAtttt                                                                                                                                                                         >  1:463384/1‑251 (MQ=255)                  aaCCCAACGTGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCa                                                                                                                                                        <  1:368654/251‑1 (MQ=255)                         cGTGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAacg                                                                                                                                                 >  2:47334/1‑251 (MQ=255)                         cGTGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAacg                                                                                                                                                 >  2:37006/1‑251 (MQ=255)                                               aGCGGCGACGTTTAGGCAGCTTCCTGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGt                                                                                                                           >  2:340073/1‑251 (MQ=255)                                                        gTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATccacc                                                                                                                  >  1:256727/1‑251 (MQ=255)                                                              gCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTg                                                                                                             >  1:275978/1‑250 (MQ=255)                                                                      aTGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGAc                                                                                                    >  2:121115/1‑251 (MQ=255)                                                                        gcaaGATTCCGGAGACTCTGGATCAACGGTGAGGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTGTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGTTACACAGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGAAGTGCATCCAACACATGATTGCAGACCg                                                                                                  >  2:365234/3‑251 (MQ=255)                                                                                             aTCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCa                                                                             <  1:37006/251‑1 (MQ=255)                                                                                                      gATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAAt                                                                     <  2:112782/250‑1 (MQ=255)                                                                                                         gCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCTGTACTATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCGGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCTCCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTcc                                                                  <  1:30060/250‑1 (MQ=255)                                                                                                               caccGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCt                                                           <  2:284449/251‑1 (MQ=255)                                                                                                                      ttttCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTAcc                                                    >  2:166491/1‑251 (MQ=255)                                                                                                                                      tttGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACg                                     <  2:49863/250‑1 (MQ=255)                                                                                                                                       ttgttgCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACgg                                    <  2:256727/250‑1 (MQ=255)                                                                                                                                                                        aTTTTTTGTGCGTTGTTTCAAGATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGATTttgtt  <  2:15474/251‑1 (MQ=255)                                                                                                                                                                                                             |                                                                                                                                                                                                                      GCACACGGCTGTGTCTTAACCCAACGTGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGATTTTGTT  >  NC_000913/1643475‑1643892

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