Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_bme_pgi 4,188,282 G→A R978H (CGT→CAT)  rpoC → RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi4,188,2820GA97.9% 136.1 / NA 47R978H (CGT→CAT) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (25/21);  minor base C (1/0);  total (26/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.99e-01

CTGACCATGCGTACGTTCCACATCGGTGGTGCGGCATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAGCAACGTGAAGTCGGTTGTGAACTCCAGCGGTAAACTGGTTATCACTTCCCGTAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCGGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACACCATGC  >  NC_000913_3_bme_pgi/4188137‑4188425
                                                                                                                                                 |                                                                                                                                               
cTGACCATGCGTACGTTCCACATCGGTGGTGCGGCATCGCGTGCCGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAGCAACGTGAAGTCGGTTGTGAACTCCAGCGGTAAACTGGTTATCACTTCCCATaa                                                                                                                                              <  2:630882/149‑1 (MQ=255)
  taCCATGCGGACGTTCCACATCGGTGGTGCGGCATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAGCAACGTGAAGTCGGTTGTGAACTCCAGCGGTAAACTGGTTATCACTTCCCataata                                                                                                                                            <  1:395952/148‑1 (MQ=255)
       tGCGTACGTTCCACATCGGTGGTGCGGCATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAGCAACGTTAAGTCGGTTGTGAACTCCAGCGGTAAACTGGTTATCACTTCCCATAATactgaa                                                                                                                                       <  2:1005471/149‑1 (MQ=255)
       tGCGTACGTTCCACATCGGTGGTGCGGCATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAGCAACGTGAAGTCGGTTGTGAACTCCAGTTTTAAACTGGTTATCACTTCCCATAATactgac                                                                                                                                       >  2:465215/1‑148 (MQ=255)
         cGTACGTTCCACATCGGTGGTGCGGCATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAGCAACGTGAAGTCGGTTGTGAACTCCAGCGGCAAACTGGTTATCACTTCCCACAATactgaact                                                                                                                                     >  2:97149/1‑149 (MQ=255)
          gTACGTTCCACATCGGTGGCACGGCATCACGTGCGGCTGCTAAATCCAGCATCCAAGTGCAAAACAAAGGTAGCATCAAGCTCAGCACCGTGAAGTCGGTTGTGAACTCCAGCGGTAAACTGGTTATCTCTCCCCACAATactgaactg                                                                                                                                    >  2:688158/1‑149 (MQ=255)
           tACGTTCCACATCGGTGGTGCGGCATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTTTGAAAAGTGAAGACGTTTGTGAACTCCAGCGGTAAACTGGTTATCACTTTCCATAATactgaactga                                                                                                                                   >  1:234616/1‑149 (MQ=255)
               ttCCACATCGGTGGTGCGGCATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAGCAACGTGAAGTCGGTTGTGAACTCCAGCGGTAAACTGGTTATCACTTCCCACAATCGCGcactgaaact                                                                                                                               >  2:143809/1‑149 (MQ=255)
                  cacaTCGGGGGTGTGGCATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAGCAACGTGAAGTCGGTTGTGAACTCCAGCGGTAAACTGGTTATCACTTCCCATAATACTGaactgaaactgcc                                                                                                                            <  2:244100/149‑3 (MQ=255)
                   acaTCGGTGGTGCGGCATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTCGCATCAAGCTCAGCAACGTGAAGTCGGTTGTGAACTCCAGCGGTAAACTGGTTATCACTTCCCATAATACTGAACTGAAACTGATc                                                                                                                           <  1:989121/149‑1 (MQ=255)
                     aTCGGTGGTGCGGCATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAGCAACGTGAAGTCGGTTGTGAACTCCAGCGGTAAACTGGTTATCACTTCCCATAATACTGAACTGAAACTGATcga                                                                                                                         <  1:521050/149‑1 (MQ=255)
                           ggtgCGGCATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAGCAACGTGAAGTCGGTTGTGAACTCCAGCGGTAAACTGGTTATCACTTCCCATAATACTGAACTGAAACTGATCGACGAAtt                                                                                                                   <  1:779309/149‑1 (MQ=255)
                            gtgCGGCATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAGCAACGTGAAGTCGGTTGTGAACTCCAGCGGTAAACTGGTTATCACTTCCCATAATACTGAACTGAAACTGATCGACGAATTc                                                                                                                  >  1:572466/1‑149 (MQ=255)
                                ggCAACTCGTACCGCCGCCTAATCAAGCATACTAGTCAAAAACAAAGGTAGCGTTAAGCTCAGCAACGTGAAGTCTGCTGGGAACTCTAGCGGTAAACTGGTGATCACTTCCCATAATGCTGAACTGAAACTGATCGACGAATTCGGTc                                                                                                              <  2:541620/149‑1 (MQ=255)
                                   aTCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAGCAACGTGAAGTCGGTTGTGAACTCCAGCGGTAAACTGGTTATCACTTCCCATAATACTGAACTGAAACTGATCGACGAATTCGGTCGTa                                                                                                           >  1:421376/1‑149 (MQ=255)
                                    tctcGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAGCAACGTGAAGTCGGTTGTGAACTCCAGCGGTAAACTGGTTATCACTTCCCATAATACTGAACTGAAACTGATCAACGAATTCGGTCGTAc                                                                                                          >  1:56508/1‑149 (MQ=255)
                                                      tCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAGCAACGTGAAGTCGGTTGTGAACTCCAGCGGTAAACTGGTTATCACTTCCCATAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGt                                                                                        >  1:599240/1‑149 (MQ=255)
                                                       ccAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAGCAACGTGAAGTCGGTTGTGAACTCCAGCGGTTAACTGGTTATCACTTCCCATAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTAAAGAAAACTACAAAGTa                                                                                       >  1:553044/1‑149 (MQ=255)
                                                        cAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAGCAACGTGAAGTCGGTTGTGAACTCCAGCGGTAAACTGGTTATCACTTCCCATAATACTGAACTGAAACTGATCGACGAATTCGGTCCTACTAAAGAAAGCTACAAAGTAc                                                                                      >  2:685238/1‑149 (MQ=255)
                                                              ccAAGTGAAAAACAAAGGTAGCATCAAGCTCAGCAACGTTAAGTCGGTTGTGAACTCCAGCGGTAAACTGGTTATCACTTCCCATAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACg                                                                                <  2:922970/149‑1 (MQ=255)
                                                                aaGTGAAAAACAAAGGTAGCATCAAGCTCAGCAACGTGAAGTCGGTTGTGAACTCCAGCGGTAAACTGGTTATCACTTCCCATAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGt                                                                              >  1:966246/1‑149 (MQ=255)
                                                                    gAAAAACAAAGGTAGCATCAAGCTCAGCAACGTGAAGTCGGTTGTGAACTCCAGCGGTAAACTGGTTATCACTTCCCATAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCgg                                                                          <  2:845106/149‑1 (MQ=255)
                                                                          ccaaGGTAGGATCAAGCCCAGCAACGGCACGTCGGTTGTGAACTCCAGCGGTAAACTGGTTATCACTTCCCATAATACTGAACTGAAACTGATCGACGAATTCTGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCGGTACTgg                                                                    <  1:60730/147‑1 (MQ=255)
                                                                                tatgctCAAGGTCAGCAAGGTAAAGCGGGTTGTGAAATCCATCGGTAACCTGGGTATCACTTCCCATAATACGGAACTGAACCTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCGGTACTGGCGAAAg                                                              <  2:139298/144‑1 (MQ=255)
                                                                                      cAAGCTCAGCAACGTGAAGTCGGTTGTGAACTCCAGCGGTAAACTGGTTATCACTTCCCATAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCGGTACTGGCGAAAGGCGATg                                                        <  1:718061/149‑1 (MQ=255)
                                                                                      cAAGCTCAGCAACGTGAAGTCGGTTGTGAACTCCAGCGGTAAACTGGTTATCACTTCCCATAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCGGTACTGGCGAAAGGCGATg                                                        <  1:477146/149‑1 (MQ=255)
                                                                                           tCAGCAACGTGAAGTCGGTTGTGAACTCCAGCGGTAAACTGGTTATCCCTTCCCCTAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACGAAGTACCTTACGGTGCGGTACTGGCGAGAGGCGATGACGac                                                   >  2:633445/1‑148 (MQ=255)
                                                                                            cAGCAACGTGAAGTCGGTTGTGAACTCCAGCGGTAAACTGGTTATCACTTCCCATAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCGGTACTGGCGAAAGGCGATgcccacc                                                  >  1:730986/1‑143 (MQ=255)
                                                                                               cAACGTGAAGTCGGTTGTGAACTCCAGCGGTAAACTGGTTATCACTTCCCATAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCGGTACTGGCGAAAGGCGATGGCGAACAgg                                               >  2:798028/1‑149 (MQ=255)
                                                                                                aaCGTGAAGTCGGTTGTGAACTCCAGCGGTAAACTGGTTATCACTTCCCATAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACGAAGTACCTTAGGGTGCGGTACTGGCGAAAGGCGATGGCGAACAGGt                                              <  2:880892/149‑1 (MQ=255)
                                                                                                        gTCGGTTGTGAACTCCAGCGGTAAACTGGTTATCACTTCCCATAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCGGTACTGGCGAAAGGCGATGGTGAACAGGTTGCTggcg                                      <  1:878617/149‑1 (MQ=255)
                                                                                                        gTCGGTTGTGAACTCCAGCGGTAAACTGGTTATCACTTCCCATAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCGGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTggcg                                      >  1:145334/1‑149 (MQ=255)
                                                                                                         tCGGTTGTGAACTCCAGCGGTAAACTGGTTATCACTTCCCATAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCGTACGGTGCGGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTggcgg                                     <  2:421376/149‑1 (MQ=255)
                                                                                                            gatgtaAACACCCGCGGTAAACTGGTTATCAATTCCCATAATACTGAACTGAAACTGATCGACGAATGCGGTCGTACTAAAGAAAGCTACAAAGTACCGTACGGTGCGGTACTGCCGAACGAGCATGGCCACCATGTTGCGggcgcgga                                  >  2:448259/3‑145 (MQ=255)
                                                                                                            gTTGTGAACTCCAGCGGTAAACTGGTTATCACTTCCCATAATACTGAACTGAAACTGATCGACGAATTCGGTCGTATTAAAGCAAGCTACAAAGTACCTTACGGTGCGGTACTGGCGAAAGGGGATGGTGAACAGGTTGCTggcggggg                                  >  2:8191/1‑146 (MQ=255)
                                                                                                             ttGTGAACTCCAGCGGTAAACTGGTTATCACTTCCCATAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCGGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGaa                                 >  2:981478/1‑149 (MQ=255)
                                                                                                             ttGTGAACTCCAGCGGTAAACTGGTTATCACTTCCCATAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCGGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGaa                                 >  2:981552/1‑149 (MQ=255)
                                                                                                                   acaccAGCGGTAAAGTGGTGATCACTTCCCATAATACTGAACTGAAACGGATCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCGGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGtt                           <  2:60769/146‑1 (MQ=255)
                                                                                                                   aatCCAGCGGTAAACTGGTTATCACTTCCCATAAGACTGAACTGAAACTTATCGACGAATGCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCGGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGtt                           <  2:60837/147‑1 (MQ=255)
                                                                                                                           ggTAAACTGGTTATCCCTTCCCATAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCGGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTg                   <  1:171060/149‑1 (MQ=255)
                                                                                                                            gTAAACTGGTTATCACTTCCCATAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTCCCTTACGGTGCGGTACTGGCGAAAGGCGATGGGGAACAGGTTGCGGGCGGGGAAACCGTTGCGAGCTgg                  >  1:76441/1‑149 (MQ=255)
                                                                                                                            gTAAACTGGTTATCACTTCCCATAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCGGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGTAACCGTTGCAAACTgg                  <  1:879661/149‑1 (MQ=255)
                                                                                                                                  tGGTTATCACTTCCCATAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTCCCTTACGGTGCGGTACTGGCGAAAGGCGATGGCCAACAGGTTGCTGGCGGCGAAACCGCTGCAAACTGGGACCCg            >  2:613295/1‑149 (MQ=255)
                                                                                                                                    gTTATCACTTCCCATAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCGGTACTGGCGAAAGGCGATGGCGAACAGGTTCCTCGCGGCGAAACCGTTGCAAACTGGGACCCGcc          >  2:436655/1‑148 (MQ=255)
                                                                                                                                     ttATCACTCCCCATAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACAAACAACCTAACGGGACAGTACCGGCGAAAGACGATGGCGAACAGGTTACTCGCGCCGAGACCATGCCCAACGGGGACCCgacc         >  2:418011/1‑146 (MQ=255)
                                                                                                                                           cTTCCCATAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTCCGCTGCGGTACTGGCGAAAGCCGATGGTGAACAGGTTGCTGGCGGCGAAACCGTTGGCAACTGGCACCCGCccaccacg   >  2:983452/1‑147 (MQ=255)
                                                                                                                                           cTTCCCATAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCGGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACACCATg   >  1:857834/1‑149 (MQ=255)
                                                                                                                                            ttCCCATAATACTGAACTGAAACTGATCGATGAATTCGGTCTTACTAAAGAAAGCTACAAAGTACCTTACGGTGCGGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACACCATGc  <  2:234616/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
CTGACCATGCGTACGTTCCACATCGGTGGTGCGGCATCTCGTGCGGCTGCTGAATCCAGCATCCAAGTGAAAAACAAAGGTAGCATCAAGCTCAGCAACGTGAAGTCGGTTGTGAACTCCAGCGGTAAACTGGTTATCACTTCCCGTAATACTGAACTGAAACTGATCGACGAATTCGGTCGTACTAAAGAAAGCTACAAAGTACCTTACGGTGCGGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACACCATGC  >  NC_000913_3_bme_pgi/4188137‑4188425

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: ‑