New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 1474796 =31 (0.630)6 (0.120) 4/56 NT 18.5% coding (1030/1209 nt) ydcM predicted transposase
?REL606 = 2574503 22 (0.360)intergenic (+345/+473) sseA/sseB 3‑mercaptopyruvate sulfurtransferase/rhodanase‑like enzyme, sulfur transfer from thiosulfate

AAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTC  >  REL606/2574477‑2574538
                          |                                   
aaaTTCACTTTCATTAAAAAAAAAAAAAACCGCCaa                            <  1:2056995/36‑1 (MQ=25)
  aTTCACTTTCCTTTAAAAAACAAAAAACCGCCaaaa                          >  1:1247007/1‑36 (MQ=255)
       cTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAgg                     >  1:813171/1‑36 (MQ=255)
        tttCCTTTAAAAAACAAAAAACCGCCAAAATCAGGc                    >  1:3844853/1‑36 (MQ=255)
        tttCCTTTAAAAAACAAAAAACCGCCAAAATCAGGc                    >  1:715051/1‑36 (MQ=255)
        tttCCTTTAAAAAACAAAAAACCGCCAAAATCAGGc                    >  1:2059516/1‑36 (MQ=255)
         ttCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCg                   >  1:7527996/1‑36 (MQ=255)
              aaaaaaaaaaaaaaacCGCCAAAATCAGGCGGtttt              <  1:1926217/34‑1 (MQ=25)
                agaaaaaaaaaaacCGCC.GAATTTGGCGGtttttta           <  1:1454297/34‑2 (MQ=17)
                aaaaaaaaaaaaaccgcctgtcacaaagtaaattca            >  1:3252829‑M2/1‑11 (MQ=255)
                aaaaaaaaaaaaacCGCCAAATTC.GGCGGtttttta           <  1:4168538/36‑2 (MQ=17)
                aaaaaaaaaaaaacCGCCAAAATCAGGCGGtttttt            <  1:915082/36‑1 (MQ=255)
                  aaaaaaaaaaaccgcctgtcacaaagtaaattcaga          >  1:4724606‑M2/1‑9 (MQ=255)
                  aaaaaaaaaaaccgcctgtcacaaagtaaattcaga          >  1:1455998‑M2/1‑9 (MQ=255)
                  aaaaaaaaaaacCGCCAAAATCAGGCGGTTTTttgt          <  1:94583/36‑1 (MQ=255)
                   aaaaaaaaaacCGCCAAAATCAGGCGGTTTTttgtt         <  1:1494628/36‑1 (MQ=255)
                    aaaaaaaaacCGCCAAAATCAGGCGGTTTTttgttg        <  1:6957349/33‑1 (MQ=255)
                    aaaaaaaaacCGCCAAAATCAGGCGGTTTTttgttg        <  1:1202871/33‑1 (MQ=255)
                     aaaaaaaacCGCCAAAATCAGGCGGTTTTTTGTTGc       <  1:5115937/34‑1 (MQ=255)
                      caaaaaacCGCCAAAATCAGGCGGGTTTTTGTTTCt      >  1:4017061/1‑36 (MQ=25)
                       aaaaaacCGCCAAATTCGGCGGGTTTTTTTTTTCtt     >  1:3197719/1‑35 (MQ=17)
                        aaaaaccgcctgtcacaaagtaaattcagatgccag    >  1:6445006‑M2/1‑3 (MQ=255)
                        aaaaacCGCCAAAATCAGGCGGTTTTTTTTTGCTgg    >  1:2700602/1‑36 (MQ=255)
                         aaacccgcctgtcacaaagtaaattcagatgccagg   >  1:2486346‑M2/1‑2 (MQ=255)
                         aaaaccgcctgtcacaaagtaaattcagatgccagg   >  1:3425720‑M2/1‑2 (MQ=255)
                         aaaacCGCCAAAATCAGGCGGTTTTTTGTTGCTGGt   <  1:2526855/36‑1 (MQ=255)
                         aaaacCGCCAAAATCAGGCGGTTTTTTGTTGCTGGt   <  1:1577395/36‑1 (MQ=255)
                          aaacCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTc  <  1:7334547/36‑1 (MQ=255)
                          |                                   
AAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTC  >  REL606/2574477‑2574538

Base quality scores:  ATCG  < 3 ≤  ATCG  < 11 ≤  ATCG  < 26 ≤  ATCG  < 33 ≤  ATCG  < 37 ≤  ATCG