Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,168,196	C→T	noncoding (1538/1542 nt)	rrsB →	16S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,168,196	0	C	T	100.0%	106.6 / NA	31	noncoding (1538/1542 nt)	rrsB	16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/10);  total (21/10)

ACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAG  >  NC_000913/4168158‑4168257                                       |                                                              aCAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTTCTTACCTTAAAGAAGCGTACTTTGTAGTg                                  >  1:3093904/1‑68 (MQ=25) aCAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTTCTTACCTTAAAGAAGCGTACTTTGTAGTg                                  <  1:3069697/68‑1 (MQ=25)  cAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTTCTTACCTTAAAGAAGCGTACTTTGTAGTGc                                 >  1:1665295/1‑68 (MQ=25)     ggTAACCGTAGGGGAACCTGCGGTTGGATCACCTTCTTACCTTAAAGAAGCGTACTTTGTAGTGCTca                              <  1:1098218/68‑1 (MQ=25)         gccGTAGGGGAACCTGCGGTTGGATCACCTTCTTACCTTAAAGAAGCGTACTTTGTAGTGCTcacaca                          <  1:3919516/67‑1 (MQ=255)         aCCGTAGGGGAACCTGCGGTTGGATCACCTTCTTACCTTAAAGAAGCGTACTTTGTAGTGCTcacaca                          <  1:1392998/68‑1 (MQ=255)         aCCGTAGGGGAACCTGCGGTTGGATCACCTTCTTACCTTAAAGAAGCGTACTTTGTAGTGCTcacaca                          <  1:170690/68‑1 (MQ=255)            gTAGGGGAACCTGCGGTTGGATCACCTTCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGAt                       >  1:2676671/1‑68 (MQ=255)               ggggAACCTGCGGTTGGATTACCTTCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGt                    <  1:3601892/68‑1 (MQ=255)               ggggAACCTGCGGTTGGATCACCTTCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGt                    <  1:3956497/68‑1 (MQ=255)                gggAACCTGCGGTTGGATCACCTTCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTc                   >  1:1130220/1‑68 (MQ=255)                 ggAACCTGCGGTTGGATGACCTTCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCt                  >  1:2326511/1‑68 (MQ=255)                 ggAACCTGCGGTTGGATCACCTTCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCt                  >  1:164342/1‑68 (MQ=255)                 ggAACCTGCGGTTGGATCACCTTCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCt                  >  1:3526215/1‑68 (MQ=255)                 ggAACCTGCGGTTGGATCACCTTCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCt                  >  1:3361182/1‑68 (MQ=255)                 ggAACCTGCGGTTGGATCACCTTCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCt                  >  1:674978/1‑68 (MQ=255)                 ggAACCTGCGGTTGGATCACCTTCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCt                  >  1:3169547/1‑68 (MQ=255)                 ggAACCTGCGGTTGGATCACCTTCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCt                  >  1:692940/1‑68 (MQ=255)                 ggAACCTGCGGTTGGATCACCTTCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCt                  >  1:2976812/1‑68 (MQ=255)                 ggAACCTGCGGTTGGATCACCTTCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCt                  >  1:2696574/1‑68 (MQ=255)                 ggAACCTGCGGTTGGATCACCTTCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCt                  >  1:2475786/1‑68 (MQ=255)                 ggAACCTGCGGTTGGATCACCTTCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCt                  >  1:2316802/1‑68 (MQ=255)                 ggAACCTGCGGTTGGATCACCTTCTTACCTTAAAGAAGCGCACTTTGTAGTGCTCACACAGATTGTCt                  >  1:3415720/1‑68 (MQ=11)                 ggAACCTGCGGTTGGAACACCTTCTTACCTTAAAGAAGCGTCCTTTGTAGTGCTCACACAGATTGTCt                  >  1:2194604/1‑68 (MQ=11)                   aaCCTGCGGTTGGATCACCTTCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGa                >  1:1800562/1‑68 (MQ=255)                    aCCTGCGGTTGGATCACCTTCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGAt               >  1:703400/1‑68 (MQ=255)                      cTGCGGTTGGATCACCTTCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAg             >  1:763182/1‑68 (MQ=255)                           gTTGGATCACCTTCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGt        <  1:3133987/68‑1 (MQ=255)                             tGGATCACCTTCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGACAGAAAGTGa      <  1:179930/68‑1 (MQ=255)                               gATCACCTTCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGaaa    <  1:112628/68‑1 (MQ=255)                                 tCACCTTCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAg  <  1:597025/68‑1 (MQ=255)                                 tCACCTTCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAg  >  1:3217250/1‑68 (MQ=255)                                       |                                                              ACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAG  >  NC_000913/4168158‑4168257

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑