Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,188,171	C→T	A941V (GCA→GTA)	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,171	0	C	T	100.0%	130.8 / NA	38	A941V (GCA→GTA)	rpoC	RNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (20/18);  total (20/18)

CACTGCTACGGTCGTGACCTGGCGCGTGGCCACATCATCAACAAGGGTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAACCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGCATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAGCAACGTGAAGTCGGTTGTGAACTCCAGCGGTAAACTGGTTATCACTTCCCGTAATACTGAACTGAAACTG  >  NC_000913/4188038‑4188301                                                                                                                                      |                                                                                                                                   cACTGCTACGGTCGTGACCTGGCGCGTGGCCACATCATCAACAAGGGTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAACCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGTAtctc                                                                                                                               >  2:284949/1‑139 (MQ=255) cACTGCTACGGTCGTGACCTGGCGCGTGGCCACATCATCAACAAGGGTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAACCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGTAtctc                                                                                                                               >  1:192502/1‑139 (MQ=255) cACTGCTACGGTCGTGACCAGGCGCGTGGTCACATCATCAACAAGGGTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAACCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGTAtctc                                                                                                                               >  2:407368/1‑139 (MQ=255)   cTGCTACGGTCGTGACCTGGCGCGTGGCCACATCATCAACAAGGGTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAACCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGTATCTCGt                                                                                                                             >  1:296581/1‑139 (MQ=255)    tGCTACGGTCGTGACCTGGCGCGTGGCCACATCATCAACAAGGGTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAACCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGTATCTCGTg                                                                                                                            >  2:65938/1‑139 (MQ=255)             cGTGACCTGGCGCGTGGCCACATCATCAACAAGGGTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAACCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGTATCTCGTGCGGctgctg                                                                                                                   <  2:598521/139‑1 (MQ=255)              gTGACCTGGCGCGTGGCCACATCATCAACAAGGGTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAACCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGTATCTCGTGCGGCTGCTGa                                                                                                                  >  1:174783/1‑139 (MQ=255)              gTGACCTGGCGCGTGGCCACATCATCAACAAGGGTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAACCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGTATCTCGTGCGGCTGCTGa                                                                                                                  <  1:65938/139‑1 (MQ=255)              gTGACCTGGCGCGTGGCCACATCATCAACAAGGGTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAACCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGTATCTCGTGCGGCTGCTGa                                                                                                                  >  2:248689/1‑139 (MQ=255)              gTGACCTGGCGCGTGGCCACATCATCAACAAGGGTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAACCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGTATCTCGTGCGGCTGCTGa                                                                                                                  >  1:87142/1‑139 (MQ=255)                     ggCGCGTGGCCACATCATCAACAAGGGTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAACCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGTATCTCGTGCGGCTGCTGAATCCAGc                                                                                                           >  2:84153/1‑139 (MQ=255)                                             gggTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAACCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGTATCTTGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGc                                                                                   >  1:254695/1‑139 (MQ=255)                                             gggTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAACCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGTATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGc                                                                                   >  1:448430/1‑139 (MQ=255)                                                        tCGGTGTTATCGCGGCACAGTCCATCGGTGAACCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGTATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAg                                                                        <  1:84153/139‑1 (MQ=255)                                                        tCGGTGTTATCGCGGCACAGTCCATCGGTGAACCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGTATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAg                                                                        <  2:448430/139‑1 (MQ=255)                                                          ggTGTTATCGCGGCACAGTCCATCGGTGAACCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGTATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAGCa                                                                      >  1:281909/1‑139 (MQ=255)                                                                 tCGCGGCACAGTCCATCGGTGAACCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGTATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAGCAACGTGaa                                                               <  1:108741/139‑1 (MQ=255)                                                                      gCACAGTCCATCGGTGAACCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGTATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAGCAACGTGAAGTCgg                                                          >  2:286696/1‑139 (MQ=255)                                                                           gTCCATCGGTGAACCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGTATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAGCAACGTGAAGTCGGTtgtg                                                     >  2:563307/1‑139 (MQ=255)                                                                           gTCCATCGGTGAACCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGTATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAGCAACGTGAAGTCGGTtgtg                                                     <  1:286696/139‑1 (MQ=255)                                                                            tCCATCGGTGAACCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGTATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAGCAACGTGAAGTCGGTTGTGa                                                    <  1:192323/139‑1 (MQ=255)                                                                                tCGGTGAACCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGTATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAGCAACGTGAAGTCGGTTGTGAACTc                                                <  1:457025/139‑1 (MQ=255)                                                                                tCGGTGAACCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGTATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAGCAACGTGAAGTCGGTTGTGAACTc                                                <  1:325188/139‑1 (MQ=255)                                                                                tCGGTGAACCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGTATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAGCAACGTGAAGTCGGTTGTGAACTc                                                <  2:296581/139‑1 (MQ=255)                                                                                      aaCCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGTATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATc                                                                                >  2:223334/1‑101 (MQ=255)                                                                                      aaCCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGTATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATc                                                                                <  1:223334/101‑1 (MQ=255)                                                                                      aaCCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGTATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAGCAACGTGAAGTCGGTTGTGAACTCCAGCgg                                          >  1:439618/1‑139 (MQ=255)                                                                                        ccGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGTATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAGCAACGTGAAGTCGGTTGTGAACTCCAGCGGTa                                        <  1:283398/139‑1 (MQ=255)                                                                                               cacaGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGTATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAGCAACGTGAAGTCGGTTGTGAACTCCAGCGGTAAACTGGt                                 <  2:174783/139‑1 (MQ=255)                                                                                                  aGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGTATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAGCAACGTGAAGTCGGTTGTGAACTCCAGCGGTAAACTGGTTAt                              <  1:284949/139‑1 (MQ=255)                                                                                                      gACCATGCGTACGTTCCACATCGGTGGTGCGGTATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAGCAACGTGAAGTCGGTTGTGAACTCCAGCGGTAAACTGGTTATCACt                          <  1:13789/139‑1 (MQ=255)                                                                                                       aCCATGCGTACGTTCCACATCGGTGGTGCGGTATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATc                                                                                >  2:481621/1‑84 (MQ=255)                                                                                                       aCCATGCGTACGTTCCACATCGGTGGTGCGGTATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATc                                                                                <  1:481621/84‑1 (MQ=255)                                                                                                            gCGTACGTTCCACATCGGTGGTGCGGTATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAGCAACGTGAAGTCGGTTGTGAACTCCAGCGGTAAACTGGTTATCACTTCCCGt                    >  2:579954/1‑139 (MQ=255)                                                                                                                 cGTTCCACATCGGTGGTGCGGTATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAGCAACGTGAAGTCGGTTGTGAACTCCAGCGGTAAACTGGTTATCACTTCCCGTAATac               <  1:212607/139‑1 (MQ=255)                                                                                                                      cacaTCGGTGGTGCGGTATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAGCAACGTGAAGTCGGTTGTGAACTc                                                <  2:249442/101‑1 (MQ=255)                                                                                                                      cacaTCGGTGGTGCGGTATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAGCAACGTGAAGTCGGTTGTGAACTc                                                >  1:249442/1‑101 (MQ=255)                                                                                                                              tggtgCGGTATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAGCAACGTGAAGTCGGTTGTGAACTCCAGCGGTAAACTGGTTATCACTTCCCGTAATACTGaactgaaactg  >  1:519312/1‑139 (MQ=255)                                                                                                                                      |                                                                                                                                   CACTGCTACGGTCGTGACCTGGCGCGTGGCCACATCATCAACAAGGGTGAAGCAATCGGTGTTATCGCGGCACAGTCCATCGGTGAACCGGGTACACAGCTGACCATGCGTACGTTCCACATCGGTGGTGCGGCATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAGCAACGTGAAGTCGGTTGTGAACTCCAGCGGTAAACTGGTTATCACTTCCCGTAATACTGAACTGAAACTG  >  NC_000913/4188038‑4188301

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