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

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913_3_bme_pgi 1978503–1979268 1979278 11–776 62 [0] [0] 62 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)55 (1.250) 39/256 NT 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                                                                                                                                             <  1:631660‑M1/149‑4 (MQ=255)
  aTGGGTCTGTGGACTTATCGAAGTTGTTTCCTATATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACATAACACTCAATTACTTAACATAAATGtataa                                                                                                                                           <  2:93063‑M1/149‑6 (MQ=255)
    gggTCTTTACACTTATGTAAGATTTGACCTAAATCGACGCAACTGTACTCGCAACTACACGCACATCCACCGGAGGGTGGCTGCTATCTTCAATAATGCCTGATCCAGATCCCACAAAACTCTCAATTACTTAACATAAATgtacaagt                                                                                                                                         <  2:575296/149‑8 (MQ=255)
         tttACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCGGACAAAACACTCAATTACTTAACATAAATGtataagtcatac                                                                                                                                    >  1:744489‑M1/1‑137 (MQ=255)
         tttACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatac                                                                                                                                    >  2:866935‑M1/1‑137 (MQ=255)
         tttACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatac                                                                                                                                    >  1:450347‑M1/1‑137 (MQ=255)
         tttACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGAGGGCTGCGATTTTCAATACTCCCCCATGCACATCACACAAAACACACAATTACTTAACATAAATGtataagtcatac                                                                                                                                    >  2:317639‑M1/1‑137 (MQ=255)
                    cTAATATTTTTCTTAAATCGGCGCAACTGTACTCGTCACTACACGCACATACAAAGAAGGATGGATGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgt                                                                                                                         <  2:680457‑M1/149‑24 (MQ=255)
                    cTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACAGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                                         <  2:199684‑M1/149‑24 (MQ=255)
                     tAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGCGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTAAATTACTTAACATAAATGtataagtcatacttttgttttggg                                                                                                                        >  1:612784‑M1/1‑125 (MQ=255)
                     tAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtaatacttttgttttggg                                                                                                                        >  1:463454‑M1/1‑125 (MQ=255)
                     tAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGAGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtttttgg                                                                                                                        >  2:232872‑M1/1‑125 (MQ=255)
                         aTTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGAGTGATGCAAATCACACAAAAAACTCAATTACTTAACATAAATGaataagccataattttgtttcgggtata                                                                                                                    >  2:655016‑M1/1‑121 (MQ=255)
                         aTTTTTCCTAAATCGACGCAACTGCACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgta                                                                                                                    >  1:779255‑M1/1‑121 (MQ=255)
                         aTTTTTCCTAAATCAACGCAACTGCACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgta                                                                                                                    >  1:779004‑M1/1‑121 (MQ=255)
                                 tAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatc                                                                                                            <  1:441556‑M1/149‑37 (MQ=255)
                                 tAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgatttgggtgtatttcaatc                                                                                                            >  1:632816‑M1/1‑113 (MQ=255)
                                     tCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgtt                                                                                                        >  2:754777‑M1/1‑109 (MQ=255)
                                        aCGCAACTGTACTCGTCACTACACTCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaa                                                                                                     <  1:674125‑M1/149‑44 (MQ=255)
                                          tcaACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagttacactttctttttgggtgtatttcaatctgttaaaaa                                                                                                   <  1:892238‑M1/148‑46 (MQ=255)
                                          ggaactGTACTCGACACTAAACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaa                                                                                                   <  2:626734‑M1/147‑46 (MQ=255)
                                            aaCTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCAGACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtattagatgtatttcaatatgataaaaagt                                                                                                 >  2:284672‑M1/1‑102 (MQ=255)
                                            aaCTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagt                                                                                                 >  1:756477‑M1/1‑102 (MQ=255)
                                            aaCTGTACTCGTCACTACACGCACATACAACGGAGGAGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtatatttcaatctgttaccaacg                                                                                                 >  2:444072‑M1/1‑102 (MQ=255)
                                                gTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttttt                                                                                             <  1:819274‑M1/149‑52 (MQ=255)
                                                gTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttttt                                                                                             <  1:218656‑M1/149‑52 (MQ=255)
                                                   cTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgc                                                                                          <  2:868084‑M1/149‑55 (MQ=255)
                                                    tCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttctggctgtattacaatatattaaaaagcttttctca                                                                                         >  2:806669‑M1/1‑94 (MQ=255)
                                                    tCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATAACACAAAAGCCTCAATTACTTAACATAAATGtataaatcatacgtttattttgggtgtatgtcaatctgttaaaaaatttttcgct                                                                                         >  1:525097‑M1/1‑94 (MQ=255)
                                                       tccctacACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACTCTCAATTACTTAACATTAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacg                                                                                      <  1:270452‑M1/146‑59 (MQ=255)
                                                                     tACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatatgttaaaaagtttttcgctccgctagcaggcacaaa                                                                        >  2:510132‑M1/1‑77 (MQ=255)
                                                                     tACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgctttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaa                                                                        >  2:679316‑M1/1‑77 (MQ=255)
                                                                      aCAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcg                                                                                           <  1:621065‑M1/129‑54 (MQ=255)
                                                                      aCAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcg                                                                                           >  2:621065‑M1/1‑76 (MQ=255)
                                                                         aCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaa                                                                    >  2:124175‑M1/1‑73 (MQ=255)
                                                                          cGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaa                                                                   >  1:446426‑M1/1‑72 (MQ=255)
                                                                             aGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacag                                                                <  1:113132‑M1/149‑81 (MQ=255)
                                                                                       cGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagattttcgctactctagccagccaaaatggaaaagcaataatcca                                                      >  2:226552‑M1/1‑59 (MQ=255)
                                                                                         aTTTTCAATAATGCGTGATGCAGATCAAAAAAATCATTCAATTACTTAACATAAATGtataagttatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaa                                                    <  2:721474‑M1/149‑93 (MQ=255)
                                                                                           tttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttaaatatgttaaaaagtttttagctacgatagcaagcaaaaatgaaacaggaataatagaaatg                                                  >  1:742661‑M1/1‑55 (MQ=255)
                                                                                              cAATAATGCATGATGAAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgtgtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaacgaaacaggaataatcgaaatggga                                               <  2:553761‑M1/149‑98 (MQ=255)
                                                                                                       gTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttggatgtattttaatctgttaaaaactttttcgctacgctagcaagcaaaaataacacaggaataatcgaaatggtatgttgcgac                                      >  2:337074‑M1/1‑43 (MQ=255)
                                                                                                         gATGCAGATTATACAAAATACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaacagtttctcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcaca                                    <  1:849432‑M1/149‑109 (MQ=255)
                                                                                                         gATGCAGATCACACAAAATACTTAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaattttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcaca                                    <  1:393741‑M1/149‑109 (MQ=255)
                                                                                                         gATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcaca                                    <  2:549031‑M1/149‑109 (MQ=255)
                                                                                                             cAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtca                                <  1:820522‑M1/149‑113 (MQ=255)
                                                                                                               gATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtaaaa                              >  1:86839‑M1/1‑35 (MQ=255)
                                                                                                                aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggaagttacgcacagtcaaaa                             >  2:890852‑M1/1‑34 (MQ=255)
                                                                                                                 ttacacaaaACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaat                            <  1:32511‑M1/147‑117 (MQ=255)
                                                                                                                     acaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataact                        >  2:785851‑M1/1‑29 (MQ=255)
                                                                                                                     aaaaaacaCTCAATTACTTAACATAAATGtataagttatacttttgttttaggtgtatttcaatctgttaaaaagaatttcgcgacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataact                        <  1:123227‑M1/147‑121 (MQ=255)
                                                                                                                        aaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctccgataggaagcaaaaattaaaaaggaataatcgaaatgggatgttgcgcacagtcaaaataacttca                     >  2:802244‑M1/1‑26 (MQ=255)
                                                                                                                              gcaATTACATAAAATAAATGaaatagtcatagttttgttttcggtgtatttcaatctgttaaaaagtttttcgctaggctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaa               <  2:93153‑M1/148‑130 (MQ=255)
                                                                                                                                aaTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctactctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataaatcaccgtaaata             >  1:514366‑M1/1‑18 (MQ=255)
                                                                                                                                 aTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataa            >  1:483245‑M1/1‑17 (MQ=255)
                                                                                                                                  ttACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgc                                         <  2:168327‑M1/119‑104 (MQ=255)
                                                                                                                                  ttACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgc                                         >  1:168327‑M1/1‑16 (MQ=255)
                                                                                                                                    aCTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaaaaatca         >  1:762385‑M1/1‑14 (MQ=255)
                                                                                                                                     ttgacaataaatgaataagttatacttttgttttggttgtatttctatcttttaaaaatttttttgctacgctagaaagccaaaatgaaacaggaaaaatcgaaaggggatgttgcgcacagtcaaaataactcaccgtaaataatcat        <  1:450200‑M1/143‑137 (MQ=255)
                                                                                                                                     cTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactacccgtaaacaaacat        >  2:668448‑M1/1‑13 (MQ=255)
                                                                                                                                     cTTAACATAAATGtaaaagttacacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcatcgtaaataatcat        <  2:446426‑M1/149‑137 (MQ=255)
                                                                                                                                           aTAAATGtataagtcatactttttttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagaaaaaatgaaaaaggaataatcgaaatgggatgttgcgcacagtcaaaataaattacagtaaataatcatctgctc  >  2:534165‑M1/1‑7 (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 < 27 ≤ 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.