New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 257451122 (0.480)4 (0.090) 4/56 NT 10.1% intergenic (+353/+465) sseA/sseB 3‑mercaptopyruvate sulfurtransferase/rhodanase‑like enzyme, sulfur transfer from thiosulfate
?REL606 3715185 = 49 (0.850)coding (2195/4851 nt) ECB_03460 putative adhesin

TGTCTTCACTGCAGCTCATGGCAACAATACCGCCAGCAAAATCACCAATATCCTGGACGGCACAGTCACTG  >  REL606/3715150‑3715220
                                   |                                   
tgtCTTCACTGCAGCTCATGGCAACAATACCGCCAg                                     >  1:4093287/1‑36 (MQ=255)
 gtCTTCACTGCAGCTCATGGCAACAATACCGCCAGc                                    >  1:3044624/1‑36 (MQ=255)
  tCTTCACTGCAGCTCATGGCAACAATAAAGCCAGCa                                   >  1:2392410/1‑36 (MQ=37)
   cTTCACTGCAGCTCATGGCAACAATACCGCCAGCaa                                  >  1:6352105/1‑36 (MQ=255)
    ttCACTGCAGCTCATGGCAACAATACCGCCAGCaaa                                 >  1:4397692/1‑36 (MQ=255)
    ttCACTGCAGCTCATGGCAACAATACCGCCAGCaaa                                 >  1:3168620/1‑36 (MQ=255)
    ttCACTGCAGCTCATGGCAACAATACCGCCAGCaaa                                 >  1:5218751/1‑36 (MQ=255)
    ttCACTGCAGCTCATGGCAACAATACCGCCAGCaaa                                 >  1:2036782/1‑36 (MQ=255)
     tCACTGCAGCTCATGGCAACCATACCGCCAGCaaaa                                >  1:124272/1‑36 (MQ=255)
     tCACTGCAGCTCATGGCAACAATACCGCCAGCaaaa                                >  1:5375278/1‑36 (MQ=255)
      cACTGCAGCTCATGGCAACAATACCGCCAGCAAAAt                               >  1:1750584/1‑36 (MQ=255)
      cACTGCAGCTCATGGCAACAATACCGCCAGCAAAAt                               >  1:6346015/1‑36 (MQ=255)
      cACTGCAGCTCAGGAAAAAAATACCGCCAGCAAAAt                               <  1:1530204/36‑1 (MQ=16)
       acggCAGCTAAGGGAAAAAATACCGCCAGCAAAATc                              <  1:2003826/33‑1 (MQ=12)
       aCTGCAGCTCATGGCAACAATACCGCCAGCAAAATc                              >  1:6133440/1‑36 (MQ=255)
       aCTGCAGCTCATGGCAACAATACCGCCACCAAAATc                              >  1:6527808/1‑36 (MQ=255)
        cTGCAGCTCATGGNAACAATACCGCCAGCAAAATCa                             >  1:2892870/1‑36 (MQ=255)
         tGCAGCTCATGGCAACAATACCGCCAGCAAAATCAc                            >  1:1646602/1‑36 (MQ=255)
          gCAGCTCATGGCAACAATACCGCCAGCAAAATCAcc                           >  1:1357003/1‑36 (MQ=255)
            aGCTCATGGCAACAATACCGCCAGCAAAATCACCaa                         >  1:5169890/1‑36 (MQ=255)
            aGCTCATGGCAACAATACCGCCAGCAAAATCAACaa                         >  1:4285060/1‑36 (MQ=255)
             gCTCATGGCAACAATACCGCCAGCAAAATCACCAat                        >  1:5201737/1‑36 (MQ=255)
              cTCATGGCAACAATACCGACTGCAAAATCACCAata                       >  1:4905761/1‑36 (MQ=37)
               tCATGGCAACAATACCGCCAGCAAAATCACCAatat                      >  1:4535695/1‑36 (MQ=255)
                aaaaaaaaaaaaaccgccaGCAAAATCACCAATATc                     <  1:4110355‑M2/17‑1 (MQ=255)
                 aTGGCAACAATACCGCCAGCAAAATCACCAATATcc                    >  1:2350883/1‑36 (MQ=255)
                 aTGGCAACAATACCGCCAGCAAAATCACCAATATcc                    >  1:3334201/1‑36 (MQ=255)
                 aTGGCAACAATACCGCCAGCAAAATCACCAATATcc                    >  1:2561975/1‑36 (MQ=255)
                  tGGCAACAATACCGCCAGCAAAATCACCAATATCCt                   >  1:7350043/1‑36 (MQ=255)
                  tGGCAACAATACCGCCAGCAAAATCACCAATATCCt                   >  1:2904043/1‑36 (MQ=255)
                  tGGCAACAATACCGCCAGCAAAATCACCAATATCCt                   >  1:7491699/1‑36 (MQ=255)
                  tGGCAACAATACAGCCAGCAAAATCACCAATATCCt                   >  1:7115325/1‑36 (MQ=255)
                   ggCAACCATACCGCCAGCAAAATCACCAATATCCTg                  >  1:2515255/1‑36 (MQ=255)
                   ggCAACAATACCGCCAGCAAAATCACCAATATCCTg                  >  1:4148922/1‑36 (MQ=255)
                   aaaaaaaaaaccgccaGCAAAATCACCAATATCCTg                  <  1:4945465‑M2/20‑1 (MQ=255)
                    gCAACAATACCGCCAGCAAAATCACCAATATCCTgg                 >  1:3547113/1‑36 (MQ=255)
                     caacaaTACCGCCAGCAAAATCACCAATATCCTTGa                >  1:4756425/1‑36 (MQ=255)
                     caacaaTACCGCCAGCAAAATCACCAATATCCTGGa                >  1:5308887/1‑36 (MQ=255)
                     caacaaTACCGCCAGCAAAATCACCAATATCCTGGa                >  1:268169/1‑36 (MQ=255)
                     caacaaTACCGCCAGCAAAATCACCAATATCCTGGa                >  1:2249875/1‑36 (MQ=255)
                     caacaaTACCGCCAGCAAAATCACCAATATCCTGGa                >  1:4586007/1‑36 (MQ=255)
                      aacaaTACCGCCAGCAAAATCACCAATATCCTGGAc               >  1:3888846/1‑36 (MQ=255)
                      aacaaTACCGCCAGCAAAATCACCAATATCCTGGAc               >  1:5415325/1‑36 (MQ=255)
                      aacaaTACCGCCAGCAAAATCACCAATATCCTGGAc               >  1:3146223/1‑36 (MQ=255)
                      aacaaTACCGCCAACAAAATCACCAATATCATGGAc               >  1:7143962/1‑36 (MQ=37)
                        aaacaccgccaGCAAAATCACCAATATCCTGGACgg             <  1:3389681‑M2/25‑1 (MQ=255)
                         aaTACCGCCAGCAAAATCACCAATATCCTGGACGGc            >  1:5820595/1‑36 (MQ=255)
                         aaTACCGCCAGCAAAATCACCAATATCCTGGACGGc            >  1:286606/1‑36 (MQ=255)
                         aaTACCGCCAGCAAAATCACCAATATCCTGGACGGc            >  1:2839367/1‑36 (MQ=255)
                         aaTACCGCCAGCAAAATCACCAATATCCTGGACGGc            >  1:3576238/1‑36 (MQ=255)
                         aaTACCGCCAGCAAAATCACCAATATCCTGGACGGc            >  1:3973699/1‑36 (MQ=255)
                          aTACCGCCAGCAAAATCACCAATATCCTGGACGGca           >  1:3607351/1‑36 (MQ=255)
                           tACCGCCAGCAAAATCACCAATATCCTGGACGGcac          >  1:3844331/1‑36 (MQ=255)
                           aaccgccaGCAAAATCACCAATATCCTGGACGGcac          <  1:4617993‑M2/28‑1 (MQ=255)
                            accgccaGCAAAATCACCAATATCCTGGACGGcaca         <  1:1513545‑M2/29‑1 (MQ=255)
                             ccgccaGCAAAATCACCAATATCCTGGACGGCACAg        >  1:6437419‑M2/7‑36 (MQ=255)
                             ccgccaGCAAAATCACCAATATCCTGGACGGCACAg        <  1:4510608‑M2/30‑1 (MQ=255)
                             ccgccaGCAAAATCACCAATATCCTGGACGGCACAg        <  1:886271‑M2/30‑1 (MQ=255)
                              cgccaGCAAAATCACCAATATCCTGGACGGCACAGt       <  1:209875‑M2/31‑1 (MQ=255)
                              cgccaGCAAAATCACCAATATCCTGGACGGCACAGt       >  1:949467‑M2/6‑36 (MQ=255)
                              cgccaGCAAAATCACCAATATCCTGGACGGCACAGt       >  1:3440792‑M2/6‑36 (MQ=255)
                              cgccaGCAAAATCACCAATATCCTGGACGGCACAGt       <  1:3497867‑M2/31‑1 (MQ=255)
                               gccaGCAAAATCACCAATATCCTGGACGGCACAGTc      >  1:5668774‑M2/5‑36 (MQ=255)
                               gccaGCAAAATCACCAATATCCTGGACGGCACAGTc      >  1:5429663‑M2/5‑36 (MQ=255)
                               gccaGCAAAATCACCAATATCCTGGACGGCACAGTc      <  1:4469476‑M2/32‑1 (MQ=255)
                               gccaGCAAAATCACCAATATCCTGGACGGCACAGTc      >  1:3627214‑M2/5‑36 (MQ=255)
                                ccaGCAAAATCACCAATATCCTGGACGGCACAGTCa     <  1:4346070‑M2/33‑1 (MQ=255)
                                ccaGCAAAATCACCAATATCCTGGACGGCACAGTCa     <  1:948256‑M2/33‑1 (MQ=255)
                                 caGCAAAATCACCAATATCCTGGACGGCACAGTCAc    >  1:5155554‑M2/3‑36 (MQ=255)
                                  aGCAAAATCACCAATATCCTGGACGGCACAGTCACt   <  1:5579443‑M2/35‑1 (MQ=255)
                                  aGCAAAATCACCAATATCCTGGACGGCACAGTCACt   >  1:5899783‑M2/2‑36 (MQ=255)
                                  aGCAAAATCACCAATATCCTGGACGGCACAGTCACt   >  1:2915379‑M2/2‑36 (MQ=255)
                                  aGCAAAATCACCAATATCCTGGACGGCACAGTCACt   <  1:1540444‑M2/35‑1 (MQ=255)
                                   gCAAAATCACCAATATCCTGGACGGCACAGTCACTg  <  1:5064834/36‑1 (MQ=255)
                                   gCAAAATCACCAATATCCTGGACGGCACAGTCACTg  >  1:4537474/1‑36 (MQ=255)
                                   gCAAAATCACCAATAGCCTGGACGGCACAGTCACTg  <  1:1318244/36‑1 (MQ=255)
                                   |                                   
TGTCTTCACTGCAGCTCATGGCAACAATACCGCCAGCAAAATCACCAATATCCTGGACGGCACAGTCACTG  >  REL606/3715150‑3715220

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