Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	4,029,029	A→T	noncoding (139/1542 nt)	BW25113_RS19980 →	16S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	4,029,029	0	A	T	100.0%	28.0 / NA	15	noncoding (139/1542 nt)	BW25113_RS19980	16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (6/9);  total (6/9)

AGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCAC  >  NZ_CP009273/4028971‑4029161                                                           |                                                                                                                                     agCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGTTGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGTACCTTCGGGCCTCtt                                                      <  2:432692/139‑1 (MQ=17)             tGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGTTGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGTACCTTCGGGCCTCTTGCCATCGGAt                                            <  2:39990/137‑1 (MQ=17)                    gAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGTTGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGa                                   <  1:370926/139‑1 (MQ=17)                    gAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGTTGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGa                                   <  2:6859/139‑1 (MQ=17)                         gCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGTTGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGa                              >  2:120280/1‑139 (MQ=17)                               gggTGAGTAATGTCTGGGAAACTGCCTGTTGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCt                        <  1:296670/139‑1 (MQ=17)                                 gTGAGTAATGTCTGGGAAACTGCCTGTTGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGa                                   <  1:153403/126‑1 (MQ=17)                                 gTGAGTAATGTCTGGGAAACTGCCTGTTGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGa                                   >  2:153403/1‑126 (MQ=17)                                     gTAATGTCTGGGAAACTGCCTGTTGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAgg                  <  1:120280/139‑1 (MQ=17)                                     gTAATGTCTGGGAAACTGCCTGTTGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAgg                  <  2:43515/139‑1 (MQ=17)                                           tCTGGGAAACTGCCTGTTGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGcc                                       <  1:482223/112‑1 (MQ=17)                                           tCTGGGAAACTGCCTGTTGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGcc                                       >  2:482223/1‑112 (MQ=17)                                           tCTGGGAAACTGCCTGTTGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTgggg             >  2:216098/1‑138 (MQ=17)                                           tCTGGGAAACTGCCTGTTGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTgggg             >  2:80468/1‑138 (MQ=17)                                                     tGCCTGTTGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCAc  >  2:474707/1‑139 (MQ=17)                                                           |                                                                                                                                     AGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCAC  >  NZ_CP009273/4028971‑4029161

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