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 33 [0] [0] 33 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)30 (0.710) 22/270 1.2 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

AAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGC  >  NC_000913_3_bme_pgi/1978357‑1978644
                                                                                                                                                 |                                                                                                                                              
aaaTGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtat                                                                                                                                             <  2:193399‑M1/149‑4 (MQ=255)
aaaTGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtat                                                                                                                                             <  1:307376‑M1/149‑4 (MQ=255)
aaaTGGGTCTTTACACTTATCTAACAGTTTTCCTGGATCGAAGCAACTGTACTCGTCACTACACGCACATACAACGTAGGGTGGCTGCGATTTTAAATAATGTGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAAGGtat                                                                                                                                             <  2:51164‑M1/149‑4 (MQ=255)
             cacTTATCTAAGTTTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatactttt                                                                                                                                <  1:439066‑M1/149‑17 (MQ=255)
             cacTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatactttt                                                                                                                                <  1:290173‑M1/149‑17 (MQ=255)
                  aTCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTTATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtttt                                                                                                                           <  1:725442‑M1/149‑22 (MQ=255)
                    cTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                                         <  2:579015‑M1/149‑24 (MQ=255)
                    cTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                                         <  1:474395‑M1/149‑24 (MQ=255)
                    cTAAGATTTTTCCTAAAACGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCGCACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                                         <  1:508288‑M1/149‑24 (MQ=255)
                                cTAAATCGACGCAACTGTACTCGTGACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATTCGTGATGCAGATCACACAAAACACTCAATTACTTAACATCAATGtgtaagtgataattttgttttgggtgtatttcaat                                                                                                             >  2:87628‑M1/1‑114 (MQ=255)
                                            aaCTGTACTCGTCACTACACGCACTTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATTCAGATCACACAAAACACTCAATTACTTAACATAAATGtataaatcatactttttttttgggtgtttttcattctggttaaaagc                                                                                                 >  2:159263‑M1/1‑102 (MQ=255)
                                                    tCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcattctgttaaaaagtttttcgct                                                                                         >  1:643947‑M1/1‑94 (MQ=255)
                                                    tCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgct                                                                                         >  2:51471‑M1/1‑94 (MQ=255)
                                                        cACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgc                                                                                     <  2:510832‑M1/149‑60 (MQ=255)
                                                                     tACAACGGAGGGGGGCTGCGATTTGCAATAAGGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaa                                                                        <  1:546629‑M1/149‑73 (MQ=255)
                                                                     tACAACGGAGAGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtttttggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaa                                                                        <  1:679093‑M1/149‑73 (MQ=255)
                                                                        aaCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatga                                                                     <  1:477467‑M1/149‑76 (MQ=255)
                                                                         aCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtca                                                                                                                                       >  2:163454‑M1/1‑73 (MQ=255)
                                                                         aCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtca                                                                                                                                       <  1:163454‑M1/82‑10 (MQ=255)
                                                                                               aataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgcaacaggaataatcgaaatgggat                                              <  1:741706‑M1/149‑99 (MQ=255)
                                                                                                ataatGCGTGATGCAGATCAAACAAAAAACTCAATTACTTAACATAAATGgataagtaatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaaacaaaaatgaaacaggaataacagaaacaggatt                                             >  2:265350‑M1/1‑50 (MQ=255)
                                                                                                   atGCGCGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttg                                          <  2:247272‑M1/149‑103 (MQ=255)
                                                                                                                aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaa                             <  1:762302‑M1/149‑116 (MQ=255)
                                                                                                                aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaa                             >  2:424076‑M1/1‑34 (MQ=255)
                                                                                                                 tCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcattctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgatatgggatgttgcgcacagtcacaat                            >  1:118096‑M1/1‑33 (MQ=255)
                                                                                                                    cacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataac                         >  2:179758‑M1/1‑30 (MQ=255)
                                                                                                                              tCAATTACTTAACATAAATGtataagtcatactcttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgtgatgttgcgcacaatcaaaataactcagcgtaaa               >  1:619747‑M1/1‑20 (MQ=255)
                                                                                                                                        aaCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctg     <  1:424076‑M1/149‑140 (MQ=255)
                                                                                                                                        aaCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataagcctctg     >  1:274884‑M1/1‑10 (MQ=255)
                                                                                                                                        aaCATAAATGtataagtcatacttttgttttgggtgtatttcaatctggtaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctg     >  1:184193‑M1/1‑10 (MQ=255)
                                                                                                                                           aTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgcta  >  1:340016‑M1/1‑7 (MQ=255)
                                                                                                                                             aaaTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgc                                       <  1:684342‑M1/110‑106 (MQ=255)
                                                                                                                                             aaaTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgc                                       >  2:684342‑M1/1‑5 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
AAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGC  >  NC_000913_3_bme_pgi/1978357‑1978644

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.