Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,941,874	C→T	noncoding (67/1542 nt)	rrsC →	16S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,941,874	0	C	T	100.0%	47.8 / NA	21	noncoding (67/1542 nt)	rrsC	16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (16/5);  total (16/5)

CGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCC  >  NC_000913/3941635‑3942119                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      cgcgCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGa                                                                                                                                                                                                                                             >  2:163878/1‑250 (MQ=14)      cgcTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGc                                                                                                                                                                                                                                       >  1:171014/1‑251 (MQ=18)         tGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTg                                                                                                                                                                                                                                    >  2:180507/1‑251 (MQ=18)            gaAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTg                                                                                                                                                                                                                                 >  1:189767/1‑251 (MQ=18)                      agcGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGAc                                                                                                                                                                                                                       <  1:180507/251‑1 (MQ=17)                                                                                                                                       aTTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCgg                                                                                                      <  2:231438/251‑1 (MQ=14)                                                                                                                                            ttACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCctc                                                                                                 >  1:291584/1‑251 (MQ=14)                                                                                                                                                                  cATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGAt                                                                           <  2:274894/251‑1 (MQ=18)                                                                                                                                                                    tCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGTCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTATTGCCATCCGATGGGCCCAGATgg                                                                         >  2:249438/1‑251 (MQ=9)                                                                                                                                                                                   gAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGt                                                          >  1:164416/1‑251 (MQ=25)                                                                                                                                                                                             aTCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACgg                                                >  1:281622/1‑251 (MQ=37)                                                                                                                                                                                             aTCATGGCCCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGGAACgg                                                >  2:116120/1‑251 (MQ=255)                                                                                                                                                                                                       aGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGAGGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAgg                                      >  1:115536/1‑251 (MQ=255)                                                                                                                                                                                                             aaCGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAAAAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGAGACCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGcgacg                                 >  1:180079/1‑250 (MQ=255)                                                                                                                                                                                                                              tAACACATGCAAGGCGAATGGTAACAGGAACCAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGAGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAAAGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCCGATGGGATTAGATAGTAGGGGGGGTCACGGCACACCTAGGCGACGATCCACAGCTGGTCTg                 >  2:196059/1‑249 (MQ=14)                                                                                                                                                                                                                               aaCACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTga                <  2:171014/249‑1 (MQ=39)                                                                                                                                                                                                                                  acaTGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGAGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAgg            >  1:12308/1‑250 (MQ=255)                                                                                                                                                                                                                                     tGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGGAGGTGGGGTAACGGCTCACCTAGGCGACGAACCCTAACCGGGCTGAGAGGaga         >  2:228640/1‑248 (MQ=14)                                                                                                                                                                                                                                       cAAGGCGAATGGTAACAGGAAACAGCGTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGAcc      <  2:87688/251‑1 (MQ=255)                                                                                                                                                                                                                                          gTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGc   >  1:265030/1‑251 (MQ=38)                                                                                                                                                                                                                                            cGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGGGAGTAATGTCTGGGACACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCACAGAGGGGGACCTTCCGGACTCTTGCCATCCGATGCGCCCAGATGGGGTTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTCACTGGGCTGGGAGGATGACCCGcc  >  2:229595/1‑250 (MQ=11)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      CGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCC  >  NC_000913/3941635‑3942119

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