Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,213,803	T→C	S253P (TCT→CCT)	rpoD →	RNA polymerase, sigma 70 (sigma D) factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,213,803	0	T	C	100.0%	114.4 / NA	31	S253P (TCT→CCT)	rpoD	RNA polymerase, sigma 70 (sigma D) factor
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (16/15);  total (16/15)

CGAAAGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGAC  >  NC_000913/3213756‑3213847                                                |                                             cGAAAGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCt                                            >  1:934328/1‑50 (MQ=255)  gAAAGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCt                                            >  1:1940825/1‑49 (MQ=255)      ggTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGt                                       <  1:2576339/50‑1 (MQ=255)       gTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTa                                      >  1:1488274/1‑50 (MQ=255)       gTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTa                                      <  1:2071382/50‑1 (MQ=255)          gCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTAtt                                    <  1:567256/49‑1 (MQ=255)           cAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCa                                  <  1:1878018/50‑1 (MQ=255)           cAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCa                                  <  1:2561477/50‑1 (MQ=255)               cACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCaaa                                >  1:1743344/1‑48 (MQ=255)                aCGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACa                              <  1:1398357/49‑1 (MQ=255)                  gCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGtt                           >  1:2026487/1‑50 (MQ=255)                   cTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTc                          >  1:1119638/1‑50 (MQ=255)                    tACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTcc                         >  1:1629921/1‑50 (MQ=255)                     aCCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCg                        >  1:2203463/1‑50 (MQ=255)                      ccGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTcc                         <  1:1547832/48‑1 (MQ=255)                      ccGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTcc                         <  1:196109/48‑1 (MQ=255)                       cGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGcc                      <  1:2399240/50‑1 (MQ=255)                        gCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCt                     <  1:988140/50‑1 (MQ=255)                           cAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGt                  >  1:831529/1‑50 (MQ=255)                           cAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGt                  >  1:288799/1‑50 (MQ=255)                             ggAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGc                >  1:107033/1‑50 (MQ=255)                              gAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGcc               >  1:1305861/1‑50 (MQ=255)                                agagATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCg              >  1:1554168/1‑49 (MQ=255)                                   gaTCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGc          <  1:571154/50‑1 (MQ=255)                                   gaTCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGc          <  1:928578/50‑1 (MQ=255)                                     tCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAg        <  1:2433508/50‑1 (MQ=255)                                     tCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAg        >  1:2311954/1‑50 (MQ=255)                                     tCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAg        >  1:2142253/1‑50 (MQ=255)                                       cTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGtt      <  1:2309132/50‑1 (MQ=255)                                          aaaCTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTg    <  1:2447477/49‑1 (MQ=255)                                            aCTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGAc  >  1:863072/1‑49 (MQ=255)                                                |                                             CGAAAGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGAC  >  NC_000913/3213756‑3213847

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