Predicted mutation
evidence seq id position mutation freq annotation gene description
MC JC NC_000913 1,978,503 Δ776 bp 100% insB5insA5 insB5, insA5

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 1978503–1979270 1979278 9–776 20 [0] [1] 21 insB5–insA5 insB5,insA5

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 19785020 (0.000)19 (1.430) 17/384 NT 100% intergenic (‑305/+16) flhD/insB5 DNA‑binding transcriptional dual regulator FlhD/IS1 family protein InsB
?NC_000913 1979279 = 0 (0.000)intergenic (‑64/‑474) insA5/uspC IS1 family protein InsA/universal stress protein C

ATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCATACA  >  NC_000913/1978253‑1978708
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                              
tgggccccgggaagCCCCGGTAAAAAATAATTAGCATTAGAAGAATTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTGAATATGTCCTTACAAAGAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAATTCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGTGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACATAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggt                                                                                                                                                                                       <  2:168248‑M1/269‑26 (MQ=255)
                                                       cTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtt                                                                                                                                                                                             >  2:141124‑M1/1‑195 (MQ=255)
                                                           aaTAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                                                                                                         <  1:141124‑M1/214‑24 (MQ=255)
                                                                                                aCAAATAGAAATGGGTCTTTACACTTATCTCAGATTTTTCCTAAATCGACGCAACTGTACTAGTCCCTACACGCAAATACACCGGCGGGGGGCTGCGATCTTCAATAATGCGTGCTGCAGATCACACAAAACACCCAATTACTTACCCTACATGtataagtcatacttttgttttgggtgtatttcaatctgttaacaagtttttcgc                                                                                                                                                          >  1:326‑M1/1‑154 (MQ=255)
                                                                                                                                                                cACTACACGCACATACCACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtttg                                                                                                                                                                                           >  1:130100‑M1/1‑90 (MQ=255)
                                                                                                                                                                        gCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataattcatactttttttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagaaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcacttt                                                    >  1:31378‑M1/1‑82 (MQ=255)
                                                                                                                                                                        gCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacagcaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaaaaatcatctgctataaataatcacttt                                                    >  1:19330‑M1/1‑82 (MQ=255)
                                                                                                                                                                          acaTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctcgcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaata                                                                                           >  1:46258‑M1/1‑80 (MQ=255)
                                                                                                                                                                             tACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcacttt                                                    >  1:98730‑M1/1‑77 (MQ=255)
                                                                                                                                                                                 aCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcactttcatg                                                <  2:98730‑M1/233‑161 (MQ=255)
                                                                                                                                                                                     ggggggtttttatttttaaatattgcgtaaTACAGATCACAAAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaaggttttcgctacgctagcaagcaaaaatgacacaggaataatcgaaatgggatgttgcgaacagtcaaaagaactcaccgtaaataatcatctgctataaacaatcactttcatgcacgaccagctaagctatttttncacagaaacttaccgcacaanaa  <  2:261880‑M1/262‑207 (MQ=255)
                                                                                                                                                                                      ggggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgg                                                                                                                 >  1:224066‑M1/1‑68 (MQ=255)
                                                                                                                                                                                          ggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggctgttgcgcacagtcaaaa                                                                                             >  1:49176‑M1/1‑64 (MQ=255)
                                                                                                                                                                                          ggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgagcacagtcaaaataactcaccttaaataatcatctgctataaataatcactttcatgcaataccagataagctatttttaaacagacacttaccgcaca      >  1:261880‑M1/1‑64 (MQ=255)
                                                                                                                                                                                          ggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatctcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatg                                                                                                             <  2:224066‑M1/163‑100 (MQ=255)
                                                                                                                                                                                              gCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataac                                                                                         <  2:49176‑M1/179‑120 (MQ=255)
                                                                                                                                                                                                                             acaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcactttcatgcaataccagataagctatttt                           >  1:201827‑M1/1‑29 (MQ=255)
                                                                                                                                                                                                                                 aaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaacaactcaccgtaaataatcatctgctataaataatcactttcatgcaataccagataagctatttttaaa                       <  2:201827‑M1/210‑186 (MQ=255)
                                                                                                                                                                                                                                    acTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcacagtaaataatcatctgctataaataatcactttcatgcactaccagataagctatttttaaacagacacttaccg          >  1:75554‑M1/1‑22 (MQ=255)
                                                                                                                                                                                                                                                       aTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcacttt                                                    >  1:169316‑M1/1‑3 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                              
ATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCATACA  >  NC_000913/1978253‑1978708

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