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 1979278 776 60 [0] [0] 61 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)56 (1.220) 46/266 0.0 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

CTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATG  >  NC_000913/1978365‑1978640
                                                                                                                                         |                                                                                                                                          
cTTTACACTCATCTAAGATTTCTCCCAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCCTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataa                                                                                                                                       <  1:849327‑M1/143‑6 (MQ=255)
    acacTTATCTAAGATTTTTCCTAAATCGACGCCACTGTACTCGTCACTACACGCACATAAAACGGACGGGGGCTGCGATTTTCAATAATGAGCGATGCAGACCACACAAAACACTCAACTACTTAACATAAATGtatacgtcc                                                                                                                                   >  2:400389‑M1/1‑134 (MQ=255)
    acacTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACCAAATTAATGAACATAAATGtataactca                                                                                                                                   >  2:124871‑M1/1‑134 (MQ=255)
              aaGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATCAATGtataagtcatacttttgtt                                                                                                                         >  1:137202‑M1/1‑124 (MQ=255)
              aaGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtt                                                                                                                         >  2:678188‑M1/1‑124 (MQ=255)
                  tttttCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtatacgtcatacttttatcttgg                                                                                                                     >  2:40855‑M1/1‑120 (MQ=255)
                    tttCCTAACTCTACGCAACTGTACTCCTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAACAATGCGTGATGCAGATCACACAAAACTCTCAATTACTTAACATAAATGtataagtcatacttttgttttgggt                                                                                                                   <  1:264479‑M1/143‑26 (MQ=255)
                    tttCCTAAATCGACGGAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggt                                                                                                                   >  2:750568‑M1/1‑118 (MQ=255)
                        cTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGCTGCAGATCACACAACACACTCAATTACTTAACATAAATGtataagccatacttttgtttttggtgtat                                                                                                               >  2:706852‑M1/1‑114 (MQ=255)
                         tAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatt                                                                                                              <  1:698311‑M1/143‑31 (MQ=255)
                              cGACGCGACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaat                                                                                                         <  2:201025‑M1/143‑36 (MQ=255)
                                aCGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGGGTGATGCAGATCACACAAAACACTCAATTACTTTACATCAATGtataagtcatacttttgttttgggtgtatttcaatct                                                                                                       <  1:386066‑M1/143‑38 (MQ=255)
                                aCGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtctaagtcatacttttgttttgggtgtatttcaatct                                                                                                       >  2:298868‑M1/1‑106 (MQ=255)
                                  gCAACTGTGCTCGTCACTCCACGCACATCCAGCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgt                                                                                                     <  2:661427‑M1/143‑40 (MQ=255)
                                      cTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATCAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaa                                                                                                 >  2:21531‑M1/1‑100 (MQ=255)
                                       tGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcattctgttagaa                                                                                                >  2:433473‑M1/1‑99 (MQ=255)
                                             cGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttt                                                                                          <  1:754755‑M1/143‑51 (MQ=255)
                                             cGTCACTACACGCACATACAACGGAGGGAGACTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttt                                                                                          <  2:697935‑M1/143‑51 (MQ=255)
                                              gtgggttcacccacatgCTTCGGAGGATGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttttt                                                                                         <  1:641997‑M1/136‑52 (MQ=255)
                                              gTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCTGATCACACATAACACTCAAGTACTTAACATAAATGtagaagtcatacttctgttttgggtgtatttcaatctgttaaaaagttttt                                                                                         <  2:587872‑M1/143‑52 (MQ=255)
                                                cACTACACGCACATACACCGGAGGGGGGCTGCTATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcg                                                                                       <  1:551761‑M1/143‑54 (MQ=255)
                                                      acGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgc                                                                                 <  1:850288‑M1/143‑60 (MQ=255)
                                                        gCACATACAACGGAGGGGGGCTGCGATTTTCAATGCTGCGTGATGCAGATCACACAAAACACTCAAGTACTTAACATAAATGtatacgtcatacttttgttttcggtgtatttcaatctggtcaaacggtgttcgctacgctc                                                                               >  2:683705‑M1/1‑82 (MQ=255)
                                                          acaTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagc                                                                             >  2:525342‑M1/1‑80 (MQ=255)
                                                            aTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaa                                                                           <  1:649128‑M1/143‑66 (MQ=255)
                                                              aCAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagc                                                                         <  1:818903‑M1/143‑68 (MQ=255)
                                                              aCAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagc                                                                         >  1:394126‑M1/1‑76 (MQ=255)
                                                                 aCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaa                                                                      <  1:458103‑M1/143‑71 (MQ=255)
                                                                  cGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaa                                                                     >  1:63542‑M1/1‑72 (MQ=255)
                                                                  cGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaa                                                                     <  1:750568‑M1/143‑72 (MQ=255)
                                                                   ggggtggggctgcgCTTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaa                                                                    <  2:74150‑M1/138‑73 (MQ=255)
                                                                        ggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaa                                                               <  1:632453‑M1/143‑78 (MQ=255)
                                                                        ggggCAGTGTTTTTTAATAATACGTGATGCAGATCACAGAAAACACTGGATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttgtcgctacgctagcaagcaaaaatgaaa                                                               <  2:680490‑M1/143‑78 (MQ=255)
                                                                            cTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaacatgaaacagg                                                           >  2:781424‑M1/1‑62 (MQ=255)
                                                                                      cAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaa                                                 >  1:43729‑M1/1‑52 (MQ=255)
                                                                                      cAATAATACGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGgataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaa                                                 <  1:573618‑M1/143‑92 (MQ=255)
                                                                                       aataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaa                                                >  1:84045‑M1/1‑51 (MQ=255)
                                                                                       aataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttctcgctacgctagcaagccaaaatgaaacaggaataatcgaaa                                                >  1:373774‑M1/1‑51 (MQ=255)
                                                                                            tGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggga                                           <  1:771690‑M1/143‑98 (MQ=255)
                                                                                              cGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatg                                         >  1:612113‑M1/1‑44 (MQ=255)
                                                                                                tGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgtt                                       <  2:63542‑M1/143‑102 (MQ=255)
                                                                                                  aTGAAGATCACACAAAACACCCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgc                                     <  1:864461‑M1/143‑104 (MQ=255)
                                                                                                    gcggatcACACATAACCCTCATCTACTTAACCTAAATGtataagtcatacttttgttttgggtgtatttcaagctgttaaaaagtttttcgttacgctagcaagcaaaaatgaaacaggaataatcggaatgggatgttgcgc                                   <  2:749018‑M1/140‑106 (MQ=255)
                                                                                                     cAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatactcttgttttgggtgtatttaaatctgttaaaaagtttttcgctacgccagcaagaaaaaatgcaacaggaataatcaacacgggacgtcccgca                                  >  2:743172‑M1/1‑37 (MQ=255)
                                                                                                       gATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcaca                                <  1:626493‑M1/143‑109 (MQ=255)
                                                                                                        aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacag                               >  1:826456‑M1/1‑34 (MQ=255)
                                                                                                             acaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaa                          >  1:28461‑M1/1‑29 (MQ=255)
                                                                                                             acaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaa                          >  1:90121‑M1/1‑29 (MQ=255)
                                                                                                               aaaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgcgacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaat                        <  2:537616‑M1/143‑117 (MQ=255)
                                                                                                               aaaaCACTCAATTACTCACCATAAATGtataaggtatagttttgttttgagtgcagttcaatctgttaaaatgtttttagccacgcgagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaat                        <  2:312046‑M1/143‑117 (MQ=255)
                                                                                                                aaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagaaaaactgaaacaggaataatcgaaattcgatgttgcgcacagtacaaaaa                       >  2:832033‑M1/1‑26 (MQ=255)
                                                                                                                 aaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataa                      <  1:438421‑M1/143‑119 (MQ=255)
                                                                                                                        aaTAACTTAACATAAATGtataagtcatccttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatggaactgggatgttgccgaccgtcaaaataccacaccg               >  2:775781‑M1/1‑18 (MQ=255)
                                                                                                                           tACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaa            <  1:652949‑M1/143‑129 (MQ=255)
                                                                                                                              ttAACATTAATGtataagtcatacttttgtgttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaata         <  2:439240‑M1/143‑132 (MQ=255)
                                                                                                                                 aCATAAATGtataagtcatactgttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatc      <  2:453320‑M1/143‑135 (MQ=255)
                                                                                                                                   aTAAATGtataagtcatacttttgtttttggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaccaggaataatcgaaatgggatgttgcggacagtcaaaataacccaccgtaaataatcat    <  2:489087‑M1/143‑137 (MQ=255)
                                                                                                                                    tAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatc   >  1:204030‑M1/1‑6 (MQ=255)
                                                                                                                                     acatgtagaagtcagacttttgttttgggtgtatttcaatctgttaaaaagcttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatct  <  1:683705‑M1/141‑139 (MQ=255)
                                                                                                                                     aaaTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatct  >  2:91317‑M1/1‑5 (MQ=255)
                                                                                                                                         |                                                                                                                                          
CTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATG  >  NC_000913/1978365‑1978640

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.