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

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 1978503–1979269 1979278 10–776 36 [0] [0] 36 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)32 (0.660) 27/260 NT 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

ATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCG  >  NC_000913/1978359‑1978643
                                                                                                                                               |                                                                                                                                             
aTGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGTGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataa                                                                                                                                          <  2:757281‑M1/149‑6 (MQ=255)
              ttATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCTCATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACCTAAATGtataagtcatacttttgtt                                                                                                                            >  2:103996‑M1/1‑130 (MQ=255)
                aTCTAAGATTTTTCCTAAATCGACGCAACTGTAGTCGTCACTACACGCACATACAACGGAAGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTCATTAACATGAATGtataagtcatactttattttt                                                                                                                          >  2:314143‑M1/1‑128 (MQ=255)
                  cTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGCGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                                        <  2:798132‑M1/149‑24 (MQ=255)
                    aaGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggt                                                                                                                      >  1:215190‑M1/1‑124 (MQ=255)
                     aGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtg                                                                                                                     >  2:89435‑M1/1‑123 (MQ=255)
                            tCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGTGGGCGGCGTCGAGTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttca                                                                                                              <  1:623873‑M1/149‑34 (MQ=255)
                               tAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaagc                                                                                                           >  2:188544‑M1/1‑113 (MQ=255)
                                aaaTCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATCCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtttttggtgtattttaatct                                                                                                          >  2:476390‑M1/1‑112 (MQ=255)
                                 aaTCGACGCAACCGTACCCGTCACCACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctg                                                                                                         <  1:184636‑M1/149‑39 (MQ=255)
                                              gTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttttt                                                                                            >  1:289921‑M1/1‑98 (MQ=255)
                                                  tCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttggca                                                                                        >  2:637062‑M1/1‑94 (MQ=255)
                                                        cTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgcta                                                                                  >  2:523905‑M1/1‑88 (MQ=255)
                                                          acacGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagc                                                                                >  1:534145‑M1/1‑86 (MQ=255)
                                                                        cGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaa                                                                  <  1:530921‑M1/149‑78 (MQ=255)
                                                                              ggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaa                                                            >  2:315735‑M1/1‑66 (MQ=255)
                                                                              ggggCTGCGATTTTCAATAATACGTGATGCAGAGCACACAAAACACTCAATTACTTAACATAAATGtataagtgatactttagttttgggtgtatttcaatctgttgaaaagttttatgctacgctagcaagaaaaaaggaaaaaggat                                                            >  2:512013‑M1/1‑66 (MQ=255)
                                                                                ggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaata                                                          <  2:115512‑M1/149‑86 (MQ=255)
                                                                                ggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaata                                                          >  2:118187‑M1/1‑64 (MQ=255)
                                                                                     cGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcga                                                     >  1:620482‑M1/1‑59 (MQ=255)
                                                                                       aTTTTCAATAATGCCTGATGCAGATCACACAAAACACTCAATTACTTAAGATAAATGtataagtcatacttttgtattgggtgtatttcaatccgttaaaacgtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaa                                                   >  1:673297‑M1/1‑57 (MQ=255)
                                                                                       aTTTTCAATAATGCATGATGCAGATCACACAAAACACTCAATTACTTAAAATAAATGtataagtcatacttttgtgttgggtgtatttcaatctgttaaaacgtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaa                                                   >  1:673264‑M1/1‑57 (MQ=255)
                                                                                          ttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtctaagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcactacgcgagcaatcaaaaatggaacaggaatgatcgatatag                                                >  2:129697‑M1/1‑54 (MQ=255)
                                                                                             aataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtatgggacagactgctgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggat                                             <  1:664260‑M1/149‑99 (MQ=255)
                                                                                                     gTGATGCAGATCACACAAAACACTTAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgca                                     <  2:253920‑M1/149‑107 (MQ=255)
                                                                                                     gTGATGCAGATCACAAAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgca                                     <  1:597272‑M1/149‑107 (MQ=255)
                                                                                                      tGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttctgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggattttgcgcac                                    >  2:592260‑M1/1‑42 (MQ=255)
                                                                                                         tGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagt                                 >  1:767552‑M1/1‑39 (MQ=255)
                                                                                                         tGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagt                                 >  1:720264‑M1/1‑39 (MQ=255)
                                                                                                           cAGATCACACAAAACACTCCATTCCTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtca                               <  1:197187‑M1/149‑113 (MQ=255)
                                                                                                              aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaa                            <  2:363352‑M1/149‑116 (MQ=255)
                                                                                                              aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaa                            >  2:488366‑M1/1‑34 (MQ=255)
                                                                                                              aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaa                            >  1:753176‑M1/1‑34 (MQ=255)
                                                                                                                                   cTTAACATAAATGtataagtcatacttttcttttgggtgtatttcaatctgtttaaaagtttttcgctacgctagcaagcaaaaatgaatcaggaataatcgaaatgggatgttgcgcacagtcaaaataattcaccgtaaataatcat       <  2:673264‑M1/149‑137 (MQ=255)
                                                                                                                                   cTTAACATAAATGtataagtcatacttttcttttgggtgtatttcaatctgttcaaaagtttttcgctacgctagcaagcaaaaatgaagcaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcat       <  2:673297‑M1/149‑137 (MQ=255)
                                                                                                                                        aatAAATGtataggtcctacttttgttttgggtgtatttcaatgtgttaaaaagtttgtcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgct  <  2:656817‑M1/148‑142 (MQ=255)
                                                                                                                                               |                                                                                                                                             
ATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCG  >  NC_000913/1978359‑1978643

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ 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.