New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 809428 =43 (0.910)8 (0.170) 6/60 4.4 16.0% intergenic (‑155/‑242) ybhK/moaA predicted transferase with NAD(P)‑binding Rossmann‑fold domain/molybdenum cofactor biosynthesis protein A
?REL606 = 4141016 41 (0.740)coding (89/705 nt) yijC DNA‑binding transcriptional repressor
Rejected: Position hash score below cutoff.

TGATCGTGCTATTTCGCTACAGAATAAGGGAAAATGGCGG..................................  <  REL606/809467‑809428
...................................GGCGGGTTTTTTCTTTTTGTTGCGCTCTTACGCCCATCA  <  REL606/4141016‑4140983
                                                                          
     GTGCTATTTCGCTACAGAATAAGGGAAAATGGCGGG                                   >  2:3339525/1‑36 (MQ=255)
     GTGCTATTTCGCTACAGAATAAGGGAAAATGCCGGG                                   >  1:2249962/1‑36 (MQ=255)
      TGCTATTTCGCTACAGAATAAGGGAAAATGGCGGGT                                  >  1:2678671/1‑36 (MQ=255)
         TATTTCGCTACAGAATAAGGGAAAATGGCGGGTTTT                               >  2:2258734/1‑36 (MQ=255)
         TATTTCGCTACAGAATAAGGGAAAATGGCGGGGTTT                               >  1:1318725/1‑36 (MQ=255)
          ATTTCGCTACAGAATAAGGGAAAATGGCGGGTTTTT                              >  2:2321664/1‑36 (MQ=255)
                CTACAGAATAAGGGAAAATGGCGGGTTTTTTTTTTT                        >  1:1913099/1‑36 (MQ=255)
                 TACAGAATAAGGGAAAATGGCGGGTTTTTTTTTTTT                       >  1:358746/1‑36 (MQ=255)
                                                                          
TGATCGTGCTATTTCGCTACAGAATAAGGGAAAATGGCGG..................................  <  REL606/809467‑809428
...................................GGCGGGTTTTTTCTTTTTGTTGCGCTCTTACGCCCATCA  <  REL606/4141016‑4140983

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