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.

GCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACG  >  REL606/809393‑809463
                                   |                                   
gCTATGCCTCCCTGATGATGTATATCAAAGTAAAAc                                     >  1:1655807/1‑36 (MQ=255)
  tATGCCTCCCTGATGATGTATATCAAAGTAAAACCg                                   <  2:1814134/36‑1 (MQ=255)
   aTGCCTCCCTGATGATGTATATCAAAGTAAAACCGc                                  >  2:990484/1‑36 (MQ=255)
   aTGCCTCCCTGATGATGTATATCAAAGTAAAACCGc                                  <  2:3209058/36‑1 (MQ=255)
   aTGCCTCCCTGATGATGTATATCAAAGTAAAACCGc                                  >  2:2316370/1‑36 (MQ=255)
   aTGCCTCCCTGATGATGTATATCAAAGTAAAACCGc                                  >  1:2238380/1‑36 (MQ=255)
   aTGCCTCCCTGATGATGTATATCAAAGTAAAAACGc                                  >  2:871796/1‑36 (MQ=255)
    tGCCTCCCTGATGATGTATATCAAAGTAAAACCGcc                                 >  1:969221/1‑36 (MQ=255)
            tgatgaTGTATATCAAAGTAAAACCGCCATTTTccc                         >  1:2461136/1‑36 (MQ=255)
             gatgatGTATATCAAAGTAAAACCGCCATTTTCTCt                        >  1:1118904/1‑36 (MQ=255)
              atgatgTATATCAAAGTAAAACCGGCATTTTCCCtt                       >  1:1962422/1‑36 (MQ=255)
              atgatgTATATCAAAGTAAAACCGCCATTTTCCCtt                       >  2:2250139/1‑36 (MQ=255)
              atgatgTATATCAAAGTAAAACCGCCATCATCCCtt                       >  2:3386482/1‑36 (MQ=37)
               tgataTATAAAAAAATAACCCCGCCATTTTCCCTTa                      <  2:581669/36‑1 (MQ=9)
                gatgTATATCAAAGTAAAACCGCCATTTTCCCTTAt                     >  1:2085510/1‑36 (MQ=255)
                 atgTATATCAAAGTAAAACCGCCATTTTCCCTTAtt                    >  1:2604827/1‑36 (MQ=255)
                  tgTATATCAAAGTAAAACCGCCATTTTCCCTTATTc                   >  2:211550/1‑36 (MQ=255)
                   gTATATCAAAGTAAAACCGCCATTTTCCCTTATTCt                  >  1:2527744/1‑36 (MQ=255)
                   gTATATCAAAGTAAAACCGCCATTTTCCCTTATTCt                  >  1:1163536/1‑36 (MQ=255)
                   gTATATCAAAGTAAAACCGCCATTTTCCCTTATTCt                  >  2:248969/1‑36 (MQ=255)
                   gTATATCAAAGTAAAACCGCCATTTTCCCTTATTCt                  >  1:948535/1‑36 (MQ=255)
                    tataTCAAAGTAAAACCGCCATTTTCCCTTATTCTg                 >  2:2204860/1‑36 (MQ=255)
                    tataTCAAAGTAAAACCGCCATTTTCCCTTATTCTg                 >  1:2615837/1‑36 (MQ=255)
                     atatCAAAGTAAAACCGCCATTTTCCCTCATTCTGt                >  2:1724023/1‑36 (MQ=255)
                     atatCAAAGTAAAACCGCCATCTTCCCTTATTTTGt                >  2:1876189/1‑36 (MQ=37)
                      tatCAAAGTAAAACCGCCATTTTCCCTTATTCTGTa               >  2:1249944/1‑36 (MQ=255)
                       atCAAAGTAAAACCGCCATTTTCCCTTATTTTGTAg              >  2:2890609/1‑36 (MQ=255)
                       atCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAg              >  1:1747494/1‑36 (MQ=255)
                       atCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAg              >  1:3609828/1‑36 (MQ=255)
                         cAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCg            >  1:1137524/1‑36 (MQ=255)
                         cAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCg            >  1:2070906/1‑36 (MQ=255)
                          aaaGTAAAACCGCCATTTTCCCTTATTCTGTAGCGa           >  2:271707/1‑36 (MQ=255)
                            aGTAAAACCGCCATTTTCCCTTATTCTGTAGCGaaa         >  2:103046/1‑36 (MQ=255)
                            aGTAAAACCGCCATTTTCCCTTATTCTGTAGCGaaa         >  2:1336943/1‑36 (MQ=255)
                            aGTAAAACCGCCATTTTCCATTATTCTGTAGCGaaa         >  1:2541060/1‑36 (MQ=255)
                            aGTAAAACCGCCATCTTCCCTTATTCTGTAGCGaaa         >  1:1235161/1‑36 (MQ=255)
                            aGTAAAACCGCAATTTTCCCTTATTCTGAAGCGaaa         >  2:1127828/1‑36 (MQ=37)
                             gTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAAt        >  2:1231646/1‑36 (MQ=255)
                             gTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAAt        >  1:857190/1‑36 (MQ=255)
                             aaaaacCCGCCATTTTCCCTTATTCTGTAGCGAAAt        <  2:2321664/34‑1 (MQ=37)
                              tAAAACCGCCCTTTTCCCTTATTCTGTAGCGAAATa       >  2:949968/1‑36 (MQ=255)
                              aaaacCCGCCATTTTCCCTTATTCTGTAGCGAAATa       <  2:2258734/35‑1 (MQ=255)
                               aaaaCCGCCATTTTCCCTTATTCTGTAGCGAAATAg      >  2:2752145/1‑36 (MQ=255)
                                 acccGCCATTTTCCCTTATTCTGTAGCGAAATAGCa    <  1:2678671/34‑1 (MQ=255)
                                 aaCCGCCATTTTCCCTTATTCTGTAGCGAAATAGCa    >  1:69025/1‑36 (MQ=255)
                                  cccGGCATTTTCCCTTATTCTGTAGCGAAATAGCAc   <  1:2249962/35‑1 (MQ=255)
                                  cccGCCATTTTCCCTTATTCTGTAGCGAAATAGCAc   <  2:3339525/35‑1 (MQ=255)
                                  aCCGCCATTTTCCCTTATTCTGTAGCGAAATAGCAc   >  1:1548725/1‑36 (MQ=255)
                                  aCCGCCATTTTCCCTTATTCTGTAGCGAAATAGCAc   >  2:3347956/1‑36 (MQ=255)
                                   ccGCCATTTTCCCTTATTCTGTAGCGAAATAGCACg  >  2:1746270/1‑36 (MQ=255)
                                   ccGCCATTTTCCCTTATTCTGTAGCGAAATAGCACg  <  2:2963442/36‑1 (MQ=255)
                                   ccGCCATTTTCCCTTATTCTGTAGCGAAATAGCACg  <  1:2224857/36‑1 (MQ=255)
                                   |                                   
GCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACG  >  REL606/809393‑809463

Base quality scores:  ATCG  < 3 ≤  ATCG  < 11 ≤  ATCG  < 19 ≤  ATCG  < 30 ≤  ATCG  < 32 ≤  ATCG