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–1979270 1979278 9–776 68 [2] [0] 64 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.030)60 (1.080) 48/282 0.1 98.5% 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

ATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTC  >  NC_000913_3_bme_pgi/1978353‑1978651
                                                                                                                                                     |                                                                                                                                                     
atattaatgggtCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACAAGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTTATGCAGATCACACAAAACACTCAATTACTTAACATAAATG                                                                                                                                                       <  2:918139‑M1/145‑1 (MQ=255)
aTAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAATGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATG                                                                                                                                                       <  1:908742‑M1/150‑1 (MQ=255)
 tAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGt                                                                                                                                                      >  1:200112‑M1/1‑149 (MQ=255)
      aTGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataag                                                                                                                                                 <  2:494483‑M1/150‑7 (MQ=255)
        gggTCTTTACACTTATCTAAGATTTTTCCTACATCGACGCAACTGTCCTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtc                                                                                                                                               <  2:809307‑M1/150‑9 (MQ=255)
         ggTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtca                                                                                                                                              >  1:413412‑M1/1‑141 (MQ=255)
          gTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcat                                                                                                                                             >  2:372651‑M1/1‑140 (MQ=255)
                 cacTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGAGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttg                                                                                                                                      >  1:681160‑M1/1‑133 (MQ=255)
                       tCTAAGATTTTTCCTAAATCGACGCAACAGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtat                                                                                                                                                    <  2:818756‑M1/130‑4 (MQ=255)
                       tCTAAGATTTTTCCTAAATCGACGCAACAGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtat                                                                                                                                                    >  1:818756‑M1/1‑127 (MQ=255)
                         tAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggt                                                                                                                              >  1:174788‑M1/1‑125 (MQ=255)
                             aTTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtat                                                                                                                          >  1:754099‑M1/1‑121 (MQ=255)
                             aTTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtat                                                                                                                          <  2:200112‑M1/150‑30 (MQ=255)
                                    cTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatc                                                                                                                   >  2:787972‑M1/1‑114 (MQ=255)
                                      aaaTCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctg                                                                                                                 <  1:441455‑M1/150‑39 (MQ=255)
                                          cgccgcaACTGTACTCGTCAATACACGCACATACAACGTAGCGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATTtataagtcaaagttttgttttgggtgtatttcaatcagtaaa                                                                                                             >  2:27513‑M1/4‑108 (MQ=255)
                                          cGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaa                                                                                                             >  2:92628‑M1/1‑108 (MQ=255)
                                             cGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaa                                                                                                          <  1:897660‑M1/150‑46 (MQ=255)
                                             agcAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaa                                                                                                          <  1:722005‑M1/149‑46 (MQ=255)
                                                 aCTGTACTCGTCACTACACGCACATACAACGGATTGGGTCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttt                                                                                                      <  1:691838‑M1/150‑50 (MQ=255)
                                                   tGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttttt                                                                                                    >  1:818362‑M1/1‑99 (MQ=255)
                                                    gTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttc                                                                                                   <  1:92628‑M1/150‑53 (MQ=255)
                                                         cGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTAAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctaa                                                                                              >  1:925086‑M1/1‑93 (MQ=255)
                                                         cGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTACAATAAAGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagacatacttttgttttgggtgtatttcaatatgttaaaaagtttttcgctac                                                                                              >  1:477761‑M1/1‑93 (MQ=255)
                                                              cTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctag                                                                                         <  1:599596‑M1/150‑63 (MQ=255)
                                                               tACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagc                                                                                        <  1:93954‑M1/150‑64 (MQ=255)
                                                                   cGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagc                                                                                    <  2:709512‑M1/150‑68 (MQ=255)
                                                                         tACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaat                                                                              >  2:809252‑M1/1‑77 (MQ=255)
                                                                          aCAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatt                                                                             >  2:604746‑M1/1‑76 (MQ=255)
                                                                            aaCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACAATCAATTACTTAACATAAATGtattagtcatacttttgttttgcgtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaa                                                                           >  1:396951‑M1/1‑74 (MQ=255)
                                                                             aCGGAGGGGGGCTGCGTTTTTCAATAATGCTTGATGCAGATCACACAAAACAATCAATTACTTAACATATATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaattttttcgctacgctagcaagcaaaaatgaaa                                                                          <  2:477761‑M1/150‑78 (MQ=255)
                                                                                gAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacag                                                                       >  1:471450‑M1/1‑70 (MQ=255)
                                                                                gAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacag                                                                       >  1:278690‑M1/1‑70 (MQ=255)
                                                                                 aGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacagg                                                                      <  1:860200‑M1/150‑82 (MQ=255)
                                                                                        cTGCGATTTTCAAAAATGCGTGATGCAGATCACACATAACACTCAATTACTTAACATAAATGtataagtcatacttttgttatgggtgtatttcaatctgttaaaaagtttttcgctacgctagctagcaaaaatgaaacaggaataatc                                                               <  2:733168‑M1/150‑89 (MQ=255)
                                                                                         tGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcg                                                              <  2:681160‑M1/150‑90 (MQ=255)
                                                                                         tGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcg                                                              >  2:309264‑M1/1‑61 (MQ=255)
                                                                                         tGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcg                                                              <  2:389591‑M1/150‑90 (MQ=255)
                                                                                            gATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaa                                                           >  1:91534‑M1/1‑58 (MQ=255)
                                                                                            gATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaa                                                           >  2:903271‑M1/1‑58 (MQ=255)
                                                                                            gATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatagaaa                                                           >  2:548791‑M1/1‑58 (MQ=255)
                                                                                                ttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggg                                                       <  2:174788‑M1/150‑97 (MQ=255)
                                                                                                 tCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggga                                                      >  2:61987‑M1/1‑53 (MQ=255)
                                                                                                      aatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttatgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttg                                                 >  1:757743‑M1/1‑48 (MQ=255)
                                                                                                         gCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgtgtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgc                                              <  1:489469‑M1/150‑106 (MQ=255)
                                                                                                              aTGCAGATCACACAAACCACTCAATTACTTAACATAAATGgataagtcatactgttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaaggggatgttgcgcacagt                                         >  2:69274‑M1/1‑40 (MQ=255)
                                                                                                                    aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaat                                   <  1:117778‑M1/150‑117 (MQ=255)
                                                                                                                    aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaat                                   <  2:91534‑M1/150‑117 (MQ=255)
                                                                                                                    aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaat                                   <  1:431088‑M1/150‑117 (MQ=255)
                                                                                                                    aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaat                                   >  1:109459‑M1/1‑34 (MQ=255)
                                                                                                                    aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaaa                                   >  1:297175‑M1/1‑34 (MQ=255)
                                                                                                                    aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaaacaaaaatgaaacaggaataatcgaaatgggatgttgcgcacaggcaaaat                                   <  1:807140‑M1/150‑117 (MQ=255)
                                                                                                                      cacacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttggtttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggattaatcgaaatgggatgttgcgcacagtcaaaataa                                 >  2:46736‑M1/1‑32 (MQ=255)
                                                                                                                        cacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatctaaatgggatgttgcgcacagtcaaaataact                               >  1:757351‑M1/1‑30 (MQ=255)
                                                                                                                            aaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcacc                           >  2:790516‑M1/1‑26 (MQ=255)
                                                                                                                                acTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaa                       >  2:346032‑M1/1‑22 (MQ=255)
                                                                                                                                      ttACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatc                 >  1:94107‑M1/1‑16 (MQ=255)
                                                                                                                                       tACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatca                <  1:372651‑M1/150‑136 (MQ=255)
                                                                                                                                        aCTTAACATAATTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaattgaaacaggaataatcgaaatgggatgttgcgcacagtcaaattaactcaccgtaaataatcat               >  2:867043‑M1/1‑14 (MQ=255)
                                                                                                                                          ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatct             <  1:346032‑M1/150‑139 (MQ=255)
                                                                                                                                          ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctaagctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatct             <  2:925086‑M1/150‑139 (MQ=255)
                                                                                                                                           tAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctg            >  1:225223‑M1/1‑11 (MQ=255)
                                                                                                                                              cATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttctctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctt         >  2:192567‑M1/1‑8 (MQ=255)
                                                                                                                                               cgtattGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCt        <  1:900462/145‑1 (MQ=14)
                                                                                                                                                tAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctata       >  1:237634‑M1/1‑6 (MQ=255)
                                                                                                                                                tAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctata       <  1:809252‑M1/150‑145 (MQ=255)
                                                                                                                                                     gtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataa  >  2:196206‑M1/1‑1 (MQ=255)
                                                                                                                                                     gggTAATGACTCCAACTTATTGATAGTGATTTATGTTCAGATAATGCACGATGACTTTGTCATGCAGCTCGCCCGATTTTGAGAACGACAGCGACTTCCGACACAGCAGTGCCAGGTGCTGCCTCAGATTCTGGTTATGCCTGTCAAtta  >  1:648172/1‑149 (MQ=11)
                                                                                                                                                     |                                                                                                                                                     
ATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTC  >  NC_000913_3_bme_pgi/1978353‑1978651

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