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.

ACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGC  >  REL606/4157595‑4157665
                                   |                                   
aCCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCa                                     >  1:3037922/1‑36 (MQ=255)
aCCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCa                                     >  1:2814052/1‑36 (MQ=255)
 ccTGGTGCGGGGCGTACCGCGTGTGATGGGCTTCAt                                    <  1:1502251/36‑1 (MQ=37)
 ccTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCAt                                    <  1:826601/36‑1 (MQ=255)
 ccTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCAt                                    >  2:2487183/1‑36 (MQ=255)
  ttGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATc                                   >  1:1431629/2‑36 (MQ=255)
  cTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATc                                   <  1:3691591/36‑1 (MQ=255)
   tGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATcg                                  <  2:3691640/36‑1 (MQ=255)
   tGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATcg                                  <  2:1850802/36‑1 (MQ=255)
    ggTGCGCAGCGTACCGCGTGTGATGGGCTTCATcgg                                 <  1:1439329/36‑1 (MQ=255)
    ggTCCGCAGCGTACCGCGTGTGATGGGCTTCATcgg                                 <  2:268689/36‑1 (MQ=255)
     gTGCGCAGCGTCCCGCGGGTGATGGGCTTCATcggc                                <  2:3285704/36‑1 (MQ=37)
      tGCTCAGCATACCGCGTGTGAGGGGCTTCATcggcg                               <  1:263605/36‑1 (MQ=25)
        cgcAGCGTACCGCGTGTGATGGGCTTCATcggcggg                             >  1:241148/1‑35 (MQ=255)
         gcAGCGTACCGCGTGTGATGGGCTTCATcggcgggg                            >  2:3145682/1‑34 (MQ=255)
         gcAGCGTACCGCGTGTGATGGGCTTCATcggcgggg                            >  2:2501921/1‑34 (MQ=255)
         gcAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTa                            <  2:809211/36‑1 (MQ=255)
          cAGCGTACCGCGTGTGATGGGCTTCATcggcgggag                           >  1:2425187/1‑33 (MQ=255)
          cAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTAc                           <  1:3261049/36‑1 (MQ=255)
           aGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACt                          <  2:494431/36‑1 (MQ=255)
            gCGTACCGCGTGTGATGGGCTTCATCGGCGGTACtt                         <  2:3415879/36‑1 (MQ=255)
              ttACCGCGTGTGATGGGCTTCATCGGCGGTACTTcc                       <  2:2278758/35‑1 (MQ=255)
               tACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCg                      <  2:2846775/36‑1 (MQ=255)
                tccGCGTGTGATGGGCTTCATCGGCGGTACTTCCGa                     <  2:502784/35‑1 (MQ=255)
                aCCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGa                     <  1:402234/36‑1 (MQ=255)
                  cgcgTGTGATGGGCTTCATCGGCGGTACTTCCGATc                   <  1:2448522/36‑1 (MQ=255)
                  cgcgTGTGATGGGCTTCATCGGCGGTACTTCCGATc                   <  2:2325288/36‑1 (MQ=255)
                  cgcgTGTGATGGGCTTCATCGGCGGTACTTCCGATc                   <  1:1487504/36‑1 (MQ=255)
                   gggtgtGATGGGCTTCATCGGCGGTACTTCCGATCg                  <  1:492168/34‑1 (MQ=255)
                    cgTGTGATGGGCTTCATCGGCGGTACTGCCGATCGt                 <  2:853277/36‑1 (MQ=255)
                      tgtgATGGGCTTCATCGGCGGTACTTCCGATCGTcc               <  2:1456311/36‑1 (MQ=255)
                      tgtgATGGGCTTCATCGGCGGTACTTCCGATCGTcc               <  2:2124747/36‑1 (MQ=255)
                      tgtgATGGGCTTCATCGGCGGGGGTTCCGTTCGttc               >  1:2416163/1‑34 (MQ=16)
                       gtgATGGGCTTCATCGGCGGGGATTGCGATCGTCCt              >  1:2853273/1‑36 (MQ=19)
                        tgAGGGGCTTCATCGGCGGTACTTCCGATCGTCCTg             <  1:3732744/36‑1 (MQ=255)
                         gATGGGCTTCATCGGCGGTACTTCCGATCGTCCTgc            <  2:557402/36‑1 (MQ=255)
                         gATGGGCTGCATCGGCGGTACTTCCGATCGTCCTgc            <  2:3064921/36‑1 (MQ=255)
                          aTGGTCTTCATCGGCGGTACTTCCGATCGTCCTgcg           <  2:1191329/36‑1 (MQ=255)
                            gggTTTCATCGGCGGTACTTCCGATCGTCCTGCGcc         <  2:2318373/36‑1 (MQ=255)
                            gggCTTCATCGGCGGTACTTCCGATCGTCCTGCGcc         <  1:2399836/36‑1 (MQ=255)
                             ggCTTCATCGGCGGCACTTTCCGTCGTCCTGCCCCa        >  1:2516364/1‑36 (MQ=12)
                              gCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCaa       <  2:336922/36‑1 (MQ=255)
                              gCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCaa       <  1:2725992/36‑1 (MQ=255)
                               cTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAAt      <  1:423584/36‑1 (MQ=255)
                               cTTCATCGGCGGTACTCCCCATCGTCCTTCGCCact      >  1:2028682/1‑34 (MQ=21)
                                ttCATCGGCGGTACTTCCGATCGTCCTGCGCCAATc     <  2:3236058/36‑1 (MQ=255)
                                 tCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCa    <  2:2204208/36‑1 (MQ=255)
                                  cATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAg   <  1:1277124/36‑1 (MQ=255)
                                   aTCGGCGGTACTTCCGATCGTCCTGCGCCAGTCAGc  <  2:2707049/36‑1 (MQ=255)
                                   aTCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGc  <  2:1115154/36‑1 (MQ=255)
                                   |                                   
ACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGC  >  REL606/4157595‑4157665

Base quality scores:  ATCG  < 3 ≤  ATCG  < 5 ≤  ATCG  < 13 ≤  ATCG  < 27 ≤  ATCG  < 32 ≤  ATCG