Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 4,187,816 Δ15 bp coding (2467‑2481/4224 nt) rpoC → RNA polymerase, beta prime subunit

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 4187816 4187830 15 35 [0] [0] 34 rpoC RNA polymerase, beta prime subunit

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 41878150 (0.000)33 (0.640) 31/296 1.1 100% coding (2466/4224 nt) rpoC RNA polymerase, beta prime subunit
?NC_000913 4187831 = 0 (0.000)coding (2482/4224 nt) rpoC RNA polymerase, beta prime subunit

CTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGACGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACGAACAGTGGTGTGACCTGCTG  >  NC_000913/4187665‑4187965
                                                                                                                                                      |                                                                                                                                                      
cTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATG                                                                                                                                                        >  2:810307‑M1/1‑151 (MQ=255)
  tCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGg                                                                                                                                                       <  2:566423‑M1/150‑2 (MQ=255)
   cATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGggt                                                                                                                                                     >  1:793712‑M1/1‑148 (MQ=255)
             cacGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGggtggtgacgtta                                                                                                                                           <  1:77363‑M1/151‑14 (MQ=255)
                   gCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGggtggtgacgttaaagag                                                                                                                                      <  2:808723‑M1/150‑19 (MQ=255)
                                tGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGggtggtgacgttaaagagccgctgcgcgatcg                                                                                                                        >  1:635188‑M1/1‑119 (MQ=255)
                                 ggCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGggtggtgacgttaaagagccgctgcgcgat                                                                                                                          <  1:792560‑M1/148‑31 (MQ=255)
                                 ggCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGggtggtgacgttaaagagccgctgcgcgat                                                                                                                          <  1:609635‑M1/148‑31 (MQ=255)
                                       tACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGggtggtgacgttaaagagccgctgcgcgatcgcgtactg                                                                                                                 >  1:756762‑M1/1‑112 (MQ=255)
                                        aCCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGggtggtgacgttaaagagccgctgcgcgatcgcgtactgg                                                                                                                <  1:221837‑M1/151‑41 (MQ=255)
                                          cGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGggtggtgacgttaaagagccgctgcgcgatcgcgtactgggt                                                                                                              <  1:690259‑M1/151‑43 (MQ=255)
                                            cACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGggtggtgacgttaaagagcggctgcgcgatcgcgtactgggtcg                                                                                                            >  1:184491‑M1/1‑107 (MQ=255)
                                            cACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGggtggtgacgttaaagagccgctgcgcgatcgcgtactgggtcg                                                                                                            <  1:202825‑M1/151‑45 (MQ=255)
                                                gAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGggtggtgacgttaaagagccgctgcgcgatcgcgtactgggtcgtgt                                                                                                         >  1:194687‑M1/1‑103 (MQ=255)
                                                  aaaCTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGggtggtgacgttaaagagccgctgcgcgatcgcgtactgggtcgtgtaac                                                                                                      <  1:611472‑M1/151‑51 (MQ=255)
                                                               cGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGggtggtgacgttaaagagccgctgcgcgatcgcgtactgggtcgtgtaactgctgaagacgtt                                                                                         >  1:99359‑M1/1‑88 (MQ=255)
                                                                   tACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGggtggtgacgttaaagagccgctgcgcgatcgcgtactgggtcgtgtaactgctgaagacgttctg                                                                                      >  1:371226‑M1/1‑84 (MQ=255)
                                                                               cgtcTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGggtggtgacgttaaagagccgctgcgcgatcgcgtactgggtcgtgtaactgctgaagacgttctgaagccgggtactg                                                                         >  2:783866‑M1/1‑72 (MQ=255)
                                                                                               cgcAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGggtggtgacgttaaagagccgctgcgcgatcgcgtactgggtcgtgtaactgctgaagacgttctgaagccgggtactgctgatatcctcgttcc                                                         <  2:39230‑M1/151‑96 (MQ=255)
                                                                                                  aGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGggtggtgacgttaaagagccgctgcgcgatcgcgtactgggtcgtgtaactgctgaagacgttctgaagccgggtactgctgatatcctcgttccgcg                                                      >  1:28574‑M1/1‑53 (MQ=255)
                                                                                                     aCCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGggtggtgacgttaaagagccgctgcgcgatcgcgtactgggtcgtgtaactgctgaagacgttctgaagccgggtactgctgatatcctcgttccgcgcaa                                                   <  2:217633‑M1/151‑102 (MQ=255)
                                                                                                       cTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGggtggtgacgttaaagagccgctgcgcgatcgcgtactgggtcgtgtaactgctgaagacgttctgaagccgggtactgctgatatcctcgttccgcgcaac                                                  >  2:356394‑M1/1‑48 (MQ=255)
                                                                                                            ggtTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGggtggtgacgttaaagagccgctgcgcgatcgcgtactgggtcgtgtaactgctgaagacgttctgaagccgggtactgctgatatcctcgttccgcgcaacacgctg                                            >  1:61405‑M1/1‑43 (MQ=255)
                                                                                                                   gAAGACGATTGTGGTACCCATGAAGGTATCATGATGggtggtgacgttaaagagccgctgcgcgatcgcgtactgggtcgtgtaactgctgaagacgttctgaagccgggtactgctgatatcctcgttccgcgcaacacgctgctgcacg                                     >  1:49608‑M1/1‑36 (MQ=255)
                                                                                                                     aGACGATTGTGGTACCCATGAAGGTATCATGATGggtggtgacgttaaagagccgctgcgcgatcgcgtactgggtcgtgtaactgctgaagacgttctgaagccgggtactgctgatatcctcgttccgcgcaacacgctgctgcacgaa                                   <  1:748971‑M1/151‑118 (MQ=255)
                                                                                                                       aCGATTGTGGTACCCATGAAGGTATCATGATGggtggtgacgttaaagagccgctgcgcgatcgcgtactgggtcgtgtaactgctgaagacgttctgaagccgggtactgctgatatcctcgttccgcgcaacacgctgctgcacgaaca                                 <  1:779529‑M1/151‑120 (MQ=255)
                                                                                                                             gtgGTACCCATGAAGGTATCATGATGggtggtgacgttaaagagccgctgcgcgatcgcgtactgggtcgtgtaactgctgaagacgttctgaagccgggtactgctgatatcctcgttccgcgcaacacgctgctgcacgaacagtggtg                           <  1:813280‑M1/151‑126 (MQ=255)
                                                                                                                                  aCCCATGAAGGTATCATGATGggtggtgacgttaaagagccgctgcgcgatcgcgtactgggtcgtgtaactgctgaagacgttctgaagccgggtactgctgatatcctcgttccgcgcaacacgctgctgcacgaacagtggtgtgacc                      >  2:81426‑M1/1‑21 (MQ=255)
                                                                                                                                      aTGAAGGTATCATGATGggtggtgacgttaaagagccgctgcgcgatcgcgtactgggtcgtgtaactgctgaagacgttctgaagccgggtactgctgatatcctcgttccgcgcaacacgctgctgcacgaacagtggtgtgacctgct                  >  2:489644‑M1/1‑17 (MQ=255)
                                                                                                                                       tGAAGGTATCATGATGggtggtgacgttaaagagccgctgcgcgatcgcgtactgggtcgtgtaactgctgaagacgttctgaagccgggtactgctgatatcctcgttccgcgcaacacgctgctgcacgaacagtggtgtgacctgctg                 >  2:351345‑M1/1‑16 (MQ=255)
                                                                                                                                          aGGTATCATGATGggtggtgacgttaaagagccgctgcgcgatcgcgtactgggtcgtgtaactgctgaagacgttctgaagccgggtactgctgatatcctcgttccgcgcaacacgctgctgcacgaacagtggtgtgacctgctgga               >  2:394581‑M1/1‑13 (MQ=255)
                                                                                                                                                 atgatgggtggtgacgttaaagagccgctgcgcgatcgcgtactgggtcgtgtaactgctgaagacgttctgaagccgggtactgctgatatcctcgttccgcgcaacacgctgctgcacgaacagtggtgtgacctgctggaagagaact       >  1:561466‑M1/1‑6 (MQ=255)
                                                                                                                                                   gatgggtggtgacgttaaagagccgctgcgcgatcgcgtactgggtcgtgtaactgctgaagacgttctgaagccgggtactgctgatatcctcgttccgcgcaacacgctgctgcacgaacagtggtgtgacctgctggaagagaactct     <  2:251995‑M1/151‑148 (MQ=255)
                                                                                                                                                   gatgggtggtgacgttaaagagccgctgcgcgatcgcgtactgggtcgtgtaactgctgaagacgttctgaagccgggtactgctgatatcctcgttccgcgcaacacgctgctgcacgaacagtggtgtgacctgctggaagagaactct     <  1:505692‑M1/151‑148 (MQ=255)
                                                                                                                                                      gggtggtgacgttaaagagccgctgcgcgatcgcgtactgggtcgtgtaactgctgaagacgttctgaagccgggtactgctgatatcctcgttccgcgcaacacgctgctgcacgaacagtggtgtgacctgctggaagagaactctgtc  <  1:9904‑M1/151‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
CTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGACGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACGAACAGTGGTGTGACCTGCTG  >  NC_000913/4187665‑4187965

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: