Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,189,512 +C coding (4164/4224 nt) rpoC → RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,189,5111.C95.3% 125.9 / 0.4 43coding (4163/4224 nt)rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (21/20);  total (22/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

ACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTG  >  NC_000913/4189373‑4189648
                                                                                                                                            |                                                                                                                                         
tccaaCTCCTCGGCCTGGCAGCGATCGTTATTCTGTGTCGGCTGTGCCCGGCCGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCTGCTGCACCGCAGGTGACTGCAGAAGAccgc                                                                                                                                       <  2:34783/140‑4 (MQ=255)
        cgcgGCCTGAAAGGGAACGTTATCGTGGGGCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACCGCATCTGCCa                                                                                                                                <  2:489861/143‑1 (MQ=255)
         gcgGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACCGCATCTGCCAg                                                                                                                               >  1:610072/1‑143 (MQ=255)
          cgGCCTGAAAGAGAACGTTATCGTGGGTCGTCGGATCCCGGCAGGTACCGGGTAGGCGTACCACCTGGTTCGTATGCGTCGCCGTGCGGCGGGTGAAGCTCCGGCTGCACCGCGGGTGACTGCCGAAGACCGCATCTGccggg                                                                                                                              >  1:404972/1‑140 (MQ=255)
             tcgggcggggaACGGTAGCGGGGGTCGCCTGATCCCGGCAGGTACCGGTTAGGCGTACCACCAGGATCGTCAGCGGCGCCGTGCGGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACCGCAGCTGCCAGCCTg                                                                                                                           <  2:360990/134‑1 (MQ=255)
              cagAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGCGAAGCTCCGTTTTTACCGCAGGTTACTGCAGAAGACCGCATCTGCCAGCCTgg                                                                                                                          <  1:862894/141‑1 (MQ=255)
               tGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACCGCATCTGCCAGCCTGGc                                                                                                                         >  1:158835/1‑143 (MQ=255)
               tGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCGCCGGCTGCACCAGAGGTGACCGACGCAGACCGCAACCGCCAGCATGGc                                                                                                                         >  2:327261/1‑143 (MQ=255)
                   agagAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACCGCATCTGCCAGCCTGGCGGaa                                                                                                                     >  1:306374/1‑143 (MQ=255)
                        aCGTTATCGGGGGTCGTCTGGTCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACCGCATCTGCCAGCCTGGCAGAActgct                                                                                                                <  1:526970/143‑1 (MQ=255)
                           ttATCCGGGGTCGTCTGATGCGGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACCGCATCTGCCAGCCTGGCAGAACTGCTGaa                                                                                                             <  2:748040/143‑1 (MQ=255)
                              tCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACCGCATCTGCCAGCCTGGCAGACCTGCTGAACGc                                                                                                          >  1:497678/1‑143 (MQ=255)
                                gTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCCGGTGAAGCTCCGGCGGCACCGCAGGTGACTGCAGAAGACCTCATCTGCCAACCTGGCAGAAATGCTGAACGCAg                                                                                                        >  1:58934/1‑143 (MQ=255)
                                      cgtcTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGCATCGTATGCGTCGCCGTGCTGCGCGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACCGCATCTGCCAGCCCGGCAGAACTGCTGAACGCAGGTCTgg                                                                                                  <  2:863589/143‑1 (MQ=255)
                                       ctcTGCTCCGGGCAGGGACCGGTCACGCGTACCACCAGGAGCGTATGCGTCGCCGTGCTGCGGGTGACGGCCCGGCTGCTCCGCAGGTGACTGCAGACGCCCGCAGCTGCCAGCCTGGCAGAACTGCTGAGCGCAGGTCTggg                                                                                                 <  2:468216/142‑1 (MQ=255)
                                        tcTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCCCCGTGCTGCAGGAGAAGCTCCGGATGCACCGCAGGTGATTGCAGAAGAACGCACCCGCCAGCCAAGCACACCTGCTGAACGCaactagggga                                                                                                >  2:711655/1‑134 (MQ=255)
                                          tGATCCCGGCGGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCgg                                                                                              <  2:548246/143‑1 (MQ=255)
                                            agcccGGCGGGGACCGGTTACGCGGACCACGAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGCGACTGCAGAAGACCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGtt                                                                                            <  1:270044/141‑1 (MQ=255)
                                              ccggggaGGTACCGGGTACGCGTACCACCAGGAGCGTATGCGTCGCTGTGGTGCGGACGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCt                                                                                          <  1:416376/137‑1 (MQ=255)
                                                ccccgcgcaCCGGTGACGCGTACCACGAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGa                                                                                        <  2:581887/135‑1 (MQ=255)
                                                cGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCCGCACCGCAGGTGACTGCAGAAGACCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCCGa                                                                                        >  2:545821/1‑143 (MQ=255)
                                                   cAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATaa                                                                                     >  2:768630/1‑143 (MQ=255)
                                                      ttACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGTGTGAAGCTCCGGTTGTACCGCAGGTGACTGCAGAAGACCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGa                                                                                  <  1:662498/142‑1 (MQ=255)
                                                         ccGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCCACGGCTGCGCAGGCGGTGACCGCAGAAGACCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGAGGTTCTGATaaagcata                                                                               >  2:789639/1‑137 (MQ=255)
                                                                          caGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGAATGGAAAGTGTTTAGGACGAGTAATCGATAATgcgcccat                                                              >  1:164847/1‑139 (MQ=255)
                                                                          caGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCGCCGGCTGCACCGCAGGTGACTGCAGAAGACCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAAt                                                              <  1:836874/143‑1 (MQ=255)
                                                                           aGGATCGTAGGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATa                                                             <  2:523221/143‑1 (MQ=255)
                                                                                       gcctccGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAAccc                                                 <  2:704609/139‑1 (MQ=255)
                                                                                       gTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGACAACGAGTAATCGTTAATCCGCAACTAACGTAAAAAccc                                                 >  1:439365/1‑143 (MQ=255)
                                                                                          gCCGTTCTGCGGGTCAAGCTTCGCCTGCACCCCAGGTGACTGCAGAAGACCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCt                                              <  1:615342/143‑1 (MQ=255)
                                                                                           ccGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTGAAAGCCCGCtt                                             >  1:433527/1‑143 (MQ=255)
                                                                                                gtgcGGGTTAAGCTCCGGCTGCACAGCAGGTGACTGCAGAAGACCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTcggcg                                        <  1:301669/142‑1 (MQ=255)
                                                                                                   cGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTCAAAACCCGCTTCGGCgggg                                     >  2:483346/1‑142 (MQ=255)
                                                                                                   cGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGt                                     <  1:888890/143‑1 (MQ=255)
                                                                                                         aagcctcGGCTGCATCGCAGGTGACTGCAGAAGACCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGCGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGttttttt                               <  1:239325/137‑1 (MQ=255)
                                                                                                           gCTCCGGCTGCACCGCAGGTGACTGCAGAAGACCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTAt                             >  2:521512/1‑143 (MQ=255)
                                                                                                              ccGGCTGCACCGCAGGTGACTGCAGAAGACCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATggg                          >  1:91664/1‑143 (MQ=255)
                                                                                                              ccGGCTGCACCGCAGGTGACTGCAGAAGACCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGCTTTTTTATggg                          >  1:389035/1‑143 (MQ=255)
                                                                                                                ggCTGCACCGCAGGTGACTGCAGAAGACCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGCTCCTGATAACGAGTACTCGTTAATCCGCAAATAACGTAAACACCCGCTTCGGACGGTTTTTTTATgggac                        >  2:680132/1‑141 (MQ=255)
                                                                                                                   tGCACCGCAGGTGACTGCAGAAGACCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGACTTTTTATGGGGGGAg                     >  2:272403/1‑143 (MQ=255)
                                                                                                                      aCCGCAGGTGACTGCAGAAGACCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAACAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGttt                  >  2:267239/1‑143 (MQ=255)
                                                                                                                       ccGCAGGTGACTGCAGAAGACCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTa                 >  1:525419/1‑143 (MQ=255)
                                                                                                                        cGCAGGTGACGGCAGAAGACCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATTACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAg                <  1:789639/143‑1 (MQ=255)
                                                                                                                                  tGCAGAAGACCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCa      <  1:768630/143‑1 (MQ=255)
                                                                                                                                      gaagacCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGCGCATTTg  >  1:455170/1‑143 (MQ=255)
                                                                                                                                            |                                                                                                                                         
ACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTG  >  NC_000913/4189373‑4189648

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: