New junction evidence
  seq id position reads (cov) reads (cov) score skew annotation gene product
* ? REL606 72407 =24 (0.510)5 (0.110) 4/60 5.3 coding (446/1701 nt) araB ribulokinase
?REL606 1286712 = 49 (0.890)intergenic (+36/‑147) narI/ECB_01206 nitrate reductase 1, gamma (cytochrome b(NR)) subunit/hypothetical protein
Rejected: Position hash score below cutoff.

ACGCGCCGGGCAATGTTGACTACTCCCGCTATATTGGCGG..................................  <  REL606/72446‑72407
...................................GGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGT  >  REL606/1286712‑1286745
                                                                          
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:4669691/36‑1 (MQ=1)
   CGCCGGGCAATGTTGACTACTCCCGCTATATTGGCG                                     >  1:5598957/1‑36 (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‑36 (MQ=32)
                                     CGGGGTTTTTTTATGGGCACGGTGCGGGGTGCGTTG   >  1:4697124/1‑36 (MQ=0)
                                      GGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGT  <  1:753203/36‑1 (MQ=1)
                                                                          
ACGCGCCGGGCAATGTTGACTACTCCCGCTATATTGGCGG..................................  <  REL606/72446‑72407
...................................GGCGGGGTTTTTTTATGGGCACGGTGCGGGGTGAGTTGT  >  REL606/1286712‑1286745

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