Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,171,788	G→T	intergenic (+32/‑269)	rrfB → / → murB	5S ribosomal RNA/UDP‑N‑acetylenolpyruvoylglucosamine reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,171,788	0	G	T	100.0%	18.7 / NA	13	intergenic (+32/‑269)	rrfB/murB	5S ribosomal RNA/UDP‑N‑acetylenolpyruvoylglucosamine reductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (9/4);  total (9/4)

TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATAC  >  NC_000913/4171570‑4172019                                                                                                                                                                                                                           |                                                                                                                                                                                                                                        tGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTTGGCCTTTCGTTTTATCTGTTGTTTGTCGGTg                                                                                                                                                                                                          >  1:15577/1‑250 (MQ=18)       aagaTTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTTGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGctc                                                                                                                                                                                                   >  2:81209/1‑251 (MQ=18)         gaTTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTTGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTTGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCtctc                                                                                                                                                                                                 >  2:175377/1‑251 (MQ=18)                     gagaCAGATTAAAGCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCGGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTTGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACa                                                                                                                                                                                     <  2:221312/248‑1 (MQ=16)                              tAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTTGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATccgccg                                                                                                                                                                            >  2:72081/1‑251 (MQ=18)                                 aTCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTTGGCCTTTCGTTTTATCTGTTGTTGTTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCggg                                                                                                                                                                          >  1:210216/1‑250 (MQ=18)                                  tCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTTGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCggg                                                                                                                                                                          >  1:212322/1‑249 (MQ=18)                                                                                                                      aaaggaaaccgccgaagcgcggggggTAGTGGGGGGTCCCCCCATGCGATAGTAGGAAACTGCCAGCCATAAAAAAAAACAAAAGGCTAAGTCAAAAGATTTGGCTTTTTGTTTTTTGTGTTGTTTTTCGGTGAACCCTCTCGGGAGTAGAACAACCCCGCGGGGGGGGTATTTAACCTTTGCGAAGAAACGGCCTGGTGGGTGGGGGGCAGAACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGa                                                                                    <  1:246060/227‑1 (MQ=2)                                                                                                                        aGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTTGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAATTCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGa                                                                                    >  1:4571/1‑249 (MQ=18)                                                                                                                               cgcggtggcgcAGATGGTAGGGGGGGGCCCCCCCATGCGAAAGAAGGGAACTGCCAGGCAAAAAAAAAAAAGAAAGGCTCAGTCGAAAGACTTGGCTTTTCGTTTTATCTGTTGTTTGGCGGTGAAGGCTCTCGGGAGTAGAAAAAATCCGCGGGGAGCGGATTTGAACGTGGCGAAGCAGCGGCCGGGGGGGTGGGGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATcc                                                                           <  1:156975/244‑1 (MQ=2)                                                                                                                                    gaGCGCGGAGGGTAGTGGGGGGTCCCCCCATGGGAGAGGAGGAAACTGCCAGGCAAAAAATAAAACGAAGGGCTCAGTCAAAAGACTTGGCTTTTCTTTTTATTTTTTGTTTGTGGGTGAACGCTCTCTGGAGGAGGAAAAACCCGCGGGGAGGGGATTTGAGGGTTGCAAAGAACCGGCGGGGGGGGTGGGGGGCAGGCCGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACg                                                                      <  1:66480/250‑1 (MQ=2)                                                                                                                                                                                                      gAAAGGCTCAGTCGAAAGACTTGGCGTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGCGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGCGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCACCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATcccccaaccgat    >  2:239793/1‑244 (MQ=17)                                                                                                                                                                                                        aaGGCTCAGTCGAAAGACTTGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATAc  >  1:68531/1‑251 (MQ=255)                                                                                                                                                                                                                           |                                                                                                                                                                                                                                        TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATAC  >  NC_000913/4171570‑4172019

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