New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_bme_pgi = 87242946 (0.990)3 (0.070) 3/262 NT 6.3% coding (540/912 nt) gsiD glutathione ABC transporter permease
?NC_000913_3_bme_pgi = 872441 45 (1.000)coding (552/912 nt) gsiD glutathione ABC transporter permease

GCTGATCATGCGCATTTGCGATGTGCTGTTTGCCTTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTGGTGCCAGCGATATGACCGTTTTGTTGCGTCATATCCTGCCTGGGACCGTCTCTTCTATCGTGGTGTTTTTCACCATGCGCATTGGTACCTCG  >  NC_000913_3_bme_pgi/872294‑872588
                                                                                                                                                   |                                                                                                                                                   
gCTGATCATGCGCATTTGCGATGTGCTGTTTGCCTTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCg                                                                                                                                                    <  1:625944‑M2/149‑2 (MQ=255)
      cATGCGCATTTGGGATGTGTTGTTTGCCTTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTCCCCGCCTGGTGCGCGGCAaca                                                                                                                                              <  2:792476/149‑1 (MQ=255)
      cATGCGCATTTGCGATGTGCTGTTTGCCTTCCCGGGTATTTTACTGGCGAGCGCTTTTGTTGCGGTGTTGGGAAGCGGCAGTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAaca                                                                                                                                              <  1:353724/149‑1 (MQ=255)
             aTTTGGGATGTGCTGTTTGCCTTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGgt                                                                                                                                       >  2:152372/1‑149 (MQ=255)
                       tgCTGTTTGCCTTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAg                                                                                                                             >  1:888900/1‑149 (MQ=255)
                                gctgcctgggTAGTTTTCTGGTGGTTGCTGTTGTTGCTGTGTTGGGAAGCGGCATTGCTAACGTTATTATTGCAGTCGCCATTTTTTCCATCACCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCttt                                                                                                                    <  2:63890/142‑1 (MQ=255)
                                 cTTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTACCGTGATTATTGCAGTCGCCATTTATTCCATCACCGCGTTTACCCGCATGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAACCCTTTa                                                                                                                   >  1:18198/1‑149 (MQ=255)
                                  ttCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAac                                                                                                                                               >  2:424435/1‑120 (MQ=255)
                                  ttCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAac                                                                                                                                               <  1:424435/120‑1 (MQ=255)
                                  ttCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAAc                                                                                                                        <  2:261711/143‑1 (MQ=255)
                                  ttCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAAc                                                                                                                        >  1:261711/1‑143 (MQ=255)
                                   tgtcGGGTATTTTACTGGCGATCGCTTTTGTTGGGGTGTTTGGAAGCGGTATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTAtt                                                                                                                 <  2:575740/146‑1 (MQ=255)
                                     ccGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGa                                                                                                               >  1:816351/1‑149 (MQ=255)
                                       gggTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGCGATTATTGCAGTCGCCATTTTTTCCATCCCCGCTTTTGCCCGCCTGTTGCCCGGCACCACGCTGGTTCTTAAACACCAAACCTTCCTTGaac                                                                                                             >  1:656312/1‑147 (MQ=255)
                                            ttttACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCa                                                                                                        >  1:461318/1‑149 (MQ=255)
                                            ttttACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCACGGAGACACTCGTTTGTTGAAACAGCAAACCTTTATTGAGTCAGCa                                                                                                        >  1:468753/1‑149 (MQ=255)
                                                 cTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAg                                                                                                   >  1:315039/1‑149 (MQ=255)
                                                    gCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTAt                                                                                                <  1:588407/149‑1 (MQ=255)
                                                      gATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTg                                                                                              >  1:819362/1‑149 (MQ=255)
                                                        ggggtgttgtGTCGGGGTTGGGGCGCGGAATTGGTAATGTGTTTATTGTAGTCGGCATTTTTGCCATCCCCGCGTTTGCCCGCCTGTGGCGCGGCAACACGGTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGTAGTATTGGt                                                                                            <  2:402319/145‑1 (MQ=255)
                                                             gttgttGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTAGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTGGTGCCAg                                                                                       >  1:842665/1‑149 (MQ=255)
                                                                gttgCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTAGTGCCAGCGa                                                                                    <  2:686312/149‑1 (MQ=255)
                                                                   gCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTGGTGCCAGCGatat                                                                                 >  1:637392/1‑149 (MQ=255)
                                                                     ggTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTGGTGCCAGCGATATGa                                                                               >  1:612705/1‑149 (MQ=255)
                                                                      gtgtTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTGGTGCTAGCGATATGAc                                                                              <  2:479053/149‑1 (MQ=255)
                                                                         gtgagAAGCGGCATTGATAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCGCGCCTGCTGCGCGGCAACACGCGCGTGTTTAAACATCAAACCTTTATTGAGTCACCACGCAGTATTGATGCCATCGCTACTACCgc                                                                           >  2:616597/5‑148 (MQ=255)
                                                                              aaGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGATGTTGAAACGGCAAACCTTTATTGAGTCAGCACGCAGTATTGGTGCCAGCGATATAACCGTTtagt                                                                      >  1:431151/1‑149 (MQ=255)
                                                                               aGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTGGTGCCAGCGATATGACCGTTttgtt                                                                     <  2:625944/149‑1 (MQ=255)
                                                                                 ggtcATTGCTAAAGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTGGTGCCAGCGATATGACCGTTTTGTTGc                                                                   <  1:398568/146‑1 (MQ=255)
                                                                                  ggCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGCGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTGGTGCCAGCGATATGACCGTTTTGTTGCg                                                                  <  2:736937/149‑1 (MQ=255)
                                                                                   gCATTGCTAACTTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTGGTGCCAGCGATATGACCGTTTTGTTGCGt                                                                 <  1:682743/149‑1 (MQ=255)
                                                                                     aTTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTGGTGCCAGCGATATGACCGTTTTGTTGCGTCa                                                               <  1:509517/149‑1 (MQ=255)
                                                                                         cTAACGTGATTATTGCAGTCGCCATTTTTTTCATCGCCGCGTTTGCCCGCCTGGTGCGCggcaaacgcggggctggaaaaaatggcgactgcaataatcacgttagcaattccgcttcccaccaccgcaacaacagcgatcgccaggaa                                                           <  1:438327‑M2/149‑91 (MQ=255)
                                                                                             cGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTGGTGCCAGCGATATGACCGTTTTGTTGCGTCATATcctgc                                                       <  1:152372/149‑1 (MQ=255)
                                                                                                 attattGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTGGTGCCAGCGATATGACCGTTTTGTTGCGTCATATCTTGCCTgg                                                   >  2:234260/1‑149 (MQ=255)
                                                                                                   tattGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTGGTGCCAGCGATATGACCGTTTTGTTGCGTCATATCCTGTCTGGGa                                                 <  2:630871/149‑1 (MQ=255)
                                                                                                   tattGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTGGTGCCAGCGATATGACCGTTTTGTTGCGTCATATCCTGCCTGGGa                                                 <  2:888900/149‑1 (MQ=255)
                                                                                                     ttGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTGGTGCCAGCGATATAACCGTTTTGTTGCGTCATATCCTGCCTGGGAcc                                               >  1:598136/1‑149 (MQ=255)
                                                                                                        tggTCGGCTTTTTTTCCTTCCCTGCTTTTGGCCGCCTGGTGCTCGGCAACACGCTGGTGTTGATACAGCAAACCTTTATTGAGTCTGCACGCAGTATTGGTGCCAGCGATATGTCCGTTTTGTTGCTTCATATCCTGCCTGGGACCGtc                                            <  1:221425/147‑1 (MQ=255)
                                                                                                           tCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTGGTGCCAGCGATATGACCGTTTTGTTGCGTCATATCCTGCCTGGGACCGTctct                                         >  1:45082/1‑149 (MQ=255)
                                                                                                              ccATTTTTTCCGTCCCCGCGTTTGCCCGCCTGGTGGGCggcaaacgcggggatggaaaaaatggcgactgcaataatcacgttagcaatgccgcttcccaacaccgcaacaacagcgatcgccagtaaaatacccgggaaggcaaacag                                      >  2:438327‑M2/1‑38 (MQ=255)
                                                                                                                 ttttttCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTGGTGCCAGCGATATGACCGTTTTGTTGCGTCATATCCTGCCTGGGACCGTCTCTTCTATc                                   >  1:899254/1‑149 (MQ=255)
                                                                                                                 tttttaCCACCCCCTGGTTTGCCCGCCTGGTGCGCGTCAACCCGCTGGTGTTGAAGCAGCAAACCTTTATTGAGTCAGCACGCATTATTGGTGCCAGCGATATGACCGTTTTGTTGCGTCATATCCTGTCTGGGACCGTCTCGTCTATc                                   <  2:60192/149‑1 (MQ=255)
                                                                                                                    tttCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGAAGTATTGGTGCCAGCGATATGACCGTTTTGTTGCGTCATATCCTGCCTGGGACCGTCTCTTCTATCgtg                                <  1:20246/149‑1 (MQ=255)
                                                                                                                         aTCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTGGTGCCAGCGATATGACCGTTTTGTTGCCTCATATCCTCCCTGGGACCGTCTCTTCTATCGTGGCGtt                           >  2:876735/1‑149 (MQ=255)
                                                                                                                         aTCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTGGTGCCAGCGATATGACCGATTTGTTGCGTCATATCCTGCCTGGGACCGTCTCTTCTCTCGTGGTGtt                           >  2:193453/1‑149 (MQ=255)
                                                                                                                            cccggcgttGTCCGCCTGGTGCGCGGCAAGACGTGTGTTTTGAACCAGCAAACCTTTATTGAGTCCGCATGCAGTATTGGTGCCAGCGATATGACCGCTTTGTTGCTTCTTATCCTGCCTTGGACCGTCTCTTCTATCGTGGTGttttt                        <  2:51504/142‑1 (MQ=255)
                                                                                                                            cccGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTGGTGCCAGCGATATGACCGTTTTGTTGCGTCATATCCTGCCTGGGACCGTCTCTTCTATCGTGGTGttttt                        >  1:42301/1‑149 (MQ=255)
                                                                                                                             ccGCGTTTGCCCGCCTGGTGCGGggcaacacgcggggttggaaaaaatggcgactgcaaaaatcacgttagcaatgccgcttcccaacaccgcaacaacagcgatcgccagtaaaataccccggaaggcaaacagcacaccgcaagtgc                       <  2:795697‑M2/149‑127 (MQ=255)
                                                                                                                             ccGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTGGTGCCAGCGATATGACCGTTTTGTTGCGTCATATCCTGCCTGGGACCGTCTCTTCTATCGTGGTGTTTTTc                       >  1:88221/1‑149 (MQ=255)
                                                                                                                                  tttGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTGGTGCCAGCGATATGACCGTTTTGTTGCGTCATATCCTGCCTGGGACCGTCTCTTCTATCGTGGTGTTTTTCACCAt                  >  1:352422/1‑149 (MQ=255)
                                                                                                                                      cccGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTGGTGCCAGCGATATGACCGTTTTGTTGCGTCATATCCTGCCTGGGCCCGTCTCTTCTATCGTGGTGTTTTTCACCATgcgc              >  2:603739/1‑149 (MQ=255)
                                                                                                                                              gTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTGGTGCCAGCGATATGACCGTTTTGTTGCGTCATATCCTGCCTGGGACCGTCTCTTCTATCGTGGTGTTTTTCACCATGCGCATTGGTAc      <  1:509221/149‑1 (MQ=255)
                                                                                                                                                gcgcGGTAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTGGTGCCAGCGATATGACCGTTTTGTTGCGTCATATCCTGCCTGGGACCGTCTCTTCTATCGTGGTGTTTTTCACCATGCGCATTGGTACCt    <  2:431151/149‑1 (MQ=255)
                                                                                                                                                  gcgGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTGGTGCCAGCGATATGACCGTTTTGTTGCGTCATATCCTGCCTGGGACCGTCTCTTCTATCGTGGTGTTTTTCACCATGCGCATTGGTACCTCg  >  2:748451/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCTGATCATGCGCATTTGCGATGTGCTGTTTGCCTTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTGGTGCCAGCGATATGACCGTTTTGTTGCGTCATATCCTGCCTGGGACCGTCTCTTCTATCGTGGTGTTTTTCACCATGCGCATTGGTACCTCG  >  NC_000913_3_bme_pgi/872294‑872588

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