New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 2753982 =41 (0.800)4 (0.080) 4/60 NT 9.3% coding (2030/2562 nt) mutS DNA mismatch repair protein
?REL606 = 4141016 37 (0.620)coding (89/705 nt) yijC DNA‑binding transcriptional repressor

CGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAGTACAG  >  REL606/2753947‑2754015
                                   |                                 
cGGGCGCTCAACCTTTATGGTGGAGATGACTGAAAc                                   <  1:7447318/36‑1 (MQ=255)
 gggCGCTCAACCTTTATGGTGGAGATGACTGAAAcc                                  <  1:5057806/36‑1 (MQ=255)
 gggCGCTCAACCTTTATGGTGGAGATGACTGAAAcc                                  >  1:1940300/1‑36 (MQ=255)
 gggCGCTCAACCTTTATGGTGGAGATGACTGAAAcc                                  >  1:3419955/1‑36 (MQ=255)
  ggCGCTCAACCTTTATGGTGGAGATGACTGAAACCg                                 <  1:3147202/36‑1 (MQ=255)
   gcgcTCAACCTTTATGGTGGAGATGACTGAAACCGc                                <  1:4237782/36‑1 (MQ=255)
    cgcTCAACCTTTATGGTGGAGATGACTGAAACCGcc                               >  1:790457/1‑36 (MQ=255)
     gcTCAACCTTTATGGTGGAGATGACTGAAACCGCCa                              >  1:1303119/1‑36 (MQ=255)
     gcTCAACCTTTATGGTGGAGATGACTGAAACCGCCa                              >  1:1456004/1‑36 (MQ=255)
     gcTCAACCTTTATGGTGGAGATGACTGAAACCGCCa                              >  1:2960185/1‑36 (MQ=255)
     gcTCAACCTTTATGGCGGAGATGACTGAAACCGcta                              >  1:6212750/1‑34 (MQ=38)
      cTCAACCTTTATGGTGGAGATGACTGAAACCGCCaa                             >  1:7706465/1‑36 (MQ=255)
       tCAACCTTTATGGTGGAGATGACTGAAACCGCCAat                            >  1:3025634/1‑36 (MQ=255)
         aaCCTTTATGGTGGAGATGACTGAAACCGCCAatat                          >  1:6607149/1‑36 (MQ=255)
         aaCCTTTATGGTGGAGATGACTGAAACCGCCAatat                          >  1:2147924/1‑36 (MQ=255)
         aaCCTTTATGGTGGAGATGACTGAAACCGCCAatat                          >  1:4707733/1‑36 (MQ=255)
          aCCTTTATGGTGGAGATGACTGAAACCGCCAATAtt                         >  1:2624150/1‑36 (MQ=255)
          aCCTTTATGGTGGAGATGACTGAAACCGCCAATAtt                         >  1:90174/1‑36 (MQ=255)
           ccTTTATGGTGGAGATGACTGAAACCGCCAATCttt                        >  1:2662879/1‑36 (MQ=255)
           ccTTTATGGTGGAGATGACTGAAACCGCCAATAttt                        >  1:3210081/1‑36 (MQ=255)
           ccTTTATGGTGGAGATGACTGAAACCGCCAATAttt                        >  1:4012617/1‑36 (MQ=255)
            cTTTATGGTGGAGATGACTGAAACCGCCAATAtttt                       >  1:7318658/1‑36 (MQ=255)
              ttATGGTGGAGATGACTGAAACCGCCAATATTTTAc                     >  1:5752486/1‑36 (MQ=255)
              ttATGGTGGAGATAACTGAAACCGCCAATATTTTAc                     >  1:2670777/1‑36 (MQ=255)
               tATGGTGGAGATGACTGAAACCGCCAATATTTTACa                    >  1:3297280/1‑36 (MQ=255)
               tATGGTGGAGATGACTGAAACCGCCAATATTTTACa                    >  1:5987085/1‑36 (MQ=255)
                aTGGTGGAGATGACTGAAACCGCCAATATTTTACAt                   >  1:1367069/1‑36 (MQ=255)
                aTGGTGGAGATGACTGAAACCGCCAATATTTTACAt                   >  1:2313031/1‑36 (MQ=255)
                  ggtggAGATGACTGAAACCGCCAATATTTTACATaa                 >  1:7268693/1‑36 (MQ=255)
                  ggtggAGATGACTGAAACCGCCAATATTTTACATaa                 >  1:760516/1‑36 (MQ=255)
                  ggtggAGATGACTGAAACCGCCAATATTTTACATaa                 >  1:3486593/1‑36 (MQ=255)
                  ggtggAGATGACTGAAACCGCCAATATTTTACATaa                 >  1:5719946/1‑36 (MQ=255)
                   gtggAGATGACTGAAACCGCCAATATTTTACATAAc                >  1:5216165/1‑36 (MQ=255)
                   gtggAGATGACTGAAACCGCCAATATTTTACATAAc                >  1:7437747/1‑36 (MQ=255)
                     ggAGATGACTGAAACCGCCTATATTTTACATAACGc              >  1:5635470/1‑36 (MQ=255)
                     ggAGATGACTGAAACCGCCAATATTTTACATAACGc              >  1:1071457/1‑36 (MQ=255)
                     ggAGATGACTGAAACCGCCAATATTTTACATAACGc              >  1:2015190/1‑36 (MQ=255)
                       agaTGACTGAAACCGCCAATATTTTACATAACGCCa            >  1:5861136/1‑36 (MQ=255)
                       aaaaaaaaaaacCCGCCAATATTTTACATAACGCCa            >  1:4587770‑M1/13‑36 (MQ=255)
                        gacgACTGAAACCGCCAATATTTTACATAACGCCAc           >  1:7362325/4‑36 (MQ=255)
                        aaaaaaaaaacCCGCCAATATTTTACATAACGCCAc           >  1:5391464‑M1/12‑36 (MQ=255)
                          tGACTGAAACCGCCAATATTTTACATAACGCCACCg         >  1:5093462/1‑36 (MQ=255)
                          tGACTGAAACCGCCAATATTTTACATAACGCCACCg         >  1:137971/1‑36 (MQ=255)
                          tGACTGAAACCGCCAATATTTTACATAACGCCACCg         >  1:5670718/1‑36 (MQ=255)
                           gACTGAAACCGCCAATATTTTACATAACGCCACCGa        >  1:1697876/1‑36 (MQ=255)
                           gACTGAAACCGCAAATATTTTACATAACGCCACCGa        >  1:2796343/1‑36 (MQ=255)
                            aCTGAAACCGCCAATATTTTACATAACGCCACCGAg       >  1:6365479/1‑36 (MQ=255)
                             acaaaaCCGCCAATATTTTACATAACGCCACCGAGt      <  1:6852069/33‑1 (MQ=255)
                               gAAACCGCCAATATTTTACATAACGCCACCGAGTAc    >  1:7511596/1‑36 (MQ=255)
                                aacCCGCCAATATTTTACATAACGCCACCGAGTACa   >  1:5249996‑M1/4‑36 (MQ=255)
                                 acCCGCCAATATTTTACATAACGCCACCGAGTACAg  >  1:7163836‑M1/3‑36 (MQ=255)
                                 aaCCGCCAATATTTTACATAACGCCACCGAGTACAg  <  1:1647287/36‑1 (MQ=255)
                                 aaCCGCCAATATTTTACATAACGCCACCGAGTACAg  >  1:98926/1‑36 (MQ=255)
                                   |                                 
CGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAGTACAG  >  REL606/2753947‑2754015

Base quality scores:  ATCG  < 3 ≤  ATCG  < 13 ≤  ATCG  < 25 ≤  ATCG  < 33 ≤  ATCG  < 37 ≤  ATCG