New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 277140951 (0.960)4 (0.080) 3/60 NT 8.9% intergenic (‑219/‑33) cysD/iap sulfate adenylyltransferase subunit 2/aminopeptidase in alkaline phosphatase isozyme conversion
?REL606 4412609 = 31 (0.500)intergenic (‑376/‑158) ECB_04076/ytfA putative oxoacyl reductase/predicted transcriptional regulator

TAAACAGACGGTTTATCATCACTTTGGTAATAAAGAAAATTTATTTATTGAAGTACTGGAATATACCTGGA  >  REL606/4412574‑4412644
                                   |                                   
taAACAGACGGTTTATCATCACTTTGGTAATAAAGa                                     <  1:2246021/36‑1 (MQ=255)
agcgcagcggcagtttcatactatggcgggtaaaaa                                     >  1:2188333‑M2/36‑36 (MQ=255)
 aaaCAGACGGTTTATCATCACTTTGGTAATAAAGaa                                    <  1:4434091/36‑1 (MQ=255)
  aaCAGACGGTTTATCATCACTTTGGTAATAAAGaaa                                   <  1:6426571/36‑1 (MQ=255)
  aaCAGACGGGTTATCATCACTTTGGTAATAAAGaaa                                   >  1:6132411/1‑36 (MQ=255)
     aGACGGTTTATCATCACTTTGGTAATAAAGAAAatt                                <  1:2371958/36‑1 (MQ=255)
     aGACGGTTTATCATCACTTTGGTAATAAAGAAAatt                                >  1:7072906/1‑36 (MQ=255)
      gACGGTTTATCATCACTTTGGTAATAAAGAAAattt                               <  1:1575250/36‑1 (MQ=255)
      gACGGTTTATCATCACTTTGGTAATAAAGAAAattt                               >  1:4603919/1‑36 (MQ=255)
      gACGGTTTATCATCACTTTGGTAATAAAGAAAattt                               >  1:2282949/1‑36 (MQ=255)
      gACGGTTTATCATCACTTTGGTAATAAAGAAAattt                               <  1:3203852/36‑1 (MQ=255)
       aCGGTTTATCATCACTTTGGTAATAAAGAAAattta                              >  1:2019645/1‑36 (MQ=255)
        cGGTTTATCATCACTTTGGTAATAAAGAAAatttat                             <  1:1405435/36‑1 (MQ=255)
         ggTTTATCATCACTTTGGTAATAAAGAAAatttatt                            <  1:3625621/36‑1 (MQ=255)
         ggTTTATCATCACTTTGGTAATAAAGAAAatttatt                            <  1:6977481/36‑1 (MQ=255)
          gTTTATCATCACTTTGGTAATAAAGAAAatttattt                           <  1:3285213/36‑1 (MQ=255)
          gTTTATCATCACTTTGGTAATAAAGAAAatttattt                           <  1:3289848/36‑1 (MQ=255)
          cagtttcatactatggcggttaaaaAAAatttattt                           >  1:4897116‑M2/26‑36 (MQ=255)
          cagtttcatactatggcggttaaaaAAAatttattt                           >  1:1213105‑M2/26‑36 (MQ=255)
             tttcatactatggcggttaaaaAAAATTtatttttt                        >  1:4957704‑M2/23‑36 (MQ=255)
              atcatcACTTTGGTAATAAAGAAAATTtatttatta                       <  1:6106535/36‑2 (MQ=255)
               tcatcaCTTTGGTAATAAAGAAAATTTATTTATTGa                      >  1:6188360/1‑36 (MQ=255)
               tcatcaCTTTGGTAATAAAGAAAATTTATTTATTGa                      <  1:1951480/36‑1 (MQ=255)
                catcaCTTTGGTAATAAAGAAAATTTATTTATTGaa                     <  1:3317888/36‑1 (MQ=255)
                catcaCTTTGGTAATAAAGAAAATTTATTTATTGaa                     >  1:7205634/1‑36 (MQ=255)
                 atcaCTTTGGTAATAAAGAAAATTTATTTATTGAAg                    <  1:4438334/36‑1 (MQ=255)
                   caCTTTGGTAATAAAGAAAATTTATTTATTGAAGTa                  <  1:6346438/36‑1 (MQ=255)
                   caCTTTGGTAATAAAGAAAATTTATTTATTGAAGTa                  <  1:380380/36‑1 (MQ=255)
                    aCTTTGGTAATAAAGAAAATTTATTTATTGAAGTAc                 >  1:3894293/1‑36 (MQ=255)
                     cTTTGGTAATAAAGAAAATTTATTTATTGAAGTACt                >  1:4379006/1‑36 (MQ=255)
                     cTTTGGTAATAAAGAAAATTTATTTATTGAAGTACt                >  1:2833532/1‑36 (MQ=255)
                      tttGGTAATAAAGAAAATTTATTTATTGAAGTACTg               <  1:7198580/36‑1 (MQ=255)
                          gTAATAAAGAAAATTTATTTATTGAAGTACTGGAat           >  1:3363609/1‑36 (MQ=255)
                           taataaAGAAAATTTATTTATTGAAGTACTGGAata          >  1:3005077/1‑36 (MQ=255)
                           taataaAGAAAATTTATTTATTGAAGTACTGGAata          <  1:8206916/36‑1 (MQ=255)
                               aaaGAAAATTTATTTATTGAAGTACTGGAATATAcc      >  1:8166889/1‑36 (MQ=255)
                                aaGAAAATTTATTTATTGAAGTACTGGAATATACCt     <  1:2613273/36‑1 (MQ=255)
                                  gAAAATTTATTTATTGAAGTACTGGAATATACCTgg   <  1:5750249/36‑1 (MQ=255)
                                  gAAAATTTATTTATTGAAGTACTGGAATATACCTgg   >  1:2621870/1‑36 (MQ=255)
                                   aaaaTTTATTTATTGAAGTACTGGAATATACTTGGa  >  1:7328591/1‑36 (MQ=255)
                                   |                                   
TAAACAGACGGTTTATCATCACTTTGGTAATAAAGAAAATTTATTTATTGAAGTACTGGAATATACCTGGA  >  REL606/4412574‑4412644

Base quality scores:  ATCG  < 3 ≤  ATCG  < 16 ≤  ATCG  < 26 ≤  ATCG  < 34 ≤  ATCG  < 37 ≤  ATCG