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

AAAAAACAAAAAACCGCCA...................................  >  REL606/2574493‑2574511
...........AACCGCCATAGAAAGGACCAAAGAATGACACTTTCAGCGCAGC  >  REL606/4356096‑4356130
                                                      
AAAAAAAAAAAAACCGCCATAGAAAGGACCAAAGAA                    <  1:3458943/36‑1 (MQ=255)
AAAAAAAAAAAAACCGCCATAGAAAGGACCAAAGAA                    <  1:5602230/36‑1 (MQ=255)
 AAAAAAAAAAAACCGCCATAGAAAGGACCAAAGAAT                   <  1:4823930/36‑1 (MQ=255)
     aaAAAAAACCGCCATAGAAAGGACCAAAGAATGACA               <  1:2424060/34‑1 (MQ=17)
        AAAAACCGCCATAGAAAGGACCAAAGAATGACACTT            <  1:2627217/36‑1 (MQ=1)
         AAAACCGCCATAGAAAGGACCAAAGAATGACACTTT           <  1:3416513/36‑1 (MQ=1)
          AAACCGCCATAGAAAGGACCAAAGAATGACACTTTC          <  1:379095/36‑1 (MQ=1)
          AAACCGCCATAGAAAGGACCAAAGAATGACACTTTC          <  1:995655/36‑1 (MQ=1)
                                                      
AAAAAACAAAAAACCGCCA...................................  >  REL606/2574493‑2574511
...........AACCGCCATAGAAAGGACCAAAGAATGACACTTTCAGCGCAGC  >  REL606/4356096‑4356130

Base quality scores:  ATCG  < 3 ≤  ATCG  < 6 ≤  ATCG  < 15 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG