New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 3944628 =29 (0.610)12 (0.250) 6/58 NT 25.8% intergenic (+17/‑130) proM/aslB tRNA‑Pro/predicted regulator of arylsulfatase activity
?REL606 = 4258591 40 (0.700)coding (1180/1350 nt) yjcD predicted permease

TTGCTGATTTTGTTGATGCGATGGCGGGTCTGGTTACGGCGGTATTCATCGTGCTGACCTGTAACATCGTA  >  REL606/4258556‑4258626
                                   |                                   
ttGCTGATTTTGTTGATGCGATGGCGGGTCTGGTTa                                     <  1:2385716/36‑1 (MQ=255)
ttGCTGATTTTGTTGATGCGATGGATGGTTTGGTTa                                     >  1:5786238/1‑36 (MQ=25)
 tGCTGATTTTGTTGATGCGATGGCGGGTCTTGTTAc                                    <  1:5540883‑M2/36‑2 (MQ=255)
 tGCTGATTTTGTTGATGCGATGGCGGGTCTGGTTAc                                    <  1:5862650‑M2/36‑2 (MQ=255)
 tGCTGATTTTGTTGATGCGATGGCGGGTCTGGTTAc                                    >  1:1630300‑M2/1‑35 (MQ=255)
  gCTGATTTTGTTGATGCGATGGCGTGTCTGGTTAcg                                   <  1:1347015‑M2/36‑3 (MQ=255)
  gCTGATTTTGTTGATGCGATGGCGGGTCTGGTTAcg                                   >  1:4500875‑M2/1‑34 (MQ=255)
   cTGATTTTGTTGATGCGATGGCGGGTCTGGTTAcgg                                  <  1:5960592‑M2/36‑4 (MQ=255)
   cTGATTTTGTGGATTGGATGGCGGGTCTGGTTAcgg                                  <  1:4022175/36‑1 (MQ=37)
    tGATTTTGTTGATGCGATGGCGGGTCTGGTTAtggc                                 <  1:523203‑M2/36‑5 (MQ=255)
    tGATTTTGTTGATGCGATGGCGGGTCTGGTTAcggc                                 >  1:3222682‑M2/1‑32 (MQ=255)
    tGATTTTGTTGATGCGATGGCGGGTCTGGTTAcggc                                 >  1:642950‑M2/1‑32 (MQ=255)
    tGATTTTGTTGATGCGATGGCGGGTCTGGTTAcggc                                 >  1:6475518‑M2/1‑32 (MQ=255)
     gATTTTGTTGATGCGATGGCGGGTTTGGTTAcggcg                                <  1:230604‑M2/36‑6 (MQ=255)
     gATTTTGTTGATGCGATGGCGGGTCTGGTTAcggcg                                >  1:7503151‑M2/1‑31 (MQ=255)
     gATTTTGTTGATGCGATGGCGGGTCTGGTTAcggcg                                <  1:4388956‑M2/36‑6 (MQ=255)
      aTTTTGTTGATGCGATGGCGGGTCTGGTTAcggcgg                               <  1:6053996‑M2/36‑7 (MQ=255)
      aTTTTGTTGATGCGATGGCGGGTCTGGTTAcggcgg                               >  1:7436157‑M2/1‑30 (MQ=255)
       ttttGTTTATGAGATGGCGGGTCTGGTTACGGCGGt                              <  1:4450654/36‑1 (MQ=38)
       ttttGTTGATGCGATGGCGGGTCTGGTTAcggcggg                              <  1:7121795‑M2/36‑8 (MQ=255)
       ttttGTTGATGCGATGGCGGGTCTGGTTACGGCGGt                              <  1:1453641/36‑1 (MQ=255)
        tttGTTGATGCGATGGCGGGTCTGGTTAcggcgggg                             <  1:255326‑M2/36‑9 (MQ=255)
        tttGTTGATGCGATGGCGGGTCTGGTTACGGCGGTa                             <  1:3358331/36‑1 (MQ=255)
         ttgttgATGCGATGGCGGGTCTGGTTAcggcgggtt                            <  1:6656470‑M2/36‑10 (MQ=255)
         ttgttgATGCGATGGCGGGTCTGGTTAcggcggggt                            <  1:4005323‑M2/36‑10 (MQ=255)
         ttgttgATGCGATGGCGGGTCTGGTTACGGCGGTAt                            <  1:6605459/36‑1 (MQ=255)
         ttgttgATGCGATGGCGGGTCTGGTTACGGCGGTAt                            <  1:1977787/36‑1 (MQ=255)
          tgttgATGCGATGGCGGGTCTGGTTAcggcggggtt                           <  1:5694489‑M2/36‑11 (MQ=255)
           gttgATGCGATGGCGGGTCTGGTTAcggcgggggtt                          <  1:3741564‑M2/36‑12 (MQ=255)
           gttgATGCGATGGCGGGTCTGGTTACGGCGGTATTc                          <  1:356706/36‑1 (MQ=255)
           gttgATGCGATGGCGGGTCTGGTTACGGCGGTATTc                          <  1:36809/36‑1 (MQ=255)
            ttgATGCGATGGCGAGTCTGGTTACGGCGGTATTCa                         <  1:806851/36‑1 (MQ=255)
             tgATGCGATGGCGGGTCTGGTTAcggcggtgttttt                        <  1:609812‑M2/36‑14 (MQ=255)
             tgATGCGATGGCGGGTCTGGTTACGGCGGTATTCAt                        <  1:2649325/36‑1 (MQ=255)
              gATGGGATGGAGGCTCTGGTTACCGCGGTATTCATc                       <  1:29225/36‑1 (MQ=19)
              gATGCGATGGCGGGTCTGGTTACGGCGGTATTCATc                       <  1:5843945/36‑1 (MQ=255)
              gATGCGATGGCGGGTCTGGTTACGGCGGTATTCATc                       <  1:2876404/36‑1 (MQ=255)
               aTGCGATGGCGGGTCTGGTTAcggcggggttttttt                      <  1:1206388‑M2/36‑16 (MQ=255)
               aTGCGATGGCGGGTCTGGTTAcggcggggtttttct                      <  1:7513764‑M2/36‑16 (MQ=255)
               aTGCGATGGCGGGTCTGGTTACGGCGGTATTCATCg                      <  1:1596457/36‑1 (MQ=255)
                tGCGATGGCGGGTCTGGTTAcggcggggttttttgt                     <  1:6108207‑M2/36‑17 (MQ=255)
                tGCGATGGCGGGTCTGGTTAcggcggggtttttgtt                     <  1:2537912‑M2/36‑17 (MQ=255)
                tGCGATGGCGGGTCTGGTTAcggcgggggtttttgg                     <  1:5373226‑M2/36‑17 (MQ=255)
                tGCGATGGCGGGTCTGGTTACGGCGGTATTCATCGt                     <  1:2280541/36‑1 (MQ=255)
                tGCGATGGCGGGTCTGGTTACGGCGGTATTCATCGt                     <  1:4943377/36‑1 (MQ=255)
                tGCGATGGCGGGTCTGGTTACGGCGGTATTCATCGt                     <  1:5647589/36‑1 (MQ=255)
                  cGATGGCGTGTCTGGTTACGGCGGTATTCATCGTGc                   <  1:1793594/36‑1 (MQ=255)
                  cGATGGCGGGTCTGGTTACGGCGGTATTCATCGTGc                   <  1:1135025/36‑1 (MQ=255)
                  cGATGGCGGGTCTGGTTACGGCGGTATTCATCGTGc                   <  1:3737951/36‑1 (MQ=255)
                  cGATGGCGGGTCTGGTTACGGCGGTATTCATCGTGc                   <  1:4949669/36‑1 (MQ=255)
                   gATGGCGGGTCTGGTTACGGCGGTATTCATCGTGCt                  <  1:4632790/36‑1 (MQ=255)
                    aTGGCGGGTCTGGTTACGGCGGTATTCATCGTGCTg                 <  1:1240991/36‑1 (MQ=255)
                      ggCGGGTCTGGTTACGGCGGTATTCATCGTGCTGAc               <  1:3141649/36‑1 (MQ=255)
                       gCGGGTCTGGTTACGGCGGTATTCATCGTGCTGAcc              <  1:5542860/36‑1 (MQ=255)
                       gCGGGTCTGGTTACGGCGGTATTCATCGTGCTGAcc              <  1:5512903/36‑1 (MQ=255)
                       gCGGGTCTGGTTACGGCGGTATTCATCGTGCTGAcc              <  1:4059565/36‑1 (MQ=255)
                       gCGGGTCTGGTTACGGCGGTATTCATCGTGCTGAcc              <  1:5992407/36‑1 (MQ=255)
                       ccGGGTCTGGTTACGGCTTTATTCATCGTGCTGAcc              <  1:4803523/35‑1 (MQ=37)
                        cGGGTCTGGTTACGGCGGTATTCATCGTGCTGACCt             <  1:4215230/36‑1 (MQ=255)
                        cGGGTCTGGTTACGGCGGTATTCATCGTGCTGACCt             <  1:7071069/36‑1 (MQ=255)
                        cGGGTCTGGTTACGGCGGTATTCATCGTGCTGACCt             <  1:5581079/36‑1 (MQ=255)
                          ggTCTGGTTACGGCGGTATTCATCGTGCTGACCTGt           <  1:479782/36‑1 (MQ=255)
                          ggTCTGGTTACGGCGGTATTCATCGTGCTGACCTGt           <  1:6683776/36‑1 (MQ=255)
                           gTCTGGTTACGGCGGTATTCATCGTGCTGACCTGTa          <  1:2307473/36‑1 (MQ=255)
                             cTGGTTACGGCGGTATTCATCGTGCTGACCTGTAAc        <  1:3976608/36‑1 (MQ=255)
                              tGGTTACGGCGGTATTCATCGTGCTGACCTGTAACa       <  1:4968560/36‑1 (MQ=255)
                                gTTACGGCGGTATTCATCGTGCTGACCTGTAACATc     <  1:2508439/36‑1 (MQ=255)
                                 ttACGGCGGTATTCATCGTGCTGACCTGTAACATCg    <  1:2469444/36‑1 (MQ=255)
                                   aCGGCGGTATTCATCGTGCTGACCTGTAACATCGTa  <  1:7278149/36‑1 (MQ=255)
                                   aCGGCGGTATTCATCGTGCTGACCTGTAACATCGTa  <  1:3596193/36‑1 (MQ=255)
                                   |                                   
TTGCTGATTTTGTTGATGCGATGGCGGGTCTGGTTACGGCGGTATTCATCGTGCTGACCTGTAACATCGTA  >  REL606/4258556‑4258626

Base quality scores:  ATCG  < 3 ≤  ATCG  < 9 ≤  ATCG  < 19 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG