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 40 [0] [3] 43 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.410) 27/408 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

CCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCATACAGCGGCCAG  >  NC_000913/1978268‑1978716
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                      
cccGGTAAAAAATAATTAGCATTAGAATAGGTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtat                                                                                                                                                                                                                     <  1:6711‑M1/238‑4 (MQ=255)
      aaaaaaTAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtca                                                                                                                                                                                                               <  2:968602‑M1/238‑10 (MQ=255)
      aaaaaaTAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtca                                                                                                                                                                                                               <  1:714365‑M1/238‑10 (MQ=255)
           ataatTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCAACTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttt                                                                                                                                                                                                         <  1:764809‑M1/239‑16 (MQ=255)
                 aGCATTAGAATAGTTGCGATAAGCTTCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtttt                                                                                                                                                                                                   <  1:491973‑M1/239‑22 (MQ=255)
                                                     aaCCACCTTTTTGGTTTAATATGCCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACTCACTTACAACGAAGGGGGGCTCCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACGCAATTACTTAACATAAATGtttaggccatacttttattttgtatttattttaatcttttaaaaaatatttgcct                                                                                                                                                                 >  2:945055‑M1/1‑182 (MQ=255)
                                                                                  aaaatAGAAATGGGTCTTTACCCTGAACTAAGATTTTTCCTAAATCGACGCAAATGTACTCGTCACTACAAGCACATACAACGGAGGGGGGGTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataa                                                                                                                                  <  1:945055‑M1/238‑87 (MQ=255)
                                                                                                     acccttttttaagatttttCCTAAATCGACCCAATTTTACTCGTCCCTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgc                                                                                                               <  2:378342‑M1/230‑106 (MQ=255)
                                                                                                     acacTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgc                                                                                                               <  1:130077‑M1/239‑106 (MQ=255)
                                                                                                      catttatcTAAGATTTTTCCCAAACCGACGTTACTGTACTCGGCACTACACGAACATACAACGGAGGGGGGCTGCCACTTGCAATAATGCGTGATGCAGATCACAAAAACCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgca                                                                                                              <  1:1106572‑M1/236‑107 (MQ=255)
                                                                                                          tATCTAAGTTTTTTCCTAAATGGACGCAAAGGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagt                                                                                                          <  1:97993‑M1/239‑111 (MQ=255)
                                                                                                           aTCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTACCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttccctacgctagcaagcaaaaatgaaacaggaataatcgaaaaggaatattggccacagca                                                                                                         >  2:97993‑M1/1‑128 (MQ=255)
                                                                                                                 gATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatact                                                                                                                                                                                                           <  1:362072‑M1/135‑14 (MQ=255)
                                                                                                                 gATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatact                                                                                                                                                                                                           >  2:362072‑M1/1‑122 (MQ=255)
                                                                                                                                 aCGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagaaagcaaaattgaaacaggaaaaatcgaaatgggatgttggcaacagccaaaaaaaataaccgaaaaaaaa                                                                                   >  2:980579‑M1/1‑106 (MQ=255)
                                                                                                                                          tACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgcta                                                                          <  1:980579‑M1/239‑143 (MQ=255)
                                                                                                                                               gTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataataatctgctataaat                                                                     >  1:484171‑M1/1‑92 (MQ=255)
                                                                                                                                                         gCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagt                                                                                                          >  1:20359‑M1/1‑82 (MQ=255)
                                                                                                                                                         gCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagt                                                                                                          <  2:20359‑M1/192‑111 (MQ=255)
                                                                                                                                                                cAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcactttcatgcaat                                                    >  1:123152‑M1/1‑75 (MQ=255)
                                                                                                                                                                   cGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataataatctcctataaataatcactttcatgcaaaacc                                                 >  2:873188‑M1/1‑72 (MQ=255)
                                                                                                                                                                       gggggtCTTCGATTTTCAATTATGGCTGATGCAATTCACACAAAACCCTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtttatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactccccgtaaataatcatctgctataaataatcactttcatgcaataccaga                                              <  2:896344‑M1/238‑171 (MQ=255)
                                                                                                                                                                           ggatgcgATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcactttcatgcaataccagataag                                          <  1:967309‑M1/235‑175 (MQ=255)
                                                                                                                                                                           ggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtaaaaataactcaccgaaaataatcatctgctataaataatcactttcatgcaaaaccagaaaaga                                         >  2:967309‑M1/1‑64 (MQ=255)
                                                                                                                                                                                         ataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcactttcatgcaataccagataagctatttttaaacaga                           <  1:873188‑M1/239‑190 (MQ=255)
                                                                                                                                                                                                              acaAAACACTCAATTACTTAACATAAATGtataagtcatccttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctttaaataatcactttcatgcaataccagataagccatttttaaaaaaacacttaccccacaacaaacaa      >  2:269166‑M1/1‑29 (MQ=255)
                                                                                                                                                                                                                   tacagcaattaCTTAACATAAATTgaaaagtcatacatgttttttgggtgtatttcaatctgttaaaaaggttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcactttcatgcaataccagataagctagttttaaacagacacttaccgcacaacaaactaagga  <  2:539698‑M1/233‑215 (MQ=255)
                                                                                                                                                                                                                       tCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcactttcatgcaataccagataagc                                         >  1:822308‑M1/1‑20 (MQ=255)
                                                                                                                                                                                                                       tCAATTACTTAACAAAAATGtaaaagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcactttcatgcaataccagataagc                                         <  2:822308‑M1/195‑176 (MQ=255)
                                                                                                                                                                                                                        cAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataact                                                                                                >  1:748476‑M1/1‑19 (MQ=255)
                                                                                                                                                                                                                        cAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataact                                                                                                <  2:748476‑M1/139‑121 (MQ=255)
                                                                                                                                                                                                                        cAATTACTTAACATAAATGtataagtcatacttttgttttgggt                                                                                                                                                                                               >  2:1051801‑M1/1‑19 (MQ=255)
                                                                                                                                                                                                                        cAATTACTTAACATAAATGtataagtcatacttttgttttgggt                                                                                                                                                                                               <  1:1051801‑M1/44‑26 (MQ=255)
                                                                                                                                                                                                                                 aaCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcacc                                                                                            >  2:456779‑M1/1‑10 (MQ=255)
                                                                                                                                                                                                                                 aaCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcacc                                                                                            <  1:456779‑M1/134‑125 (MQ=255)
                                                                                                                                                                                                                                 aaCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagt                                                                                                          <  2:81723‑M1/120‑111 (MQ=255)
                                                                                                                                                                                                                                 aaCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagt                                                                                                          >  1:81723‑M1/1‑10 (MQ=255)
                                                                                                                                                                                                                                        aTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaa                                                                                                                                           >  2:377994‑M1/1‑3 (MQ=255)
                                                                                                                                                                                                                                        aTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaa                                                                                                                                           <  1:377994‑M1/80‑78 (MQ=255)
                                                                                                                                                                                                                                         tGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaatcaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcactttcatgcaataccagataagctatttttaaa                               <  1:452068‑M1/187‑186 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                      
CCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCATACAGCGGCCAG  >  NC_000913/1978268‑1978716

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