New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 885072 =22 (0.420)9 (0.170) 5/60 NT 29.0% coding (35/858 nt) ECB_00824 protein similar to retron EC67 DNA adenine methylase
?REL606 3944634 = 22 (0.360)intergenic (+23/‑124) proM/aslB tRNA‑Pro/predicted regulator of arylsulfatase activity

CTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACTTTCCCT  >  REL606/3944599‑3944669
                                   |                                   
ctctcGCCGACCAATTTTGAACCCCGCTTCGGCggg                                     >  1:3613881/1‑36 (MQ=255)
aggaaggtattttttcagttcggacataatggcggg                                     >  1:7832785‑M2/36‑36 (MQ=255)
  ctcGCCGACCAATTTTGAACCCCGCTTCGGCGGGGt                                   >  1:5922069/1‑36 (MQ=255)
  ctcGCCGACCAATTTTGAACCCCGCTTCGGCGGGGt                                   <  1:3639156/36‑1 (MQ=255)
   tcGCCGACCAATTTTGAACCCCGCTTCGGCGGGGtt                                  <  1:5354007/36‑1 (MQ=255)
   aaggtattttttcagttcggacataatggcggGGtt                                  >  1:7409053‑M2/33‑36 (MQ=255)
   aaggtattttttcagttcggacataatggcggGGtt                                  >  1:624798‑M2/33‑36 (MQ=255)
   aaggtattttttcagttcggacataatggcggGAtt                                  >  1:6032381‑M2/33‑36 (MQ=255)
      gtattttttcagttcggacataatggcggGTttttt                               >  1:2597505‑M2/30‑36 (MQ=255)
      gtattttttcagttcggacataatggcggGTttttt                               >  1:4566477‑M2/30‑36 (MQ=255)
      ccGACCAATTTTGAACCCCGCTTCGGCGGGGttttt                               <  1:986829/36‑1 (MQ=255)
       tattttttcagttcggacataatggcggTGtttttt                              >  1:2067614‑M2/29‑36 (MQ=255)
        gACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTg                             >  1:1341728/1‑36 (MQ=255)
          tttttcagttcggacataatggcggGGTTTTTTTtt                           >  1:5484827‑M2/26‑36 (MQ=255)
           cAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGttt                          <  1:2529378/36‑1 (MQ=255)
             aTTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTc                        >  1:2100994/1‑36 (MQ=255)
              ttttGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCt                       <  1:6209174/36‑1 (MQ=255)
              tcagttcggacataatggcggGGTTTTTTTTTTtcc                       >  1:3772690‑M2/22‑35 (MQ=255)
                ttGACCCCCGCTTCGGCGGGGTTTTTTGTTTTCtgt                     <  1:1674004/36‑1 (MQ=255)
                ttGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCtgt                     >  1:3442591/1‑36 (MQ=255)
                 tGAACCCTGCTTCGGCTGGGTTTTTTGTTTTCtgtg                    <  1:1315657/36‑1 (MQ=21)
                 tGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCtgtg                    <  1:3423358/36‑1 (MQ=255)
                  gAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGc                   >  1:2458964/1‑36 (MQ=255)
                  gAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGc                   <  1:7098413/36‑1 (MQ=255)
                   aaCCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCa                  >  1:5863046/1‑36 (MQ=255)
                    aCCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTTCAt                 >  1:1236316/1‑36 (MQ=255)
                      cccGCTTCGGTGGGGTTTTTTGTTTTCTTTGCAttt               >  1:2152492/1‑36 (MQ=25)
                        cGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCg             >  1:4626986/1‑36 (MQ=255)
                         gCTTCGGCGGGGTTTTTTGTTTTTTGTGCAttttat            >  1:1441501/1‑33 (MQ=25)
                         gCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGt            <  1:7083326/36‑1 (MQ=255)
                         gCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGt            <  1:6034705/36‑1 (MQ=255)
                         gCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGt            >  1:1191925/1‑36 (MQ=255)
                          cTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTc           <  1:7506534/36‑1 (MQ=255)
                           ttCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCa          <  1:2397328/36‑1 (MQ=255)
                             cggcggGGTTTTTTGTTTTCTGTGCATTTCGTCACt        <  1:4837421/36‑1 (MQ=255)
                               gcggGGTTTTTTGTTTTCTGTGCATTTCGTCACttt      >  1:1088977‑M2/5‑36 (MQ=255)
                                cggGGTTTTTTGTTTTCTGTGCATTTCTTCACTTTc     >  1:6882411‑M2/4‑36 (MQ=255)
                                cggGGTTTTTTGTTTTCTGTGCATTTCGTCACTTTc     >  1:7702630‑M2/4‑36 (MQ=255)
                                 ggGGTTTTTTGTTTTCTGTGCATTTCGTCACTTTcc    >  1:6203280‑M2/3‑36 (MQ=255)
                                 ggGGTTTTTTGTTTTCTGTGCATTTCGTCACTTTcc    >  1:6144908‑M2/3‑36 (MQ=255)
                                 ggGGTTTTTTGTTTTCTGTGCATTTCGTCACTTTcc    >  1:4390377‑M2/3‑36 (MQ=255)
                                 ggGGTTTTTTGTTTTCTGTGCATTTCGTCACTTTcc    >  1:1509820‑M2/3‑36 (MQ=255)
                                 ggGGTTTTTTGTTTTCTGTGCATTTCGTCACTTTcc    >  1:1337399‑M2/3‑36 (MQ=255)
                                  gGGTTTTTTGTTTTCTGTGCATTTCGTCACTTTccc   <  1:439797‑M2/35‑1 (MQ=255)
                                   ggTTTTTTGTTTTCTGTGCATTTCGTCACTTTCCCt  <  1:423164/36‑1 (MQ=255)
                                   ggTTTTTTGTTTTCTGTGCATTTCGTCACTTTCCCt  >  1:3758408/1‑36 (MQ=255)
                                   |                                   
CTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACTTTCCCT  >  REL606/3944599‑3944669

Base quality scores:  ATCG  < 3 ≤  ATCG  < 11 ≤  ATCG  < 24 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG