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

AATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGC  >  REL606/2574478‑2574546
                                 |                                   
aaTTCACTTTCCTTTAAAAAACAAAAAACCGCCaaa                                   >  1:1686856/1‑36 (MQ=255)
  ttCACTTTCCTTTAAAAAACAAAAAACCGCCAAAAt                                 >  1:122810/1‑36 (MQ=255)
   tCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATc                                >  1:7419680/1‑36 (MQ=255)
    cACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCa                               >  1:3027627/1‑36 (MQ=255)
      cTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAgg                             >  1:5308574/1‑36 (MQ=255)
       tttCCTTTAAAAAACAAAAAACCGCCAAAATCAGGc                            >  1:4407696/1‑36 (MQ=255)
         tCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCgg                          >  1:124636/1‑36 (MQ=255)
            tttAAAAAACAAAAAACCGCCAAAATCAGGCGGttt                       >  1:1297190/1‑36 (MQ=255)
            aaaaaaaaaaaaaaaccCGCCAAAATCAGGCGGttt                       <  1:1983213/33‑1 (MQ=25)
             ttAAAAAACAAAAAACCGCCAAAATCAGGCGGtttt                      >  1:7273318/1‑36 (MQ=255)
              tAAAAAACAAAAAACCGCCAAAATCAGGCGGttttt                     >  1:7042608/1‑36 (MQ=255)
               agaaaataaaaaccCGCC.GAATTTGGCGGtttttta                   <  1:4660203/34‑2 (MQ=17)
               aaaaaaaaaaaaacCGCCAgcaaaatcaccaatatc                    <  1:4110355‑M1/36‑18 (MQ=255)
               aaaaaaaaaaaaacCGCCAAAATCAGGCGGtttttt                    <  1:1250294/36‑1 (MQ=255)
               aaaaaaaaaaaaacCGCC.GAATGTGGCGGtttttta                   <  1:5784263/36‑2 (MQ=17)
                 aaaacaaaaaacCGCCAAAATCAGGCGGtttttttt                  >  1:6836717/1‑34 (MQ=255)
                  aaaaaaaaaacCGCCAgcaaaatcaccaatatcctg                 <  1:4945465‑M1/36‑21 (MQ=255)
                   aacaaaaaacCGCCAAAATCAGGCGGTTTTtttttt                >  1:2502716/1‑35 (MQ=255)
                    aaaaaaaacCGCCAAAATCAGGCGGTTTTTTGTTGc               <  1:3624828/34‑1 (MQ=255)
                      aaaaaacCGCCAAATTCGGCGGGTTTTTTTTTTCtt             >  1:3874100/1‑35 (MQ=17)
                       aaacacCGCCAgcaaaatcaccaatatcctggacgg            <  1:3389681‑M1/36‑26 (MQ=255)
                       aaaaacCGCCAAAATCAGGCGGTTTTTTTTTTCTgt            >  1:1254888/1‑35 (MQ=21)
                       aaaaacCGCCAAAATCAGGCGGTTTTTTGTTGCttg            >  1:5338627/1‑34 (MQ=255)
                          aacCGCCAgcaaaatcaccaatatcctggacggcac         <  1:4617993‑M1/36‑29 (MQ=255)
                           acCGCCAgcaaaatcaccaatatcctggacggcaca        <  1:1513545‑M1/36‑30 (MQ=255)
                            ccGCCAgcaaaatcaccaatatcctggacggcacag       >  1:6437419‑M1/1‑6 (MQ=255)
                            ccGCCAgcaaaatcaccaatatcctggacggcacag       <  1:4510608‑M1/36‑31 (MQ=255)
                            ccGCCAgcaaaatcaccaatatcctggacggcacag       <  1:886271‑M1/36‑31 (MQ=255)
                             cGCCAgcaaaatcaccaatatcctggacggcacagt      >  1:949467‑M1/1‑5 (MQ=255)
                             cGCCAgcaaaatcaccaatatcctggacggcacagt      <  1:3497867‑M1/36‑32 (MQ=255)
                             cGCCAgcaaaatcaccaatatcctggacggcacagt      >  1:3440792‑M1/1‑5 (MQ=255)
                             cGCCAgcaaaatcaccaatatcctggacggcacagt      <  1:209875‑M1/36‑32 (MQ=255)
                             cGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGt      <  1:4375358/36‑1 (MQ=255)
                              gCCAgcaaaatcaccaatatcctggacggcacagtc     <  1:4469476‑M1/36‑33 (MQ=255)
                              gCCAgcaaaatcaccaatatcctggacggcacagtc     >  1:5429663‑M1/1‑4 (MQ=255)
                              gCCAgcaaaatcaccaatatcctggacggcacagtc     >  1:5668774‑M1/1‑4 (MQ=255)
                              gCCAgcaaaatcaccaatatcctggacggcacagtc     >  1:3627214‑M1/1‑4 (MQ=255)
                              gCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGtt     <  1:3696225/36‑1 (MQ=255)
                               ccAgcaaaatcaccaatatcctggacggcacagtca    <  1:4346070‑M1/36‑34 (MQ=255)
                               ccAgcaaaatcaccaatatcctggacggcacagtca    <  1:948256‑M1/36‑34 (MQ=255)
                                cAgcaaaatcaccaatatcctggacggcacagtcac   >  1:5155554‑M1/1‑2 (MQ=255)
                                 agcaaaatcaccaatatcctggacggcacagtcact  <  1:5579443‑M1/36‑36 (MQ=255)
                                 agcaaaatcaccaatatcctggacggcacagtcact  >  1:5899783‑M1/1‑1 (MQ=255)
                                 agcaaaatcaccaatatcctggacggcacagtcact  >  1:2915379‑M1/1‑1 (MQ=255)
                                 agcaaaatcaccaatatcctggacggcacagtcact  <  1:1540444‑M1/36‑36 (MQ=255)
                                 |                                   
AATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGC  >  REL606/2574478‑2574546

Base quality scores:  ATCG  < 3 ≤  ATCG  < 11 ≤  ATCG  < 23 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG