New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 257451122 (0.470)8 (0.170) 5/54 NT 27.1% intergenic (+353/+465) sseA/sseB 3‑mercaptopyruvate sulfurtransferase/rhodanase‑like enzyme, sulfur transfer from thiosulfate
?REL606 4356096 = 21 (0.350)coding (338/537 nt) blc outer membrane lipoprotein (lipocalin)

TCCCGATCAAGTGCAATGACGTTATAACCGCCATAGAAAGGACCAAAGAATGACACTTTCAGCGCAGC  >  REL606/4356063‑4356130
                                 |                                  
tCCCGATCAAGTGCAATGACGTTATAACCGCCATAg                                  >  1:5611500/1‑36 (MQ=255)
tCCCGATCAAGTGCAATGACGTTATAACCGCCATAg                                  >  1:6269860/1‑36 (MQ=255)
 cccGATCAAGTGCAATGACGTTATAACCGCCATAGa                                 >  1:5873648/1‑36 (MQ=255)
   cGATCAAGTGCAATGACGTTATAACCGCCATAGaaa                               >  1:3633312/1‑36 (MQ=255)
    gATCAAGTGCAATGACGTTATAACCGCCATAGAAAg                              >  1:1528275/1‑36 (MQ=255)
        aaGTGCAATGACGTTATAACCGCCATAGAAAGGAcc                          >  1:4661389/1‑36 (MQ=255)
         aGTGCAATGACGTTATAACCGCCATAGAAAGGACCa                         >  1:1367280/1‑36 (MQ=255)
         aGTGCAATGACGTTATAACCGCCATAGAAAGGACCa                         >  1:1827679/1‑36 (MQ=255)
          gTGCAATGACGTTATAACCGCCATAGAAAGGACCaa                        >  1:3730885/1‑36 (MQ=255)
           tGCAATGACGTTATAACCGCCATAGAAAGGACCaaa                       >  1:6718597/1‑36 (MQ=255)
           tGCAATGACGTTATAACCGCCATAGAAAGGACCaaa                       >  1:7774217/1‑36 (MQ=255)
              aaaaaaaaaaaaaccgccaTAGAAAGGACCAAAGaa                    <  1:5602230‑M2/17‑1 (MQ=255)
              aaaaaaaaaaaaaccgccaTAGAAAGGACCAAAGaa                    <  1:3458943‑M2/17‑1 (MQ=255)
              aaTGACGTTATAACCGCCATAGAAAGGACCCAAGaa                    >  1:3985756/1‑36 (MQ=255)
              aaTGACGTTATAACCGCCATAGAAAGGACCAAAGaa                    >  1:6228988/1‑36 (MQ=255)
              aaTGACGTTATAACCGCCATAGAAAGGACCAAAGaa                    >  1:7108217/1‑36 (MQ=255)
              aaTGACGTTATAACCGCCATAGAAAGGACCAAAGaa                    >  1:2301882/1‑36 (MQ=255)
              aaTGACGTTATAACCGCCATAGAAAGGACCAAAGaa                    >  1:1837845/1‑36 (MQ=255)
               aaaaaaaaaaaaccgccaTAGAAAGGACCAAAGAAt                   <  1:4823930‑M2/18‑1 (MQ=255)
                  aCGTTATAACCGCCATAGAAAGGACCAAAGAATGac                >  1:7407225/1‑36 (MQ=255)
                  aCGTTATAACCGCCATAGAAAGGACCAAAGAATGac                >  1:2297884/1‑36 (MQ=255)
                   cGTTATAACCGCCATAGAAAGGACCAAAGAATGaca               >  1:5541375/1‑36 (MQ=255)
                   aaaaaaaaccgccaTAGAAAGGACCAAAGAATGaca               <  1:2424060‑M2/22‑1 (MQ=255)
                      aaaaaccgccaTAGAAAGGACCAAAGAATGACACtt            <  1:2627217‑M2/25‑1 (MQ=255)
                       ataACCGCCATAGAAAGGACCAAAGAATGACACttt           >  1:892671/1‑36 (MQ=255)
                       aaaaccgccaTAGAAAGGACCAAAGAATGACACttt           <  1:3416513‑M2/26‑1 (MQ=255)
                        taACCGCCATAGAAAGGACCAAAGAATGACACTTTc          >  1:7299351/1‑36 (MQ=255)
                        aacccgccaTAGAAAGGACCAAAGAATGACACTTTc          <  1:1611615‑M2/27‑1 (MQ=255)
                        aaaccgccaTAGAAAGGACCAAAGAATGACACTTTc          <  1:995655‑M2/27‑1 (MQ=255)
                        aaaccgccaTAGAAAGGACCAAAGAATGACACTTTc          <  1:379095‑M2/27‑1 (MQ=255)
                         aaccgccaTAGAAAGGACCAAAGAATGACACTTTca         <  1:5700996‑M2/28‑1 (MQ=255)
                         aaccgccaTAGAAAGGACCAAAGAATGACACTTTca         >  1:5310401‑M2/9‑36 (MQ=255)
                         aaccgccaTAGAAAGGACCAAAGAATGACACTTTca         <  1:3221591‑M2/28‑1 (MQ=255)
                          accgccaTAGAAAGGACCAAAGAATGACACTTTcag        >  1:2057214‑M2/8‑36 (MQ=255)
                           ccgccaTAGAAAGGACCAAAGAATGACACTTTcagc       <  1:4601609‑M2/30‑1 (MQ=255)
                             gccaTAGAAAGGACCAAAGAATGACACTTTcagcgc     >  1:7046702‑M2/5‑36 (MQ=255)
                               caTAGAAAGGACCAAAGAATGACACTTTcagcgcag   <  1:5621854‑M2/34‑1 (MQ=255)
                               caTAGAAAGGACCAAAGAATGACACTTTcagcgcag   <  1:5284100‑M2/34‑1 (MQ=255)
                               caTAGAAAGGACCAAAGAATGACACTTTcagcgcag   >  1:3314373‑M2/3‑36 (MQ=255)
                                aTAGAAAGGACCAAAGAATGACACTTTcagcgcagc  <  1:4689121‑M2/35‑1 (MQ=255)
                                aTAGAAAGGACCAAAGAATGACACTTTcagcgcagc  >  1:7177607‑M2/2‑36 (MQ=255)
                                aTAGAAAGGACCAAAGAATGACACTTTcagcgcagc  <  1:4230952‑M2/35‑1 (MQ=255)
                                 |                                  
TCCCGATCAAGTGCAATGACGTTATAACCGCCATAGAAAGGACCAAAGAATGACACTTTCAGCGCAGC  >  REL606/4356063‑4356130

Base quality scores:  ATCG  < 3 ≤  ATCG  < 14 ≤  ATCG  < 24 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG