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

AATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACT  >  REL606/3944593‑3944663
                                   |                                   
tgctgattttgttgatgcgatggcgggtcttgttac                                     <  1:5540883‑M1/1‑1 (MQ=255)
tgctgattttgttgatgcgatggcgggtctggttac                                     >  1:1630300‑M1/36‑36 (MQ=255)
tgctgattttgttgatgcgatggcgggtctggttac                                     <  1:5862650‑M1/1‑1 (MQ=255)
aaTCCTCTCTCGCCGACCAATTTTGAACCCCGCTTc                                     <  1:16276/36‑1 (MQ=255)
aaTCCTCTCTCGCCGACCAATTTTGAACCCCGCTTc                                     >  1:6361671/1‑36 (MQ=255)
aaTCCTCTCTCGCCGACCAATTTTGAACCCCGCTTc                                     <  1:723578/36‑1 (MQ=255)
 gctgattttgttgatgcgatggcgtgtctggttacg                                    <  1:1347015‑M1/2‑1 (MQ=255)
 gctgattttgttgatgcgatggcgggtctggttacg                                    >  1:4500875‑M1/35‑36 (MQ=255)
 aTCCTCTCTCGCCGACCAATTTTGAACCCCGCTTcg                                    >  1:4710496/1‑36 (MQ=255)
  tCCTCTCTCGCCGACCAATTTTGAACCCCGCTTcgg                                   <  1:5629084/36‑1 (MQ=255)
  ctgattttgttgatgcgatggcgggtctggttacgg                                   <  1:5960592‑M1/3‑1 (MQ=255)
   tgattttgttgatgcgatggcgggtctggttatggc                                  <  1:523203‑M1/4‑1 (MQ=255)
   tgattttgttgatgcgatggcgggtctggttacggc                                  >  1:6475518‑M1/33‑36 (MQ=255)
   tgattttgttgatgcgatggcgggtctggttacggc                                  >  1:3222682‑M1/33‑36 (MQ=255)
   tgattttgttgatgcgatggcgggtctggttacggc                                  >  1:642950‑M1/33‑36 (MQ=255)
   ccTCTCTCGCCGACCAATTTTGAACCCCGCTTcggc                                  <  1:6692231/36‑1 (MQ=255)
    gattttgttgatgcgatggcgggtttggttacggcg                                 <  1:230604‑M1/5‑1 (MQ=255)
    gattttgttgatgcgatggcgggtctggttacggcg                                 <  1:4388956‑M1/5‑1 (MQ=255)
    gattttgttgatgcgatggcgggtctggttacggcg                                 >  1:7503151‑M1/32‑36 (MQ=255)
    ctctctCGCCGACCAATTTTGAACCCCGCTTcggcg                                 <  1:5959250/36‑1 (MQ=255)
     tctctcGCCGACCAATTTTGAACCCCGCTTcggcgg                                <  1:4007609/36‑1 (MQ=255)
     tctctcGCCGACCAATTTTGAACCCCGCTTcggcgg                                >  1:6342004/1‑36 (MQ=255)
     attttgttgatgcgatggcgggtctggttacggcgg                                >  1:7436157‑M1/31‑36 (MQ=255)
     attttgttgatgcgatggcgggtctggttacggcgg                                <  1:6053996‑M1/6‑1 (MQ=255)
      ttttgttgatgcgatggcgggtctggttaCGGCggg                               <  1:7121795‑M1/7‑1 (MQ=255)
       tttgttgatgcgatggcgggtctggttaCGGCgggg                              <  1:255326‑M1/8‑1 (MQ=255)
       tctcGCCGACCAATTTTGAACCCCGCTTCGGCgggg                              <  1:4838317/36‑1 (MQ=255)
        ttgttgatgcgatggcgggtctggttaCGgcgggtt                             <  1:6656470‑M1/9‑3 (MQ=255)
        ttgttgatgcgatggcgggtctggttaCGGCGGGGt                             <  1:4005323‑M1/9‑1 (MQ=255)
         tgttgatgcgatggcgggtctggttaCGGCGGGGtt                            <  1:5694489‑M1/10‑1 (MQ=255)
         tcGCCGACCAATTTTGAACCCCGCTTCGGCGTGGtt                            <  1:6053237/36‑1 (MQ=255)
          gttgatgcgatggcgggtctggttaCggcgggggtt                           <  1:3741564‑M1/11‑4 (MQ=255)
           gCCGACCAATTTTGAACCCCGCTTCGGCGGGGtttt                          <  1:4380037/36‑1 (MQ=255)
            tgatgcgatggcgggtctggttaCGGCGGTGttttt                         <  1:609812‑M1/13‑1 (MQ=255)
            ccGACCAATTTTGAACCCCGCTTCGGCGGGGttttt                         <  1:2571178/36‑1 (MQ=255)
            ccGACCAATTTTGAACCCCGCTTCGGCGGGGttttt                         >  1:5953668/1‑36 (MQ=255)
              atgcgatggcgggtctggttaCGGCGGggtttttct                       <  1:7513764‑M1/15‑3 (MQ=255)
              atgcgatggcgggtctggttaCGGCGGGgttttttt                       <  1:1206388‑M1/15‑2 (MQ=255)
               tgcgatggcgggtctggttaCGGCGgggtttttgtt                      <  1:2537912‑M1/16‑4 (MQ=255)
               tgcgatggcgggtctggttaCGGCGGGGTTTTTTGt                      <  1:6108207‑M1/16‑1 (MQ=255)
               tgcgatggcgggtctggttaCGGCGGGGGTTTTtgg                      <  1:5373226‑M1/16‑2 (MQ=255)
               aCCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGt                      <  1:6209124/36‑1 (MQ=255)
                ccAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGtt                     <  1:4852079/36‑1 (MQ=255)
                 cAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGttt                    >  1:1885232/1‑36 (MQ=255)
                  aaTTTTGAACCCCGCTTCGGCGGGGTTTTTTGtttt                   <  1:5329211/36‑1 (MQ=255)
                  aaTTTTGAACCCCGCTTCGGCGGGGTTTTTTGtttt                   >  1:6546714/1‑36 (MQ=255)
                   aTTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTc                  <  1:1925614/36‑1 (MQ=255)
                    ttttGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCt                 <  1:2127596/36‑1 (MQ=255)
                     tttGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCtg                >  1:5280268/1‑36 (MQ=255)
                      ttGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCtgt               >  1:868371/1‑36 (MQ=255)
                      ttGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCtgt               <  1:3569897/36‑1 (MQ=255)
                        gAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGc             >  1:2507447/1‑36 (MQ=255)
                        gAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGc             <  1:104786/36‑1 (MQ=255)
                         aaCCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCa            <  1:6495594/36‑1 (MQ=255)
                          aCCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCAt           <  1:5420968/36‑1 (MQ=255)
                          aCCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCAt           <  1:6761726/36‑1 (MQ=255)
                          aCCCCGCTTCGGCGGGGTTTTTTGTTTTCTCTGCAt           >  1:6839529/1‑36 (MQ=255)
                           ccccGCTTCGGCGGGGTTTTTTGTTTTCTGTGCAtt          >  1:2834290/1‑36 (MQ=255)
                            cccGCTTCGGCGGGGTTTTTTGTTTTCTGTGCAttt         <  1:4450689/36‑1 (MQ=255)
                            cccGCTTCGGCGGGGTTTTTTGTTTTCTGTGCAttt         <  1:7152762/36‑1 (MQ=255)
                             ccGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTc        <  1:3472315/36‑1 (MQ=255)
                               gCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGt      >  1:427206/1‑36 (MQ=255)
                                cTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTc     <  1:2588267/36‑1 (MQ=255)
                                 ttCGGCGGGGTTTTTTGTTTTTTTTTCATTTTGTCa    >  1:5608089/1‑36 (MQ=16)
                                 ttCGGCGGGGTTTTTTGTTTTCTGTTCATTTCGTCa    >  1:3445305/1‑36 (MQ=255)
                                   CGGCGGGGTTTTTTGTTTTCTTTGCATTTCGTCACt  <  1:1310374‑M1/36‑1 (MQ=255)
                                   |                                   
AATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACT  >  REL606/3944593‑3944663

Base quality scores:  ATCG  < 3 ≤  ATCG  < 11 ≤  ATCG  < 22 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG