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 41 [2] [1] 42 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 = 19785022 (0.040)36 (0.860) 32/270 0.3 97.4% 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

GTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAAT  >  NC_000913_3_bme_pgi/1978363‑1978649
                                                                                                                                           |                                                                                                                                                   
gTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATCAATGtataagtca                                                                                                                                            >  2:422120‑M1/1‑140 (MQ=255)
       cacTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatactttt                                                                                                                                     <  1:494976‑M1/149‑17 (MQ=255)
              cTAAGATTGTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAATCACTTAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                                              <  2:303636‑M1/149‑24 (MQ=255)
               tAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatactttagatttggg                                                                                                                             >  1:559341‑M1/1‑125 (MQ=255)
                        tCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttca                                                                                                                    <  2:699485‑M1/149‑34 (MQ=255)
                          cTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATCACGCGTGAGGCAGAGCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaat                                                                                                                  >  2:510034‑M1/1‑114 (MQ=255)
                          cTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGCTCAAACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaat                                                                                                                  >  1:372055‑M1/1‑114 (MQ=255)
                          cTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTACAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAAATAAATGtataagtcatagttttgttttggtggtatttcaat                                                                                                                  >  2:320260‑M1/1‑114 (MQ=255)
                           tAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatc                                                                                                                 >  2:635684‑M1/1‑113 (MQ=255)
                            aaaTCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTAAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatct                                                                                                                <  1:411896‑M1/149‑38 (MQ=255)
                              ctcGACGCAACTCTGCTCTCCACAACACGCACATACAACTGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgt                                                                                                              <  2:411545‑M1/148‑40 (MQ=255)
                               tCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataaatcatacttttgttttgggtgtatctcaatctgtt                                                                                                             >  2:727059‑M1/1‑109 (MQ=255)
                                           tACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGAGTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttc                                                                                                 >  2:709219‑M1/1‑97 (MQ=255)
                                              tCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgct                                                                                              >  1:554377‑M1/1‑94 (MQ=255)
                                               cGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgcta                                                                                             <  1:24803‑M1/149‑57 (MQ=255)
                                                gTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctac                                                                                            >  1:732335‑M1/1‑92 (MQ=255)
                                                 tCACTACACGCCCATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttaaatctgttcaaacgtttttcgctccg                                                                                           >  2:271989‑M1/1‑91 (MQ=255)
                                                     ttcacgcACATACAACGGAGGCGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTGAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctag                                                                                       <  2:207723‑M1/147‑63 (MQ=255)
                                                     tACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcactctgttgtaatgtttttcgctacgctag                                                                                       >  1:561036‑M1/1‑87 (MQ=255)
                                                            acaTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaa                                                                                <  2:732335‑M1/149‑70 (MQ=255)
                                                               tACAGCGGAGGGGGGCTGCTATTTACAATAATGCGTGCTGCAGATCACATAAAACACTTTATTACTTAACATAAATGtataagtcatatttttgttttgagtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaa                                                                             <  2:118075‑M1/149‑73 (MQ=255)
                                                                   aCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaa                                                                         >  2:721164‑M1/1‑73 (MQ=255)
                                                                                       tCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggg                                                     >  1:268570‑M1/1‑53 (MQ=255)
                                                                                           taatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgt                                                 >  1:633877‑M1/1‑49 (MQ=255)
                                                                                             atGCGTGAAGCAGATCACAAAAAACGCTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttg                                               <  1:422120‑M1/149‑103 (MQ=255)
                                                                                                  tGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtctttcaatctgctgaaaagttgttcgctacgctagcaagcaaaaataaaacagtaataatcgaaatgggatgttgccgac                                          >  2:277675‑M1/1‑42 (MQ=255)
                                                                                                   gATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcaca                                         <  1:267839‑M1/149‑109 (MQ=255)
                                                                                                     tGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagt                                       >  2:108923‑M1/1‑39 (MQ=255)
                                                                                                          aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaa                                  >  1:46961‑M1/1‑34 (MQ=255)
                                                                                                          aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaa                                  >  2:686564‑M1/1‑34 (MQ=255)
                                                                                                              cacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtattttaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataac                              >  2:741470‑M1/1‑30 (MQ=255)
                                                                                                              cacaAAACACTCAATGACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataac                              <  2:414041‑M1/149‑120 (MQ=255)
                                                                                                                      acTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggattttgcgcacagtcaaaataactccccttc                      >  2:308938‑M1/1‑22 (MQ=255)
                                                                                                                          aaTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcagtctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatggaaatgggatgttgcgcacagtcaaaataactcaccgtaaata                  >  1:156913‑M1/1‑18 (MQ=255)
                                                                                                                          aaTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaacataactcaccgtaaata                  >  1:425628‑M1/1‑18 (MQ=255)
                                                                                                                                  aaCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactaaccgtaaataatcatctg          >  2:237358‑M1/1‑10 (MQ=255)
                                                                                                                                     aTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatg                                                       >  2:388923‑M1/1‑7 (MQ=255)
                                                                                                                                     aTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatg                                                       <  1:388923‑M1/101‑95 (MQ=255)
                                                                                                                                          tGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaacaaggaataatcaaaatgggatgttgcgcacagtcaaaatgactcaccgtaaataatcatctgctactaat  >  2:519595‑M1/1‑2 (MQ=255)
                                                                                                                                          cgggTAATGCCTGCAACTTATGGATAGTGTTTTATGTCGAGATAGTGTGTGATGTCTTTGTCATGCAGCTCCACTGATTGTGAGAGCGACAGGGACATCCGGATCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTCTGCCGCTCaaa  <  2:448619/148‑2 (MQ=2)
                                                                                                                                          cgggTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTAATCCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCCGGTTATGCCGCTCAAt  >  1:215926/2‑149 (MQ=2)
                                                                                                                                           |                                                                                                                                                   
GTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAAT  >  NC_000913_3_bme_pgi/1978363‑1978649

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