New junction evidence
  seq id position reads (cov) reads (cov) score skew annotation gene product
* ? REL606 = 33947552 (1.510)5 (0.150) 4/44 4.0 coding (592/1083 nt) lacI lac repressor
?REL606 = 2022999 50 (0.910)coding (600/1071 nt) hisC histidinol‑phosphate aminotransferase
Rejected: Position hash score below cutoff.

AGTTTTGCCCACAGGCATCGCTGGCTGGCTGGCTGACGGAATATCCGCACCTGGCTATTTTGCGCACAC  >  REL606/2022965‑2023033
                                  |                                  
aattttGCCCACAGGCATCGCTGGCTGGCTGGCTGa                                   <  1:3822815/34‑1 (MQ=255)
aGTTTTGCCCACAGGCATCGCTGGCTGGGTGGCTGa                                   >  1:6938428/1‑36 (MQ=255)
aGTTTTGCCCACAGGCATCGCTGGCTGGCTGGGTGa                                   >  1:2424225/1‑36 (MQ=255)
aGTTTTGCCCACAGGCATCGCTGGCTGGCTGGCTGa                                   <  1:1124105/36‑1 (MQ=255)
aGTTTTGCCCACAGGCATCGCTGGCTGGCTGGCTGa                                   <  1:6703609/36‑1 (MQ=255)
 gTTTTGCCCACAGGCATCGCTGGCTGGCTGGCTGAc                                  <  1:7401414/36‑1 (MQ=255)
 gTTTTGCCCACAGGCATCGCTGGCTGGCTGGCTGAc                                  <  1:7051232/36‑1 (MQ=255)
 gTTTTGCCCACAGGCATCGCTGGCTGGCTGGCTGAc                                  >  1:5556354/1‑36 (MQ=255)
 gTTTTGCCCACAGGCATCGCTGGCTGGCTGGCTGAc                                  <  1:4775905/36‑1 (MQ=255)
   tttGCCCACAGGCATCGCTGGCTGGCTGGCTGACgg                                >  1:6172636/1‑36 (MQ=255)
     tGCCCACAGGCATCGCTGGCTGGCTGGGTGACggga                              >  1:1512962/1‑34 (MQ=25)
      gCCCACAGGCATCGCTGGCTGGCTGGCCGACGGGat                             >  1:3573540/1‑36 (MQ=25)
       cccACAGGCATCGCTGGCTGGCTGGCTGACGGAata                            <  1:4072617/36‑1 (MQ=255)
        ccACAGGCATCGCTGGCTGGCTGGCTGACGGAatat                           <  1:4265296/36‑1 (MQ=255)
        ccACAGGCATCGCTGGCTGGCTGGCTGACGGAatat                           >  1:2424515/1‑36 (MQ=255)
         cacaGGCATCGCTGGCTGGCTGGCTGACGGGatatt                          >  1:330574/1‑35 (MQ=255)
          acaGGCATCGCTGGCTGGGTGGGTGGGGGGatattg                         >  1:1273434/1‑34 (MQ=9)
          acaGGCATCGCTGGCTGGCTGGCTGACGGAATATcc                         <  1:1484638/36‑1 (MQ=255)
          acaGGCATCGCTGGCTGGCTGGCTGACGGAAAATcc                         >  1:2556798/1‑36 (MQ=255)
           caGGCATCGCTGGCTGGCTGGGTGACGGAATATCCg                        >  1:5819991/1‑36 (MQ=255)
           caGGCATCGCTGGCTGGCTGGCTGACGGAATATCCg                        >  1:3757666/1‑36 (MQ=255)
            aGGCATCGCTGGCTGGCTGGCTGACGGAATATCCGc                       <  1:664497/36‑1 (MQ=255)
             ggCATCGCTGGCTGGCTGGCTGACGGGATATCCGCa                      >  1:6388502/1‑36 (MQ=255)
             ggCATCGCTGGCTGGCTGGCTGACGGAATATCCGCa                      >  1:5739235/1‑36 (MQ=255)
               cATCGGTGGCTGGCTGGCTGACGGAATATCCGCAcc                    <  1:3483494/36‑1 (MQ=255)
               cATCGCTGGCTGGCTGGCTGACGGAATATCCGCAcc                    >  1:816325/1‑36 (MQ=255)
               cATCGCTGGCTGGCTGGCTGACGGAATATCCGCAcc                    <  1:5354959/36‑1 (MQ=255)
                aTCGCTGGCTGGCTGGCTGACGGAATATCCGCACCt                   <  1:5750483/36‑1 (MQ=255)
                 tCGCTGGCTGGCTGGCTGACGGAATATCCGCACCTg                  <  1:243655/36‑1 (MQ=255)
                   gctggctggctggctgACGGAATATCCGCACCTGGc                <  1:6003203/36‑1 (MQ=255)
                   gctggccgcctgtctgACGGAATATCCGCACCTGGc                <  1:6864435/36‑1 (MQ=21)
                    cgggctggcgggctgACGGAATATCCGCACCTGGCt               <  1:1637550/34‑1 (MQ=25)
                     tggctgggtggctgACGGAATATCCGCACCTGGCTa              <  1:5485849/36‑1 (MQ=255)
                     tggctggctggctgACGGAATATCCGCACCTGGCTa              <  1:1602291/36‑1 (MQ=255)
                     tggctggctggctgACGGAATATCCGCACCTGGCTa              <  1:2846262/36‑1 (MQ=255)
                     tggctggctggctgACCGAAAATCCGCCCCCGCCTa              >  1:3158692/1‑36 (MQ=12)
                       ggtggctggctgACGGAATATCCGCACCTGGCTAtt            <  1:2379079/34‑1 (MQ=255)
                       gctggctggctgATGGAATATCCGCACCTGGCTAtt            <  1:7021108/36‑1 (MQ=255)
                        ctgggtggctgACGGATTATCCGCACCTGGCTAttt           <  1:1383497/36‑1 (MQ=25)
                           gctggctgACGGAATATCCGCACCTGGCTATTTTgc        <  1:7192222/36‑1 (MQ=255)
                           gctggctgACGGAATATCCGCACCTGGCTATTTTgc        <  1:1422545/36‑1 (MQ=255)
                           gctggctgACGGAATATCCGCACCTGGCTATTTTgc        >  1:6640141/1‑36 (MQ=255)
                           gctggcagACGGAATATCCGCACCTGGCTATTTTgc        <  1:178556/36‑1 (MQ=255)
                            ctggctgGCGGAATATCCGCACCTGGCTATTTTgcg       <  1:4851638/36‑1 (MQ=255)
                              ggctgACGGAATATCCGCACCTGGCTATTTTGCGca     <  1:1418415/36‑1 (MQ=255)
                              ggctgACGGAATATCCGCACCTGGCTATTTTGCGca     <  1:6995263/36‑1 (MQ=255)
                              ggctgACGGAATATCCGCACCTGGCTATTTTGCGca     >  1:2277487/1‑36 (MQ=255)
                                ctgACGGAATATCCGCACCTGGCTATTTTGCGcaca   >  1:55902/1‑36 (MQ=255)
                                 tgACGGAATATCCGCACCTGGCTATTTTGCGCacac  >  1:3595717/1‑36 (MQ=255)
                                 tgACGGAATATCCGCACCTGGCTATCTTTCGCacac  >  1:748051/1‑36 (MQ=25)
                                  |                                  
AGTTTTGCCCACAGGCATCGCTGGCTGGCTGGCTGACGGAATATCCGCACCTGGCTATTTTGCGCACAC  >  REL606/2022965‑2023033

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 14 ≤  ATCG  < 27 ≤  ATCG  < 32 ≤  ATCG