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

AGCAGGAAGGTATTTTTTCAGTTCGGACATAATGGCGG...................................  <  REL606/885109‑885072
.................................GGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACTTTCCC  >  REL606/3944634‑3944668
                                                                         
   AGGAAGGTATTTTTTCAGTTCGGACATAATGGCGGG                                    >  1:7832785/1‑36 (MQ=255)
      AAGGTATTTTTTCAGTTCGGACATAATGGCGGGGTT                                 >  1:624798/1‑36 (MQ=255)
      AAGGTATTTTTTCAGTTCGGACATAATGGCGGGGTT                                 >  1:7409053/1‑36 (MQ=255)
      AAGGTATTTTTTCAGTTCGGACATAATGGCGGGATT                                 >  1:6032381/1‑36 (MQ=255)
         GTATTTTTTCAGTTCGGACATAATGGCGGGTTTTTT                              >  1:2597505/1‑36 (MQ=255)
         GTATTTTTTCAGTTCGGACATAATGGCGGGTTTTTT                              >  1:4566477/1‑36 (MQ=255)
          TATTTTTTCAGTTCGGACATAATGGCGGTGTTTTTT                             >  1:2067614/1‑36 (MQ=255)
             TTTTTCAGTTCGGACATAATGGCGGGGTTTTTTTTT                          >  1:5484827/1‑36 (MQ=255)
                 TCAGTTCGGACATAATGGCGGGGTTTTTTTTTTTCc                      >  1:3772690/1‑35 (MQ=255)
                                                                         
AGCAGGAAGGTATTTTTTCAGTTCGGACATAATGGCGG...................................  <  REL606/885109‑885072
.................................GGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACTTTCCC  >  REL606/3944634‑3944668

Base quality scores:  ATCG  < 3 ≤  ATCG  < 6 ≤  ATCG  < 17 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG