Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,188,352	G→A	A1001A (GCG→GCA)	rpoC →	RNA polymerase subunit beta'

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,188,352	0	G	A	100.0%	152.2 / NA	44	A1001A (GCG→GCA)	rpoC	RNA polymerase subunit beta'
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (27/17);  total (27/17)

TAGCATCAAGCTCAGCAACGTGAAGTCGGTTGTGAACTCCAGCGGTAAACTGGTTATCACTTCCCGTAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCGGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACACCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTA  >  NC_000913/4188217‑4188488                                                                                                                                        |                                                                                                                                         tAGCATCAAGCTCAGCAACGTGAAGTCGGTTGTGAACTCCAGCGGTAAACTGGTTATCACTTCCCGTAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCAGTAc                                                                                                                                      >  2:412993/1‑140 (MQ=255)               gCAACGTGAAGTCGGTTGTGAACTCCAGCGGTAAACTGGTTATCACTTCCCGTAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTGCCTTACGGTGCAGTACTGGCGAAAGGCGAt                                                                                                                        >  2:560285/1‑140 (MQ=255)                 aaCGTGAAGTCGGTTGTGAACTCCAGCGGTAAACTGGTTATCACTTCCCGTAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCAGTACTGGCGAAAGGCGATgg                                                                                                                      <  2:558475/140‑1 (MQ=255)                 aaCGTGAAGTCGGTTGTGAACTCCAGCGGTAAACTGGTTATCACTTCCCGTAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCAGTACTGGCGAAAGGCGATgg                                                                                                                      >  1:558475/1‑140 (MQ=255)                        aGTCGGTTGTGAACTCCAGCGGTAAACTGGTTATCACTTCCCGTAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCAGTACTGGCGAAAGGCGATGGCGAACAg                                                                                                               >  1:452443/1‑140 (MQ=255)                             gTTGTGAACTCCAGCGGTAAACTGGTTATCACTTCCCGTAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCAGTACTGGCGAAAGGCGATGGCGAACAGGTTGc                                                                                                          >  2:411689/1‑140 (MQ=255)                                        cAGCGGTAAACTGGTTATCACTTCCCGTAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCAGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGaaa                                                                                               >  2:910488/1‑140 (MQ=255)                                        cAGCGGTAAACTGGTTATCACTTCCCGTAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCAGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGaaa                                                                                               <  1:917521/140‑1 (MQ=255)                                               aaaCTGGTTATCACTTCCCGTAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCAGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGc                                                                                        <  2:523529/140‑1 (MQ=255)                                                aaCTGGTTATCACTTCCCGTAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCAGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCa                                                                                       >  2:660373/1‑140 (MQ=255)                                                     gTTATCACTTCCCGTAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCAGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTg                                                                                  >  2:702829/1‑140 (MQ=255)                                                      ttATCACTTCCCGTAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCAGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTgg                                                                                 <  1:702829/140‑1 (MQ=255)                                                      ttATCACTTCCCGTAATACTGAACTGAAACTGATCGACGAATCCGGTCGTACTAAAGAAAGCTACAAAGTCCCTTACGGTGCAGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTgg                                                                                 <  2:452443/140‑1 (MQ=255)                                                        aTCACTTCCCGTAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCAGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGa                                                                               >  1:121657/1‑140 (MQ=255)                                                        aTCACTTCCCGTAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCAGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGa                                                                               >  1:635613/1‑140 (MQ=255)                                                              tCCCGTAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCAGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGca                                                                         <  2:752564/140‑1 (MQ=255)                                                               cccGTAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCAGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGcac                                                                        >  1:462384/1‑140 (MQ=255)                                                                ccGTAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCAGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGcaca                                                                       >  2:324613/1‑140 (MQ=255)                                                                       actgaactgaAACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCAGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACACCATGcc                                                                >  1:114002/1‑140 (MQ=255)                                                                          gaactgaaACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCAGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACACCATGCCGGt                                                             >  1:151252/1‑140 (MQ=255)                                                                                atactgaTCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCAGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACACCATGCCGGTTATCAc                                                       >  2:403205/3‑140 (MQ=255)                                                                                aaactgaTCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCAGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACACCATGCCGGTTATCAc                                                       >  2:41260/1‑140 (MQ=255)                                                                                 aactgaTCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCAGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACACCATGCCGGTTATCAcc                                                      >  1:522761/1‑140 (MQ=255)                                                                                  actgaTCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCAGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACACCATGCCGGTTATCACCg                                                     >  1:509046/1‑140 (MQ=255)                                                                                  actgaTCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCAGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACACCATGCCGGTTATCACCg                                                     >  1:637337/1‑140 (MQ=255)                                                                                    tgaTCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCAGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACACCATGCCGGTTATCACCGaa                                                   <  2:635613/140‑1 (MQ=255)                                                                                     gaTCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCAGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACACCATGCCGGTTATCACCGAAg                                                  >  1:559643/1‑140 (MQ=255)                                                                                         gacgaATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCAGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACACCATGGCGGTTATCACCGAAGTAAg                                              <  1:324613/140‑1 (MQ=255)                                                                                         gacgaATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCAGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACACCATGCCGGTTATCACCGAAGTAAg                                              <  2:121657/140‑1 (MQ=255)                                                                                          acgaATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCAGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACACCATGCCGGTTATCACCGAAGTAAGc                                             >  2:526504/1‑140 (MQ=255)                                                                                                   gTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCAGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACACCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTa                                    <  2:462384/140‑1 (MQ=255)                                                                                                     cGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCAGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACACCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACg                                  <  2:522761/140‑1 (MQ=255)                                                                                                     cGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCAGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACACCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACg                                  <  1:403205/140‑1 (MQ=255)                                                                                                       tACTAAAGAAAGCTACAAAGTACCTTACGGTGCAGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACACCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCt                                <  1:412993/140‑1 (MQ=255)                                                                                                          tAAAGAAAGCTACAAAGTACCTTACGGTGCAGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACACCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTa                             <  1:910488/140‑1 (MQ=255)                                                                                                            aagaaagCTACAAAGTACCTTACGGTGCAGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACACCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACt                           <  1:560285/140‑1 (MQ=255)                                                                                                              gaaagCTACAAAGTACCTTACGGTGCAGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACACCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGa                         >  1:709607/1‑140 (MQ=255)                                                                                                                      cAAAGTACCTTACGGTGCAGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACACCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCg                 >  2:204951/1‑140 (MQ=255)                                                                                                                         aGTACCTTACGGTGCAGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACACCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACg              >  1:536751/1‑140 (MQ=255)                                                                                                                           tACCTTACGGTGCAGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACACCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGc            >  1:202787/1‑140 (MQ=255)                                                                                                                             ccTTACGGTGCAGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACACCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCa          >  1:902467/1‑140 (MQ=255)                                                                                                                                tACGGTGCAGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACACCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGAc       <  2:151252/140‑1 (MQ=255)                                                                                                                                tACGGTGCAGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACACCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGAc       <  1:204951/140‑1 (MQ=255)                                                                                                                                     tGCAGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACACCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTa  >  1:313416/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                         TAGCATCAAGCTCAGCAACGTGAAGTCGGTTGTGAACTCCAGCGGTAAACTGGTTATCACTTCCCGTAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCGGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACACCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTA  >  NC_000913/4188217‑4188488

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