New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 1474796 =36 (0.740)9 (0.190) 6/56 NT 25.7% coding (1030/1209 nt) ydcM predicted transposase
?REL606 = 2574503 16 (0.260)intergenic (+345/+473) sseA/sseB 3‑mercaptopyruvate sulfurtransferase/rhodanase‑like enzyme, sulfur transfer from thiosulfate

AGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTA  >  REL606/1474761‑1474831
                                   |                                   
aGCGTTGCGCGTACTGTGGTCATACAGCGAAAgaga                                     <  1:4435646/36‑1 (MQ=255)
aGCGTTGCGCGTACTGTGGTCATACAGCGAAAgaga                                     >  1:662367/1‑36 (MQ=255)
 gCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGaa                                    <  1:7458153/36‑1 (MQ=255)
 gCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGaa                                    >  1:3126782/1‑36 (MQ=255)
  cGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAAc                                   <  1:4597019/36‑1 (MQ=255)
   gTTGCGCGTACTGTGGTCATACAGCGAAAGAGAAcc                                  <  1:2892694/36‑1 (MQ=255)
   catgtaaattcactttcctttaaaaaacaaaaAAcc                                  <  1:1566169‑M1/4‑1 (MQ=255)
   catgtaaattcactttcctttaaaaaacaaaaAAcc                                  >  1:2554732‑M1/33‑36 (MQ=255)
   catgtaaattcactttcctttaaaaaacaaaaAAcc                                  >  1:2870012‑M1/33‑36 (MQ=255)
    ttGCGCGTACTGTGGTCATACAGCGAAAGAGAACCg                                 <  1:3884326/36‑1 (MQ=255)
    atgtaaattcactttcctttaaaaaacaaaaAACCg                                 <  1:1055674‑M1/5‑1 (MQ=255)
    atgtaaattcactttcctttaaaaaacaaaaAACCg                                 >  1:2131311‑M1/32‑36 (MQ=255)
     tgtaaattcactttcctttaaaaaacaaaaAACCGc                                <  1:3382430‑M1/6‑1 (MQ=255)
       cgcgTACTGTGGTCATACAGCGAAAGAGAACCGCCt                              >  1:4632158/1‑36 (MQ=255)
       cgcgTACTGTGGTCATACAGCGAAAGAGAACCGCCt                              >  1:5044994/1‑36 (MQ=255)
         cgTACTGTGGTCATACAGCGAAAGAGAACCGCCTGt                            >  1:6113032/1‑36 (MQ=255)
         cgTACTGTGGTCATACAGCGAAAGAGAACCGCCTGt                            >  1:5214875/1‑36 (MQ=255)
          gTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTc                           >  1:4945089/1‑36 (MQ=255)
           tACTGTGGTCATACAGCGAAAGAGAACCGCCTGTca                          >  1:5282349/1‑36 (MQ=255)
           tACTGTGGTCATACAGCGAAAGAGAACCGCCTGTca                          >  1:4868502/1‑36 (MQ=255)
           tACTGTGGTCATACAGCGAAAGAGAACCGCCTGTca                          >  1:7694459/1‑36 (MQ=255)
             cTGTGGTCATACAGCGAAAGAGAACCGCCTGTcaca                        >  1:3309410/1‑36 (MQ=255)
              tgtgGTCATACAGCGAAAGAGAACCGCCTGTCACaa                       >  1:5533635/1‑36 (MQ=255)
               gtgGTCATACAGCGAAAGAGAACCGCCTGTCACaaa                      >  1:1313422/1‑36 (MQ=255)
                tgGTCATACAGCGAAAGAGAACCGCCTGTCACAAAg                     >  1:1319742/1‑36 (MQ=255)
                tgGTCATACAGCGAAAGAGAACCGCCTGTCACAAAg                     >  1:5945224/1‑36 (MQ=255)
                 ggTCATACAGCGAAAGAGAACCGCCTGTCACAAAGt                    >  1:2167439/1‑36 (MQ=255)
                  gTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTa                   >  1:2721724/1‑36 (MQ=255)
                   tCATACAGCGAAAGAGAACCGCCTGTCACAAAGTaa                  >  1:3218985/1‑36 (MQ=255)
                     aTACAGCGAAAGAGAACCGCCTGTCACAAAGTAAAt                >  1:5564018/1‑36 (MQ=255)
                     aTACAGCGAAAGAGAACCGCCTGTCACAAAGTAAAt                >  1:7253456/1‑36 (MQ=255)
                     aTACAGCGAAAGAGAACCGCCTGTCACAAAGTAAAt                >  1:3490466/1‑36 (MQ=255)
                     aTACAGCGAAAGAGAACCGCCTGTCACAAAGTAAAt                >  1:2939075/1‑36 (MQ=255)
                      tACAGCGAAAGAGAACCGCCTGTCACAAAGTaaagt               >  1:2744444/1‑34 (MQ=255)
                      tACAGCGAAAGAGAACCGCCTGTCACAAAGTAAAtt               >  1:354640/1‑36 (MQ=255)
                       aCAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTc              >  1:5898785/1‑36 (MQ=255)
                        aaaaaaaaaaaAACCGCCTGTCACAAAGTAAATTCa             >  1:5694973‑M1/12‑36 (MQ=255)
                         aGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAg            >  1:7875483/1‑36 (MQ=255)
                          gCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGa           >  1:5172280/1‑36 (MQ=255)
                          aaaaaaaaaAACCGCCTGTCACAAAGTAAATTCAGa           >  1:937353‑M1/10‑36 (MQ=255)
                           aaaaaaaaAACCGCCTGTCACAAAGTAAATTCAGAt          >  1:2622782‑M1/9‑36 (MQ=255)
                           aaaaaAAGAACCGCCTGTCACAAAGTAAATTCAGAt          <  1:980487/34‑1 (MQ=37)
                            gaaaaaaAACCGCCTGTCACAAAGTAAATTCAGATg         >  1:7772580‑M1/8‑36 (MQ=255)
                            aaaaaaaAACCGCCTGTCACAAAGTAAATTCAGATg         >  1:1717926‑M1/8‑36 (MQ=255)
                             aaaGAGAACCGCCTGTCACAAAGTAAATTCAGATGc        >  1:2788200/1‑36 (MQ=255)
                              aaaaaAACCGCCTGTCACAAAGTAAATTCAGATGcc       >  1:4483207‑M1/6‑36 (MQ=255)
                              aaGAGAACCGCCTGTAACAAAGTAAATTCAGATGcc       >  1:4241695/1‑36 (MQ=255)
                               agagAACCGCCTGTCACAAAGTAAATTCGGATGCCa      >  1:7335700/1‑36 (MQ=255)
                               agagAACCGCCTGTCACAAAGTAAATTCAGATGCCa      >  1:4166078/1‑36 (MQ=255)
                               agagAACCGCCTGTCACAAAGTAAATTCAGATGCCa      >  1:223224/1‑36 (MQ=255)
                                aacAACCGCCTGTCACAAAGTAAATTCAGATGCCAg     >  1:1734542‑M1/4‑36 (MQ=255)
                                aaaAACCGCCTGTCACAAAGTAAATTCAGATGCCAg     >  1:4374847‑M1/4‑36 (MQ=255)
                                  aAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGt   >  1:3841867‑M1/2‑36 (MQ=255)
                                   AACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTa  <  1:3903245‑M1/36‑1 (MQ=255)
                                   |                                   
AGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTA  >  REL606/1474761‑1474831

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 26 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG