Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 1,978,503 Δ776 bp insB1insA insB1, insA

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 1978503–1979270 1979278 9–776 56 [0] [0] 57 insB1–insA insB1,insA

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 19785020 (0.000)53 (1.140) 44/268 0.1 100% intergenic (‑305/+16) flhD/insB1 flagellar class II regulon transcriptional activator, with FlhC/IS1 transposase B
?NC_000913 1979279 = 0 (0.000)intergenic (‑64/‑474) insA/uspC IS1 repressor TnpA/universal stress protein

GGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCG  >  NC_000913/1978361‑1978643
                                                                                                                                             |                                                                                                                                             
gggTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACCAACCACTCAATTCCTTAACATAAATGt                                                                                                                                              >  1:792233‑M1/1‑142 (MQ=255)
     tttACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATCAATGtataag                                                                                                                                         >  1:770143‑M1/1‑137 (MQ=255)
     tttACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataag                                                                                                                                         >  2:631420‑M1/1‑137 (MQ=255)
     tttACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataag                                                                                                                                         >  1:180635‑M1/1‑137 (MQ=255)
          acTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACCCTCACATACAACGGAGGGGGGCTGCGATTTTCAAAACTGCGCCATACAGATCGCACAAGACACTCAATTACTTAACATAAATGcataactccac                                                                                                                                    >  2:592349‑M1/1‑132 (MQ=255)
           cTTATCTAAGAATTTTCCTAAAGCGACGCAACTGTACGCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatac                                                                                                                                   <  1:670223‑M1/143‑13 (MQ=255)
               tCTAAGTTTCTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGTGCTGCGATTTTCAATAATTCGTGATGCAGATCACACAAAACTCTCAATTACTTAACATAAATGtataagtcatactttt                                                                                                                               <  2:479111‑M1/143‑17 (MQ=255)
                 tAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttga                                                                                                                             >  2:37508‑M1/1‑125 (MQ=255)
                  aagtcttttcgcaaATCGACGCAACGGGAGTCGTCACCACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACGCTCAATTACTTAACATAAATGtataagtcatacttttgtt                                                                                                                            <  2:476286‑M1/138‑20 (MQ=255)
                    gATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtttt                                                                                                                          <  1:52947‑M1/143‑22 (MQ=255)
                     aTTTTGTCTAAGTCGACGCAACTGTACTCGTCACTACATGCACATACAACTGAGGGTGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttg                                                                                                                         <  2:247399‑M1/143‑23 (MQ=255)
                      tttttCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                                        <  2:437300‑M1/143‑24 (MQ=255)
                        tttCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAACTGtataattcatacttttgttttggga                                                                                                                      >  1:752136‑M1/1‑118 (MQ=255)
                            cTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtat                                                                                                                  >  1:830833‑M1/1‑114 (MQ=255)
                             tAAATCGACGCAACTGTACTCGTCACTACCCGCACATACGACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACGCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatt                                                                                                                 >  1:38738‑M1/1‑113 (MQ=255)
                             tAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACAAAAGGATGAGAAGTACTTAACATAAATGtataagtcatacttttgttttgggtgtatt                                                                                                                 >  1:659607‑M1/1‑113 (MQ=255)
                              aaaTCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtattt                                                                                                                >  1:340552‑M1/1‑112 (MQ=255)
                              aaaTCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtattt                                                                                                                <  2:827375‑M1/143‑32 (MQ=255)
                                aTCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttca                                                                                                              >  1:890834‑M1/1‑110 (MQ=255)
                                  cGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaat                                                                                                            <  1:156527‑M1/143‑36 (MQ=255)
                                   gACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatc                                                                                                           >  1:624808‑M1/1‑107 (MQ=255)
                                    aCGCAACTGTACTCGTCACTACACGCACATACCACGGCGCGACCCTGCGATTTTCAATAATACGTGATGCAGATCACACAAAACACTCAATTACTTAACATAACTGtataagtcctacttctgctttcggtctatttcaatct                                                                                                          >  2:404781‑M1/1‑106 (MQ=255)
                                       cAACTGTACCCGGCACTACACGCACACACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgtt                                                                                                       <  1:93863‑M1/143‑41 (MQ=255)
                                        aaCTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACCCTCAATTCCTTAACATAAATGtataagtcatacttttgtttttgttctatttccatctcttc                                                                                                      >  1:454541‑M1/1‑102 (MQ=255)
                                           tGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaa                                                                                                   >  1:304315‑M1/1‑99 (MQ=255)
                                            gTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaa                                                                                                  <  2:1510‑M1/143‑46 (MQ=255)
                                                tCGTCACGACACGAACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttt                                                                                              <  1:86543‑M1/143‑50 (MQ=255)
                                                       tccacgcACATAGCACTGAGGGTGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgcta                                                                                       <  2:771060‑M1/141‑57 (MQ=255)
                                                       tACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTCCTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctc                                                                                       >  2:269955‑M1/1‑87 (MQ=255)
                                                       tACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgcta                                                                                       >  1:65340‑M1/1‑87 (MQ=255)
                                                          acGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgc                                                                                    <  2:431983‑M1/143‑60 (MQ=255)
                                                              acaTACAGCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACGCTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagc                                                                                <  2:289857‑M1/143‑64 (MQ=255)
                                                                  aCAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtaatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagc                                                                            <  2:850917‑M1/143‑68 (MQ=255)
                                                                    aaCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaa                                                                          >  1:526829‑M1/1‑74 (MQ=255)
                                                                        gAGGGGGGCTGTGATTTTCAATAATGCGTGATGCAGATCACACAAACCACTCCATTACTTATCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaat                                                                      <  2:121088‑M1/143‑74 (MQ=255)
                                                                        gAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACCCTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaat                                                                      <  1:512862‑M1/143‑74 (MQ=255)
                                                                                   cGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaat                                                           >  2:130510‑M1/1‑59 (MQ=255)
                                                                                   cGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaat                                                           >  1:158780‑M1/1‑59 (MQ=255)
                                                                                     aTTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataa                                                         >  1:274349‑M1/1‑57 (MQ=255)
                                                                                      ttttCAATAATGCGTGATGCAGATCACACTAATCACTCCAGTACCTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataat                                                        <  1:165425‑M1/143‑88 (MQ=255)
                                                                                         tCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcga                                                     <  1:695457‑M1/143‑91 (MQ=255)
                                                                                          cAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaa                                                    >  1:346399‑M1/1‑52 (MQ=255)
                                                                                                ttggtgaTGGAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggga                                              <  1:727726‑M1/140‑98 (MQ=255)
                                                                                                tGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgtgtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaacaatgaaacaggaataatcgaaatggga                                              >  2:439666‑M1/1‑46 (MQ=255)
                                                                                                            aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacag                                  >  1:152866‑M1/1‑34 (MQ=255)
                                                                                                             ccacacaaAACACTCAATTACTTAACATAAACGtataagtcatacttttgtttggggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagt                                 <  1:244337‑M1/142‑111 (MQ=255)
                                                                                                                     aaCACTCAATTACTTAGCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataa                         >  2:744372‑M1/1‑25 (MQ=255)
                                                                                                                            aaTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgagcacagtcaaaataactcaccg                  >  2:3935‑M1/1‑18 (MQ=255)
                                                                                                                                  ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaata            >  1:437141‑M1/1‑12 (MQ=255)
                                                                                                                                  ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatattgcgcacagtcaaaataagtcaccgtcaata            >  1:461129‑M1/1‑12 (MQ=255)
                                                                                                                                  ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataaccaaaatgggatcgcgcacacagtcaaaataactcaccgtaaata            >  2:434275‑M1/1‑12 (MQ=255)
                                                                                                                                    aaCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataat          >  1:15189‑M1/1‑10 (MQ=255)
                                                                                                                                     aCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatc         >  1:185482‑M1/1‑9 (MQ=255)
                                                                                                                                       aTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgtgatgttgcgcacagtcaaaataactcaccgtccataatcat       >  2:152189‑M1/1‑7 (MQ=255)
                                                                                                                                          aaTGtataagtcatccttgtgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggattttgcgcacagtcaaaataactcaccgtaaataatcatctg    <  1:305980‑M1/143‑140 (MQ=255)
                                                                                                                                            tGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgct  <  1:137680‑M1/143‑142 (MQ=255)
                                                                                                                                             |                                                                                                                                             
GGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCG  >  NC_000913/1978361‑1978643

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.