Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	3,937,289	A→G	noncoding (226/2904 nt)	BW25113_RS19530 →	23S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	3,937,289	0	A	G	100.0%	30.1 / NA	16	noncoding (226/2904 nt)	BW25113_RS19530	23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (9/7);  total (9/7)

CTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAA  >  NZ_CP009273/3937210‑3937409                                                                                |                                                                                                                         cTATCATTAACTGAATCCATAGGTTAAGGAGGCGAACC‑GGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCGACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTAt                                                              <  2:325637/139‑1 (MQ=9)              aaTCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCGACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTg                                                   >  1:69351/1‑138 (MQ=17)              aaTCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCGACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTg                                                   >  2:465576/1‑138 (MQ=17)                          aaTGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCGACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTCTGGa                                       >  1:98793/1‑138 (MQ=17)                                 cGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCGACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGgcgcg                               >  1:217421/1‑139 (MQ=17)                                 cGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCGACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGgcgcg                               <  2:98793/139‑1 (MQ=17)                                 cGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCGACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGgcgc                                <  1:293083/138‑1 (MQ=17)                                 cGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCGACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGgcgc                                >  2:293083/1‑138 (MQ=17)                                    aCCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCGACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGa                            <  1:263021/139‑1 (MQ=17)                                       gggggAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCGACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATAc                         <  2:217421/139‑1 (MQ=17)                                          ggAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCGACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAgg                      >  1:27373/1‑139 (MQ=17)                                                   aCATCTAAGTACCCCGAGGAAAAGAAATCGACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGcc             <  2:27373/139‑1 (MQ=17)                                                         aaGTACCCCGAGGAAAAGAAATCGACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTa        >  2:334220/1‑138 (MQ=17)                                                               cccGAGGAAAAGAAATCGACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACaaa  >  1:109188/1‑138 (MQ=17)                                                                    ggAAAAGAAATCGACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACa                <  1:97498/119‑1 (MQ=17)                                                                    ggAAAAGAAATCGACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACa                >  2:97498/1‑119 (MQ=17)                                                                                |                                                                                                                         CTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAA  >  NZ_CP009273/3937210‑3937409

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