Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,039,385	C→T	noncoding (1867/2905 nt)	rrlA →	23S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,039,385	0	C	T	100.0%	16.2 / NA	11	noncoding (1867/2905 nt)	rrlA	23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (6/5);  total (6/5)

CTTGAGAGAACTCGGGTGAAGGAACTAGGCAAAATGGTGCCGTAACTTCGGGAGAAGGCACGCTGATATGTAGGTGAAGCGACTTGCTCGTGGAGCTGAAATCAGTCGAAGATACCAGCTGGCTGCAACTGTTTATTAAAAACACAGCACTGTGCAAACACGAAAGTGGACGTATACGGTGTGACGCCTGCCCGGTGCCGGAAGGTTAATTGATGGGGTTAGCCGCAAGGCGAAGCTCTTGATCGAAGCCCCGGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCGGGTAAGTTCCGACCTGCACGAATGGCGTAATGATGGCCAGGCTGTCTCCACCCGAGACTCAGTGAAATTGAACTCGCTGTGAAGATGCAGTGTACCCGCGGCAAGACGGAAAGACCCCGTGAACCTTTACTATAGC  >  NC_000913/4039163‑4039603                                                                                                                                                                                                                               |                                                                                                                                                                                                                           cTTGAGAGAACTCGGGTGAAGGAACTAGGCAAAATGGCGCCGTAACTTAGGGAGAAGGCACGCTGATATGCAGGCGAAGCGACTTGCTCGTGGAGCTGAACTCAGTCGAAGATACCAGCTGGCTGCAACTGTTTATTAAAAACACAGCACTGTGCAAACACGCACGAGGAAGTATAAGGTGTGACGCCTGCCCGGTGCCGGACGGATAATTGATGGGGTTAGTCGCAAGGCGAAGCTCTTGATCGAAGccc                                                                                                                                                                                                >  1:102689/1‑251 (MQ=11)                                                                 gATATGTAGGTGAAGCGACTTGCTCGTGGAGCTGAAAGCAGTCGAAGATACCAGCTGTCTGCAACTGTTTATTAAAAACACAGCACTGTGCAAACGCGAAACTGGACGTATACGGTGTGACGCCTGCCCGGTGCCGGAAGGTTAATTGATGGGGTTAGTCGCAAGGCGAAGCTCTTGATCGAAGCCCCGGTAAACGGCGGCCGTAACTAGAACGGTCCTAAGGTAGCGAAATTCCTTGTCGGGTAAGTTcc                                                                                                                                <  1:126986/251‑1 (MQ=255)                                                                      gTAGGTGAAGCGACTTGCTCGTGGAGCTGAAATCAGTCGAAGATACCAGCTGGCTGCAACTGTTTATTAAAAACACAGCACTGTGCAAACACGAAAGTGGACGTATACGGTGTGACGCCTGCCCGGTGCCGGAAGGTTAATAGATGGGGTTAGTCGCAAGGCGAAGCTCTTGATCGAAGCCCCGGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCGGGTAAGTTCCGACCt                                                                                                                           >  1:67155/1‑251 (MQ=18)                                                                            gAAGCGACTTGCTCGTGGAGCTGAAATCAGTCGAAGATACCAGCTGGCTGCAACTGTTTATTAAAAACACAGCACTGTGCAAACACGAAAGTGGACGTATACGGTGTGACGCCTGCCCGGTGCCGGGAGGTTAATTGATGGGGTTAGTCGCAAGGCGAAGCTCTTGATCGAAGCCCCGGTAAACGGCGGCCGTAAc                                                                                                                                                                            >  1:8460/1‑196 (MQ=21)                                                                                                         gTCGAAGATACCAGCTGGCTGCAACTGTTTATTAAAAACACAGCACTGTGCAAACACGAAAGTGGACGTATACGGTGTGACGCCTGCCCGGTGCCGGAAGGTTAATTGATGGGGTTAGTCGCAAGGCAAAGCTCTTGATCGAAGCCCCGGTAAACGGCGGCCGTAACTATAAc                                                                                                                                                                      <  1:14267/173‑1 (MQ=21)                                                                                                                   ncAGCTGGCTGCAACTGTTTATTAAAAACACAGCACTGTGCAAACACGAAAGTGGACGTATACGGTGTGACGCCTGCCCGGTGCCGGAAGGTTAATTGATGGGGTTAGTCGCAAGGCGAAGCTCTTGATCGAAGCCCCGGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCGGGTAAGTTCCGACCTGCACGAAAGGCGTAATGATGGCCAGGCTGTCTCCACCCGAGACTc                                                                              >  1:39044/2‑251 (MQ=18)                                                                                                                                                      cTGTGCAAACACGAAAGTGGACGTATACGGGGTGACGCCTGCCCGGTGCCGGAAGGTTAATTGATGGGGTTAGTCGCAAGGCGAAg                                                                                                                                                                                                                <  1:67525/86‑1 (MQ=255)                                                                                                                                                                            gTATACGGTGTGACGCCTGCCCGGTGCCGGAAGGTTAATTGATGGGGTTAGTCGCAAGGCGAAGCTCTTGATCGAAGCCCCGGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCGGGTAAGTTCCGACCTGCACGAATGGCGTAATGATgg                                                                                                      <  1:125055/170‑1 (MQ=25)                                                                                                                                                                                     gtgACGCCTGCCCGGTGCCGGAAGGTTAATTGATGGGGTTAGTCGCAAGGCGAAGCTCTTGATCGAAGCCCCGGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCGGGTAAGTTCCGACCTGCACGAATGGCGTAATGATGGCCAGGCTGTCTcc                                                                                         >  1:88650/1‑174 (MQ=21)                                                                                                                                                                                          nccTGACCGGTGCCGGAAGGTTAATTGATGGGGTTAGTCGCAAGGCGAAGCTCTTGATCGAAGCCCCGGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCGGGTAAGTTCCGACCTGCACGAATGGCGTAATGATGGCCAGGCTGTCTCCACCCGAGACTCAGTGAAATTGAACTCGCTGTGAAGATGCAGTGTACCCGCGGCAAGACGGAAAGACCCCGTGAACCTTTACt       >  1:176500/2‑251 (MQ=18)                                                                                                                                                                                               cccGGTGCCGGAAGGTTAATTGATGGGGTTAGTCGCAAGGCGAAGCTCTTGATCGAAGCCCCGGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAAGTCCTTGTCGGGTAAGGGCCGACCTGCACGAATGGCGTAATGATGGCCAGGCTGTCTCCACCCGATACGCAGTGAGATTGAGCGCGCTGTGAAGAGGCAGTGTACCCGCGGCAAGACGGAAAGACCCCGTGAACCTTTACTATAGc  <  1:208311/251‑1 (MQ=11)                                                                                                                                                                                                                               |                                                                                                                                                                                                                           CTTGAGAGAACTCGGGTGAAGGAACTAGGCAAAATGGTGCCGTAACTTCGGGAGAAGGCACGCTGATATGTAGGTGAAGCGACTTGCTCGTGGAGCTGAAATCAGTCGAAGATACCAGCTGGCTGCAACTGTTTATTAAAAACACAGCACTGTGCAAACACGAAAGTGGACGTATACGGTGTGACGCCTGCCCGGTGCCGGAAGGTTAATTGATGGGGTTAGCCGCAAGGCGAAGCTCTTGATCGAAGCCCCGGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCGGGTAAGTTCCGACCTGCACGAATGGCGTAATGATGGCCAGGCTGTCTCCACCCGAGACTCAGTGAAATTGAACTCGCTGTGAAGATGCAGTGTACCCGCGGCAAGACGGAAAGACCCCGTGAACCTTTACTATAGC  >  NC_000913/4039163‑4039603

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