Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,718,414	G→A	P39S (CCA→TCA)	yfiF ←	putative methyltransferase YfiF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,718,414	0	G	A	100.0%	40.0 / NA	16	P39S (CCA→TCA)	yfiF	putative methyltransferase YfiF
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (9/7);  total (9/7)

GGGTTTCTTCCGCACGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCTTTCATTTCATCGTTCATGCTTTTCTCCACCAGCGCTGCGCGAAGCGCGCAGATTACCCGAAGTCCGCGCGGTTCGCCATGATTTCGTACCAAAGCCTGCGACTATCATACCTATTGAATAAAACAGATTGTTGTCTGGAAC  >  NC_000913/2718204‑2718652                                                                                                                                                                                                                   |                                                                                                                                                                                                                                               gggTTTCTTCCGCACGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGAACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGt                                                                                                                                                                                                        >  1:171288/1‑251 (MQ=255)        tgccGCACGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGAACGCCCGCCCCCTTTCCCGGTAcgcn                                                                                                                                                                                                                      <  1:134071/228‑2 (MQ=255)         tcCGCACGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGAACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCan                                                                                                                                                                                                 <  1:174675/250‑2 (MQ=255)                    acgacgCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCCCCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGGTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTGGCCTGAACGCCCGCCCTCTTTCCCGGTACGCGGGTTGTGGGGACGTTTATCAGAATCATCATCAc                                                                                                                                                                                     <  1:177255/251‑1 (MQ=255)                         gCAACNNTTNNNGANNAATAAAACTTTTACCACTGATGCNACCGTGATCGGCCTTTTCCGGCGTCTCATNANCNGGCGNNCGGNAAACCGTGCGNCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCANCNCGGGCGGGACGACGGCCACCGTCAGCACGAGNTTTTCCTGACCGCCCGCCCCCTTTNCCGGNNccngnnnnnnn                                                                                                                                                                                                               >  2:160734/1‑210 (MQ=255)                                         aaTAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTACCAGGCGTCTCATCCACCGGAGACCGGAAAACAGTGCGCCACAGCGAATCATCCCACTTGCGGTCACGGGGATGACTCTGATTATCATCGCGGGCGGGGAGACGGCCAGCGTCAGCACGAGATTTTCATGAACGCCCGCACCCTTGCCAGGTCCGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCgg                                                                                                                                                                                >  1:33672/1‑235 (MQ=255)                                                                                                           ngAAACCGTGCGCAACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGAACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCTTTCATTTCATCGTTCATGCTTTTCTCCACCAGCGCTGCGCGAAgcgcgc                                                                                              >  1:61424/2‑251 (MQ=255)                                                                                                                    gcgcCNNGGNNNATNNTCCCACTTGCGGTCACGGGGCTGNCTCTGTTTGTCATCGCGGGCGGGGCGAAGNCNANCGTCNNCACNAGATTTTCCTNAACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGNANCATCATCACTGCGGACATACATCACTTTGACNTTGCCGCTTTTACCTTTCATTTCANCGTTNNTGNtnnnnnnnn                                                                                                                   >  2:162947/1‑213 (MQ=255)                                                                                                                    gcgcCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCAGCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGAACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTTCCGCTTTTACCTTTCATTTCATCGTTCATGCTTTTCTCCACCAGCGCTgcg                                                                                                        <  1:162947/232‑1 (MQ=255)                                                                                                                    gcgcCACGGCGAATCGTCCCACTTGCGGGCACGGGGCTGACTCTGTTTGTCAGCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGAGTTTCCTGAACGCCCGCCCCCGTTCCCGGTACGCGGGTGGTGGGTACGTTTATCAGAATCAGGATCACGGCGGACATACATCACTTTGACCTTGTCGCGTTTACCATTCATGGCATCGTTCATGCGTTTCTCCACCAGCGCTGCGCGAAGCGCGCAGAttccc                                                                                      <  1:108800/250‑4 (MQ=255)                                                                                                                                                                              gcggggcgACGGCCACCGTCAGCACGAGATTTTCCTGAACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCTTTCATTTCATCGTTCATGCTTTTCTCCACCAGCGCTGCGCGAAGCGCGCAGATTACCCGAAGTCCGCGCGGTTCGCCATGATTTCGTACCAAAGCCTGCGACTATCATACCTAtt                            >  1:64588/1‑250 (MQ=255)                                                                                                                                                                                         nccACCGTCAGCACGAGATTTTCCTGAACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCTTTCATTTCATCGTTCATGCTTTTCTCCACCAGCGCTGCGCGAAGCGCGCAGATTACCCGAAGTCCGCGCGGTTCGCCATGATTTCGTACCAAAGCCTGCGCCTATCATACCTATTGAATaaaaa                   >  1:127612/2‑247 (MQ=255)                                                                                                                                                                                         gCCACCGTCAGCACGAGATTTTCCTGAACGCCCGCCCCCTTTCCCGGTACGGGGGTTGTGGGCACGTTTATCAGCATCCTCATCACTGCGGCCATACATCACCTTGACCTTGCAGCTTTTATCTTTCATTTAACAGTTCATGCTTTTTTCCACCAGCGCTGCGCGAAGCGAGCAGATTACCCGCAGTCCGCGCGGATCGCCAt                                                                >  1:194558/1‑203 (MQ=255)                                                                                                                                                                                              nnnnnnnnnnagnNNNNCCTNAACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTTGGTACGTTTCTCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCTTTCATTTCATCTTTCATGCTTTTCTCCACCAgcgc                                                                                                            <  2:137024/144‑1 (MQ=255)                                                                                                                                                                                                   ncACGAGATTTTCCTGAACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCTTTCATTCCATCGTTCATGCTTTTCTCCACCAGCGCTGCGCGAAGCGCGCAGATTACCCGAAGTCCGCGCGGTTCGCCATGATTTCGTACCAAAGCCTGCGACTATCATACCTATTGAATAAAACAGATTGTTGTCTg      >  1:24868/2‑251 (MQ=255)                                                                                                                                                                                                       gagaTTTTCCTGAACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCTTTCATTTCATCGTTCATGCTTTTCTCCACCAGCGCTGCGCGAAGCGCGCAGATTACCCGAAGTCCGCGCGGTTCGCCATGATTTCGTACCAAAGCCTGCGACTATCATACCTATTGAATAAAACAGATTGTTGTCTGGAAc  <  1:160734/251‑1 (MQ=255)                                                                                                                                                                                                                   |                                                                                                                                                                                                                                               GGGTTTCTTCCGCACGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCTTTCATTTCATCGTTCATGCTTTTCTCCACCAGCGCTGCGCGAAGCGCGCAGATTACCCGAAGTCCGCGCGGTTCGCCATGATTTCGTACCAAAGCCTGCGACTATCATACCTATTGAATAAAACAGATTGTTGTCTGGAAC  >  NC_000913/2718204‑2718652

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