Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	221,335	T→C	noncoding (1078/1542 nt)	BW25113_RS01010 →	16S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	221,335	0	T	C	100.0%	63.1 / NA	22	noncoding (1078/1542 nt)	BW25113_RS01010	16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (14/8);  total (14/8)

TTGACATCCACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGAC  >  NZ_CP009273/221248‑221415                                                                                        |                                                                                 ttGACATCCACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGa                                                                                <  1:456092/90‑1 (MQ=255)      aTCCACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGt                                                                           <  2:306760/90‑1 (MQ=255)        ccACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCGTGGTTGTCGTCAGCTCGTGGTGCGCAATGGTg                                                                         >  2:505322/1‑90 (MQ=18)               aCTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAg                                                                  >  2:506874/1‑90 (MQ=255)                         gaTGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACg                                                        >  1:263666/1‑90 (MQ=255)                         gaTGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATAGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACg                                                        >  1:205825/1‑90 (MQ=255)                               ttGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCa                                                  <  2:675230/90‑1 (MQ=255)                               ttGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCa                                                  <  2:348439/90‑1 (MQ=255)                               ttGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGCTGGGTTAAGTCCCGCAACGAGCGCa                                                  <  1:220527/90‑1 (MQ=255)                                 ggTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAAc                                                >  2:372629/1‑90 (MQ=37)                                   tGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAAccc                                              >  1:119252/1‑90 (MQ=25)                                       ttCGGGAACTGTGAGACAGGGGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTAt                                          >  1:346517/1‑90 (MQ=255)                                              aCTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTttgt                                   >  2:755547/1‑90 (MQ=255)                                                           tgctgcATGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTccg                      >  2:426571/1‑90 (MQ=18)                                                                 aTGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCggg                <  2:732958/90‑1 (MQ=25)                                                                 aTGGCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCggg                <  1:107737/90‑1 (MQ=25)                                                                    gCTGTCGTCAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAAc             >  2:12914/1‑90 (MQ=18)                                                                       gtcgtcAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCa          >  1:536689/1‑90 (MQ=18)                                                                         cgtcAGCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCaaa        >  1:488844/1‑90 (MQ=18)                                                                              gCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGgaga   >  1:650466/1‑90 (MQ=18)                                                                              gCTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGgaga   >  1:450932/1‑90 (MQ=18)                                                                               cTCGTGTTGCGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGAc  <  2:307602/90‑1 (MQ=18)                                                                                        |                                                                                 TTGACATCCACAGAACTTTCCAGAGATGGATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGAC  >  NZ_CP009273/221248‑221415

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