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 57 [1] [1] 58 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 = 19785021 (0.020)56 (1.070) 45/270 0.2 99.2% 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

AAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCG  >  NC_000913_3_bme_pgi/1978357‑1978643
                                                                                                                                                 |                                                                                                                                             
aaaTGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtat                                                                                                                                            <  2:32471‑M1/149‑4 (MQ=255)
aaaTGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtat                                                                                                                                            <  1:646423‑M1/149‑4 (MQ=255)
        cTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcata                                                                                                                                    <  2:44367‑M1/149‑12 (MQ=255)
         tttACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTCCTTAACATAAATGtataagtcatac                                                                                                                                   >  1:618802‑M1/1‑137 (MQ=255)
         tttACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcacac                                                                                                                                   >  1:57103‑M1/1‑137 (MQ=255)
             cacTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACGCTCAATTACTTAACATAAATGtataagtcatactttt                                                                                                                               <  2:770787‑M1/149‑17 (MQ=255)
                ttATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACTACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCCATTACTTAACATAAATGtgtacgtcatacttttgtt                                                                                                                            >  1:700602‑M1/1‑130 (MQ=255)
                 tATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttt                                                                                                                           <  1:707587‑M1/149‑21 (MQ=255)
                  aTCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtttt                                                                                                                          <  1:412036‑M1/149‑22 (MQ=255)
                    cTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                                        <  1:300906‑M1/149‑24 (MQ=255)
                    cTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                                        <  2:661309‑M1/149‑24 (MQ=255)
                      aaGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatccttttgttttgggt                                                                                                                      >  1:36109‑M1/1‑124 (MQ=255)
                       aGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtg                                                                                                                     >  1:657053‑M1/1‑123 (MQ=255)
                         aTTTTTCCTAAATCGACGTAACTGTACTCGTGACTACACGTACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgta                                                                                                                   <  1:510192‑M1/149‑29 (MQ=255)
                                          gCAACTGTACTCGTCACTACACGCACATAGAACGGAGGGGGGCTGCGATTTTCAGTACTGCATGATGCAGGTCACACAAAACACTCGATTACTTAACATAAATGtataagtcatacttttgttttaggtgtatttcaatctgttcaaaa                                                                                                  >  1:851033‑M1/1‑104 (MQ=255)
                                          gCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaa                                                                                                  <  2:566355‑M1/149‑46 (MQ=255)
                                          gCAACTGTACTCGTCACGACACGCACATACAACGGAGGGGGGCTGGGATTTTCACTGATGCGTGATGCAGTTGACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtatatttcaatctgttaaaaa                                                                                                  <  2:851033‑M1/149‑46 (MQ=255)
                                            aaCTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagt                                                                                                >  1:760225‑M1/1‑102 (MQ=255)
                                            aaCTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtcttgggtgtatttcaatctgttaaaaatt                                                                                                >  1:844810‑M1/1‑102 (MQ=255)
                                             aCTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCTAATACTTAACATAAATGtataagtcatacttttgttttggttgtatttcaatctgttaaaaagtt                                                                                               <  1:991311‑M1/149‑49 (MQ=255)
                                                gTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATATATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttttt                                                                                            <  2:1002039‑M1/149‑52 (MQ=255)
                                                 ttctcGGCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTTAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttc                                                                                           <  2:36109‑M1/147‑53 (MQ=255)
                                                 tACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttc                                                                                           >  1:257517‑M1/1‑97 (MQ=255)
                                                 tACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttc                                                                                           >  1:641443‑M1/1‑97 (MQ=255)
                                                      gTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctac                                                                                      >  2:878784‑M1/1‑92 (MQ=255)
                                                      gTCACTACACGCACATACAACGGAGGGGGGCTGCGATTATCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctac                                                                                      >  2:883688‑M1/1‑92 (MQ=255)
                                                          cTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgcta                                                                                  <  1:859634‑M1/149‑62 (MQ=255)
                                                               cGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAACTAtatcgagcaaactatgggtttgggtgtatttcaatctgttaaaaagtttttcgctacgctcgcagg                                                                             >  1:234859‑M1/1‑83 (MQ=255)
                                                                        aaCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctcgaaagcaaaaatga                                                                    >  2:780333‑M1/1‑74 (MQ=255)
                                                                             aGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacag                                                               >  1:244546‑M1/1‑69 (MQ=255)
                                                                             aGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacag                                                               <  2:257517‑M1/149‑81 (MQ=255)
                                                                                 gggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttgaatctgttaaaaagtttttccctacgctagcaacacaaaatgaaaacggaat                                                           >  2:581784‑M1/1‑65 (MQ=255)
                                                                                             tCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAAATAAATGtataagtcataatttcgttttgggtgtatttcaatctgttaacaagtttttcgctacgctagcaagcaaaaattaaacagaaataatcgaaactga                                               >  2:238278‑M1/1‑53 (MQ=255)
                                                                                             ccaATAAGGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggg                                               <  1:198586‑M1/148‑97 (MQ=255)
                                                                                              cAATAATGCGTGATGCAGATCACACAAAACACGCAAATACTGAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggga                                              <  1:91038‑M1/149‑98 (MQ=255)
                                                                                                 taatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgt                                           <  2:641443‑M1/149‑101 (MQ=255)
                                                                                                   atGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatccttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttg                                         <  1:954422‑M1/149‑103 (MQ=255)
                                                                                                   atGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaatactcgaaatgggatgttg                                         <  1:878784‑M1/149‑103 (MQ=255)
                                                                                                   atGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttg                                         <  1:883688‑M1/149‑103 (MQ=255)
                                                                                                          aTGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacag                                  <  1:151307‑M1/149‑110 (MQ=255)
                                                                                                              aGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaa                              >  2:157186‑M1/1‑36 (MQ=255)
                                                                                                              aGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatg                                                 <  2:661971‑M1/130‑95 (MQ=255)
                                                                                                              aGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatg                                                 >  1:661971‑M1/1‑36 (MQ=255)
                                                                                                                   acacaAAACAGTGGGTTACTTAACATAGAGGtataagtcgtactgttgtgttgggtggatttcaatctgttaaaaaggtttgcgctacgctagcaagcaaaaaggaagcaggaataatcgaaatgggatgttgcgcacagtcaaaataa                         <  2:157709‑M1/149‑119 (MQ=255)
                                                                                                                   acacaAAACACTTAATTCTTTAACTTATATTgataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttgttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataa                         <  2:430190‑M1/149‑119 (MQ=255)
                                                                                                                       aaaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactca                     <  1:157186‑M1/149‑123 (MQ=255)
                                                                                                                       aaaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgtt                                          <  2:34061‑M1/128‑102 (MQ=255)
                                                                                                                       aaaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgtt                                          >  1:34061‑M1/1‑27 (MQ=255)
                                                                                                                         aaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcacc                   >  1:57543‑M1/1‑25 (MQ=255)
                                                                                                                            acTCAATTACTTAACATAAATGtataagtcatacttttgtattgggtgtatttcaatatgttagaaagtttttcgctacgcaagcaagcaaaaatgaaacaggaataatcgaactgggatgttgctaacagcccaaataacgccccgta                >  2:606675‑M1/1‑22 (MQ=255)
                                                                                                                              tCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaa              >  1:551551‑M1/1‑20 (MQ=255)
                                                                                                                               cAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaat             >  1:607334‑M1/1‑19 (MQ=255)
                                                                                                                                 aTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataa           <  1:253626‑M1/149‑133 (MQ=255)
                                                                                                                                     cTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcat       <  1:449300‑M1/149‑137 (MQ=255)
                                                                                                                                      ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcccagtcaaaataactcaccgtaaataatcatc      >  1:474502‑M1/1‑12 (MQ=255)
                                                                                                                                        aaCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgagcaagatgaaattaactgaccgtaaataatcatctg    >  2:963337‑M1/1‑10 (MQ=255)
                                                                                                                                          gcgtattGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGATGCTGCCTCAGGTTCAGGTTATGCCg  >  2:489678/7‑149 (MQ=9)
                                                                                                                                                 |                                                                                                                                             
AAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCG  >  NC_000913_3_bme_pgi/1978357‑1978643

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