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.

CGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAC.CCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAA  >  REL606/4140981‑4141050
                                   |                                   
cGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAc                                     <  2:1812273/36‑1 (MQ=255)
cGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAc                                     <  2:1300945/36‑1 (MQ=255)
cGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAc                                     <  2:3576722/36‑1 (MQ=255)
  tGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAc.cc                                  <  1:615051/36‑1 (MQ=255)
  tGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAc.cc                                  >  1:380763/1‑36 (MQ=255)
  tGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAc.cc                                  <  2:3738721/36‑1 (MQ=255)
   gATGGGCGTAAGAGCGCAACAAAAAGAAAAAAC.ccg                                 >  2:580504/1‑36 (MQ=255)
   gATGGGCGTAAGAGCGCAACAAAAAGAAAAAAC.ccg                                 >  2:464498/1‑36 (MQ=255)
    aTGGGCGTAAGAGCGCAACAAAAAGAAAAAAC.ccgc                                <  2:1529374/36‑1 (MQ=255)
    aTGGGCGTAAGAGAGCAACAAAAAGAAAAAAC.ccgc                                <  2:2360322/36‑1 (MQ=255)
    aTGAGCGTAAGAGCGCAACAAAAAGAAAAAAC.ccgc                                <  2:2353459/36‑1 (MQ=255)
     tGGGCGTAAGAGCGCAACAAAAAGAAAAAAC.ccgcc                               >  2:2082902/1‑36 (MQ=255)
       ggCGTAAGAGCGCAACAAAAAGAAAAAAC.CCGCCGt                             >  1:2662804/1‑36 (MQ=255)
            aaGAGCGCAACAAAAAGAAAAAAC.CCGCCGTTCGCt                        >  2:2041119/1‑36 (MQ=255)
             agagCGCAACAAAAAGAAAAAAC.CCGCCGTTCGCtg                       >  1:2621355/1‑36 (MQ=255)
             agagCGCAACAAAAAGAAAAAAC.CCGCCGTTCGCtg                       >  1:3068826/1‑36 (MQ=255)
               agCGCAACAAAAAGAAAAAAC.CCGCCGTTCGCtggt                     >  2:1518151/1‑36 (MQ=255)
                gggcAAAAAAAAGAAAAACC.CCGCCGTTCGCtggtg                    <  1:2321715/34‑1 (MQ=37)
                  gcAACAAAAAGAAAAAAC.CCGCCGTTCGCTGGTGGa                  >  1:2604602/1‑36 (MQ=255)
                   caacaaAAAGAAAAAAC.CCGCCGTTCGCTGGTGGaa                 >  2:1955152/1‑36 (MQ=255)
                   caacaaAAAGAAAAAAC.CCGCCGTTCGCTGGTGGaa                 >  2:159588/1‑36 (MQ=255)
                     acaaAAAGAAAAACC.CCGCCGTTCGCTGGTGGAAGc               <  1:2232866/36‑1 (MQ=255)
                     acaaAAAGAAAAAAC.CCGCCGTTCGCTGGTGGAAGc               >  1:1700814/1‑36 (MQ=255)
                     acaaAAAGAAAAAAC.CCGCCGTTCGCTGGTGGAAGc               >  1:1321065/1‑36 (MQ=255)
                      caaAAAGAAAAAAC.CCGCCGTTCGCTGGTGGAAGcc              >  2:865780/1‑36 (MQ=255)
                      caaAAAGAAAAAAC.CCGCCGTTCGCTGGTGGAAGcc              >  1:1276951/1‑36 (MQ=255)
                      aaaaaaGAAAAACC.CCGCCGTTCGCTGGTGGAAGcc              <  1:524132/35‑1 (MQ=255)
                      aaaaaaAAAAAACC.CCGCCGTTCGCTGGTGGAAGcc              <  2:1013802/35‑1 (MQ=25)
                       aaaaaGAAAAAAC.CCGCCGTTCGCTGGTGGAAGcgt             >  2:2959273/1‑34 (MQ=255)
                       aaaaaGAAAAAAC.CCGCCGTTCGCTGGTGGAAGCCg             <  2:2654366/36‑1 (MQ=255)
                       aaaaaAAAAAAACACCGCCGTTCGCTGGTGGAAGcc              <  2:2318109/36‑1 (MQ=21)
                        aaaaGAAAAACC.CCTCCGTTCGCTGGTGGAAGCCGc            <  1:2661837/36‑1 (MQ=37)
                        aaaaGAAAAACC.CCGCCGTTCGCTGGTGGAAGCCGc            <  2:3663772/36‑1 (MQ=255)
                         aaaGAAAAAAC.CCGCCGTTCGCTGGTGGAAGCCGCa           >  1:3520853/1‑36 (MQ=255)
                         aaaAAAAAACC.CCGCCGTTCGCTGGTGGAAGCCGCa           <  2:3602655/36‑1 (MQ=37)
                          aaGAAAAAAC.CCGCCGTTCGCTGGTGGAAGGCGCAt          >  1:3022736/1‑36 (MQ=255)
                          aaGAAAAAAC.CCGCCGTTCGCTGGTGGAAGCCGCAt          >  2:89873/1‑36 (MQ=255)
                          aaGAAAAAAC.CCGCCGTTCGCTGGTGGAAGCCGCAt          >  1:1824492/1‑36 (MQ=255)
                          aaGAAAAAAA.ACGACGTTAGCTGGTGGAAGACGcca          >  2:700979/1‑34 (MQ=12)
                           aGAAAAAAC.CCGCCGTTCGCTGGTGGAAGCCGCAtt         >  1:2047123/1‑36 (MQ=255)
                            gAAAAAAC.CCGCCGTTCGCTGGTGGAAGCCGCAttt        >  2:2516415/1‑36 (MQ=255)
                             aaaaaaC.CCGCCGTTCGCTGGTGGAAGCCGCATTTa       <  1:2363223/36‑1 (MQ=255)
                              aaaacC.CCGCCGTTCGCTGGTGGAAGCCGCATTTAg      <  2:575608/36‑1 (MQ=255)
                               aaaaC.CCGCCGTTCGCTGGTGGAAGCCGCATTTAGc     >  1:2608689/1‑36 (MQ=255)
                                 acc.ccGCCGTTCGCTGGTGGAAGCCGCATTTAGCCa   <  2:3099180/34‑1 (MQ=255)
                                  cc.ccGCGGTTGGCTGGTGGAAGCCGCATTTAGCCaa  <  1:1105259/35‑1 (MQ=37)
                                   |                                   
CGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAAC.CCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAA  >  REL606/4140981‑4141050

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