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.

AATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGAT  >  REL606/4081010‑4081079
                                  |                                   
aaTAAAAACAATGATCCGATAAAAATAAAACAGCGt                                    <  2:2534240/36‑1 (MQ=255)
 aTAAAAACAATGATCCGATAAAAATAAAACAGCGtt                                   >  2:1214803/1‑36 (MQ=255)
 aTAAAAACAATGATCCGATAAAAATAAAACAGCGtt                                   >  2:2763301/1‑36 (MQ=255)
 aTAAAAACAATGATCCGACAAAAATAAACCAGCGtt                                   >  1:3564401/1‑36 (MQ=37)
  tAAAAACAATGATCCGATAAAAATAAAACAGCGttt                                  <  1:3566381/36‑1 (MQ=255)
  tAAAAACAATGATCCGATAAAAATAAAACAGCGttt                                  >  2:2561006/1‑36 (MQ=255)
  tAAAAACAATGATCCGATAAAAATAAAACAGCGttt                                  >  2:2051113/1‑36 (MQ=255)
   aaaaaCAATGATCCGA.AAAAATAAAACAGCGTTTCa                                <  2:3331707/36‑1 (MQ=37)
    aaaaCAATGATCCGATAAAAATAAAACAGCGTTTCa                                <  2:462751/36‑1 (MQ=255)
    aaaaCAATGATCCGATAAAAATAAAACAGCGTTTCa                                <  2:1356030/36‑1 (MQ=255)
    aaaaCAATGATCCGATAAAAATAAAACAGCGTTTCa                                <  2:1495177/36‑1 (MQ=255)
    aaaaCAATGATCCGATAAAAATAAAACAGCGTTTCa                                <  2:1702475/36‑1 (MQ=255)
     aaaCAATGATCCGATAAAAATAAAACAGCGTTTCaa                               <  1:436219/36‑1 (MQ=255)
      aaCAATGATCCGATAAAAATAAAACAGCGTTTCAAt                              >  1:2366451/1‑36 (MQ=255)
      aaCAATGATCCGATAAAAATAAAACAGCGTTTCAAt                              <  1:1135583/36‑1 (MQ=255)
       aCAATGATCCGATAAAAATAAAACAGCGTTTCAAtt                             >  1:2634700/1‑36 (MQ=255)
       aCAATGATCCAATAAAAATAAAACAGCGTTTCAAtt                             <  1:3048878/36‑1 (MQ=255)
         aaTGATCCGATAAAAATAAAACAGCGTTTCAATTGa                           >  1:3382245/1‑36 (MQ=255)
          aTGATCCGATAAAAATAAAACAGCGTTTCAATTGAt                          <  1:1570541/36‑1 (MQ=255)
             aTCCGATTAAAATAAGACAGCGTTTCAATTGAtgtg                       <  2:3451303/36‑1 (MQ=37)
             aTCCGATAAAAATAAAACAGCGTTTCAATTGAtgtg                       >  2:906319/1‑36 (MQ=255)
             aTCCGATAAAAATAAAACAGCGTTTCAATTGAtgtg                       >  1:2487656/1‑36 (MQ=255)
               ccGATAAAAATAAAACAGCGTTTCAATTGATGTGGt                     <  2:1505205/36‑1 (MQ=255)
                cGATAAAAATAAAACAGCGTTTCAATTGATGTGGtt                    >  2:2061611/1‑36 (MQ=255)
                cGATAAAAAT.AAACAGCGTTTCAATTGATGTGGttt                   >  1:554703/1‑36 (MQ=38)
                 gATAAAAATAAAACAGCGTTTCAATTGATGTGGttt                   <  2:1087903/36‑1 (MQ=255)
                  ataaaaataaaaCAGCGTTTCAATTGATGTGGtttt                  <  1:2740954/36‑1 (MQ=255)
                   taaaaataaaaCAGCGTTTCAATTGATGTGGTTTTg                 <  1:866039/36‑1 (MQ=255)
                   taaaaataaaaCAGCGTTTCAATTGATGTGGTTTTg                 >  1:2035117/1‑36 (MQ=255)
                    aaaaataaaaCAGCGTTTCAATTGATGTGGTTTTGa                <  1:1423807/36‑1 (MQ=255)
                    aaaaataaaaCAGCGTTTCAATTGATGTGGTTTTGa                >  1:35360/1‑36 (MQ=255)
                     aaaataaaaCAGCGTTTCAATTGATGTGGTTTTGAt               <  1:1141504/36‑1 (MQ=255)
                     aaaataaaaCAGCGTTTCAATTGATGTGGTTTTGAt               <  1:663988/36‑1 (MQ=255)
                         taaaaCAGCGTTTCAATTGATGTGGTTTTGATCACt           >  1:2432992/1‑36 (MQ=255)
                         taaaaCAGCGTTTCAATTGATGTGGTTTTGATCACt           <  2:2454986/36‑1 (MQ=255)
                          aaaaCAGCGTTTCAATTGATGTGGTTTTGATCTCtt          >  1:3682382/1‑36 (MQ=255)
                          aaaaCAGCGTTTCAATTGATGTGGTTTTGATCACtt          >  1:3313834/1‑36 (MQ=255)
                           aaaCAGCGTTTCAATTGATGTGGTTTTGATCACttt         >  2:2342497/1‑36 (MQ=255)
                           aaaCAGCGTTTCAATTGATGTGGTTTTGATCACttt         <  2:1395013/36‑1 (MQ=255)
                           aaaCAGCGTTTCAATTGATGTGGTTTTGATCACttt         >  1:1347759/1‑36 (MQ=255)
                            aaCAGCGTTTCAGTTGATGTGGTTTTGATCACtttt        >  1:2613476/1‑36 (MQ=255)
                             aCAGCGTTTCAATTGATGTGGTTTTTATCACTTTTa       >  2:2794559/1‑36 (MQ=255)
                              cAGCGTTTCAATTGATGTGGTTTTGATCACTTTTAt      <  2:2077218/36‑1 (MQ=255)
                              cAGCGTTTCAATTGATGTGGTTTTGATCACTTTTAt      <  1:724468/36‑1 (MQ=255)
                               aGCGTTTCAATTGATGTGGTTTTGATCATTTTTAtt     >  2:2971339/1‑36 (MQ=255)
                                gCGTTTCAATTGATGTGGTTTTGATCACTTTTATTg    <  2:2595784/36‑1 (MQ=255)
                                gCGTTTCAATTGATGTGGTTTTGATCACTTTTATTg    <  2:3529025/36‑1 (MQ=255)
                                  gTTTCAATTGATGTGGTTTTGATCACTTTTATTGAt  >  1:185462/1‑36 (MQ=255)
                                  gTTTCAATTGATGTGGTTTTGATCACTTTTATTGAt  <  2:721328/36‑1 (MQ=255)
                                  |                                   
AATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGAT  >  REL606/4081010‑4081079

Base quality scores:  ATCG  < 3 ≤  ATCG  < 13 ≤  ATCG  < 23 ≤  ATCG  < 30 ≤  ATCG  < 32 ≤  ATCG