Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,186,085 C→A P246Q (CCG→CAG)  rpoC → RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,186,0850CA97.9% 135.3 / NA 47P246Q (CCG→CAG) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (21/25);  minor base G (0/1);  total (21/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATC  >  NC_000913/4185952‑4186221
                                                                                                                                     |                                                                                                                                        
cattaaGAGCTGAACGAAAACAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCAGCCAGATCt                                                                                                                                 <  2:841362/139‑1 (MQ=255)
 gTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGCATGATCCTGACCGTTCTGCCGGTACTGCAGCAAGATCTg                                                                                                                                >  1:90898/1‑143 (MQ=255)
  tGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCAGCCAGATCTGc                                                                                                                               <  1:850607/143‑1 (MQ=255)
  gcaaaaggtGAACGAAACTAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCAGCCAGATCTGc                                                                                                                               <  2:840891/135‑1 (MQ=255)
         cTGAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCAGCCAGATCTGCGTCCGCt                                                                                                                        <  2:940575/143‑1 (MQ=255)
            aaCGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTATGGTAACAAACCAGAGCGGATGATCCTGACCGTTCTGCCGGTACCGCAGCCAGATCTGCCTACGCTGGt                                                                                                                     >  1:133632/1‑143 (MQ=255)
                  aCCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCt                                                                                                               <  2:151276/143‑1 (MQ=255)
                    cgaCTCCGAAGCCAAGCGTAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGGAACAAACCCGAGTACTTTATCCTGTCCGTTCTGCCGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTgg                                                                                                             <  2:548651/141‑1 (MQ=255)
                         ccGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGAtggt                                                                                                        >  2:925564/1‑143 (MQ=255)
                         ccGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCAGCCAGATCTGCGTCCGCTGGGTCCGCTGGAtggt                                                                                                        >  2:173437/1‑143 (MQ=255)
                           gAAACCAAGCGTAAAAAGCTGACCAATCGTATCAAACTGCTAGAAGAATTCGTTAAGTCTGGTAACAAACAAGAGTGGATGATCCTGACCGTCCTGCAGGAACTGCAGCCAGATCTCCGTCCCCGCGTTCCTCTGGAtgcggg                                                                                                      >  1:67028/1‑139 (MQ=255)
                            aaaCCAAGCGTAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGATggtggt                                                                                                     >  2:674206/1‑143 (MQ=255)
                                    cGTAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTcgc                                                                                             <  2:430231/143‑1 (MQ=255)
                                      tAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGa                                                                                           <  2:799834/143‑1 (MQ=255)
                                      tAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGa                                                                                           >  2:240028/1‑143 (MQ=255)
                                      gaaaaaGCTGACCAAGCGTATCACACTGCTGGAAGCGTTCTTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCAGCCACAGCTTCTTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGa                                                                                           <  1:291918/142‑1 (MQ=255)
                                             cTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGa                                                                                    <  2:50889/143‑1 (MQ=255)
                                             cTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTGCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGa                                                                                    <  2:513946/143‑1 (MQ=255)
                                               gACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGCCCGTTCTGCCGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGAcc                                                                                  >  2:209482/1‑143 (MQ=255)
                                                aCCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGGTCTGCCGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGAGACTTGTGGCCt                                                                                 >  1:301872/1‑143 (MQ=255)
                                                  cAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCGGa                                                                               >  2:921867/1‑143 (MQ=255)
                                                   aaGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCCGACCGTTCTGCCGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGAGGGTGGTCGTTCCGCGACTTCTGGccatga                                                                              >  1:88599/1‑139 (MQ=255)
                                                          tCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGCTGGTGGTCGTTTCGCGAATTCTGACCTGAACGATCTg                                                                       >  1:750568/1‑143 (MQ=255)
                                                           cAAACTGCCGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGt                                                                      <  2:92228/143‑1 (MQ=255)
                                                                   tgGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTcgc                                                              <  2:545799/143‑1 (MQ=255)
                                                                    ggAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTcgcg                                                             <  1:407192/143‑1 (MQ=255)
                                                                       aGCGTTCTTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCAGCCGGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCa                                                          <  2:940115/143‑1 (MQ=255)
                                                                        gcgtgcggggggggtGGTGACAGACGGGAGTGGAGGATCCTGGCCGTGCGGCCGGTACTGCGGCCAGATCTGCGTCGGCTGGTTCCGCTGGATGGTGGTCGTGTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCAt                                                         <  1:28499/129‑1 (MQ=255)
                                                                             cgttcAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAAcc                                                    <  2:335798/143‑1 (MQ=255)
                                                                              gttcAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGCCCTGCACGATCTGTATCGTCGCGTCATGAACCg                                                   >  1:235127/1‑143 (MQ=255)
                                                                               ttcAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGt                                                  >  2:621534/1‑143 (MQ=255)
                                                                                 cattCTGGTAACAAGTTAGAGTGGTTGATGCTGATCGTTCTGCCGGTACTGCAGCCAGATCTGCGTTCGTTTGTTCTGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTaa                                                <  1:541076/140‑1 (MQ=255)
                                                                                   gTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGAAACAAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCATTTCGCGACTTCTGACCTGAACGATCTGTACCGGCACGTCATGAAcctgagaa                                              >  1:99017/1‑137 (MQ=255)
                                                                                    tCTGGTAACGAGCCTGAGTGGATGAGCCTGACCGTTCTGCGGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTaacaa                                             <  1:300557/143‑1 (MQ=255)
                                                                                        gTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGt                                         >  2:826416/1‑143 (MQ=255)
                                                                                               aCCGGAGTGGATGATCCGGACCGTTCTGCCGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAAc                                  <  2:492157/143‑1 (MQ=255)
                                                                                               aCCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAAc                                  <  2:722455/143‑1 (MQ=255)
                                                                                                 cAGAGTGGATGTGCTTGACCGGTCTGCCGTTACTGCAGCCAGAGCTGCGTCCGCGGGTTCCGCTGGATGGTGGTCGGTTCGCGAGTTCTGACCGGAGCGATCTGTATCGTCGCGTCATTAACCGTGACAACCGTCTGAAACGt                                <  1:6187/143‑1 (MQ=255)
                                                                                                    agTGGATGATCCTGTCCGTTCTGCCGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTctg                             <  1:393375/143‑1 (MQ=255)
                                                                                                          tgatCCTGACCGTTCTGCCGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCAACTTCTGACCTGCAGGATCTGTATCGTCGCCTCCTTAAACGTAACAACCGTCTGAAACGTCTGCAGGaa                       >  1:61804/1‑142 (MQ=255)
                                                                                                           gatCCTGACCGTTCTGCCGGTACTGCAGCCAGAACTGCGTCCGCTGGTTCCGCTGGATGGTGATCGTTTCGCGACTTCTGACCTTAACGATCTGTCTCGTCGCGTCATTAACCGTCACAAGCGTTTGCAACGTCTGCTGGATc                      >  2:631030/1‑143 (MQ=255)
                                                                                                               cTGACCGTTCTGCCGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAATGATCTGTATCGTCGCGTCATTAACCGAAACAACCGTCTGAAACGTCTGCGGGAttttgc                  >  2:671143/1‑138 (MQ=255)
                                                                                                                      ttCTGCCGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCg           <  1:842999/143‑1 (MQ=255)
                                                                                                                        cggCCGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGa         <  2:525575/141‑1 (MQ=255)
                                                                                                                          gCCGGTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCACGACTTCTGACCTGAACGATCTGTATCGTCGAGTCATTAACCGTAACAACAGTCTGAAACGTATGCTGTGTCTTAATGAGCCGGAca       >  1:727825/1‑143 (MQ=255)
                                                                                                                             ggTACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGCACGCCCTGTATCGTCGCGTCATTCACTGAAAAAATCGTCCGAACCGACTGCTGGCTCTGGCTGCGCCGGAcatca    >  1:567280/1‑143 (MQ=255)
                                                                                                                               tACTGCAGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTCTCGCTACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACatcatc  <  1:674206/143‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                        
CGTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATC  >  NC_000913/4185952‑4186221

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: