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%	123.4 / NA	32	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/19);  total (13/19)

GAAAGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGAC  >  NC_000913/3213757‑3213847                                               |                                             gAAAGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCt                                            >  1:876445/1‑49 (MQ=255)     ggTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGt                                       <  1:2075023/50‑1 (MQ=255)      gTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTa                                      <  1:2258893/50‑1 (MQ=255)      gTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTa                                      <  1:1866634/50‑1 (MQ=255)       tCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTAt                                     <  1:1153205/50‑1 (MQ=255)       tCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTAt                                     <  1:1499575/50‑1 (MQ=255)         gCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTAtt                                    <  1:546651/49‑1 (MQ=255)         gCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTc                                   <  1:119219/50‑1 (MQ=255)          cAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCa                                  <  1:528905/50‑1 (MQ=255)             tCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCaaa                                >  1:927145/1‑49 (MQ=255)               aCGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAg                             >  1:1552357/1‑50 (MQ=255)                cGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGt                            >  1:133723/1‑50 (MQ=255)                   tACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTcc                         >  1:1555913/1‑50 (MQ=255)                     ccGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTcc                         <  1:941768/48‑1 (MQ=255)                     ccGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTcc                         <  1:633911/48‑1 (MQ=255)                       gCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCt                     >  1:2265831/1‑50 (MQ=255)                          cAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGt                  >  1:1609098/1‑50 (MQ=255)                           aGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTgg                   <  1:32029/48‑1 (MQ=255)                           aGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGt                  <  1:804455/49‑1 (MQ=255)                           aGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGt                  <  1:1188305/49‑1 (MQ=255)                           aGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTg                 <  1:555759/50‑1 (MQ=255)                           aGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTg                 <  1:2232984/50‑1 (MQ=255)                            ggAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGc                >  1:2241136/1‑50 (MQ=255)                                  gaTCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAg           >  1:1129842/1‑49 (MQ=255)                                   aTCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAg           >  1:1403972/1‑48 (MQ=255)                                   aTCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCa         <  1:643184/50‑1 (MQ=255)                                    tCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAg        >  1:1065192/1‑50 (MQ=255)                                      cTGAAACTGCCTGAAGTATTCAAATAGTTCCGCCTGGTGCCGAAGCAGt       <  1:739593/49‑1 (MQ=38)                                      cTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGtt      <  1:445471/50‑1 (MQ=255)                                        gAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGttt     >  1:1652430/1‑49 (MQ=255)                                         aaaCTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGttt     <  1:1421937/48‑1 (MQ=255)                                           aCTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGAc  >  1:1594540/1‑49 (MQ=255)                                               |                                             GAAAGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGAC  >  NC_000913/3213757‑3213847

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