New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 19776584 (0.110)6 (0.170) 4/56 NT 30.8% coding (182/1005 nt) yedY hypothetical protein
?REL606 3323381 = 23 (0.510)coding (769/3801 nt) yhdP conserved membrane protein, predicted transporter

GGCACCGCCCTGTTTCAGCCAGACGTCACCGCCGGTTACATCGCCTTTGTCGATCGTCATCCAGCCTTCAA  >  REL606/3323346‑3323416
                                   |                                   
ggCACCGCCCTGTTTCAGCCAGACGTCACCGCCGGt                                     >  1:2265608/1‑36 (MQ=255)
  cACGGCCCTGTTACAGCTAGACGTCACCGCCGGTTa                                   <  1:256647/36‑1 (MQ=25)
  cACCCCCCTGTTTCAGCCAGACGTCACCGCCGGTTa                                   <  1:4565002/36‑1 (MQ=255)
     cGCCCTGTTTCAGCCAGACGTCACCGCCGGTTACAt                                <  1:1713421/36‑1 (MQ=255)
      gCCCTTTTTCAGCCGGCCGTCCCCGCCGGTTACATc                               <  1:562745/36‑1 (MQ=19)
      gCCCTGTTTCAGCCAGACGTCACCGCCGGTTACATc                               >  1:3613205/1‑36 (MQ=255)
      gCCCTGTTTCAGCCAGACGTCACCGCCGGGTACATc                               >  1:3417405/1‑36 (MQ=255)
        cttGTTTCAACCAGACGTCACCGCCGGTTACATCGc                             <  1:4554701/34‑1 (MQ=38)
        ccTGTTTCAGCCAGACGTCACCGCCGGTTACATCGc                             >  1:3140744/1‑36 (MQ=255)
           ttttCAGCCAGACGTCACCGCAGGTTACATCGCCtt                          <  1:4905520/35‑1 (MQ=255)
           gTTTCAGCCAGACGTCACCGCCGGTTACATCGCCtt                          >  1:833905/1‑36 (MQ=255)
           gTTTCAGCCAGACGTCACCGCCGGTTACATCGCCtt                          >  1:487910/1‑36 (MQ=255)
            ttacAGCCAGACGTCACCGCCGGTTACATCGCCttt                         <  1:2853848/33‑1 (MQ=255)
               cccccACACGTCCCCGCCGGTTACATCGCCTTTGTc                      <  1:3978884/33‑1 (MQ=25)
               cAGCCAGAAGTCCCCGCCGGTTACATCGCCTTTGTc                      <  1:4054422/36‑1 (MQ=37)
                aGCCAGACGTCACCGCCGGTTACATCGCCTTTGTCg                     >  1:3296433/1‑36 (MQ=255)
                aGCCAGACGTCACCGCCGGTTACATCGCCTTTGTCg                     >  1:5849793/1‑36 (MQ=255)
                 gCCAGACGTCACCGCAGGTTACATCGCCGTTGTTGa                    >  1:561023/1‑36 (MQ=21)
                  ccaaccgcccccgccggTTACATCGCCTTTGTCGAt                   <  1:1515258‑M2/19‑1 (MQ=255)
                   cacccgcccccgccggTTACATCGCCTTTGTCGATc                  <  1:5062073‑M2/20‑1 (MQ=255)
                    agccgcccccgccggTTACATCGCCTTTGTCGATCg                 <  1:2482730‑M2/21‑1 (MQ=255)
                    aGACGTCACCGCCGGTTACATCGCCTTTGTCGagct                 >  1:2653133/1‑33 (MQ=255)
                    aGACGTCACCGCCGGTTACATCGCCTTTGTCGATCg                 >  1:5721974/1‑36 (MQ=255)
                     cccgcccccgccggTTACATCGCCTTTGTCGATCGt                <  1:635301‑M2/22‑1 (MQ=255)
                      ccGTCACCGACGGTTACATCGCCTTTGTCGATCGTc               <  1:4597233/35‑1 (MQ=255)
                       cGTCACCGCCGGTTACATCGCCTTTGTCGATCGTCa              >  1:145967/1‑36 (MQ=255)
                        gTCACCGCCGGTTACATCGCCTTTGTCGATCGTCAt             >  1:2824410/1‑36 (MQ=255)
                         tCACCGCCGGTTACATCGCCTTTGTCGATCGTCATc            >  1:5887443/1‑36 (MQ=255)
                         tCACCGCCGGTTACATCGCCTATGTCGATCGTCCTc            >  1:2407755/1‑36 (MQ=37)
                          ccccgccgtTTACATCGCCTTTGTCGATCGTCATcc           <  1:1254740‑M2/27‑1 (MQ=255)
                          ccccgccggTTACATCGCCTTTGTCGACCGTCATcc           <  1:2509846‑M2/27‑1 (MQ=255)
                           aCCGCCGGTTACATCGCCTTTGTCGATCGTCATCCa          >  1:5259144/1‑36 (MQ=255)
                           aCCGCCGGTTACATCGCCTTTGTCGATCGTCATCCa          >  1:1439929/1‑36 (MQ=255)
                           aCCGCCGGTTACATCGCCTTTGTCGATCGTCATCCa          >  1:3854827/1‑36 (MQ=255)
                            ccgccggTTACATCGCCTTTGTCGATCGTCATCCAg         >  1:483546‑M2/8‑36 (MQ=255)
                             cgccgGTTACATCGCCTTTGTCGAACGTGATCCAGc        >  1:2165087/1‑36 (MQ=37)
                              gccggTTACATCGCCTTTGTCGATCGTCATCCAGcc       >  1:5490226‑M2/6‑36 (MQ=255)
                              gccggTTACATCGCCTTTGTCGATCGTCATCCAGcc       >  1:3349696‑M2/6‑36 (MQ=255)
                               ccggTTACATCGCCTTTGTCGATCGTCATCCAGCCt      >  1:5418833‑M2/5‑36 (MQ=255)
                                 ggTTACATCGCCTTTGTCGATCTTCATCCTGCCttt    >  1:1685921/1‑35 (MQ=25)
                                 ggTTACATCGCCTTTGTCGATCGTCATCCAGCCTTc    >  1:4174415‑M2/3‑36 (MQ=255)
                                 ggTTACATCGCCTTTGTCGATCGTCATCCAGCCTTc    <  1:3805170‑M2/34‑1 (MQ=255)
                                 ggTTACATCGCCTTTGTCGATCGTCATACAGCCTTc    >  1:437239‑M2/3‑36 (MQ=255)
                                   ttACATCGCCTTTGTCGATCGTCATCCAGCCTTAaa  >  1:5830847/1‑36 (MQ=255)
                                   |                                   
GGCACCGCCCTGTTTCAGCCAGACGTCACCGCCGGTTACATCGCCTTTGTCGATCGTCATCCAGCCTTCAA  >  REL606/3323346‑3323416

Base quality scores:  ATCG  < 3 ≤  ATCG  < 5 ≤  ATCG  < 12 ≤  ATCG  < 25 ≤  ATCG  < 37 ≤  ATCG