New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 144001848 (1.170)3 (0.070) 3/52 5.3 5.9% coding (894/1500 nt) feaB phenylacetaldehyde dehydrogenase
?REL606 = 4081044 47 (0.850)intergenic (+132/‑128) sodA/kdgT superoxide dismutase, Mn/2‑keto‑3‑deoxygluconate permease
Rejected: Position hash score below cutoff.

CGCAATGGGTTATTGAAGGCTTGATGACCGGAAGCTTCCTGAATCAAGGGCAAGTATGCGCCGCCAGTTC  >  REL606/1439984‑1440053
                                  |                                   
cGCAATGGGTTATTGAAGGCTTGATGACCGGAAGCt                                    <  2:3731853/36‑1 (MQ=255)
 gCAATGGGTTATTGAAGGCTTGATGACCGGAAGCtt                                   >  1:773028/1‑36 (MQ=255)
  cAATGGGTTATTGAAGGCTTGATGACCGGAAGCTTc                                  <  2:3254442/36‑1 (MQ=255)
     tGGGTTATTGAAGGCTTGATGACCGGAAGCTTTCTg                               >  2:1909942/1‑36 (MQ=255)
     tGGGTTATTGAAGGCTTGATGACCGGAAGCTTCCTg                               <  1:1377155/36‑1 (MQ=255)
     tGGGTTATTGAAGGCTTGATGACCGGAAGCTTCCTg                               <  1:572977/36‑1 (MQ=255)
      gggTTATTGAAGGCTTGATGACCGGAAGCTTCCTGa                              <  2:2180870/36‑1 (MQ=255)
      gggTTATTGAAGGCTTGATGACCGGAAGCTTCCTGa                              >  1:1051972/1‑36 (MQ=255)
       ggTTATTGAAGGCTTGATGACCGGAAGCCTCCCGaa                             >  2:2393903/1‑36 (MQ=37)
         ttATTGAAGGCTTGATGACCGGAAGCTTCCTGAATc                           <  2:815425/36‑1 (MQ=255)
         ttATTGAAGGCTTGATGACCGGAAGCTTCCTGAATc                           <  2:1481622/36‑1 (MQ=255)
          tATTGAAGGCTTGATGACCGTAAGCTTCCTGAATCa                          >  2:3213914/1‑36 (MQ=255)
          tATTGAAGGCTTGATGACCGGAAGCTTCCTGAATCa                          >  1:367232/1‑36 (MQ=255)
          tATTGAAGGCTTGATGACCGGAAGCTTCCTGAATCa                          <  1:3708738/36‑1 (MQ=255)
          tATTGAAGGCTTGATGACCGGAAGCTTCCTGAATCa                          <  2:1884566/36‑1 (MQ=255)
           aTTGAAGGCTTGATGACCGGAAGCTTTCTGAACCaa                         >  1:3699876/1‑36 (MQ=37)
           aTTGAAGGCTTGATGACCGGAAGCTTCCTGAATCaa                         <  1:3115587/36‑1 (MQ=255)
            ttGAAGGCTTGATGACCGGAAGCTTCCTGAATCAAg                        <  1:675487/36‑1 (MQ=255)
            ttGAAGGCTTGATGACCGGAAGCTTCCTGAATCAAg                        <  2:1107969/36‑1 (MQ=255)
             tGAAGGCTTGATGACCGGAAGCTTCCTGAATCAAgg                       <  2:707781/36‑1 (MQ=255)
              gAAGGCTTGATGACCGGAAGCTTCCTGAATCAAggg                      <  2:934891/36‑1 (MQ=255)
              gAAGGCTTGATGACCGGAAGCTTCCTGAATCAAggg                      <  1:1715359/36‑1 (MQ=255)
               aaGGCTTGATGACCGGAAGCTTCCTGAATCCAGGGc                     >  2:1023386/1‑36 (MQ=255)
                  gCTTGATGACCGGAAGCTTCCTGAATCAAGAGCAAg                  >  2:1933148/1‑36 (MQ=255)
                  gCTTAATAACCGGAAGCTTCCTGAATCAAGGGCAAg                  <  1:2750339/36‑1 (MQ=38)
                   cTTGATGAGCGGGAGCTTCCTGAGTCAAGGGCAAGt                 <  1:3471677/36‑1 (MQ=37)
                   cTTGATGACCGGAAGCTTCCTGAATCAAGGGCAAGt                 >  1:1550742/1‑36 (MQ=255)
                   cTTGATGACCGGAAGCTTCCTGAATCAAGGGCAAGt                 <  2:1656450/36‑1 (MQ=255)
                   cTTGATGACCGGAAGCTTCCTGAATCAAGGGCAAGt                 <  1:1853947/36‑1 (MQ=255)
                    gtgatgaCCGGAAGCTTCCTGAATCAAGGGCAAGTa                <  2:650474/35‑1 (MQ=255)
                     tgatgaCCGGAAGCTTCCTGAATCAAGGGCAAGTAt               >  1:3425704/1‑36 (MQ=255)
                       atgaCCGGAAGCTTCCTGAAGCAAGGGCAAGTATgc             >  2:2878499/1‑36 (MQ=255)
                        tgaCCGGAAGCTTCCTGATTCAAGGGCAAGTATgcg            <  1:540014/36‑1 (MQ=255)
                        tgaCCGGAAGCTTCCTGAATCAAGGGCAAGAATgcg            <  2:3596428/36‑1 (MQ=255)
                        tgaCCGGAAGCTTCCTGAATCAAGGGCAAGAACgcg            >  2:926287/1‑36 (MQ=37)
                          aCCGGAAGCTTCCTGAATCAAGGGCAAGTATGcgcc          <  1:2938453/36‑1 (MQ=255)
                          aCCGGAAGCTTCCTGAATCAAGGGCAAGTATGcgcc          <  2:1676759/36‑1 (MQ=255)
                          aCCGGAAGCTTCCTGAATCAAGGGCAAGTATGcgcc          <  2:3621697/36‑1 (MQ=255)
                           ccGGAAGCTTCCTGAATCAAGGGCAAGTATGcgccg         >  2:1693794/1‑36 (MQ=255)
                           acGGAAGCTTCCTGAATCAAGGGCAAGTATGcgccg         <  2:1175021/35‑1 (MQ=255)
                            cgaaaGCTTCCCGAATCAAGGGCAAGTATGcgccgc        <  2:2936530/33‑1 (MQ=38)
                            cGGAAGCTTCCTGAATCAAGGGCAAGTATGcgccgc        >  1:1690880/1‑36 (MQ=255)
                             ggAAGCTTCCTGAATCAAGGGCAAGTATGCgccgcc       >  2:3236420/1‑36 (MQ=255)
                             ggAAGCTTCCTGAATCAAGGGCAAGTATGCgccgcc       <  2:1975518/36‑1 (MQ=255)
                             ggAAGCTTCCTGAATCAAGGGCAAGTATGCgccgcc       <  2:2685580/36‑1 (MQ=255)
                               aaGCTTCCTGAATCAAGGGCAAGTATGCGCCGCCAg     >  2:3321862/1‑36 (MQ=255)
                               aaGCTTCCTGAATCAAGGGCAAGTATGCGCCGCCAg     <  2:2454115/36‑1 (MQ=255)
                                aGCTTCCTGAATCAAGGGAAAGTATGCTCCGCCAGt    <  1:2570271/36‑1 (MQ=38)
                                  cTTCCTGAATCAAGGGCAAGTATGCGCCGCCAGTTc  <  2:1419331/36‑1 (MQ=255)
                                  |                                   
CGCAATGGGTTATTGAAGGCTTGATGACCGGAAGCTTCCTGAATCAAGGGCAAGTATGCGCCGCCAGTTC  >  REL606/1439984‑1440053

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