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

GAAAGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGAC  >  NC_000913/3213757‑3213847                                               |                                             gAAAGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCt                                            >  1:734331/1‑49 (MQ=255)  aaaGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCt                                            >  1:1829800/1‑48 (MQ=255)  aaaGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCt                                            >  1:1630827/1‑48 (MQ=255)  aaaGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCt                                            >  1:2529299/1‑48 (MQ=255)   aaGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTg                                           >  1:99130/1‑48 (MQ=255)   aaGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTg                                           >  1:1333716/1‑48 (MQ=255)   aaGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGa                                          <  1:1019458/49‑1 (MQ=255)    aGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTg                                           >  1:2263160/1‑47 (MQ=255)     ggTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAg                                        <  1:2682210/49‑1 (MQ=255)      gTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGt                                       <  1:1743789/49‑1 (MQ=255)       tCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTa                                      <  1:1941804/49‑1 (MQ=255)       tCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTa                                      <  1:221252/49‑1 (MQ=255)        cGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTAt                                     >  1:2158512/1‑49 (MQ=255)         gCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTAtt                                    <  1:1096359/49‑1 (MQ=255)           aGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTc                                   >  1:1636171/1‑48 (MQ=255)           aGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCa                                  <  1:2345094/49‑1 (MQ=255)           aGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCa                                  <  1:1616875/49‑1 (MQ=255)            gTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCaa                                 <  1:1461408/49‑1 (MQ=255)            gTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCa                                  >  1:1736886/1‑48 (MQ=255)               aCGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCaaa                                >  1:2609929/1‑47 (MQ=255)               aCGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCaaa                                >  1:356823/1‑47 (MQ=255)               aCGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACa                              <  1:2502615/49‑1 (MQ=255)               aCGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACa                              <  1:1812511/49‑1 (MQ=255)               aCGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACa                              <  1:1906724/49‑1 (MQ=255)                cGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAg                             >  1:2590035/1‑49 (MQ=255)                  cTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGtt                           >  1:1608165/1‑49 (MQ=255)                  cTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGtt                           >  1:974923/1‑49 (MQ=255)                   tACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTc                          >  1:1109008/1‑49 (MQ=255)                    aCCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTcc                         >  1:360619/1‑49 (MQ=255)                     ccGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTc                          <  1:1980223/47‑1 (MQ=255)                     ccGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCg                        >  1:1791692/1‑49 (MQ=255)                     ccGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCg                        <  1:1985849/49‑1 (MQ=255)                       gCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGcc                      >  1:2113911/1‑49 (MQ=255)                       gCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGcc                      >  1:1196603/1‑49 (MQ=255)                       gCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGcc                      >  1:1990511/1‑49 (MQ=255)                        cTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCt                     <  1:2181724/49‑1 (MQ=255)                        cTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCt                     <  1:491419/49‑1 (MQ=255)                        cTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCt                     <  1:1753381/49‑1 (MQ=255)                          cAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTgg                   >  1:1885514/1‑49 (MQ=255)                           aGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGt                  <  1:2616391/49‑1 (MQ=255)                           aGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGt                  >  1:118105/1‑49 (MQ=255)                            ggAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTg                 <  1:2342377/49‑1 (MQ=255)                              aaGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGcc               <  1:1339514/49‑1 (MQ=255)                                 agaTCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCg              >  1:2819498/1‑47 (MQ=255)                                 agaTCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGaa            <  1:1146168/49‑1 (MQ=255)                                  gaTCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAg           <  1:1815345/49‑1 (MQ=255)                                  caTCCTGAAACTTCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAg           <  1:407184/48‑1 (MQ=255)                                   aTCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAg           >  1:1715672/1‑48 (MQ=255)                                   aTCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAg           >  1:835462/1‑48 (MQ=255)                                   aTCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAg           >  1:1315139/1‑48 (MQ=255)                                   aTCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAg           >  1:2466590/1‑48 (MQ=255)                                    tCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCa         <  1:879126/49‑1 (MQ=255)                                     ccTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAg        >  1:1423753/1‑49 (MQ=255)                                     ccTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAg        >  1:1577647/1‑49 (MQ=255)                                     ccTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAg        >  1:2523191/1‑49 (MQ=255)                                      cTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGt       <  1:1599135/49‑1 (MQ=255)                                        gAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGttt     >  1:2248821/1‑49 (MQ=255)                                         aaaCTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGttt     <  1:24567/48‑1 (MQ=255)                                         aaaCTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTg    <  1:2743095/49‑1 (MQ=255)                                          aaCTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGa   <  1:769731/49‑1 (MQ=255)                                          aaCTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGa   >  1:864350/1‑49 (MQ=255)                                          aaCTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGa   >  1:1034588/1‑49 (MQ=255)                                           aCTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGAc  <  1:2518202/49‑1 (MQ=255)                                            cTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGAc  >  1:43455/1‑48 (MQ=255)                                               |                                             GAAAGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGAC  >  NC_000913/3213757‑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: ‑