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

AGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTA  >  REL606/1474761‑1474831
                                   |                                   
aGCGTTGCGCGTACTGTGGTCATACAGCGAAAgaga                                     >  1:6927099/1‑36 (MQ=255)
 gCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGaa                                    >  1:3734879/1‑36 (MQ=255)
 gCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGaa                                    >  1:2823069/1‑36 (MQ=255)
  cGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAAc                                   <  1:4057043/36‑1 (MQ=255)
   gTTGCGCGTACTGTGGTCATACAGCGAAAGAGAAcc                                  >  1:3872546/1‑36 (MQ=255)
   gTTGCGCGTACTGTGGTCATACAGCGAAAGAGAAcc                                  >  1:6143696/1‑36 (MQ=255)
   gTTGCGCGTACTGTGGTCATACAGCGAAAGAGAAcc                                  <  1:1268537/36‑1 (MQ=255)
         cgTACTGTGGTCATACAGCGAAAGAGAACCGCCTGt                            >  1:2488526/1‑36 (MQ=255)
          gTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTc                           >  1:7922611/1‑36 (MQ=255)
          gTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTc                           >  1:8262670/1‑36 (MQ=255)
          gTACTGGGGTCATACAGCGAAAGAGAACCGCCTGTc                           >  1:4620122/1‑36 (MQ=255)
           tACTGTGGTCATACAGCGAAAGAGAACCGCCTGTca                          >  1:1974532/1‑36 (MQ=255)
           tACTGTGGTCATACAGCGAAAGAGAACCGCCTGTca                          >  1:1368827/1‑36 (MQ=255)
           tACTGTGGTCATACAGCGAAAGAGAACCGCCTGTca                          >  1:136759/1‑36 (MQ=255)
            aCTGTGGTCATACAGCGAAAGAGAACCGCCTGTcac                         >  1:7815445/1‑36 (MQ=255)
            aCTGTGGTCATACAGCGAAAGAGAACCGCCTGTcac                         >  1:1660291/1‑36 (MQ=255)
            aCTGTGGTCATACAGCGAAAGAGAACCGCCTGTcac                         >  1:2865040/1‑36 (MQ=255)
                 ggTCATACAGCGAAAGAGAACCGCCTGTCACAAAGt                    >  1:4692492/1‑36 (MQ=255)
                  gTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTa                   >  1:1013170/1‑36 (MQ=255)
                   tCATACAGCGAAAGAGAACCGCCTGTCACAAAGTaa                  >  1:3645098/1‑36 (MQ=255)
                    cATACAGCGAAAGAGAACCGCCTGTCACAAAGTaaa                 >  1:4126197/1‑36 (MQ=255)
                    cATACAGCGAAAGAGAACCGCCTGTCACAAAGTaaa                 >  1:5236709/1‑36 (MQ=255)
                     aTACAGCGAAAGAGAACCGCCTGTCACAAAGTAAAt                >  1:3560871/1‑36 (MQ=255)
                     aTACAGCGAAAGAGAACCGCCTGTCACAAAGTAAAt                >  1:4919032/1‑36 (MQ=255)
                     aTACAGCGAAAGAGAACCGCCTGTCACAAAGTAAAt                >  1:8153063/1‑36 (MQ=255)
                      tACAGCGAAAGAGAACCGCCTGTCACAAAGTAAAtt               >  1:1702807/1‑36 (MQ=255)
                      tACAGCGAAAGAGAACCGCCTGTCACAAAGTAAAtt               >  1:4975748/1‑36 (MQ=255)
                        cAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCa             >  1:808759/1‑36 (MQ=255)
                        aaaaaaaaaaaAACCGCCTGTCACAAAGTAAATTCa             >  1:3252829‑M1/12‑36 (MQ=255)
                         aGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAg            >  1:2555823/1‑36 (MQ=255)
                          gaaaaaGAAAACCGCCTGTCACAAAGTAAATTCAGa           <  1:893468/33‑1 (MQ=37)
                          gCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGa           >  1:8257296/1‑36 (MQ=255)
                          gCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGa           >  1:2880625/1‑36 (MQ=255)
                          aaaaaaaaaAACCGCCTGTCACAAAGTAAATTCAGa           >  1:1455998‑M1/10‑36 (MQ=255)
                          aaaaaaaaaAACCGCCTGTCACAAAGTAAATTCAGa           >  1:4724606‑M1/10‑36 (MQ=255)
                            gAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATg         >  1:2547628/1‑36 (MQ=255)
                               agagAACCGCCTGTCACAAAGTAAATTCAGATGcaa      >  1:8079169/1‑34 (MQ=255)
                               agagAACCGCCTGTCACAAAGTAAATTCAGATGCCa      >  1:921615/1‑36 (MQ=255)
                                aaaAACCGCCTGTCACAAAGTAAATTCAGATGCCAg     >  1:6445006‑M1/4‑36 (MQ=255)
                                 agaACCGCCTGTCACAAAGTAAATTCAGATGCCAgg    >  1:2461738/1‑36 (MQ=255)
                                 agaACCGCCTGTCACAAAGTAAATTCAGATGCCAgg    >  1:6052910/1‑36 (MQ=255)
                                 aaACCCGCCTGTCACAAAGTAAATTCAGATGCCAgg    >  1:2486346‑M1/3‑36 (MQ=255)
                                 aaAACCGCCTGTCACAAAGTAAATTCAGATGCCAgg    >  1:3425720‑M1/3‑36 (MQ=255)
                                  gaACCGCCTGTCACAAAGTAAATTCAGATGCCAGGt   >  1:851069/1‑36 (MQ=255)
                                   aaCCGCCGGTCACAAAGTAAATGCAGATGCCAGGTa  >  1:7993798/1‑36 (MQ=37)
                                   AACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTa  <  1:3040296‑M1/36‑1 (MQ=255)
                                   AACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTa  >  1:3163795‑M1/1‑36 (MQ=255)
                                   AACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTa  >  1:3336969‑M1/1‑36 (MQ=255)
                                   AACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTa  >  1:964309‑M1/1‑36 (MQ=255)
                                   |                                   
AGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTA  >  REL606/1474761‑1474831

Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 25 ≤  ATCG  < 31 ≤  ATCG  < 37 ≤  ATCG