New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 3944627 =36 (0.820)8 (0.180) 4/56 5.1 17.0% intergenic (+16/‑131) proM/aslB tRNA‑Pro/predicted regulator of arylsulfatase activity
?REL606 = 4157630 42 (0.760)coding (277/546 nt) nusG transcription antitermination protein NusG
Rejected: Position hash score below cutoff.

GAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCA  >  REL606/3944592‑3944661
                                   |                                  
gAATCCTCTCTCGCCGACCAATTTTGAACCCCGCtt                                    >  1:2920486/1‑36 (MQ=255)
    ccTTTCTCGCCGACCAATTTTGAACCCCGCTTcggc                                >  2:1696069/1‑36 (MQ=255)
    ccTCTCTCGCCGACCAATTTTGAACCCCGCTTcggc                                >  2:2998917/1‑36 (MQ=255)
     ctctctCGCCGACCAATTTTAAACCCCGCTTcggcg                               >  2:732583/1‑36 (MQ=255)
      tctctcGCCGACCAATTTTGAACCCCGCTTcggcgg                              <  1:75115/36‑1 (MQ=255)
       ctctcGCCGACCAATTTTGAACCCCGCTTCGGCggg                             <  2:500997/36‑1 (MQ=255)
       ctctcGCCGACCAATTTTGAACCCCGCTTCGGCggg                             <  1:2413784/36‑1 (MQ=255)
          tcGGCGACCAATTTTGAGCCCCGCTTCGGCGGGGtt                          <  2:1346201/36‑1 (MQ=25)
           cGCCGACCAATTTTGAACCCCGCTTCGGCGGGGttt                         <  1:582895/36‑1 (MQ=255)
            gCCGACCAATTTTGAACCCCGCTTCGGCGGGGtttt                        >  2:500882/1‑36 (MQ=255)
            gCCGACCAATTTTGAACCCCGCTTCGGCGGGGgttt                        >  2:3040757/1‑36 (MQ=255)
             ccGACCAATTTTGAACCCCGCTTCGGCGGGGttttt                       <  1:1422353/36‑1 (MQ=255)
             ccGACCAATTTTGAACCCCGCTTCGGCGGGGttttt                       <  1:1890278/36‑1 (MQ=255)
               gACCAATTTTGCACCCCGCTCCGGCGGGCTTTTTTg                     >  2:3001625/1‑36 (MQ=17)
               gACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTg                     <  2:47369/36‑1 (MQ=255)
               gACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTg                     >  2:3603343/1‑36 (MQ=255)
               gACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTg                     >  2:2305116/1‑36 (MQ=255)
               gACCAATTTTGAACCCCGCTTCGGCGGGGTTTCTTg                     >  2:1981479/1‑36 (MQ=255)
                 ccAATTTTGAACCCCGCTTCGGCGTGGTTTTTTGtt                   <  2:61179/36‑1 (MQ=255)
                  cAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGttt                  >  2:3699228/1‑36 (MQ=255)
                   aaTTTTGAACCCCGCTTCGGCGGGGTTTTTTGtttt                 <  1:906144/36‑1 (MQ=255)
                   aaTTTTGAACCCCGCTTCGGCGGGGTTTTTTGtttt                 >  2:1108485/1‑36 (MQ=255)
                      tttGAACCCCGCTTCGGCGGGGTTTTTTGttttttt              >  1:3247583/1‑33 (MQ=255)
                      tttGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCtg              >  2:536917/1‑36 (MQ=255)
                       ttGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCtgt             <  2:1085266/36‑1 (MQ=255)
                         gAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGc           <  1:2959295/36‑1 (MQ=255)
                         gAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGc           <  1:238937/36‑1 (MQ=255)
                            ccccGCTTAGGCGGGGTTTTTTGTTTTCTGCGCAtt        >  1:3190502/1‑36 (MQ=21)
                             cccGCTTCGGCGGGGTTTTTTGTTTTCTGTGCAttt       <  2:1152347/36‑1 (MQ=255)
                               cGCTTCGGCGGGGTTTTTTGTTTTCTGTGCTTTTCg     >  1:3321398/1‑36 (MQ=255)
                               cGCTTCGGCGGGGTTTTTTGTTTTCTGTGCAtttag     >  1:2719336/1‑34 (MQ=255)
                               cGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCg     <  2:1705368/36‑1 (MQ=255)
                                 cTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTc   >  2:2495634/1‑36 (MQ=255)
                                  ttCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCa  >  2:2078122/1‑36 (MQ=255)
                                  ttCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCa  >  2:1779645/1‑36 (MQ=255)
                                  ttCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCa  <  1:1837966/36‑1 (MQ=255)
                                  ttCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCa  <  1:2655250/36‑1 (MQ=255)
                                   |                                  
GAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCA  >  REL606/3944592‑3944661

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