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%	108.4 / NA	28	S253P (TCT→CCT)	rpoD	RNA polymerase, sigma 70 (sigma D) factor
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (13/15);  total (13/15)

AAAGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGAC  >  NC_000913/3213758‑3213847                                              |                                             aaaGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTg                                           <  1:242039/49‑1 (MQ=255)  aaGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGa                                          <  1:639080/49‑1 (MQ=255)     gTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAg                                        >  1:2037385/1‑48 (MQ=255)         cAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTc                                   <  1:1006703/49‑1 (MQ=255)         cAGTCACGCTACCGCTCAGGAAAAGATCCTGAAACTGCCTGAAGTAtt                                    >  1:1940662/1‑48 (MQ=38)          aGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTc                                   >  1:943970/1‑48 (MQ=255)            tCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCaaa                                <  1:1795074/49‑1 (MQ=255)             cACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCaaa                                >  1:793864/1‑48 (MQ=255)              aCGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAAc                               <  1:240223/48‑1 (MQ=255)                 cTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGtt                           <  1:1920540/49‑1 (MQ=255)                    ccGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTc                          <  1:1358976/47‑1 (MQ=255)                    ccGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTc                          <  1:1334981/47‑1 (MQ=255)                    ccGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTc                          <  1:369997/47‑1 (MQ=255)                     cGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGc                       >  1:1276595/1‑49 (MQ=255)                      gCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGcc                      >  1:1166065/1‑49 (MQ=255)                      gCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGcc                      >  1:718355/1‑49 (MQ=255)                       cTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCGAACAGTTCCGCCt                     <  1:98003/49‑1 (MQ=38)                       cTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCt                     <  1:1493561/49‑1 (MQ=255)                       cTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCt                     <  1:1720489/49‑1 (MQ=255)                        tCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTg                    >  1:1165601/1‑49 (MQ=255)                           ggAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTg                 >  1:354310/1‑49 (MQ=255)                            gAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGc                <  1:878562/49‑1 (MQ=255)                             aaGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGcc               <  1:1052325/49‑1 (MQ=255)                             aaGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGcc               >  1:781554/1‑49 (MQ=255)                             aaGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGcc               >  1:1436952/1‑49 (MQ=255)                             aaGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGcc               <  1:1307643/49‑1 (MQ=255)                                        aaaCTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGttt     >  1:317122/1‑48 (MQ=255)                                          aCTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGAc  >  1:2013862/1‑49 (MQ=255)                                              |                                             AAAGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGAC  >  NC_000913/3213758‑3213847

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