Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	21,427	A→G	I7M (ATA→ATG)	ribF →	bifunctional riboflavin kinase/FMN adenylyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	21,427	0	A	G	100.0%	121.1 / NA	35	I7M (ATA→ATG)	ribF	bifunctional riboflavin kinase/FMN adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (16/19);  total (16/19)

TTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCTTGACCGCTGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTAC  >  NC_000913/21294‑21551                                                                                                                                      |                                                                                                                             ttGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCTTGACCGCTGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATGCATAAt                                                                                                                        <  2:625778/140‑1 (MQ=255)       cATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCTTGACCGCTGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATGCATAATCTCAGc                                                                                                                  <  2:164949/140‑1 (MQ=255)         tGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCTTGACCGCTGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATGCATAATCTCAGCCa                                                                                                                >  1:19320/1‑140 (MQ=255)            cAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCTTGACCGCTGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATGCATAATCTCAGCCAGGc                                                                                                             <  1:428499/140‑1 (MQ=255)            cAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCTTGACCGCTGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATGCATAATCTCAGCCAGGc                                                                                                             >  1:759577/1‑140 (MQ=255)                      ccACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCTTGACCGCTGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATGCATAATCTCAGCCAGGCCCCGCaagaa                                                                                                   >  1:44719/1‑140 (MQ=255)                         cttcttTCTGATTTCGGTACTCAATCGCCGGTTAACCTTGACCGCTGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAACTGATACGCGGCATGCATAATCTCAGCCAGGCCCCGCAAGAAggg                                                                                                <  2:774858/140‑1 (MQ=255)                          ttctttctGATTTCGGTACTCAATCGCCGGTTAACCTTGACCGCTGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATGCATAATCTCAGCCAGGCCCCGCAAGAAGGgt                                                                                               <  2:19320/140‑1 (MQ=255)                          ttctttctGATTTCGGTACTCAATCGCCGGTTAACCTTGACCGCTGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATGCATAATCTCAGCCAGGCCCCGCAAGAAGGgt                                                                                               <  2:399285/140‑1 (MQ=255)                          ttctttctGATTTCGGTACTCAATCGCCGGTTAACCTTGACCGCTGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATGCATAATCTCAGCCAGGCCCCGCAAGAAGGgt                                                                                               <  2:44719/140‑1 (MQ=255)                                         gTACTCAATCGCCGGTTAACCTTGACCGCTGTACAAGGTATACTAGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATGCATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTg                                                                                <  1:723380/140‑1 (MQ=255)                                                      ggTTAACCTTGACCGCTGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATGCATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACgg                                                                   >  1:639900/1‑140 (MQ=255)                                                         tAACCTTGACCGCTGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATGCATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGt                                                                >  2:165799/1‑140 (MQ=255)                                                                 aCCGCTGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATGCATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCa                                                             >  1:321727/1‑135 (MQ=255)                                                                 aCCGCTGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATGCATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCa                                                             <  2:321727/135‑1 (MQ=255)                                                                 aCCGCTGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATGCATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCAGcgcg                                                        >  2:763789/1‑140 (MQ=255)                                                                      tGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATGCATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCAt                                                   >  2:161274/1‑140 (MQ=255)                                                                        tACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATGCATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATcg                                                 <  1:763789/140‑1 (MQ=255)                                                                            aGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATGCATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATcgcgcg                                             <  2:421280/140‑1 (MQ=255)                                                                            aGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATGCATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATcgcgcg                                             <  1:161274/140‑1 (MQ=255)                                                                            aGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATGCATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATcgcgcg                                             <  1:764757/140‑1 (MQ=255)                                                                             ggTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATGCATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCgcgcgc                                            <  2:759577/140‑1 (MQ=255)                                                                               tataCTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATGCATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATcgcgcg                                             >  1:683897/1‑137 (MQ=255)                                                                               tataCTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATGCATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATcgcgcg                                             <  2:683897/137‑1 (MQ=255)                                                                               tataCTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATGCATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTg                                          >  2:648289/1‑140 (MQ=255)                                                                                 taCTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATGCATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGtt                                        >  1:101124/1‑140 (MQ=255)                                                                                  aCTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATGCATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTa                                       >  2:753441/1‑140 (MQ=255)                                                                                  aCTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATGCATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTa                                       >  1:235159/1‑140 (MQ=255)                                                                                         cGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATGCATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCt                                >  2:136422/1‑140 (MQ=255)                                                                                            ttttCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATGCATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGc                             <  2:786966/140‑1 (MQ=255)                                                                                                         gAGCCAGACATGAAGCTGATACGCGGCATGCATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGgcg                <  2:80972/140‑1 (MQ=255)                                                                                                                 cATGAAGCTGATTCGCGGCATGCATAATCTCAGCCAGGGCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGgcgcaagcgca        <  2:647531/140‑1 (MQ=255)                                                                                                                      aGCTGATACGCGGCATGCATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTa   >  1:794991/1‑140 (MQ=255)                                                                                                                       gCTGATACGCGGCATGCATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTAc  <  2:794991/140‑1 (MQ=255)                                                                                                                           aTACGCGGCATGCATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGgcgcaagcgcaa       >  1:758649/1‑131 (MQ=255)                                                                                                                           aTACGCGGCATGCATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGgcgcaagcgcaa       <  2:758649/131‑1 (MQ=255)                                                                                                                                      |                                                                                                                             TTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCTTGACCGCTGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTAC  >  NC_000913/21294‑21551

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