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 38 [0] [0] 39 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)38 (0.750) 37/296 0.6 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

GCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGACGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACGAACAGTGGTGTGACCTGCTGGAAGAGAACTCTG  >  NC_000913/4187684‑4187978
                                                                                                                                                   |                                                                                                                                                   
atctccacccacggtgctcgtaaaggtctggcggataccgcactgaaaactgcgaactccggttacctgactcgtcgtctggttgacgtggcgcaggacctggtggttaccgaagacgattgtggtacccatgaaggtatcatgatgggtg                                                                                                                                                  >  1:453635‑M2/148‑151 (MQ=255)
    ccacccacggtgctcgtaaaggtctggcggataccgcactgaaaactgcgaactccggttacctgactcgtcgtctggttgacgtggcgcaggacctggtggttaccgaagacgattgtggtacccatgaaggtatcatgatgggtgg                                                                                                                                                 <  2:160457‑M2/5‑1 (MQ=255)
    ccacccacggtgctcgtaaaggtctggcggataccgcactgaaaactgcgaactccggttacctgactcgtcgtctggttgacgtggcgcaggacctggtggttaccgaagacgattgtggtacccatgaaggtatcatgatgGgtggtg                                                                                                                                               <  2:101266‑M2/7‑1 (MQ=255)
               gctcgtaaaggtctggcggataccgcactgaaaactgcgaactccggttacctgactcgtcgtctggttgacgtggcgcaggacctggtggttaccgaagacgattgtggtacccatgaaggtatcatgatgGGTGGTGACGTTAAAGAGc                                                                                                                                   >  2:140988‑M2/133‑151 (MQ=255)
                        ggtctggcggataccgcactgaaaactgcgaactccggttacctgactcgtcgtctggttgacgtggcgcaggacctggtggttaccgaagacgattgtggtacccatgaaggtatcatgatgGGTGGTGACGTTAAAGAGCCGCTGcgcg                                                                                                                          >  1:67859‑M2/124‑151 (MQ=255)
                           ctggcggataccgcactgaaaactgcgaactccggttacctgactcgtcgtctggttgacgtggcgcaggacctggtggttaccgaagacgattgtggtacccatgaaggtatcatgatgGGTGGTGACGTTAAAGAGCCGCTGCGCGATc                                                                                                                       >  1:197447‑M2/121‑151 (MQ=255)
                                ggataccgcactgaaaactgcgaactccggttacctgactcgtcgtctggttgacgtggcgcaggacctggtggttaccgaagacgattgtggtacccatgaaggtatcatgatgGGTGGTGACGTTAAAGAGCCGCTGCGCGATCGCGTa                                                                                                                  <  2:50982‑M2/36‑1 (MQ=255)
                                  ataccgcactgaaaactgcgaactccggttacctgactcgtcgtctggttgacgtggcgcaggacctggtggttaccgaagacgattgtggtacccatgaaggtatcatgatgGGTGGTGACGTTAAAGAGCCGCTGCGCGATCGCGTACt                                                                                                                <  1:682243‑M2/38‑1 (MQ=255)
                                    accgcactgaaaactgcgaactccggttacctgactcgtcgtctggttgacgtggcgcaggacctggtggttaccgaagacgattgtggtacccatgaaggtatcatgatgGGTGGTGACGTTAAAGAGCCGCTGCGCGATCGCGTACTgg                                                                                                              <  1:180818‑M2/40‑1 (MQ=255)
                                            gaaaactgcgaactccggttacctgactcgtcgtctggttgacgtggcgcaggacctggtggttaccgaagacgattgtggtacccatgaaggtatcatgatgGGTGGTGACGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCgtgt                                                                                                       >  2:734642‑M2/104‑150 (MQ=255)
                                            gaaaactgcgaactccggttacctgactcgtcgtctggttgacgtggcgcaggacctggtggttaccgaagacgattgtggtacccatgaaggtatcatgatgGGTGGTGACGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCgtgt                                                                                                       <  1:27911‑M2/47‑1 (MQ=255)
                                                  tgcgaactccggttacctgactcgtcgtctggttgacgtggcgcaggacctggtggttaccgaagacgattgtggtacccatgaaggtatcatgatgGGTGGTGACGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAActgct                                                                                                >  1:505648‑M2/98‑151 (MQ=255)
                                                   gcgaactccggttacctgactcgtcgtctggttgacgtggcgcaggacctggtggttaccgaagacgattgtggtacccatgaaggtatcatgatgGGTGGTGACGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAActgct                                                                                                <  2:259162‑M2/54‑1 (MQ=255)
                                                     gaactccggttacctgactcgtcgtctggttgacgtggcgcaggacctggtggttaccgaagacgattttggtacccatgaaggtatcatgatgGGTGGTGACGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCCGaa                                                                                             <  1:774306‑M2/57‑1 (MQ=255)
                                                           cggttacctgactcgtcgtctggttgacgtggcgcaggacctggtggttaccgaagacgattgtggtacccatgaaggtatcatgatgGGTGGTGACGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGtt                                                                                       >  2:106569‑M2/89‑151 (MQ=255)
                                                             gttacctgactcgtcgtctggttgacgtggcgcaggacctggtggttaccgaagacgattgtggtacccatgaaggtatcatgatgGGTGGTGACGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCt                                                                                     <  1:86798‑M2/65‑1 (MQ=255)
                                                                         gtcgtctggttgacgtggcgcaggacctggtggttaccgaagacgattgtggtacccatgaaggtatcatgatgGGTGGTGACGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTAc                                                                         >  2:740405‑M2/75‑151 (MQ=255)
                                                                           cgtctggttgacgtggcgcaggacctggtggttaccgaagacgattgtggtacccatgaaggtatcatgatgGGTGGTGACGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTAct                                                                        <  2:634747‑M2/78‑1 (MQ=255)
                                                                                 gttgacgtggcgcaggacctggtggttaccgaagacgattgtggtacccatgaaggtatcatgatgGGTGGTGACGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGat                                                                  >  2:717823‑M2/67‑150 (MQ=255)
                                                                                      cgtggcgcaggacctggtggttaccgaagacgattgtggtacccatgaaggtatcatgatgGGTGGTGACGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTc                                                            <  1:37236‑M2/90‑1 (MQ=255)
                                                                                         ggcgcaggacctggtggttaccgaagacgattgtggtacccatgaaggtatcatgatgGGTGGTGACGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGt                                                          <  2:505890‑M2/92‑1 (MQ=255)
                                                                                             caggacctggtggttaccgacgacgattgtggtacccatgaaggtatcatgatgGGTGGTGACGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGtt                                                         >  1:399877‑M2/55‑147 (MQ=255)
                                                                                             caggacctggtggttaccgaagacgattgtggtacccatgaaggcatcatgatgGGTGGCGACGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCcgc                                                     >  1:26942‑M2/55‑151 (MQ=255)
                                                                                              aggacctggtggttaccgaagacgattgtggtacccatgaaggcatcatgatgGGTGGTGACGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCcgcg                                                    >  1:63847‑M2/54‑151 (MQ=255)
                                                                                                  cctggtggttaccgaagacgattgtggtacccatgaaggtatcatgatgGGTGGTGACGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAac                                                >  1:350821‑M2/50‑151 (MQ=255)
                                                                                                              cgaagacgattgtggtacccatgaaggtatcatgatgGGTGGTGACGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCAc                                    <  2:713989‑M2/114‑1 (MQ=255)
                                                                                                                         gtggtacccatgaaggtatcatgatgGGTGGTGACGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACGAACAgtggtg                         <  2:428515‑M2/125‑1 (MQ=255)
                                                                                                                           ggtacccatgaaggtatcatgatgGGTGGTGACGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACGAACAGTGGtgtg                       >  2:771245‑M2/25‑151 (MQ=255)
                                                                                                                             tacccatgaaggtatcatgatgGGTGGTGACGTTAAAGAGCCGCCGCGCGATCGCGTAATGGGTCGAGTAACTGCTGAAGACGTTCTGAAGCAGGGTACAGATGATATACTCGTTCCGCGCAACACGCTGTTGCACGAACAGAGGggtg                       >  1:748612‑M2/23‑149 (MQ=255)
                                                                                                                                  atgaaggtatcatgatgGGTGGTGACGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACGAACAGTGGTGTGACctgct                <  2:61494‑M2/134‑1 (MQ=255)
                                                                                                                                  atgaaggtatcatgatgGGTGGTGACGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACGAACAGTGGTGTGACctgct                >  2:203146‑M2/18‑151 (MQ=255)
                                                                                                                                   tgaaggtatcatgatgGGTGGTGACGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACGAACAGTGGTGTGACctgctg               >  1:562817‑M2/17‑151 (MQ=255)
                                                                                                                                       ggtatcatgatgGGTGGTGACGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACGAACAGTGGTGTGACCTGCTGGa             >  2:392738‑M2/13‑149 (MQ=255)
                                                                                                                                       ggtatcatgatgGGTGGTGACGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACGAACAGTGGTGTGACCTGCTGGAag           <  2:552942‑M2/139‑1 (MQ=255)
                                                                                                                                          atcatgatgGGTGGTGACGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACGAACAGTGGTGTGACCTGCTGGAagag         <  2:449451‑M2/141‑1 (MQ=255)
                                                                                                                                            catgatgGGTGGTGACGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACGAACAGTGGTGTGACCTGCTGGAAGAGAAc      <  2:581812‑M2/144‑1 (MQ=255)
                                                                                                                                             atgatgGGTGGTGACGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACGAACAGTGGTGTGACCTGCTGGAAGAGAAct     >  1:310659‑M2/7‑151 (MQ=255)
                                                                                                                                                atgGGTGGTGACGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACGAACAGTGGTGTGACCTGCTGGAAGAGAACTCTg  >  2:427265‑M2/4‑151 (MQ=255)
                                                                                                                                                   ggtggtGAAGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACGAACAGTGGTGTGACCTGCTGGAAGAGAActct   <  2:217302/147‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGACGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACGAACAGTGGTGTGACCTGCTGGAAGAGAACTCTG  >  NC_000913/4187684‑4187978

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