New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_bme_pgi = 402740563 (1.470)6 (0.150) 5/270 5.2 9.1% noncoding (84/281 nt) REP291 (repetitive extragenic palindromic) element; contains 6 REP sequences REP291 (repetitive extragenic palindromic) element; contains 6 REP sequences
?NC_000913_3_bme_pgi = 4027473 60 (1.480)noncoding (152/281 nt) REP291 (repetitive extragenic palindromic) element; contains 6 REP sequences REP291 (repetitive extragenic palindromic) element; contains 6 REP sequences
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGC  >  NC_000913_3_bme_pgi/4027258‑4027549
                                                                                                                                                   |                                                                                                                                                 
ggAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCa                                                                                                                                                  >  2:293370/1‑149 (MQ=255)
    gaaCGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTCCAAATAGTCCGAACCGTCGGCCGGATAAGGCGTTTCCGCCGCATCCGACAAATg                                                                                                                                              >  1:169848/1‑149 (MQ=25)
    gaaCGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTg                                                                                                                                              <  1:671008/149‑1 (MQ=255)
        ggCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGc                                                                                                                                          >  1:537686/1‑149 (MQ=255)
          gaaCTAAGCTGTTTTTCAGTGGCAGTTAGATCTATGACGTATCTGGCAAACCAGGCCCGATGCGACGCTGTCGCGTCTTATCGGGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTAGGCCGCATCCGGCAATTGGTGCat                                                                                                                                        <  2:285154/148‑1 (MQ=18)
          cacgtAAACTGTTTTCGAGTGGCATGTAGATCTATGACCTATCTGGTAAACCATGCCCGATGCGACGCTGCCGCGTCCTATCGTGCCTACAAATAGTCCCGACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCat                                                                                                                                        <  2:229012/145‑1 (MQ=17)
                             tGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTg                                                                                                                     >  1:474258/1‑149 (MQ=255)
                             tGGCATTTAGATCTATGACGTATATAGCAGACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTg                                                                                                                     >  1:738578/1‑149 (MQ=255)
                              ggCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTgg                                                                                                                    >  2:382744/1‑148 (MQ=255)
                               gCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACTCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGCCGCTGTc                                                                                                                   >  2:604367/1‑149 (MQ=255)
                                       aTCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCTCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTa                                                                                                           >  1:658760/1‑149 (MQ=255)
                                       aTCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTa                                                                                                           >  2:62730/1‑149 (MQ=255)
                                          tATGACGTATCTGGCAAACCATGCCCGATGCGAAGCTGTCGCGTCTTATCGTGCCTACAGATAGTCCGAACCGTCGGCAGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCGGTTCAGAGCCTGTCGCGTCGTTTCg                                                                                                        >  1:356259/1‑149 (MQ=17)
                                              aCGTATCTGGCAAACCCTGCCCGAAGCGACGCGGTCGCGTCTTCTCGCTGATACTAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGcc                                                                                                    <  1:545978/149‑1 (MQ=17)
                                                      ggCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTGCAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTa                                                                                                  <  1:331174/143‑1 (MQ=255)
                                                      ggCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTa                                                                                                  >  2:331174/1‑143 (MQ=255)
                                                      ggCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATa                                                                                            <  1:744638/149‑1 (MQ=255)
                                                       gCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAg                                                                                           <  1:670962/149‑1 (MQ=255)
                                                        cAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGt                                                                                          <  1:291985/149‑1 (MQ=255)
                                                         aaaCCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTc                                                                                         >  2:494537/1‑149 (MQ=255)
                                                         aaaCCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGAATAAGGCTGTTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTc                                                                                         >  1:549107/1‑149 (MQ=17)
                                                          aaCCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTcc                                                                                        <  1:496771/149‑1 (MQ=255)
                                                           aCCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCg                                                                                       >  1:341769/1‑149 (MQ=255)
                                                             cATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGAAGCTGTCGCGGCTTATCGTGCCTACAAATAGTCCGaa                                                                                     <  2:537686/149‑1 (MQ=255)
                                                               tGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCCTCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTCTCGTGCCTACAAATAGCCCGGAcc                                                                                  >  2:191539/1‑149 (MQ=255)
                                                                  ccGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTc                                                                                                              >  1:347194/1‑119 (MQ=35)
                                                                       gCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCg                                                                          <  1:418152/149‑1 (MQ=255)
                                                                        cGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTCGTGCAGGATGCCTGATGCGACGCTGTCGCCTCTTATCGTGCCTCCAAAGGCTCCGAACACGTAGTCCg                                                                          >  2:660828/1‑149 (MQ=14)
                                                                          aCGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTCATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGAt                                                                       >  1:516196/1‑149 (MQ=255)
                                                                              tGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCGTGCCAACAAAGAGTCCGACCCGTCGGCCCGATAAgg                                                                   >  2:270154/1‑149 (MQ=14)
                                                                                  gcgTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGtt                                                               >  1:212045/1‑149 (MQ=255)
                                                                                        tATCGTGCCTACAAATAGTCCGAACCGTAGGTCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTAcgcc                                                         <  2:622803/149‑1 (MQ=255)
                                                                                                      aTAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTg                                           <  1:494537/149‑1 (MQ=255)
                                                                                                        aGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGt                                         <  1:196509/149‑1 (MQ=255)
                                                                                                        aGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGt                                         <  1:604072/149‑1 (MQ=255)
                                                                                                          tCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGGACCGTAGGCCGGGTAAGGCGTTTACGCCGCATCCGGCGATTGTTGc                                       >  1:475697/1‑149 (MQ=255)
                                                                                                             gAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCatg                                    >  1:506899/1‑149 (MQ=255)
                                                                                                             gAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCaga                                    >  2:206992/1‑147 (MQ=255)
                                                                                                                  gTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGcc                               >  2:313196/1‑149 (MQ=255)
                                                                                                                   tAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCt                              <  1:432085/149‑1 (MQ=255)
                                                                                                                    aGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTg                             <  1:313196/149‑1 (MQ=255)
                                                                                                                     ggCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGa                            >  2:513311/1‑149 (MQ=255)
                                                                                                                      gCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGAt                           >  1:647537/1‑149 (MQ=255)
                                                                                                                      gCCGGATAAGGCGTTTACGCCGCATCCCGAAAAAGGGGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACACATAGTCCGAACCGTAGGCCGGGTAAGGCGGTTACGCCGCATCCGGCAATTAGTGCATGCTGCCTGAt                           >  1:738445/1‑149 (MQ=17)
                                                                                                                       ccGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATg                          >  2:137195/1‑149 (MQ=255)
                                                                                                                             aaGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGGCTATTGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGc                    <  1:581172/149‑1 (MQ=255)
                                                                                                                                gCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTCGCGCTCCCTATGTTTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGAAGCTgg                 >  1:362353/1‑149 (MQ=18)
                                                                                                                                   tttACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGgcgc              <  2:425563/149‑1 (MQ=255)
                                                                                                                                   tttACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTAAAAATAGTCCGAACCGTAGGCCGGGTAAGCGGTTTAAGCCACATCCCGCCATTGGTGCATGAAGCCTGATGCGACGCTGgcgc              >  2:205/1‑149 (MQ=18)
                                                                                                                                     tACGCCCTATGCCGCATTTGGTGAACCATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGt            <  2:474258/149‑1 (MQ=11)
                                                                                                                                      aCGCCCCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTc           <  1:604367/149‑1 (MQ=255)
                                                                                                                                        gccgcATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCtt         >  1:391160/1‑149 (MQ=255)
                                                                                                                                        gccgcATCCGGCAAGTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCtt         <  2:647537/149‑1 (MQ=255)
                                                                                                                                         ccgcATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTa        >  1:326826/1‑149 (MQ=255)
                                                                                                                                         ccgcATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTa        <  2:212045/149‑1 (MQ=255)
                                                                                                                                         ccgcATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTa        >  1:659609/1‑149 (MQ=255)
                                                                                                                                         ccgcATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAAGCGTTTCCGCCGCATCCCGCAATTGGGGCAAGATGCCTGGTGCGGCGCAGgcgccgcgag        >  2:712053/1‑143 (MQ=17)
                                                                                                                                          cgcATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATAGTGCCAGATGCGACGCGGGCGCCCCTTac       >  2:402175/1‑148 (MQ=17)
                                                                                                                                             aTCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAt                                             >  2:511247/1‑108 (MQ=255)
                                                                                                                                             aTCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAt                                             <  1:511247/108‑1 (MQ=255)
                                                                                                                                               ccGGCAATTGGTGCATGACGCCTGATGCGACGCTGTCGCGTCTTATCGTGCATACACCGCTTCCTAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGc  <  1:174626/149‑1 (MQ=17)
                                                                                                                                                  gCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTc                                                                                                              <  2:333400/39‑1 (MQ=255)
                                                                                                                                                  gCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTc                                                                                                              >  1:333400/1‑39 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGC  >  NC_000913_3_bme_pgi/4027258‑4027549

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