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 1979278 776 36 [0] [1] 37 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)31 (0.730) 24/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

ATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCT  >  NC_000913/1978359‑1978645
                                                                                                                                               |                                                                                                                                               
aTGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataa                                                                                                                                            <  1:277652‑M1/149‑6 (MQ=255)
aTGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataa                                                                                                                                            <  1:560928‑M1/149‑6 (MQ=255)
  gggTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCGGCGATTTTCAATAATGCGTGATGCAAATCACCCAAAACACTCAATGACCTAACATAAATGtataagt                                                                                                                                          >  2:191377‑M1/1‑142 (MQ=255)
    gTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtca                                                                                                                                        >  1:200637‑M1/1‑140 (MQ=255)
    gTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGAGGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtca                                                                                                                                        >  1:64018‑M1/1‑140 (MQ=255)
                  cTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGATGGGGGCTGCGATTTTCAATAATGCGTGATGCACATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                                          <  2:532928‑M1/149‑24 (MQ=255)
                  cTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                                          <  2:63596‑M1/149‑24 (MQ=255)
                  cTAAGATTTGTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                                          <  2:796701‑M1/149‑24 (MQ=255)
                        tttctCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATTCAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAATCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtat                                                                                                                    <  2:322771‑M1/149‑30 (MQ=255)
                         ttttCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatt                                                                                                                   >  1:301970‑M1/1‑119 (MQ=255)
                                        gCAACTGTACGCGCCCCTACACGCACATACAACGGACGGGGGCTGCGATTTTCCATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaa                                                                                                    <  1:32874‑M1/149‑46 (MQ=255)
                                         cAACTGTACTCGTCACTACACGCACATACAACGGAGGGGTGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaag                                                                                                   <  1:448560‑M1/149‑47 (MQ=255)
                                          aaCTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagt                                                                                                  >  1:140835‑M1/1‑102 (MQ=255)
                                            cTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCCATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttataaagttt                                                                                                >  2:368814‑M1/1‑100 (MQ=255)
                                             tGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtaattcaatctgttaaaaagtttt                                                                                               >  2:752254‑M1/1‑99 (MQ=255)
                                                        cTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatactttttttttgggtgtatttcaatctgttaaaaaatttttcgctacgata                                                                                    >  1:1957‑M1/1‑88 (MQ=255)
                                                        cTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgcta                                                                                    >  1:797179‑M1/1‑88 (MQ=255)
                                                         tACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctag                                                                                   >  1:160074‑M1/1‑87 (MQ=255)
                                                            acGCACATACAACGGACGGGGGCTGCGATTTTCAATAATGCGTAACGCAGATCACACAAAACTCTCTATTACTTAACATACATGtataagacatacttttgttttgggtgtatttcattctgttaaaaagtttttcgctacgctagcaa                                                                                <  2:277543‑M1/149‑66 (MQ=255)
                                                                      aaCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACAATCAATTACTTAACATAAATGtataagtcatacttttgttttgagtgtatttcaatctgttaaaaagtttttcgctacgctaccaagcaaaaatga                                                                      >  2:442078‑M1/1‑74 (MQ=255)
                                                                             gggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacagga                                                               <  1:339504‑M1/149‑83 (MQ=255)
                                                                             gggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacagga                                                               <  1:461992‑M1/149‑83 (MQ=255)
                                                                             gggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacagga                                                               <  2:200637‑M1/149‑83 (MQ=255)
                                                                             gggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAAGTACTTAACATAGATGtataagtcattcttttgttttgggtgtatttcaatctgttaaaaagtttttcgtttcgatagaaagaaaaaatgaaacagga                                                               <  2:252830‑M1/149‑83 (MQ=255)
                                                                                  cTGCGATTTTCAATAATGCGTGATGCAGATCAGACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtcttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacgggaataat                                                          <  2:446589‑M1/149‑88 (MQ=255)
                                                                                   tGCTATTTTCAATAATGCGTAATGCATATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatc                                                         <  2:201450‑M1/149‑89 (MQ=255)
                                                                                   tGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTCCTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatc                                                         >  2:363405‑M1/1‑61 (MQ=255)
                                                                                               taatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGaataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaacatgaaacaggaataatcgaaatgggatgt                                             >  2:522460‑M1/1‑49 (MQ=255)
                                                                                                 atGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctggcaagcaaaaatgaaacaggaataatcgaaatgggatgttg                                           >  2:210909‑M1/1‑47 (MQ=255)
                                                                                                     gTGATGCAAATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctatcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgca                                       <  1:394960‑M1/149‑107 (MQ=255)
                                                                                                          gCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtc                                  <  1:34804‑M1/149‑112 (MQ=255)
                                                                                                          gCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttaggtgtctttcaatatgttaaaaaggtttttggttcgctcgcaagcaataatgaaacaggaataatagaaatgggatgttgcgcccagtg                                  >  2:685358‑M1/1‑38 (MQ=255)
                                                                                                            aGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaa                                <  2:301970‑M1/149‑114 (MQ=255)
                                                                                                                                    ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatc        >  2:537712‑M1/1‑12 (MQ=255)
                                                                                                                                     tAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatct       <  2:709062‑M1/149‑139 (MQ=255)
                                                                                                                                          tAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctat  <  2:160074‑M1/149‑144 (MQ=255)
                                                                                                                                               |                                                                                                                                               
ATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCT  >  NC_000913/1978359‑1978645

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.