Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,212,973	Δ1 bp	intergenic (+27/‑67)	rrlE → / → rrfE	23S ribosomal RNA/5S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,212,973	0	T	.	100.0%	19.7 / NA	8	intergenic (+27/‑67)	rrlE/rrfE	23S ribosomal RNA/5S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (2/6);  total (2/6)

AATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATT  >  NC_000913/4212924‑4213210                                                  |                                                                                                                                                                                                                                              aaTGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAAcc                                      <  2:15745/250‑1 (MQ=11)    gAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCg                                   <  1:49293/250‑1 (MQ=14)           gAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAGGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGGGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGAGAGAGTAGGGAACTGCCAGGCATCAAATTAGCAAAACCCCGGTCCATa                            >  2:16978/1‑250 (MQ=11)                  aaCCTTACAACGCCGAAGCTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCgg                     <  2:279619/250‑1 (MQ=255)                              ccGAAGTTGGTTTGGCGGAT‑GAGAGAAGATGTTCAGCCGGATACATAGTAAAGCAGAACGCAGAATCGGTATGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCGGACCCCATGCCGAACTCAGAAGTGAGACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGATTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTTCATAATGCCGGGGtgttcggcata        <  2:322915/251‑11 (MQ=11)                              ccGAAGCTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCat         >  1:253476/1‑250 (MQ=255)                               cGAAGCTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCTGGGTTTTTTGCatat       <  2:141060/251‑1 (MQ=255)                                      gtTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAatt  <  1:159652/249‑1 (MQ=255)                                                  |                                                                                                                                                                                                                                              AATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATT  >  NC_000913/4212924‑4213210

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