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%	36.8 / NA	19	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 (8/11);  total (8/11)

CCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATG‑CGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGAC‑TGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCC‑GCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTAC  >  NC_000913/4171541‑4171998                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                    ccTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCTGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATG‑CGAGAGTAGGGAACTGCCAGGCATCAAATAAAACTAAAGGCTCAGTCGAAAGAC‑TTGGc                                                                                                                                                                                                                  <  2:21098/251‑1 (MQ=255)         ccacGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGACCGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGTCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATG‑CGAGATTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGAC‑TTGTCCTGTCGtt                                                                                                                                                                                                          <  2:279241/248‑1 (MQ=11)            cGAAGCTGTTTTGGCTGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATG‑CGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGAC‑TTGGCCTTTCGTTTTa                                                                                                                                                                                                       <  2:255750/251‑1 (MQ=37)                 cTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGCTGGTAGTGTGGGGTCTCCCCATG‑CGAGCGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGAC‑TTGGCCTTTCGTTTTATCtgt                                                                                                                                                                                                  >  1:145290/1‑251 (MQ=255)                           ggATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGGAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATG‑CGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGAC‑TTGGCCTTTCGTTTTATCTGTTGTTTGTCgg                                                                                                                                                                                        <  1:167568/251‑1 (MQ=21)                             aTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATG‑CGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGAC‑TTGGCCTTTCGTTTTATCTGTTGTTTGTCGGTg                                                                                                                                                                                      <  2:87840/251‑1 (MQ=18)                               gagagaAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATG‑CGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGAC‑TTGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGaa                                                                                                                                                                                    <  1:63451/251‑1 (MQ=18)                                         tttCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATG‑CGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGAC‑TTGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTg                                                                                                                                                                          >  1:270733/1‑251 (MQ=18)                                                              aTCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATG‑CGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGAC‑TTGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGa                                                                                                                                                     <  2:98846/251‑1 (MQ=18)                                                                                       aaaCAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATG‑CGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGAC‑TTGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAAc                                                                                                                            >  2:88686/1‑251 (MQ=18)                                                                                                  gCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGAACCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGAGGGGTCTCCCCATG‑CGAGAGTAGGGAACTGCCAGGCATCCAATAAAACGAAAGGCTCAGTCGAAAGAC‑TTGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAggg                                                                                                                 >  1:168401/1‑251 (MQ=16)                                                                                                                           ccACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATG‑CGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGAC‑TTGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGGGGCGGGCAGGACGCCC‑GCCATaa                                                                                        >  2:216662/1‑250 (MQ=18)                                                                                                                                                             gCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATG‑CGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGAC‑TTGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGt                                                                                                                >  2:207819/1‑193 (MQ=21)                                                                                                                                                             gCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATG‑CGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGAC‑TTGGCCTTTCGTTGTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGt                                                                                                                <  1:207819/193‑1 (MQ=21)                                                                                                                                                                  ggctccaatggtatggtgGGGTC‑CCCCATGGAGAGGGTAGGGAACTGCGAGGCAAAAAAAAAAACAAAGGGCTCAATCGAA‑GACCTTGGCTTTTCTTTTTTTCTTGTGTTTGTGGGTGAACCCTCTCCGGGGTGGAACAAAGCCGCGGGGGGCGGATTTGAACGTTGCAAAAAACCGGCCGGGGGGGGGGGGGGC‑GGACCCCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACg                                                 <  1:15356/235‑1 (MQ=2)                                                                                                                                                                          gggtggtgggggttccccccAGG‑CGAGAGTAGGAAACGCCCAGACAAAAAAAAAAACAAAAGGCTCAGACGAAAGAT‑TTGCCCTTTCTTTTTCTTTGTTGTTTGGGGGTAACCGCCCGCCGGATAAGCACAAACCCGCCGGGAGGGTATTGAAACGTTACGAAGCACCGGCCGGGAGGGTGGGGGGCAGGACGCCC‑GCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCtt                                        <  1:272850/235‑1 (MQ=2)                                                                                                                                                                                        ctcCCCATG‑CGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGAC‑TTGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACTCCC‑GCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGcc                                          >  1:129418/1‑235 (MQ=18)                                                                                                                                                                                                  cGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGACAGCC‑TTGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGCGGCGGGCACGACGCAC‑CCCATACACCGCCAGGCATCGAATTAAGCAGAAGGCCACACTAACGGATGCCCTTTTTACGGTTCTACAACctctccc                 >  2:275978/1‑248 (MQ=11)                                                                                                                                                                                                                 gCCAGGCAAAAAAAAAAACGAGAGGCTCAGCCGAAAGAT‑TTGGCTTTTGGTTTTATTTTTTGTTGGTCGGTGACCGCCCTCCGGAGTAGAACAACCCCGCGGGGGGCGGTTTTGAACTTTGCGAAGCAGCGGCCGGGGGGGTGGGGGGCAGGACGCCC‑GCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTAc  <  1:216662/251‑1 (MQ=9)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                    CCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATG‑CGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGAC‑TGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCC‑GCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTAC  >  NC_000913/4171541‑4171998

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