New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 3082514 =40 (0.760)8 (0.150) 5/60 NT 15.2% coding (446/1041 nt) yghZ aldo‑keto reductase
?REL606 = 4141016 49 (0.790)coding (89/705 nt) yijC DNA‑binding transcriptional repressor

TCATCGCGTCGATGAAAATACGCCGATGGAAGAAACCGCCTCTGCGCTGGCTCATGCGGTACAAAGCGGTA  >  REL606/3082479‑3082549
                                   |                                   
tcATCGCGTCGATGAAAATACGCCGATGGAAGAAAc                                     >  1:4723779/1‑36 (MQ=255)
gtgatgggcgtaagagcgcaacaaaaagaaaaaacc                                     >  1:1767827‑M1/36‑36 (MQ=255)
gtgatgggcgtaagagcgcaacaaaaagaaaaaacc                                     <  1:4787648‑M1/1‑1 (MQ=255)
gtgatgggcgtaagagcgcaacaaaaagaaaaaacc                                     >  1:7341501‑M1/36‑36 (MQ=255)
 tgatgggcgtaagagcgcaacaaaaagaaaaaaccc                                    <  1:656620‑M1/2‑1 (MQ=255)
 tgatgggcgtaagagcgcaacaaaaagaaaaaaccc                                    <  1:6852849‑M1/2‑1 (MQ=255)
 tgatgggcgtaagagcgcaacaaaaagaaaaaaccc                                    <  1:2923902‑M1/2‑1 (MQ=255)
 cATCGCGTCGATGAAAATACGCCGATGGAAGAAAcc                                    >  1:4209781/1‑36 (MQ=255)
   tCGCGTCGATGAAAATACGCCGATGGAAGAAACCGc                                  >  1:8393652/1‑36 (MQ=255)
   tCGCGTCGATGAAAATACGCCGATGGAAGAAACCGc                                  <  1:3136805/36‑1 (MQ=255)
    tgggcgtaagagcgcaacaaaaagaaaaaacCCGcc                                 <  1:5571892‑M1/5‑1 (MQ=255)
    tgggcgtaagagcgcaacaaaaagaaaaaacCCGcc                                 >  1:7433392‑M1/32‑36 (MQ=255)
     gcgTCGATGAAAATACGCCGATGGAAGAAACCGCct                                <  1:2695519/36‑1 (MQ=255)
     gcgTCGATGAAAATACGCCGATGGAAGAAACCGCct                                >  1:8357917/1‑36 (MQ=255)
       gTCGATGAAAATACGCCGATGGAAGAAACCGCctct                              >  1:1600227/1‑36 (MQ=255)
       gTCGATGAAAATACGCCGATGGAAGAAACCGCctct                              >  1:5899029/1‑36 (MQ=255)
           ataaaaaGACCCCGAGGAAAAAAACCGCCTCTgcgc                          <  1:8005283/33‑1 (MQ=9)
            tGAAAATACGCCGATGGAAGAAACCGCCTCTGCGCt                         >  1:3005920/1‑36 (MQ=255)
              aaaaAACGCCGATGAAAAAACCCGCCTCTGCGCTgg                       <  1:655007/36‑1 (MQ=19)
               aaaTACGCCGATGGAAGAAACCGCCTCTGCGCTGGc                      >  1:1918601/1‑36 (MQ=255)
                aaTCCGCGGAAGAAAGAAACCGCCTCTGCGCTGGCt                     <  1:2531223/36‑1 (MQ=16)
                aaTACGCCGATGGAAGAAACCGCCTCTGCGCTGGCt                     >  1:5646083/1‑36 (MQ=255)
                aaTACGCCGATGGAAGAAACCGCCTCTGCGCTGGCt                     >  1:4491260/1‑36 (MQ=255)
                aaTACGCCGATGGAAGAAACCGCCTCTGCGCTGGCt                     >  1:7153062/1‑36 (MQ=255)
                aaTACGCCGATGGAAGAAACCGCCTCTGCGCTCGCt                     >  1:3985981/1‑36 (MQ=255)
                  tACGCCGATGGAAGAAACCGCCTCTGCGCTGGCTCa                   >  1:2151964/1‑36 (MQ=255)
                   aCGCCGATGGAAGAAACCGCCTCTGCGCTGGCTCAt                  >  1:8001068/1‑36 (MQ=255)
                   aCGCCGATGAAAAAAACCGCCTCTGCGCTGGCTCAt                  <  1:3382732/36‑1 (MQ=37)
                    cGCCGAAGAAACAAACCGCCTCTGCGCTGGCTCATg                 <  1:3696708/36‑1 (MQ=25)
                     gCCGATGGAAGAAACCGCCTCTGCGCTGGCTCATGc                >  1:2733717/1‑36 (MQ=255)
                     gCCGATGGAAGAAACCGCCTCTGCGCTGGCTCATGc                >  1:186014/1‑36 (MQ=255)
                      ccGAGGAAAGAAACCGCCTCTGCGCTGGCTCATGCg               <  1:5152213/36‑1 (MQ=37)
                      ccGAGGAAAGAAACCGCCTCTGCGCTGGCTCATGCg               <  1:5423105/36‑1 (MQ=37)
                      ccGAAGAAAAAAACCGCCTCTGCGCTGGCTCATGCg               <  1:947803/36‑1 (MQ=25)
                        gATGGAAGAAACCGCCTCTGCGCTGGCTCATGCGGt             >  1:8429294/1‑36 (MQ=255)
                        gATGGAAGAAACCGCCTCTGCGCTGGCTCATGCGGt             >  1:3674092/1‑36 (MQ=255)
                        gATGAAAAAAACCGCCTCTGCGCTGGCTCATGCGGt             <  1:5351066/36‑1 (MQ=37)
                        gATGAAAAAAACCGCCTCTGCGCTGGCTCATGCGGt             <  1:6645441/36‑1 (MQ=37)
                         agggAAGAAACCGCCTCTGCGCTGGCTCATGCGGTa            <  1:8449046/34‑1 (MQ=255)
                         aTGGAAGAAACCGCCTCTGCGCTGGCTCATGCGGTa            >  1:5624144/1‑36 (MQ=255)
                          tGGAAGAAACCGCCTCTGCGCTGGCTCATGCGGTAc           >  1:63836/1‑36 (MQ=255)
                          ggaaagaacCCGCCTCTGCGCTGGCTCATGCGGTAc           >  1:2980193‑M1/10‑36 (MQ=255)
                          ggaaaaaaaCCGCCTCTGCGCTGGCTCATGCGGTAc           >  1:1686106‑M1/10‑36 (MQ=255)
                          agaaaaaacCCGCCTCTGCGCTGGCTCATGCGGTAc           >  1:3421879‑M1/10‑36 (MQ=255)
                           ggAAGAAACCGCCTCTGCGCTGGCTCATGCGGTACa          >  1:8112587/1‑36 (MQ=255)
                           ggAAGAAACCGCCTCTGCGCTGGCTCATGCGGTACa          >  1:6923838/1‑36 (MQ=255)
                           gaaaaaacCCGCCTCTGCGCTGGCTCATGCGGTACa          >  1:2043722‑M1/9‑36 (MQ=255)
                            gaagaaACCGCCTCTGCGCTGGCTCATGCGGTACaa         >  1:692851/1‑36 (MQ=255)
                            gaagaaACCGCCTCTGCGCTGGCTCATGCGGTACaa         >  1:3063050/1‑36 (MQ=255)
                            gaagaaACCGCCTCTGCGCTGGCTCATGCGGTACaa         >  1:3793264/1‑36 (MQ=255)
                            aaaaaagCCGCCTCTGCGCTGGCTCATGCGGTACaa         >  1:7890023‑M1/8‑36 (MQ=255)
                             aaaaaaCCGCCTCTGCGCTGGCTCATGCGGTACaaa        >  1:3998721‑M1/7‑36 (MQ=255)
                              agaaACCGCCTCTGCGCTGGCTCATGCGGTACAAAg       >  1:5455207/1‑36 (MQ=255)
                               gaaACCGCCTCTGCGCTGGCTCATGCGGTACAAAGc      >  1:2560855/1‑36 (MQ=255)
                               gaaACCGCCTCTGCGCTGGCTCATGCGGTACAAAGc      >  1:1184532/1‑36 (MQ=255)
                               aaacCCGCCTCTGCGCTGGCTCATGCGGTACAAAGc      >  1:3501407‑M1/5‑36 (MQ=255)
                               aaaaCCGCCTCTGCGCTGGCTCATGCGGTACAAAGc      <  1:1103782/35‑1 (MQ=255)
                                aacCCGCCTCTGCGCTGGCTCATGCGGTACAAAGCg     >  1:2674592‑M1/4‑36 (MQ=255)
                                 aaCCGCCTCTGCGCTCGCTAATGCGGTACAAAGCgg    >  1:6134108/1‑36 (MQ=37)
                                  aCCGCCTCTGCGCTGGCTCATGCGGTACAAAGCGGt   <  1:7071695/36‑1 (MQ=255)
                                   CCGCCTCTGCGCTGGCTCATGCGGTCCAAAGCGGTa  >  1:310998‑M1/1‑36 (MQ=255)
                                   CCGCCTCTGCGCTGGCTCATGCGGTACAAAGCGGTa  >  1:1978682‑M1/1‑36 (MQ=255)
                                   |                                   
TCATCGCGTCGATGAAAATACGCCGATGGAAGAAACCGCCTCTGCGCTGGCTCATGCGGTACAAAGCGGTA  >  REL606/3082479‑3082549

Base quality scores:  ATCG  < 3 ≤  ATCG  < 9 ≤  ATCG  < 22 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG