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%	187.2 / NA	46	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/30);  total (16/30)

CGAAAGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTA  >  NC_000913/3213756‑3213849                                                |                                               cGAAAGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCt                                              <  1:95104/50‑1 (MQ=255) cGAAAGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCt                                              <  1:3429040/50‑1 (MQ=255)  gAAAGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCt                                              >  1:2836457/1‑49 (MQ=255)  gAAAGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTg                                             <  1:2419625/50‑1 (MQ=255)    aaGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTg                                             >  1:468987/1‑48 (MQ=255)     aGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAg                                          <  1:1207307/50‑1 (MQ=255)     aGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAg                                          <  1:3606066/50‑1 (MQ=255)      ggTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGt                                         <  1:2967020/50‑1 (MQ=255)          gCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTAtt                                      >  1:3134077/1‑49 (MQ=255)          gCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTAtt                                      <  1:433892/49‑1 (MQ=255)          gCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTAtt                                      >  1:2143685/1‑49 (MQ=255)           cAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCa                                    <  1:1834902/50‑1 (MQ=255)            aGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCaa                                   <  1:3887754/50‑1 (MQ=255)              tCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAAc                                 <  1:1862451/50‑1 (MQ=255)              tCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAAc                                 <  1:123898/50‑1 (MQ=255)              tCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAAc                                 <  1:2107813/50‑1 (MQ=255)                aCGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAg                               <  1:1725694/50‑1 (MQ=255)                aCGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAg                               >  1:768192/1‑50 (MQ=255)                 cGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGt                              <  1:913497/50‑1 (MQ=255)                   cTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTc                            <  1:2532753/50‑1 (MQ=255)                    tACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTcc                           >  1:1949500/1‑50 (MQ=255)                    tACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTcc                           >  1:2964969/1‑50 (MQ=255)                      ccGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTcc                           <  1:3685299/48‑1 (MQ=255)                      ccGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTcc                           <  1:2711913/48‑1 (MQ=255)                      ccGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCg                          <  1:3411602/49‑1 (MQ=255)                      ccGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGc                         >  1:3374636/1‑50 (MQ=255)                        gCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCt                       >  1:1209158/1‑50 (MQ=255)                        gCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCt                       <  1:1535571/50‑1 (MQ=255)                         cTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTg                      <  1:3084825/50‑1 (MQ=255)                           cAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGt                    >  1:2842825/1‑50 (MQ=255)                           cAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGt                    >  1:3278125/1‑50 (MQ=255)                            aGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGt                    <  1:1336350/49‑1 (MQ=255)                            aGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTg                   <  1:2055014/50‑1 (MQ=255)                            aGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTg                   >  1:640152/1‑50 (MQ=255)                             ggAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGc                  <  1:1947050/50‑1 (MQ=255)                              gAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGcc                 >  1:3536548/1‑50 (MQ=255)                              gAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGcc                 <  1:427102/50‑1 (MQ=255)                                  agaTCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAg             <  1:271393/50‑1 (MQ=255)                                  agaTCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAg             <  1:1931879/50‑1 (MQ=255)                                   gaTCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAg             >  1:648321/1‑49 (MQ=255)                                         gAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGttt       >  1:2105053/1‑49 (MQ=255)                                          aaaCTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGttt       <  1:182678/48‑1 (MQ=255)                                          aaaCTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGttt       <  1:843527/48‑1 (MQ=255)                                            aCTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACt   <  1:1355891/50‑1 (MQ=255)                                             cTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTa  >  1:1088316/1‑50 (MQ=255)                                             cTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTa  <  1:1327156/50‑1 (MQ=255)                                                |                                               CGAAAGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTA  >  NC_000913/3213756‑3213849

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: ‑