New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 1946057 =38 (0.780)8 (0.160) 3/56 NT 22.5% coding (151/390 nt) cheY chemotaxis regulator transmitting signal to flagellar motor component
?REL606 = 2574503 17 (0.280)intergenic (+345/+473) sseA/sseB 3‑mercaptopyruvate sulfurtransferase/rhodanase‑like enzyme, sulfur transfer from thiosulfate

ATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACCGCCTGCCTGCAACTTATTGAGAGCGTCGACGC  >  REL606/1946022‑1946092
                                   |                                   
atTGGGCATGTTCCAGTCGGAGATAACAAATCCATa                                     >  1:4084572/1‑36 (MQ=255)
 ttGGGCATGTTCCAGTCGGAGATAACAAATCCATaa                                    >  1:7292838/1‑36 (MQ=255)
  tGGGCATGTTCCAGTCGGAGATAACAAATCCATAAc                                   >  1:5141584/1‑36 (MQ=255)
  tGGGCATGTTACAGTCGGAGATAACAAATCCATAAc                                   <  1:2875048/36‑1 (MQ=255)
   gggCATGTTCCAGTCGGAGATAACAAATCCATAAcc                                  >  1:4602925/1‑36 (MQ=255)
   gggCATGTTCCAGTCGGAGATAACAAATCCATAAcc                                  <  1:1999224/36‑1 (MQ=255)
    ggCATGTTCCAGTCGGAGATAACAAATCCATAACCg                                 <  1:7155227/36‑1 (MQ=255)
    ggCATGTTCCAGTCGGAGATAACAAATCCATAACCg                                 >  1:3444938/1‑36 (MQ=255)
    ggCATGTTCCAGTCGGAGATAACAAATCCATAACCg                                 <  1:2607040/36‑1 (MQ=255)
     gCATGTTCCAGTCGGAGATAACAAATCCATAACCgc                                >  1:4097611/1‑36 (MQ=255)
      cATGTTCCAGTCGGAGATAACAAATCCATAACCgcc                               >  1:3151289/1‑36 (MQ=255)
        tGTTCCAGTCGGAGATAACAAATCCATAACCgcctg                             <  1:430728/36‑1 (MQ=255)
        tGTTCCAGTCGGAAATAACAAATCCATAACCgcctg                             <  1:2362394/36‑1 (MQ=255)
        tGTTCAAGTCGGAGATAACAAATCCATAACCgcctg                             <  1:7695970/36‑1 (MQ=255)
          ttCCAGTCGGAGATAACAAATCCATAACCgcctgcc                           >  1:340544/1‑36 (MQ=255)
           cccAGTCGGAAAAAAAAAACCAATAACCgcctgcct                          <  1:1197621/35‑1 (MQ=12)
            ccAGTCGGAAAAAAAAAATCCAAAACCGcctgcctg                         <  1:7301956/36‑1 (MQ=16)
             cAGTCGGAGATAACAAATCCATAACCGCctgcctgc                        >  1:1108114/1‑36 (MQ=255)
              aGTCGGAGATAACAAATCCATAACCGCCTGCCTGCa                       >  1:3660512/1‑36 (MQ=255)
              aGTCGGAGATAACAAATCCATAACCGCCTGCCTGCa                       >  1:175331/1‑36 (MQ=255)
              aGTCGGAGATAACAAATCCATAACCGCCTGCCTGCa                       >  1:1514171/1‑36 (MQ=255)
              aGTCGGAGATAACAAATCCATAACCGCCTGCCTGCa                       >  1:6860580/1‑36 (MQ=255)
              aGTCGGAGATAACAAATACATAACCGCCTGCCTGCa                       >  1:1049754/1‑36 (MQ=255)
               gTCGGAGATAACAAATCCATAACCGCCTGCCTGCaa                      >  1:2095030/1‑36 (MQ=255)
               gTCGGAGATAACAAATCCATAACCGCCTGCCTGCaa                      >  1:4181539/1‑36 (MQ=255)
               gTCGGAGATAACAAATCCATAACCGCCTGCCTGCaa                      >  1:6806958/1‑36 (MQ=255)
                gcGGAAAAAAAAAACCAATAACCGCCTGCCTGCAAc                     <  1:2497534/35‑1 (MQ=9)
                 cGGAGATAACAAATCCATAACCGCCTGCCTGCAACt                    >  1:1031693/1‑36 (MQ=255)
                  ggAGATAACAAATCCATAACCGCCTGCCTGCAACtt                   >  1:1920003/1‑36 (MQ=255)
                   gagaTAACAAATCCATAACCGCCTGCCTGCAACTTa                  >  1:3980633/1‑36 (MQ=255)
                   gagaTAACAAATCCACAACCCCCCCCCCGCACCTTa                  >  1:5742069/1‑36 (MQ=12)
                     gaTAACAAATCCATAACCGCCTGCCTGCAACTTAtt                >  1:4183630/1‑36 (MQ=255)
                     gaTAACAAATCCATAACCGCCTGCCTGCAACTTAtt                >  1:2316908/1‑36 (MQ=255)
                     gaTAACAAATCCATAACCGCCTGCCTGCAACTTAtt                >  1:6233211/1‑36 (MQ=255)
                      aTAAAAAATCCATAACTGCCTGCCTGCAACTTATTg               >  1:4746455/1‑36 (MQ=37)
                        aaCAAATCCATAACCGCCTGCCTGCAACTTATTgag             >  1:6714636/1‑36 (MQ=255)
                         aaaaataaaaAACCGCCTGCCTGCAACTTATTgaga            >  1:7310442‑M1/11‑36 (MQ=255)
                         aaaaaccaaaAACCGCCTGCCTGCAACTTATTgaga            >  1:7834296‑M1/11‑36 (MQ=255)
                         aaaaaTCCACCAACGCCTGCCTGCAACTTATTgaga            <  1:6340124/34‑1 (MQ=18)
                         aCAAATCCATAACCGCCTGCCTGCAACTTATTgagt            >  1:1436995/1‑35 (MQ=255)
                         aCAAATCCATAACCGCCTGCCTGCAACTTATTgaga            >  1:5638731/1‑36 (MQ=255)
                          cAAATCCATAACCGCCTGCCTGCAACTTATTGagag           >  1:349150/1‑36 (MQ=255)
                          aaaacaaaaAACCGCCTGCCTGCAACTTATTGagag           >  1:6828849‑M1/10‑36 (MQ=255)
                           aaaTCCATAACCGCCTGCCTGCAACTTATTGAGAGc          >  1:3594558/1‑36 (MQ=255)
                           aaaCCCATAACCGCCTGCCTGCAACTTATTGAGAGc          <  1:2171379/36‑1 (MQ=255)
                             aTCCATAACCGCCTGCCTGCAACTTATTGAGAGCGt        >  1:5970427/1‑36 (MQ=255)
                               caacAACCGCCTGCCTGCAACTTATTGAGAGCGTCg      >  1:7475400‑M1/5‑36 (MQ=255)
                               caaaAACCGCCTGCCTGCAACTTATTGAGAGCGTCg      >  1:6513445‑M1/5‑36 (MQ=255)
                                cATAACCGCCTGCCTGCAACTTATTGAGAGCGTCGa     >  1:3971697/1‑36 (MQ=255)
                                cATAACCGCCTGCCTGCAACTTATTGAGAGCGTCGa     >  1:4806420/1‑36 (MQ=255)
                                cATAACCGACTGCCGGCAACTTATTGAGAGCGTCGa     >  1:5855084/1‑36 (MQ=37)
                                 aaACCCGCCTGCCTGCAACTTATTGAGAGCGTCGAc    >  1:1453534‑M1/3‑36 (MQ=255)
                                 aaACCCGCCTGCCTGCAACTTATTGAGAGCGTCGAc    >  1:5471989‑M1/3‑36 (MQ=255)
                                 aaAACCGCCTGCCTGCAACTTATTGAGAGCGTCGAc    >  1:832361‑M1/3‑36 (MQ=255)
                                 aTAACCGCCTGCCTGCAACTTATTGAGAGCGTCGAc    >  1:4155738/1‑36 (MQ=255)
                                 aTAACCGCCTGCCTGCAACTTATTGAGAGCGTCGAc    >  1:7797279/1‑36 (MQ=255)
                                  tAACCGCCTGCCTGCAACTTATTGAGAGCGTCGACg   >  1:2205463/1‑36 (MQ=255)
                                   AACCGCCTGCCTGCAACTTATTGAGAGCGTCGACGc  <  1:1438505‑M1/36‑1 (MQ=255)
                                   AACCGCCTGCCTGCAACTTATTGAGAGCGCCGACGc  <  1:5821357‑M1/36‑1 (MQ=255)
                                   |                                   
ATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACCGCCTGCCTGCAACTTATTGAGAGCGTCGACGC  >  REL606/1946022‑1946092

Base quality scores:  ATCG  < 3 ≤  ATCG  < 9 ≤  ATCG  < 20 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG