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

AAAGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACT  >  NC_000913/3213758‑3213848                                              |                                              aaaGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCt                                             >  1:1915370/1‑48 (MQ=255) aaaGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTg                                            <  1:1478667/49‑1 (MQ=255)  aaGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGa                                           <  1:1906166/49‑1 (MQ=255)    ggTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAg                                         >  1:320712/1‑49 (MQ=255)     gTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAg                                         >  1:1830030/1‑48 (MQ=255)      tCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTa                                       >  1:2180329/1‑49 (MQ=255)      tCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTa                                       <  1:970175/49‑1 (MQ=255)          aGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTc                                    >  1:185250/1‑48 (MQ=255)            tCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCaaa                                 <  1:2097333/49‑1 (MQ=255)             cACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCaaa                                 >  1:898891/1‑48 (MQ=255)              aCGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAAc                                <  1:2351022/48‑1 (MQ=255)              aCGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAAc                                <  1:696051/48‑1 (MQ=255)               cGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAg                              <  1:201337/49‑1 (MQ=255)                gCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGt                             >  1:200761/1‑49 (MQ=255)                gCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGt                             <  1:1859826/49‑1 (MQ=255)                  tACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTc                           >  1:2245990/1‑49 (MQ=255)                    ccGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCg                         >  1:631953/1‑49 (MQ=255)                     cGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGc                        <  1:1268545/49‑1 (MQ=255)                      gCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGcc                       <  1:1799236/49‑1 (MQ=255)                      gCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGcc                       <  1:624926/49‑1 (MQ=255)                      gCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGcc                       <  1:2096931/49‑1 (MQ=255)                       cTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCt                      <  1:1788588/49‑1 (MQ=255)                       cTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCt                      >  1:1238391/1‑49 (MQ=255)                       cTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCt                      >  1:1118884/1‑49 (MQ=255)                        tCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTg                     >  1:1009408/1‑49 (MQ=255)                        tCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTg                     >  1:1771082/1‑49 (MQ=255)                          aGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTgg                    <  1:1678876/48‑1 (MQ=255)                          aGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGt                   <  1:876704/49‑1 (MQ=255)                          aGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGt                   >  1:1362693/1‑49 (MQ=255)                          aGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGt                   >  1:99180/1‑49 (MQ=255)                            gAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGc                 >  1:1460253/1‑49 (MQ=255)                            gAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGc                 >  1:1963227/1‑49 (MQ=255)                             aaGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGcc                >  1:838195/1‑49 (MQ=255)                              agagATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGcc                >  1:879430/1‑48 (MQ=255)                              agagATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCg               <  1:1049859/49‑1 (MQ=255)                               gagaTCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCg               >  1:1558015/1‑48 (MQ=255)                                 gaTCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAg            >  1:210362/1‑49 (MQ=255)                                  aTCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGc           <  1:1842899/49‑1 (MQ=255)                                   tCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAg            >  1:872347/1‑47 (MQ=255)                                   tCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCa          <  1:1343638/49‑1 (MQ=255)                                     cTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAg         <  1:1160467/48‑1 (MQ=255)                                     cTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGt        >  1:816378/1‑49 (MQ=255)                                     cTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGt        <  1:817180/49‑1 (MQ=255)                                      tGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGtt       >  1:742688/1‑49 (MQ=255)                                      tGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGtt       >  1:939804/1‑49 (MQ=255)                                        aaaCTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGtt       <  1:1975466/47‑1 (MQ=255)                                          aCTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGAc   <  1:811343/49‑1 (MQ=255)                                          aCTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGAc   <  1:1146312/49‑1 (MQ=255)                                           cTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGAc   >  1:1872338/1‑48 (MQ=255)                                           cTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACt  <  1:1456637/49‑1 (MQ=255)                                              |                                              AAAGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACT  >  NC_000913/3213758‑3213848

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