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%	39.8 / NA	19	noncoding (226/2904 nt)	BW25113_RS19530	23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (9/10);  total (9/10)

ACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATG  >  NZ_CP009273/3937219‑3937413                                                                       |                                                                                                                             aCTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCGACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATgtgtgt                                                            >  1:460896/1‑137 (MQ=17) aCTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCGACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATgtgtgt                                                            >  1:603936/1‑137 (MQ=17)   tGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCGACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGt                                                        >  2:452441/1‑139 (MQ=17)     aaTCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCGACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTg                                                       >  2:113215/1‑138 (MQ=17)     aaTCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACA‑‑TAAGTACCCCGAGGAAAAGAAATCGACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTgg                                                      >  2:129975/1‑137 (MQ=255)      aTCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCGACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGa                                                     <  2:433004/139‑1 (MQ=17)       tCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCGACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGaa                                                    <  1:368597/139‑1 (MQ=17)                   tGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCGACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAg                                        <  1:113215/139‑1 (MQ=17)                   tGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCGACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAg                                        >  2:108840/1‑139 (MQ=17)                        cGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCGACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGgcgcg                                   <  2:460896/139‑1 (MQ=17)                        cGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCGACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGgcgcg                                   <  1:411871/139‑1 (MQ=17)                        cGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCGACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGGCTGGAAAGgcgcg                                   >  2:24134/1‑139 (MQ=17)                                gggAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCGACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAg                           <  2:12799/139‑1 (MQ=17)                                                     ccccGAGGAAAAGAAATCGACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACaa       <  1:280051/138‑1 (MQ=17)                                                       ccGAGGAAAAGAAATCGACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACaaaaa    <  1:76041/139‑1 (MQ=17)                                                        cGAGGAAAAGAAATCGACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAAt   >  1:101964/1‑139 (MQ=17)                                                         gAGGAAAAGAAATCGACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATg  <  1:129975/139‑1 (MQ=17)                                                           ggAAAAGAAATCGACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGgcgcg                                   >  1:462562/1‑104 (MQ=17)                                                           ggAAAAGAAATCGACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGgcgcg                                   <  2:462562/104‑1 (MQ=17)                                                                       |                                                                                                                             ACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATG  >  NZ_CP009273/3937219‑3937413

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