Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913_3_bme_pgi 1,978,503 Δ776 bp insB1insA insB1, insA

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913_3_bme_pgi 1978503 1979278 776 40 [0] [0] 41 insB1–insA insB1,insA

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_bme_pgi = 19785020 (0.000)37 (0.820) 30/268 0.6 100% intergenic (‑305/+16) flhD/insB1 flagellar class II regulon transcriptional activator, with FlhC/IS1 transposase B
?NC_000913_3_bme_pgi 1979279 = 0 (0.000)intergenic (‑64/‑474) insA/uspC IS1 repressor TnpA/universal stress protein

GGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAA  >  NC_000913_3_bme_pgi/1978361‑1978648
                                                                                                                                             |                                                                                                                                                  
gggTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagt                                                                                                                                             <  1:395222‑M1/149‑8 (MQ=255)
          acTTATCTAAGATGTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTCAACATAAATGtataagtcatacttttg                                                                                                                                   >  1:363325‑M1/1‑132 (MQ=255)
                cTAAGATTTTTCCTAAATCGACGCAGCTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATACGTGATGCAGATCACACAAAGCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                                             <  1:754221‑M1/149‑24 (MQ=255)
                cTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                                             <  1:289135‑M1/149‑24 (MQ=255)
                 tAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttggt                                                                                                                            >  1:225341‑M1/1‑125 (MQ=255)
                 tAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttggg                                                                                                                            >  1:590161‑M1/1‑125 (MQ=255)
                        tttCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgcattt                                                                                                                     <  2:637032‑M1/149‑32 (MQ=255)
                        tttCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGATGGGATTTTGAATAATGAGCGAGGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttggtttgggtgtattt                                                                                                                     >  1:823392‑M1/1‑118 (MQ=255)
                            cTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaat                                                                                                                 >  2:203858‑M1/1‑114 (MQ=255)
                            cTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaat                                                                                                                 >  2:406046‑M1/1‑114 (MQ=255)
                                 tCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgtt                                                                                                            >  1:500192‑M1/1‑109 (MQ=255)
                                                tCGTCACTACACGCACATACAACGGAGGAGGGATGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgct                                                                                             >  2:825967‑M1/1‑94 (MQ=255)
                                                      cTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgcta                                                                                       <  1:751264‑M1/149‑62 (MQ=255)
                                                          acGCACATACAACGGAGGGGGGCTGCGATTTTCAACAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaa                                                                                   <  1:365409‑M1/149‑66 (MQ=255)
                                                            gCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagc                                                                                 >  2:113150‑M1/1‑82 (MQ=255)
                                                             gacatACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtttggggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagca                                                                                <  2:424952‑M1/148‑69 (MQ=255)
                                                               caTACAACGGAGGGGGGCTGCGATTTTAAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgcta                                                                                       <  1:549606‑M1/140‑62 (MQ=255)
                                                               caTACAACGGAGGGGGGCTGCGATTTTAAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgcta                                                                                       >  2:549606‑M1/1‑79 (MQ=255)
                                                                    aaCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgc                                                                                         <  2:459057‑M1/133‑60 (MQ=255)
                                                                    aaCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgc                                                                                         >  1:459057‑M1/1‑74 (MQ=255)
                                                                                cTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataat                                                             <  1:203858‑M1/149‑88 (MQ=255)
                                                                                    gATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                                             >  1:566185‑M1/1‑58 (MQ=255)
                                                                                    gATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                                             <  2:566185‑M1/81‑24 (MQ=255)
                                                                                       tttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaactg                                                      >  1:173178‑M1/1‑55 (MQ=255)
                                                                                       tttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatg                                                      >  1:576423‑M1/1‑55 (MQ=255)
                                                                                            ataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatg                                                 >  1:432978‑M1/1‑50 (MQ=255)
                                                                                              aatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataaacgaaatgggatgtt                                               >  1:494622‑M1/1‑48 (MQ=255)
                                                                                                 gCGTGATGCAGATCACACAAAACACTCAAATACTTAACGTAAATGtataactcatacgtttgttttgggtgtatttcaatctgttagaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgctgcg                                            >  1:707683‑M1/1‑45 (MQ=255)
                                                                                                   gTGATGCAGATCACACAAAACACTGAATTACTTAACATAAATGtataagtcatacttttgttttaggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttacgca                                          >  1:725816‑M1/1‑43 (MQ=255)
                                                                                                        gCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtc                                     >  1:514993‑M1/1‑38 (MQ=255)
                                                                                                        gCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatt                                                                                                                      >  2:35553‑M1/1‑38 (MQ=255)
                                                                                                        gCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatt                                                                                                                      <  1:35553‑M1/68‑31 (MQ=255)
                                                                                                            aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttggcagcatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaa                                 >  1:539106‑M1/1‑34 (MQ=255)
                                                                                                                     aaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggctgttgctcacagtcaaaatcactcacc                        >  2:337224‑M1/1‑25 (MQ=255)
                                                                                                                      acacTCAATCACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccg                       <  1:561527‑M1/149‑126 (MQ=255)
                                                                                                                             aTTACTTAACATAAATGcataactcatacttttgttttgggtgtatttcaatctgtcaaaaagtttttcgatacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcacaataactcacagtaaataa                <  2:707683‑M1/149‑133 (MQ=255)
                                                                                                                                    aaCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctg         <  2:514993‑M1/149‑140 (MQ=255)
                                                                                                                                         aaaTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagt                                                                                                     >  1:557840‑M1/1‑5 (MQ=255)
                                                                                                                                         aaaTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagt                                                                                                     <  2:557840‑M1/52‑48 (MQ=255)
                                                                                                                                           aTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaa  <  2:363325‑M1/149‑147 (MQ=255)
                                                                                                                                             |                                                                                                                                                  
GGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAA  >  NC_000913_3_bme_pgi/1978361‑1978648

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