New junction evidence
  seq id position reads (cov) reads (cov) score skew annotation gene product
* ? REL606 72407 =24 (0.510)6 (0.130) 4/60 5.3 coding (446/1701 nt) araB ribulokinase
?REL606 3944634 = 38 (0.690)intergenic (+23/‑124) proM/aslB tRNA‑Pro/predicted regulator of arylsulfatase activity
Rejected: Position hash score below cutoff.

ACGCGCCGGGCAATGTTGACTACTCCCGCTATATTGGCGG...................................  <  REL606/72446‑72407
...................................GGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACTTTCCC  >  REL606/3944634‑3944668
                                                                           
ACGCGCCGGGCAATGTTGACTACTCCCGCTATATTG                                         <  1:5540272/36‑1 (MQ=1)
 cgtGCCGGGCAATGTTGACTACTCCCGCTATATTGG                                        >  1:5183716/4‑36 (MQ=0)
 CGCGCCGGGCAATGTTGACTACTCCCGCCATATTGG                                        >  1:4901632/1‑36 (MQ=0)
  GCGCCGGGCAATGTTGACTACCCCCGCTATATTGGC                                       <  1:1608888/36‑1 (MQ=0)
   CGCCGGGCAATGTTGACTACTCCCGCTATATTGGCG                                      >  1:5598957/1‑36 (MQ=1)
   CGCCGGGCAATGTTGACTACTCCCGCTATATTGGCG                                      <  1:4669691/36‑1 (MQ=1)
    GCCGGGCAATGTTGACTACTCCCGCTATATTGGCGG                                     <  1:3842692/36‑1 (MQ=1)
    GCCGGGCAATGTTGACTACTCCCGCTATATTGGCGG                                     >  1:4550900/1‑36 (MQ=1)
      CGGGCAATGTTGACTACTCCCGCTATATTGGCGGGG                                   >  1:444221/1‑36 (MQ=1)
           AATGTTGACTACTCCCGCTATATTGGCGGGGTTTat                              >  1:1311005/1‑34 (MQ=1)
           AATGTTGACTACTCCCGCTATATTGGCGGGGTTTTT                              >  1:3187424/1‑36 (MQ=1)
            ATGTTGACTACTCCCGCTATATTGGCGGGGTTTTTT                             >  1:6514335/1‑36 (MQ=1)
             TGTTGACTACTCCCGCTATATTGGCGGGGTTTTTTt                            >  1:1097612/1‑35 (MQ=32)
                TGACTACTCCCGCTATATTGGCGGGGTTTTTTCTTT                         >  1:4710477/1‑36 (MQ=12)
                                   GGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACTT      <  1:557686/36‑1 (MQ=255)
                                     CGGGGTTTTTTGTTTTCTGTGCATTTCGTCACTTTC    >  1:6504966/1‑36 (MQ=255)
                                      GGGGTTTTTTGTTTTCTGTGCATTTCGTCACTTTtc   >  1:1580103/1‑34 (MQ=255)
                                      GGGGTTTTTTGTTTTCTGTGCATTTCGTCACTTTCC   <  1:3485908/36‑1 (MQ=255)
                                       GGGTTTTTTGTTTTATGTGCATTTCGTCACTTTCCC  >  1:6273621/1‑36 (MQ=255)
                                                                           
ACGCGCCGGGCAATGTTGACTACTCCCGCTATATTGGCGG...................................  <  REL606/72446‑72407
...................................GGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACTTTCCC  >  REL606/3944634‑3944668

Base quality scores:  ATCG  < 3 ≤  ATCG  < 8 ≤  ATCG  < 18 ≤  ATCG  < 29 ≤  ATCG  < 32 ≤  ATCG