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 32 [0] [0] 32 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)28 (0.680) 20/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

ATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAA  >  NC_000913/1978359‑1978648
                                                                                                                                               |                                                                                                                                                  
aTGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataa                                                                                                                                               <  2:454808‑M1/149‑6 (MQ=255)
aTGGGCCTTTACACTTATCTAACATTTTTCCTAAATTGACGCAACTTTACACGTCACTACACTCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCATTTACTTAACATAAAGGtataa                                                                                                                                               <  2:424160‑M1/149‑6 (MQ=255)
     tCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCAAACACAACACTCAATTACTTAACATAAATGtataagtcat                                                                                                                                          >  1:299680‑M1/1‑139 (MQ=255)
       tttACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatac                                                                                                                                        <  2:502721‑M1/149‑13 (MQ=255)
                    aaGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggt                                                                                                                           >  1:708243‑M1/1‑124 (MQ=255)
                              cTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaat                                                                                                                 >  2:763531‑M1/1‑114 (MQ=255)
                              cTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaat                                                                                                                 >  2:178606‑M1/1‑114 (MQ=255)
                               tAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatc                                                                                                                >  1:232955‑M1/1‑113 (MQ=255)
                                      aCGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggggtatttcaatctgttaaa                                                                                                         >  1:118767‑M1/1‑106 (MQ=255)
                                             tGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttt                                                                                                  >  1:722156‑M1/1‑99 (MQ=255)
                                                    tgcacTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCTATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctac                                                                                           <  1:108331‑M1/147‑58 (MQ=255)
                                                      cACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgc                                                                                         <  1:446072‑M1/149‑60 (MQ=255)
                                                        cTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgtgaaaaaatttttctctacgcta                                                                                       >  2:638507‑M1/1‑88 (MQ=255)
                                                         tACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctag                                                                                      >  2:97343‑M1/1‑87 (MQ=255)
                                                         tACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctaagctag                                                                                      >  2:275364‑M1/1‑87 (MQ=255)
                                                            acGCACATACAACGGAGTGGTGCTGCGATTTTCAATATTGCTTGATGCAGATCACACAAAACACTCAATTACTTAACATAATTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaactttttcgctacgctagcta                                                                                   <  2:404504‑M1/149‑66 (MQ=255)
                                                                           gggggggcTGCTATTTTCAATAATGCGTGATGCAGATCACACCACACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtgtttcaatctgttaaaaaggttttcgctccgctagcaagcaaaaatgaaacag                                                                    <  2:718014‑M1/148‑81 (MQ=255)
                                                                                          ttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcgagcaaaaatgaaacaggaataatcgaaatgg                                                     >  2:624356‑M1/1‑54 (MQ=255)
                                                                                          ttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgg                                                     <  1:5840‑M1/149‑96 (MQ=255)
                                                                                                        aTGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacag                                       <  2:708243‑M1/149‑110 (MQ=255)
                                                                                                              aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcggacagtcaaaa                                 >  1:629482‑M1/1‑34 (MQ=255)
                                                                                                              aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaa                                 >  1:466440‑M1/1‑34 (MQ=255)
                                                                                                                     aaaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctccgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataaccca                          >  1:367608‑M1/1‑27 (MQ=255)
                                                                                                                     aaaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtggatttcaatctgttaaaaagtttttcttttcttttgcaagcaaaaatgaaacaggaataatcgaaatgggatgttgggcccagccaaaataactcc                          >  2:419178‑M1/1‑27 (MQ=255)
                                                                                                                      aaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttttcttttttttttagcaaaaaagaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcac                         >  2:458849‑M1/1‑26 (MQ=255)
                                                                                                                      aaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcac                         >  1:325498‑M1/1‑26 (MQ=255)
                                                                                                                              aaTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaaccaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaata                 >  2:76687‑M1/1‑18 (MQ=255)
                                                                                                                              aaTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaaccaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaata                 >  2:73523‑M1/1‑18 (MQ=255)
                                                                                                                                 tACTTAACATAAATGtataattcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaacatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatc              <  2:333636‑M1/149‑135 (MQ=255)
                                                                                                                                 tACTTAACATAAACGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatc              <  2:463584‑M1/149‑135 (MQ=255)
                                                                                                                                           aaaTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctata    >  2:336240‑M1/1‑5 (MQ=255)
                                                                                                                                             aTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaa  >  1:57332‑M1/1‑3 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
ATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAA  >  NC_000913/1978359‑1978648

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.