Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	21,581	C→T	P59S (CCA→TCA)	ribF →	bifunctional riboflavin kinase/FAD synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	21,581	0	C	T	100.0%	50.1 / NA	16	P59S (CCA→TCA)	ribF	bifunctional riboflavin kinase/FAD synthetase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (6/10);  total (6/10)

AGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTC  >  NC_000913/21454‑21703                                                                                                                                |                                                                                                                           agGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAATCACTGGAACTg                                                                                                                 >  2:421409/1‑139 (MQ=255)                 aTTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAATCACTGGAACTGTTTGCTACCGATAAAg                                                                                                 >  1:401384/1‑139 (MQ=255)                        aTTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAATCACTGGAACTGTTTGCTACCGATAAAGCCCCGGc                                                                                          >  1:481273/1‑139 (MQ=255)                                                aTCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAATCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGaaaa                                                                  <  2:290532/139‑1 (MQ=255)                                                      cgcTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAATCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCg                                                            <  1:405228/139‑1 (MQ=255)                                                      cgcTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAATCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCg                                                            <  2:234300/139‑1 (MQ=255)                                                         tGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAATCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCGTTa                                                         <  2:481273/139‑1 (MQ=255)                                                              cAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAATCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCGTTACCTTg                                                    <  2:12699/139‑1 (MQ=255)                                                                   cTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAATCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCagag                                               >  2:107078/1‑139 (MQ=255)                                                                   cTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAATCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCa                                                  <  1:389092/136‑1 (MQ=255)                                                                   cTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAATCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCa                                                  >  2:389092/1‑136 (MQ=255)                                                                           aagaagGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAATCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGt                                       <  1:421409/139‑1 (MQ=255)                                                                             gaagGGCGCAAGCTCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAATCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTg                                     <  2:225290/139‑1 (MQ=255)                                                                                                        atggtgatgCTTTTTGAACCTCAATCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACa          >  2:322970/1‑139 (MQ=255)                                                                                                                gCTTTTTGAACCTCAATCACTGGAACTGTTTGCTACCGATTAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTc  <  1:322970/139‑1 (MQ=255)                                                                                                                gCTTTTTGAACCTCAATCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTc  <  1:21018/139‑1 (MQ=255)                                                                                                                                |                                                                                                                           AGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTC  >  NC_000913/21454‑21703

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