Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,187,838	A→G	D830G (GAC→GGC)	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,187,838	0	A	G	100.0%	208.7 / NA	55	D830G (GAC→GGC)	rpoC	RNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (30/25);  total (30/25)

CACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGACGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACGAACAGTGGTGTGACCTGCTGGAAGAG  >  NC_000913/4187709‑4187971                                                                                                                                  |                                                                                                                                      cACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAaga                                                                                                                              >  2:541412/1‑139 (MQ=255) cACTGAAAACTGCAAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAaga                                                                                                                              >  1:671040/1‑139 (MQ=255)      aaaaCTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGc                                                                                                                         >  2:113108/1‑139 (MQ=255)        aaCTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTg                                                                                                                       >  1:1055923/1‑139 (MQ=255)         aCTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTgc                                                                                                                      >  2:804581/1‑139 (MQ=255)                         aCCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACTggg                                                                                                      <  1:1589428/139‑1 (MQ=255)                            tGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTgcg                                                                                                                     <  1:939148/121‑1 (MQ=255)                            tGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTgcg                                                                                                                     >  2:939148/1‑121 (MQ=255)                            tGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCg                                                                                                   <  1:265940/139‑1 (MQ=255)                              aCTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCgtg                                                                                                 >  2:214176/1‑139 (MQ=255)                                 cgtcgtCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTaa                                                                                              <  1:520475/139‑1 (MQ=255)                                     gtcTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAActgc                                                                                          >  2:920963/1‑139 (MQ=255)                                           ttGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGTGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGa                                                                                    <  1:26965/139‑1 (MQ=255)                                               cGTGGCGCATGACCTGGTGGTTACCGAAGACGATCGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGtt                                                                                <  2:1055923/139‑1 (MQ=255)                                                    cgcAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGa                                                                                    <  1:1601529/130‑1 (MQ=255)                                                    cgcAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGa                                                                                    >  2:1601529/1‑130 (MQ=255)                                                       aGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTgcg                                                                                                                     <  1:837940/94‑1 (MQ=255)                                                       aGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTgcg                                                                                                                     >  2:837940/1‑94 (MQ=255)                                                        ggACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCg                                                                       >  1:88029/1‑139 (MQ=255)                                                          aCCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCggg                                                                     <  1:113108/139‑1 (MQ=255)                                                          aCCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCggg                                                                     <  1:214176/139‑1 (MQ=255)                                                           ccTGGTGGTTACCGCAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGGAGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGt                                                                    >  2:1000980/1‑139 (MQ=255)                                                                   ttACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTgcg                                                                                                                     <  2:817809/82‑1 (MQ=255)                                                                   ttACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTgcg                                                                                                                     >  1:817809/1‑82 (MQ=255)                                                                    tACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGc                                                                                                                            <  2:957221/74‑1 (MQ=255)                                                                    tACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGc                                                                                                                            >  1:957221/1‑74 (MQ=255)                                                                         aaGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCt                                                      >  1:600594/1‑139 (MQ=255)                                                                             cGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGtt                                                  <  2:671040/139‑1 (MQ=255)                                                                             cGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGATGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGtt                                                  <  1:947210/139‑1 (MQ=255)                                                                                             gAAGGTATGATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACgctgc                                  <  1:1000980/139‑1 (MQ=255)                                                                                             gAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTa                                                                   <  2:473005/106‑1 (MQ=255)                                                                                             gAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTa                                                                   >  1:473005/1‑106 (MQ=255)                                                                                             gAAGGTATCATGATGACTACGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACgctgc                                  >  2:1084519/1‑139 (MQ=255)                                                                                                   aTCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATc                                                                                                                >  2:609334/1‑55 (MQ=255)                                                                                                   aTCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATc                                                                                                                <  1:609334/55‑1 (MQ=255)                                                                                                    tCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGc                                                                                                                            >  2:739254/1‑42 (MQ=255)                                                                                                    tCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGc                                                                                                                            <  1:739254/42‑1 (MQ=255)                                                                                                    tCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTgcg                                                                                                                     <  2:371741/49‑1 (MQ=255)                                                                                                    tCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAACAGCCGCTgcg                                                                                                                     >  1:371741/1‑49 (MQ=255)                                                                                                     cATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCa                              <  2:945143/135‑1 (MQ=255)                                                                                                     cATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACGaa                          <  1:1431625/139‑1 (MQ=255)                                                                                                     cATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACTGGATCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACGaa                          >  1:1076358/1‑139 (MQ=255)                                                                                                     cATGATGACTCCGGTTAGCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAAGTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCa                              >  1:945143/1‑135 (MQ=255)                                                                                                      atgatgACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACGAAc                         >  1:46749/1‑139 (MQ=255)                                                                                                           gaCTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACGAACAgtgg                    >  1:716313/1‑139 (MQ=255)                                                                                                            aCTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCg                                                                                                                          >  2:752596/1‑36 (MQ=255)                                                                                                            aCTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCg                                                                                                                          <  1:752596/36‑1 (MQ=255)                                                                                                            aCTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACGAACAgtggt                   >  2:357633/1‑139 (MQ=255)                                                                                                                cGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACgctgc                                  >  2:753137/1‑120 (MQ=255)                                                                                                                cGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACgctgc                                  <  1:753137/120‑1 (MQ=255)                                                                                                                     aTCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTaa                                                                                              >  1:445833/1‑55 (MQ=255)                                                                                                                     aTCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTaa                                                                                              <  2:445833/55‑1 (MQ=255)                                                                                                                     aTCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACGAACAGTGGTGTGACctgc          >  2:1573794/1‑139 (MQ=255)                                                                                                                             tggtggCGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTCCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACGAACAGTGGTGTGACCTGCTGGAagag  >  1:34753/1‑139 (MQ=255)                                                                                                                               gtgGCGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTa                                                                   <  1:1248949/72‑1 (MQ=255)                                                                                                                               gtgGCGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTa                                                                   >  2:1248949/1‑72 (MQ=255)                                                                                                                                  |                                                                                                                                      CACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGACGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACGAACAGTGGTGTGACCTGCTGGAAGAG  >  NC_000913/4187709‑4187971

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