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 27 [0] [0] 27 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)27 (0.800) 22/208 0.4 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

ATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATT  >  NC_000913/1978375‑1978631
                                                                                                                               |                                                                                                                                 
aTCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcata                                                                                                                        >  1:661327‑M1/1‑128 (MQ=255)
                                       aCTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttt                                                                                 >  2:584066‑M1/1‑89 (MQ=255)
                                                     cAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtca                                                                                                                          <  2:511182‑M1/84‑10 (MQ=255)
                                                     cAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGAGCACACAAAACACTCAATTACTTAACATAAATGtataagtca                                                                                                                          >  1:511182‑M1/1‑75 (MQ=255)
                                                          gAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcttacttttgttttgggtgtatttcaatctgttaaaaagtttttcg                                                                              <  2:554087‑M1/123‑54 (MQ=255)
                                                          gAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaa                                                              >  2:79613‑M1/1‑70 (MQ=255)
                                                          gAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaa                                                              >  1:390275‑M1/1‑70 (MQ=255)
                                                          gAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcg                                                                              >  1:554087‑M1/1‑70 (MQ=255)
                                                          gAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttt                                                                                                                    >  1:202439‑M1/1‑70 (MQ=255)
                                                          gAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttt                                                                                                                    <  2:202439‑M1/85‑16 (MQ=255)
                                                            ggggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaa                                                            >  2:386578‑M1/1‑68 (MQ=255)
                                                             gggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggt                                                                                                          >  1:335491‑M1/1‑67 (MQ=255)
                                                             gggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggt                                                                                                          <  2:335491‑M1/92‑26 (MQ=255)
                                                             gggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttt                                                                                                               >  2:434545‑M1/1‑67 (MQ=255)
                                                             gggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttt                                                                                                               <  1:434545‑M1/87‑21 (MQ=255)
                                                             gggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtat                                                                                                                                <  2:474438‑M1/70‑4 (MQ=255)
                                                             gggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtat                                                                                                                                >  1:474438‑M1/1‑67 (MQ=255)
                                                                                         tGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgt                               >  1:611796‑M1/1‑39 (MQ=255)
                                                                                              aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataat                                            >  1:64373‑M1/1‑34 (MQ=255)
                                                                                              aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataat                                            <  2:64373‑M1/121‑88 (MQ=255)
                                                                                                   acaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgagtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgtgatgttgcgcacagt                     <  2:523220‑M1/139‑111 (MQ=255)
                                                                                                        acacTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaa                >  2:122029‑M1/1‑24 (MQ=255)
                                                                                                          acTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaata              >  1:532351‑M1/1‑22 (MQ=255)
                                                                                                           cTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataa             >  2:139068‑M1/1‑21 (MQ=255)
                                                                                                               aTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactca         <  2:611796‑M1/139‑123 (MQ=255)
                                                                                                                   cTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgt     <  1:386578‑M1/139‑127 (MQ=255)
                                                                                                                      aaCATAAATGtataagtaatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaa  >  2:448592‑M1/1‑10 (MQ=255)
                                                                                                                               |                                                                                                                                 
ATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATT  >  NC_000913/1978375‑1978631

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