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

GAAAGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGAC  >  NC_000913/3213757‑3213847                                               |                                             gAAAGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCt                                            >  1:1263057/1‑49 (MQ=255) gAAAGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCt                                            >  1:1512929/1‑49 (MQ=255)   aaGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTg                                           >  1:2221513/1‑48 (MQ=255)   aaGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGaa                                         <  1:1791906/50‑1 (MQ=255)    aGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAg                                        <  1:1203832/50‑1 (MQ=255)    aGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAg                                        >  1:1288405/1‑50 (MQ=255)    aGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAg                                        <  1:1857746/50‑1 (MQ=255)     ggTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAg                                        >  1:2031305/1‑49 (MQ=255)     ggTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAg                                        >  1:1197975/1‑49 (MQ=255)     ggTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGt                                       <  1:1491141/50‑1 (MQ=255)         gCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTc                                   <  1:1853519/50‑1 (MQ=255)           aGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCa                                  >  1:1545641/1‑49 (MQ=255)            gTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCaaa                                >  1:1442186/1‑50 (MQ=255)              cACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCaaa                                >  1:1009322/1‑48 (MQ=255)                  cTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTc                          <  1:719880/50‑1 (MQ=255)                  cTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTc                          <  1:465580/50‑1 (MQ=255)                  cTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTc                          <  1:1669987/50‑1 (MQ=255)                    aCCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCg                        >  1:808115/1‑50 (MQ=255)                     ccGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTcc                         <  1:1029251/48‑1 (MQ=255)                     ccGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTcc                         <  1:1571045/48‑1 (MQ=255)                     ccGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTcc                         <  1:1562205/48‑1 (MQ=255)                     ccGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCg                        <  1:106942/49‑1 (MQ=255)                          cAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGt                  >  1:285009/1‑50 (MQ=255)                           aGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTgg                   <  1:871626/48‑1 (MQ=255)                           aGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGt                  <  1:565340/49‑1 (MQ=255)                           aGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTg                 >  1:512592/1‑50 (MQ=255)                             gAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGcc               <  1:1765024/50‑1 (MQ=255)                              aaGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGcc               >  1:250551/1‑49 (MQ=255)                                gagaTCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGaa            <  1:1862914/50‑1 (MQ=255)                                   aTCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAg           >  1:414031/1‑48 (MQ=255)                                   aTCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCa         <  1:1063663/50‑1 (MQ=255)                                    tCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAg        >  1:917831/1‑50 (MQ=255)                                      cTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGt       <  1:2031503/49‑1 (MQ=255)                                        gAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGttt     >  1:1629922/1‑49 (MQ=255)                                         aaaCTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGa   >  1:1328200/1‑50 (MQ=255)                                           aCTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGAc  >  1:1510201/1‑49 (MQ=255)                                               |                                             GAAAGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGAC  >  NC_000913/3213757‑3213847

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