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_0009134,187,8380AG100.0% 158.3 / NA 42D830G (GAC→GGC) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (20/22);  total (20/22)

TACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGACGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACGAACAGTGGTGTGACC  >  NC_000913/4187704‑4187960
                                                                                                                                      |                                                                                                                          
tACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGGCGtt                                                                                                                        >  2:579381/1‑139 (MQ=255)
             aCTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTgc                                                                                                           >  2:65795/1‑139 (MQ=255)
              cTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTgcg                                                                                                          <  1:911204/139‑1 (MQ=255)
                  gAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGAt                                                                                                      >  1:415852/1‑139 (MQ=255)
                   aaCTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATc                                                                                                     <  1:579381/139‑1 (MQ=255)
                   aaCTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATc                                                                                                     <  1:807462/139‑1 (MQ=255)
                        cGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTa                                                                                                <  1:214538/139‑1 (MQ=255)
                           tttccTGGCTCGTCGTCTGGTTGGAGTTGGGGAGGACCTGGTGGGTAACGAAAACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACTg                                                                                             <  1:522477/136‑1 (MQ=255)
                           ttACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACTg                                                                                             <  1:771664/139‑1 (MQ=255)
                           ttACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACTg                                                                                             <  2:631715/139‑1 (MQ=255)
                                 gACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCgt                                                                                       >  1:824752/1‑139 (MQ=255)
                                        cgtcTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAActg                                                                                >  2:605191/1‑139 (MQ=255)
                                                              aCCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCggg                                                          <  2:599569/139‑1 (MQ=255)
                                                                    tggtTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTActgc                                                    <  1:605191/139‑1 (MQ=255)
                                                                                     tgtgGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCcgc                                   >  1:448778/1‑139 (MQ=255)
                                                                                     tgtgGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCcgc                                   >  1:518887/1‑139 (MQ=255)
                                                                                        ggTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCa                                <  2:824752/139‑1 (MQ=255)
                                                                                            cccATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCa                                <  1:728787/135‑1 (MQ=255)
                                                                                            cccATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCa                                >  2:728787/1‑135 (MQ=255)
                                                                                              cATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACgc                          >  1:229596/1‑139 (MQ=255)
                                                                                               aTGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACgct                         <  2:448778/139‑1 (MQ=255)
                                                                                                       aTCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTgcg                                                                                                          <  2:1072033/50‑1 (MQ=255)
                                                                                                       aTCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTgcg                                                                                                          >  1:1072033/1‑50 (MQ=255)
                                                                                                        tCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCgtgt                                                                                     >  1:768946/1‑70 (MQ=255)
                                                                                                        tCATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCgtgt                                                                                     <  2:768946/70‑1 (MQ=255)
                                                                                                         cATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGc                                                                                                                 >  1:634139/1‑41 (MQ=255)
                                                                                                         cATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGc                                                                                                                 <  2:634139/41‑1 (MQ=255)
                                                                                                         cATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACt                                                                                              <  1:684066/60‑1 (MQ=255)
                                                                                                         cATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACt                                                                                              >  2:684066/1‑60 (MQ=255)
                                                                                                         cATGATGACTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACGaa               <  2:788165/139‑1 (MQ=255)
                                                                                                               gaCTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACGAACAgtgg         >  1:313408/1‑139 (MQ=255)
                                                                                                                aCTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACTgg                                                                                            <  2:737820/55‑1 (MQ=255)
                                                                                                                aCTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACTgg                                                                                            >  1:737820/1‑55 (MQ=255)
                                                                                                                 cTCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACGAACAgtggtg       >  2:1047901/1‑139 (MQ=255)
                                                                                                                  tCCGGTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACGAACAGTGgtgt      >  1:560385/1‑139 (MQ=255)
                                                                                                                      gTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACGAACAGTGGTGTGAcc  >  1:133352/1‑139 (MQ=255)
                                                                                                                      gTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACGAACAGTGGTGTGAcc  <  2:652736/139‑1 (MQ=255)
                                                                                                                      gTTATCGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACGAACAGTGGTGTGAcc  <  1:512459/139‑1 (MQ=255)
                                                                                                                           cGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCa                                <  2:166192/104‑1 (MQ=255)
                                                                                                                           cGAGGGTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCa                                >  1:166192/1‑104 (MQ=255)
                                                                                                                              gggTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCcg                                    <  1:850936/97‑1 (MQ=255)
                                                                                                                              gggTGGTGGCGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCcg                                    >  2:850936/1‑97 (MQ=255)
                                                                                                                                      |                                                                                                                          
TACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGACGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACGAACAGTGGTGTGACC  >  NC_000913/4187704‑4187960

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: