Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,188,124	A→G	E925E (GAA→GAG)	rpoC →	RNA polymerase, beta prime subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,188,124	0	A	G	100.0%	237.8 / NA	65	E925E (GAA→GAG)	rpoC	RNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (32/33);  total (32/33)

AGTACGTTCTGTTGTATCTTGTGACACCGACTTTGGTGTATGTGCGCACTGCTACGGTCGTGACCTGGCGCGTGGCCACATCATCAACAAGGGTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAACCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGCATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAGCAACGTGAAGTCGGTTGTGAACTCCAG  >  NC_000913/4187992‑4188258                                                                                                                                     |                                                                                                                                       aGTACGTTCTGTTGTATCTTGTGACACCGACTTTGGTGTATGTGCGCACTGCTACGGTCGTGACCTGGCGCGTGGCCACATCATCAACAAGGGTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAGCCGGGt                                                                                                                                  <  2:82075/139‑1 (MQ=255)  gTACGTTCTGTTGTATCTTGTGACACCGACTTTGGTGTATGTGCGCACTGCTACGGTCGTGACCTGGCGCGTGGCCACATCATCAACAAGGGTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAGCCGGGTa                                                                                                                                 <  2:610725/139‑1 (MQ=255)       ttCTGTTGTATCTTGTGACACCGACTTTGGTGTATGTGCGCACTGCTACGGTCGTGACCTGGCGCGTGGCCACATCATCAACAAGGGTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAGCCGGGTACACAg                                                                                                                            <  2:1216928/139‑1 (MQ=255)       ttCTGTTGTATCTTGTGACACCGACTTTGGTGTATGTGCGCACTGCTACGGTCGTGACCTGGCGCGTGGCCACATCATCAACAAGGGTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAGCCGGGTACACAg                                                                                                                            >  2:1082370/1‑139 (MQ=255)       ttCTGTTGTATCTTGTGACACCGACTTTGGTGTATGTGCGCACTGCTACGGTCGTGACCTGGCGCGTGGCCACATCATCAACAAGGGTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAGCCGGGTACACAg                                                                                                                            >  2:1050358/1‑139 (MQ=255)                   ttGTGACACCGACTTTGGTGTATGTGCGCACTGCTACGGTCGTGACCTGGCGCGTGGCCACATCATCAACAAGGGTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAGCCGGGTACACAGCTGACCATGCGt                                                                                                                <  1:389425/139‑1 (MQ=255)                      tgACACCGACTTTGGTGTATGTGCGCACTGCTACGGTCGTGACCTGGCGCGTGGCCACATCATCAACAAGGGTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAGCCGGGTACACAGCTGACCATGCGTACg                                                                                                             >  2:109056/1‑139 (MQ=255)                       gACACCGACTTTGGTGTATGTGCGCACTGCTACGGTCGTGACCTGGCGCGTGGCCACATCATCAGCAAGGGTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGCGCCGGGTACACAGCTGACCATGCGTACGt                                                                                                            <  2:814434/139‑1 (MQ=255)                           ccGACTTTGGTGTATGTGCGCACTGCTACGGTCGTGACCTGGCGCGTGGCCACATCATCAACAAGGGTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAGCCGGGTACACAGCTGACCATGCGTACGTTCca                                                                                                        >  2:145604/1‑139 (MQ=255)                            cGACTTTGGTGTATGTGCGCACTGCTACGGTCGTGACCTGGCGCGTGGCCACATCATCAACAAGGGTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAGCCGGGTACACAGCTGACCATGCGTACGTTCcac                                                                                                       >  2:531133/1‑139 (MQ=255)                              aCTTTGGTGTATGTGCGCACTGCTACGGTCGTGACCTGGCGCGTGGCCACATCATCAACAAGGGTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAGCCGGGTACACAGCTGACCATGCGTACGTTCCACAt                                                                                                     <  1:1050358/139‑1 (MQ=255)                                   ggTGTATGTGCGCACTGCTACGGTCGTGACCTGGCGCGTGGCCACATCATCAACAAGGGTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAGCCGGGTACACAGCTGACCATGCGTACGTTCCACATCggtg                                                                                                <  1:1082370/139‑1 (MQ=255)                                         tgtgCGCACTGCTACGGTCGTGACCTGGCGCGTGGCCACATCATCAACAAGGGTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAGCCGGGTACACAGCTGACCATGCGTACGTTc                                                                                                          >  2:665097/1‑123 (MQ=255)                                         tgtgCGCACTGCTACGGTCGTGACCTGGCGCGTGGCCACATCATCAACAAGGGTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAGCCGGGTACACAGCTGACCATGCGTACGTTc                                                                                                          <  1:665097/123‑1 (MQ=255)                                              gcACTGCTACGGTCGTGACCTGGCGCGTGGCCACATCATCAACAAGGGTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAGCCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGCAtct                                                                                     >  2:953482/1‑139 (MQ=255)                                              gcACTGCTACGGTCGTGACCTGGCGCGTGGCCACATCATCAACAAGGGTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAGCCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGCAtct                                                                                     >  2:27806/1‑139 (MQ=255)                                                  tGCTACGGTCGTGACCTGGCGCGTGGCCACATCATCAACAAGGGTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAGCCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGCATCTCGTg                                                                                 <  1:953482/139‑1 (MQ=255)                                                  tGCTACGGTCGTGACCTGGCGCGTGGCCACATCATCAACAAGGGTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAGCCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGCATCTCGTg                                                                                 <  1:145604/139‑1 (MQ=255)                                                   gCTACGGTCGTGACCTGGCGCGTGGCCACATCATCAACAAGGGTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAGCCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGCATCTCGTGc                                                                                >  1:512143/1‑139 (MQ=255)                                                    cTACGGTCGTGACCTGGCGCGTTGCCACATCATCAACAAGGGTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAGCCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGCATCTCGTGCg                                                                               >  2:270897/1‑139 (MQ=255)                                                    cTACGGTCGTGACCTGGCGCGTGGCCACATCATCAACAAGGGTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAGCCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGCATCTCGTGCg                                                                               >  2:167745/1‑139 (MQ=255)                                                            gtggccTGGCGCGTGGCCACATCATCAACAAGGGTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAGCCGGGTACACAGCTGACCATGCGTACg                                                                                                             >  2:366890/5‑101 (MQ=255)                                                            gTGGCCTGGCGCGTGGCCACATCATCAACAAGGGTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAGCCGGGTACACAGCTGACCATGCGTACg                                                                                                             <  1:366890/101‑1 (MQ=255)                                                            gTGACCTGGCGCGTGGCCACATCATCAACAAGGGTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAGCCGGGTACACAGCt                                                                                                                          <  1:522040/88‑1 (MQ=255)                                                            gTGACCTGGCGCGTGGCCACATCATCAACAAGGGTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAGCCGGGTACACAGCt                                                                                                                          >  2:522040/1‑88 (MQ=255)                                                            gTGACCTGGCGCGTGGCCACATCATCAACAAGGGTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAGCCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGCATCTCGTGCGGCTGCTGa                                                                       >  2:1146732/1‑139 (MQ=255)                                                            gTGACCTGGCGCGTGGCCACATCATCAACAAGGGTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAGCCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGCATCTCGTGCGGCTGCTGa                                                                       >  2:1001759/1‑139 (MQ=255)                                                              gACCTGGCGCGTGGCCACATCATCAACAAGGGTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAGCCGGGTACACAGCTGATCATGCGTACGTTCCACATCGGTGGTGCGGCATCTCGTGCGGCTGCTGAAt                                                                     >  2:773846/1‑139 (MQ=255)                                                              gACCTGGCGCGTGGCCACATCATCAACAAGGGTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAGCCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGCATCTCGTGCGGCTGCTGAAt                                                                     >  1:407347/1‑139 (MQ=255)                                                               aCCTGGCGCGTGGCCACATCATCAACAAGGGTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAGCCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGCATCTCGTGCGGCTGCTGAATc                                                                    <  2:949908/139‑1 (MQ=255)                                                                    gcgcGTGGCCACATCATCAACAAGGGTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAGCCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGCATCTCGTGCGGCTGCTGAATc                                                                    >  2:403194/1‑134 (MQ=255)                                                                    gcgcGTGGCCACATCATCAACAAGGGTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAGCCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGCATCTCGTGCGGCTGCTGAATc                                                                    <  1:403194/134‑1 (MQ=255)                                                                         tGGCCACATCATCAACAAGGGTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAGCCGGGTACACAGCTGACCATGCGTACGTTCCACATCgg                                                                                                  <  1:124689/99‑1 (MQ=255)                                                                         tGGCCACATCATCAACAAGGGTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAGCCGGGTACACAGCTGACCATGCGTACGTTCCACATCgg                                                                                                  >  2:124689/1‑99 (MQ=255)                                                                         tGGCCACATCATCAACAAGGGTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAGCCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGCATCTCGTGCGGCTGCTGAATCCAGCATCCaa                                                          <  1:1001759/139‑1 (MQ=255)                                                                         tGGCCACATCATCAACAAGGGTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAGCCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGCATCTCGTGCGGCTGCTGAATCCAGCATCCaa                                                          <  1:1146732/139‑1 (MQ=255)                                                                               catcatCAACAAGGGTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAGCCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGCATCTCGTGCgg                                                                              >  1:1073458/1‑113 (MQ=255)                                                                               catcatCAACAAGGGTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAGCCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGCATCTCGTGCgg                                                                              <  2:1073458/113‑1 (MQ=255)                                                                                   atcaACAAGGGTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAGCCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGCATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACa                                                <  2:147669/139‑1 (MQ=255)                                                                                    tcaACAAGGGTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAGCCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGCATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACaa                                               >  2:204292/1‑139 (MQ=255)                                                                                           gggTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAGcc                                                                                                                                      >  1:1206435/1‑45 (MQ=255)                                                                                           gggTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAGCCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGCATCTCGTg                                                                                 >  1:391433/1‑98 (MQ=255)                                                                                           gggTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAGCCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGCATCTCGTg                                                                                 <  2:391433/98‑1 (MQ=255)                                                                                            ggTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAGCCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGCATCTCGTGTGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCa                                       >  2:358351/1‑139 (MQ=255)                                                                                               gAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAGCCGGGTACACAGCTGa                                                                                                                        <  1:1140174/55‑1 (MQ=255)                                                                                               gAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAGCCGGGTACACAGCTGa                                                                                                                        >  2:1140174/1‑55 (MQ=255)                                                                                                aaGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAGCCGGGTa                                                                                                                                 <  1:340225/45‑1 (MQ=255)                                                                                                aaGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAGCCGGGTa                                                                                                                                 >  2:340225/1‑45 (MQ=255)                                                                                                aaGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAGCCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGCATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAgg                                            <  2:280954/130‑1 (MQ=255)                                                                                                aaGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAGCCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGCATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAgg                                            >  1:280954/1‑130 (MQ=255)                                                                                                 aGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAGCCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGCATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAg                                  <  1:167745/139‑1 (MQ=255)                                                                                                      tCGGTGTTATCGCGGCACAGTCCATCGGTGAGCCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGCATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAg                             >  2:31247/1‑139 (MQ=255)                                                                                                      tCGGTGTTATCGCGGCACAGTCCATCGGTGAGCCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGCATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAg                             >  1:116209/1‑139 (MQ=255)                                                                                                        ggTGTTATCGCGGCACAGTCCATCGGTGAGCCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGCATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAGCa                           <  2:1111073/139‑1 (MQ=255)                                                                                                               tCGCGGCACAGTCCATCGGTGAGCCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGCATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAGCAACGTGaa                    <  1:27806/139‑1 (MQ=255)                                                                                                               tCGCGGCACAGTCCATCGGTGAGCCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGCATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAGCAACGTGaa                    <  2:407347/139‑1 (MQ=255)                                                                                                                    gCACAGTCCATCGGTGAGCCGGGTACACAGCTGACCa                                                                                                                     >  2:668339/1‑37 (MQ=255)                                                                                                                    gCACAGTCCATCGGTGAGCCGGGTACACAGCTGACCa                                                                                                                     <  1:668339/37‑1 (MQ=255)                                                                                                                    gCACAGTCCATCGGTGAGCCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGCATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCa                                       <  1:997195/115‑1 (MQ=255)                                                                                                                    gCACAGTCCATCGGTGAGCCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGCATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCa                                       >  2:997195/1‑115 (MQ=255)                                                                                                                        aGTCCATCGGTGAGCCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGCATCTCGTGCgg                                                                              >  1:1073957/1‑72 (MQ=255)                                                                                                                        aGTCCATCGGTGAGCCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGCATCTCGTGCgg                                                                              <  2:1073957/72‑1 (MQ=255)                                                                                                                              tCGGTGAGCCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGCATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAGCAACGTGAAGTCGGTTGTGAACTc     <  2:647297/139‑1 (MQ=255)                                                                                                                                 gTGAGCCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGCATCTCGTGTGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAGCAACGTGAAGTCGGTTGTGAACTCCAg  <  1:358351/139‑1 (MQ=255)                                                                                                                                 gTGAGCCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGCATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAGCAACGTGAAGTCGGTTGTGAACTCCAg  <  2:16916/139‑1 (MQ=255)                                                                                                                                     |                                                                                                                                       AGTACGTTCTGTTGTATCTTGTGACACCGACTTTGGTGTATGTGCGCACTGCTACGGTCGTGACCTGGCGCGTGGCCACATCATCAACAAGGGTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAACCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGCATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAGCAACGTGAAGTCGGTTGTGAACTCCAG  >  NC_000913/4187992‑4188258

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