New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_bme_pgi 1468596 =NA (NA)71 (0.390) 15/184 16.8 43.3% coding (231/906 nt) insD1 IS2 transposase TnpB
?NC_000913_3_bme_pgi = 1650887 99 (0.520)pseudogene (43/693 nt) insD1 IS2 transposase TnpB;IS, phage, Tn; Transposon‑related functions; extrachromosomal; transposon related
Rejected: Coverage evenness skew score above cutoff.

Only 100 of 101 total aligned reads displayed.

TGCCATTCATCGCGAATACAAGCAGACATGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCGGATCATGCGCCAGAATGCGCTGTTGCTTG  >  NC_000913_3_bme_pgi/1650788‑1650941
                                                                                                   |                                                      
tgcCATTCATCGCGAATACAAGCAGACATGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTg                                                        >  1:1856012/1‑100 (MQ=255)
 gcCATTCATCGCGAATACAAGCAGACATGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTg                                                        >  2:206622/1‑99 (MQ=255)
 gcCATTCATCGCGAATACAAGCAGACATGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGc                                                       >  1:3071158/1‑100 (MQ=255)
 gcCATTCATCGCGAATACAAGCAGACATGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGc                                                       >  1:4208414/1‑100 (MQ=255)
  ccATTCATCGCGAATACAAGCAGACATGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGc                                                       <  1:3567895/99‑1 (MQ=255)
   cATTCATCGCGAATACAAGCAGACATGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGc                                                       >  1:2613531/1‑98 (MQ=255)
   cATTCATCGCGAATACAAGCAGACATGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGa                                                     >  1:2468000/1‑100 (MQ=255)
   cATTCATCGCGAATACAAGCAGACATGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGa                                                     >  1:2368031/1‑100 (MQ=255)
   cATTCATCGCGAATACAAGCAGACATGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGa                                                     >  1:2613356/1‑100 (MQ=255)
    aTTCATCGCGAATACAAGCAGACATGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGAt                                                    <  1:992389/100‑1 (MQ=255)
    aTTCATCGCGAATACAAGCAGACATGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGAt                                                    <  2:2486187/100‑1 (MQ=255)
    aTTCATCGCGAATACAAGCAGACATGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGAt                                                    >  2:755082/1‑100 (MQ=255)
    aTTCATCGCGAATACAAGCAGACATGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGAt                                                    >  1:585764/1‑100 (MQ=255)
     ttCATCGCGAATACAAGCAGACATGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGAt                                                    >  1:347034/1‑99 (MQ=255)
     ttCATCGCGAATACAAGCAGACATGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGAt                                                    >  2:4096650/1‑99 (MQ=255)
     ttCATCGCGAATACAAGCAGACATGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATc                                                   <  2:2541847/100‑1 (MQ=255)
     tgcATCGCGAATACAAGCAGACATGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATc                                                   <  2:4215403/98‑1 (MQ=255)
       cATCGCGAATACAAGCAGACATGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATc                                                   <  2:4141255/98‑1 (MQ=255)
        aTCGCGAATACAAGCAGACATGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAAt                                                <  2:4325540/100‑1 (MQ=255)
        aTCGCGAATACAAGCAGACATGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAAt                                                >  2:2397291/1‑100 (MQ=255)
        aTCGCGAATACAAGCAGACATGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAAt                                                >  1:2886810/1‑100 (MQ=255)
         tCGCGAATACAAGCAGACATGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATg                                               <  2:2298015/100‑1 (MQ=255)
         tCGCGAATACAAGCAGACATGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATg                                               <  1:480557/100‑1 (MQ=255)
         tCGCGAATACAAGCAGACATGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATg                                               <  1:1811999/100‑1 (MQ=255)
          cgcgAATACAAGCAGACATGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGc                                              >  1:2333889/1‑100 (MQ=255)
          cgcgAATACAAGCAGACATGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGc                                              >  1:3334428/1‑100 (MQ=255)
           gcgAATACAAGCAGACATGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGcc                                             <  1:3083810/100‑1 (MQ=255)
            cgAATACAAGCAGACATGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGcc                                             >  1:445189/1‑99 (MQ=255)
                tACAAGCAGACATGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACg                                        <  1:3092927/100‑1 (MQ=255)
                tACAAGCAGACATGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACg                                        <  1:2525902/100‑1 (MQ=255)
                  cAAGCAGACATGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCCAACgtg                                      >  1:3738257/1‑100 (MQ=255)
                  cAAGCAGACATGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACgtg                                      >  1:3652968/1‑100 (MQ=255)
                  cAAGCAGACATGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACgtg                                      >  2:3027559/1‑100 (MQ=255)
                  cAAGCAGACATGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACgt                                       >  2:2326028/1‑99 (MQ=255)
                    aGCAGACATGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGtt                                    >  2:808576/1‑100 (MQ=255)
                    aGCAGACATGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGtt                                    >  1:3445581/1‑100 (MQ=255)
                    aGCAGACATGGAAAAAACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGACCTTGATGGTACGCCAGCGATCAATGCCAAAc                                         >  2:4691874/1‑95 (MQ=21)
                     gCAGAGATGGAAAGCACTGGATGACGAACTGGCCTATCGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGttt                                   <  1:3985105/100‑1 (MQ=255)
                     gCAGACATGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGtt                                    >  1:3298030/1‑99 (MQ=255)
                      cAGACATGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTa                                  <  2:237285/100‑1 (MQ=255)
                      cAGACATGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTa                                  >  1:1276537/1‑100 (MQ=255)
                          cATGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCgg                              <  1:140688/100‑1 (MQ=255)
                          cATGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCg                               >  2:1489107/1‑99 (MQ=255)
                          cATGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACACGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCg                               >  2:3219220/1‑99 (MQ=255)
                          cATGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCCTCGCCGACAGGCAGAACTTGATGGTATGCCTGCGATCACTGCCAAACGTGTTTACCg                               >  2:3784251/1‑99 (MQ=255)
                           aTGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCGGa                             >  1:4109937/1‑100 (MQ=255)
                            tGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCGGaa                            >  1:1701439/1‑99 (MQ=255)
                            tGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCGGAt                            >  1:3405406/1‑100 (MQ=255)
                            tGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCGGAt                            <  1:553739/100‑1 (MQ=255)
                            tGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCGGAt                            >  1:1701404/1‑100 (MQ=255)
                             ggAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCGGAt                            >  1:4519598/1‑99 (MQ=255)
                              gAAAACGCTGGATGACGAACTGGCCTAGGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCGGATc                           <  2:1978175/99‑1 (MQ=255)
                              gAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCGGATc                           <  1:2661689/99‑1 (MQ=255)
                              gAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCGGATc                           <  1:722375/99‑1 (MQ=255)
                              gAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCGGATc                           <  1:721208/99‑1 (MQ=255)
                               aaaaCACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCGGATCAt                         >  1:3191751/1‑100 (MQ=255)
                               aaaaCACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCGGATCAt                         >  1:3387159/1‑100 (MQ=255)
                               aaaaCACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCGGATCAt                         >  1:4443396/1‑100 (MQ=255)
                               aaaaCACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCGGATCAt                         >  1:1396743/1‑100 (MQ=255)
                                aaaCACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCGGATCATg                        <  1:2325926/100‑1 (MQ=255)
                                 aaCACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCGGATCATgc                       >  2:120819/1‑100 (MQ=255)
                                  acacTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCGGATCATgcg                      <  2:4500538/100‑1 (MQ=255)
                                    acTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCGGATCATGCGcc                    <  1:1774581/100‑1 (MQ=255)
                                      tGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCGGATCATgcgc                     >  2:2677541/1‑97 (MQ=255)
                                      tGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCGGATCATgcgc                     >  1:2827869/1‑97 (MQ=255)
                                      tGGATGACGAACTGGCCCATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCCGCGATCAATGCCAAACGTGTTTACCGGATCATccgc                     >  2:2407303/1‑97 (MQ=255)
                                       ggATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCGGATCATGCGCCAg                  >  2:4382211/1‑99 (MQ=255)
                                       ggATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCGGATCATGCGCCAg                  >  2:2827085/1‑99 (MQ=255)
                                        gATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCGGATCATGCGCCAg                  <  2:2627542/98‑1 (MQ=255)
                                        gATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCGGATCATGCGCCAg                  <  1:4072893/98‑1 (MQ=255)
                                        gATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCGGATCATGCGCCAGaa                <  1:3318920/100‑1 (MQ=255)
                                        gATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCGGATCATGCGCCAGaa                <  2:3726914/100‑1 (MQ=255)
                                        gATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAACGCCAAACGTGTTTACCGGATCATGCGCCAg                  <  1:4368815/98‑1 (MQ=255)
                                          tGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTAGGCCTGCGATCAATGCCAAACGTGTTTACCGGATCATGCGCCAGAATg              <  1:3672246/100‑1 (MQ=255)
                                          tGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTAAGCCTGCGATCAATGCCAAACGTGTTTACCGGATCATGCGCCAGAATg              <  1:3672432/100‑1 (MQ=255)
                                            aCGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCGGATCATGCGCCAGAATgcg            >  2:2492837/1‑100 (MQ=255)
                                             cGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCGGATCATGCGCCAGAATgcgn           <  2:1730158/100‑2 (MQ=255)
                                             cGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCGGATCATGCGCCAGAATgcgc           <  1:2185023/100‑1 (MQ=255)
                                             cGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCGGATCATGCGCCAGAATgcgc           <  1:98957/100‑1 (MQ=255)
                                             cGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCGGATCATGCGCCAGAATgcgc           >  1:466288/1‑100 (MQ=255)
                                              gAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCGGATCATGCGCCAGAATGCGCt          <  1:118081/100‑1 (MQ=40)
                                              gAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCGGATCATGCGCCAGAATGCGCt          <  1:1752079/100‑1 (MQ=40)
                                               aaCTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCGGATCATGCGCCAGAATGCGCTg         >  1:2935514/1‑100 (MQ=39)
                                                aCTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCGGATCATGCGCCAGAATGCGCTGt        <  2:3506813/100‑1 (MQ=38)
                                                 cTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCGGATCATGCGCCAGAATGCGCTGtt       >  1:2575640/1‑100 (MQ=37)
                                                   ggCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCGGATCATGCGCCCGAATGCGCTGTTGc     >  2:3716924/1‑100 (MQ=21)
                                                   ggCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCGGATCATGCGCCAGAATGCGCTGTTGc     >  1:4222185/1‑100 (MQ=34)
                                                   ggCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCGGATCATGCGCCAGAATGCGCTGTTGc     >  1:1104274/1‑100 (MQ=34)
                                                   ggCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCGGATCATGCGCCAGAATGCGCTGTTGc     >  1:4174016/1‑100 (MQ=34)
                                                   ggCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCGGATCATGCGCCAGAATGCGCTGTTGc     >  1:4145202/1‑100 (MQ=34)
                                                   ggCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCGGATCATGCGCCAGAATGCGCTGTTGc     >  1:160652/1‑100 (MQ=34)
                                                   ggCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCGGATCATGCGCCAGAATGCGCTGTTGc     >  2:656500/1‑100 (MQ=34)
                                                    gCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCGGATCATGCGCCAGAATGCGCt          <  2:961009/94‑1 (MQ=33)
                                                    gCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCGGATCATGCGCCAGAATGCGCTGTTGCt    >  2:2956931/1‑100 (MQ=33)
                                                     ccTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCGGATCATGCGCCAGAATGCGCTGTTGCtt   >  2:3097595/1‑100 (MQ=32)
                                                     ccTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCGGATCATGCGCCAGAATGCGCTGTTGCtt   >  2:458631/1‑100 (MQ=32)
                                                     ccTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCGGATCATGCGCCAGAATGCGCTGTTGCtt   >  2:3098463/1‑100 (MQ=32)
                                                      cTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCGGATCATGCGCCAGAATGCGCTGTTGCTTg  <  1:3059298/100‑1 (MQ=31)
                                                      cTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCGGATCATGCGCCAGAATGCGCTGTTGCTTg  <  2:4240754/100‑1 (MQ=31)
                                                      cTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCGGATCATGCGCCAGAATGCGCTGTTGCTTg  <  1:4212869/100‑1 (MQ=31)
                                                                                                   |                                                      
TGCCATTCATCGCGAATACAAGCAGACATGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCGGATCATGCGCCAGAATGCGCTGTTGCTTG  >  NC_000913_3_bme_pgi/1650788‑1650941

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