New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 1241529 =31 (0.670)7 (0.150) 4/56 NT 21.5% coding (693/1737 nt) ycgO potassium/proton antiporter
?REL606 = 2574504 20 (0.350)intergenic (+346/+472) sseA/sseB 3‑mercaptopyruvate sulfurtransferase/rhodanase‑like enzyme, sulfur transfer from thiosulfate

CCTTCCAGCGCAGTAGTTAATGAGAAAATCAGAATACCGCCACTTAATGCCAGCAATGGATATAATCCGG  >  REL606/1241494‑1241563
                                   |                                  
ccTTCCAGCGCAGTAGTTAATGAGAAAATCAGAATa                                    <  1:2020118/36‑1 (MQ=255)
 cTTCCAGCGCAGTAGTTAATGAGAAAATCAGAATAc                                   <  1:1128290/36‑1 (MQ=255)
  ttCCAGCGCAGTAGTTAATGATAAAATCAGAATAcc                                  <  1:6636922/36‑1 (MQ=255)
  ttCCAGCGCAGTAGTTAATGAGAAAATCAGAATAcc                                  >  1:6490982/1‑36 (MQ=255)
  ttCCAGCGCAGTAGTTAATGAGAAAATCAGAATAcc                                  >  1:1659459/1‑36 (MQ=255)
  ttCCAGCGCAGTAGTTAATGAGAAAATCAGAATAcc                                  <  1:6626448/36‑1 (MQ=255)
   tCCAGCGCAGTAGTTAATGAGAAAATCAGAATACCg                                 >  1:1441869/1‑36 (MQ=255)
    ccAGCGCAGTAGTTAATGAGAAAATCAGAATACCGc                                >  1:2758977/1‑36 (MQ=255)
      aGCGAAGTAGTTAATGAGAAAATCAAAATACCGCCa                              <  1:594623/36‑1 (MQ=37)
       gcgcAGTAGTTAATGAGAAAATCAGAATACCGCCAc                             >  1:7473527/1‑36 (MQ=255)
         gcAGTAGTTAATGAGAAAATCAGAATACCGCCACtt                           >  1:3066840/1‑36 (MQ=255)
         gcAGTAGTTAATGAGAAAATCAGAATACCGCCACtt                           >  1:3636109/1‑36 (MQ=255)
         gcAGAAGTTAATGAGAAAATCAGAATAACGCCACtt                           >  1:2306755/1‑36 (MQ=37)
            gtagtTAATGAGAAAATCAGAATACCGCCACTTAAt                        >  1:7159995/1‑36 (MQ=255)
            gtagtTAATGAGAAAATCAGAATACCGCCAATTAAt                        >  1:7400826/1‑36 (MQ=255)
              agtTAATGAGAAAATCAGAATACCGCCACTTAATGc                      >  1:7463263/1‑36 (MQ=255)
              agtTAATGAGAAAATCAGAATACCGCCACTTAATGc                      >  1:973813/1‑36 (MQ=255)
               gtTAATGAGAAAATCAGAATACCGCCACTTAATGcc                     >  1:369714/1‑36 (MQ=255)
                ttAATGAGAAAATCAGAATACCGCCACTTAATGCCa                    >  1:5060084/1‑36 (MQ=255)
                  aaTGAGAAAATCAGAATACCGCCACTTAATGCCAGc                  >  1:1540332/1‑36 (MQ=255)
                     gagaAAATCAGAATACCGCCACTTAATGCCAGCAAt               >  1:2737056/1‑36 (MQ=255)
                     gagaAAATCAGAATACCGCCACTTAATGCCAGCAAt               >  1:271529/1‑36 (MQ=255)
                     aaaaaaaTAAAAATACCGCCACTTAATGCCAGCAAt               <  1:4531096/33‑1 (MQ=21)
                      agaAAATCAGAATACCGCCACTTAATGCCAGCAATg              >  1:6372813/1‑36 (MQ=255)
                      aaaaaaaaaaaaaACCGCCACTTAATGCCAGCAATg              >  1:531457‑M1/14‑36 (MQ=255)
                       gaAAATCAGAATACCGCCACTTAATGCCAGCAATgg             >  1:4222741/1‑36 (MQ=255)
                         aaaaaaaaaaACCGCCACTTAATGCCAGCAATGGat           >  1:6371394‑M1/11‑36 (MQ=255)
                         aaaaaaaaaaACCGCCACTTAATGCCAGCAATGGat           >  1:416727‑M1/11‑36 (MQ=255)
                         aaaTCAGAATACCGCCACTTAATGCCAGCAATGGat           >  1:5177468/1‑36 (MQ=255)
                          aaaaaaaaaACCGCCACTTAATGCCAGCAATGGata          >  1:3225833‑M1/10‑36 (MQ=255)
                            tCAGAATACCGCCACTTAATGCCAGCAATGGAtata        >  1:3461096/1‑36 (MQ=255)
                            aaaaaaaACCGCCACTTAATGCCAGCAATGGAtata        >  1:381474‑M1/8‑36 (MQ=255)
                               gAATACCGCCACTTAATGCCAGCAATGGATATAATc     >  1:5047176/1‑36 (MQ=255)
                               aaaaACCGCCACTTAATGCCAGCAATGGATATAATc     >  1:5905322‑M1/5‑36 (MQ=255)
                                aaaCCCGCCACTTAATGCCAGCAATGGATATAATcc    >  1:1774252‑M1/4‑36 (MQ=255)
                                aaTACCGCCACTTAATGCCAGCAATGGATATAATcc    >  1:4119623/1‑36 (MQ=255)
                                 aTACCGCCACTTAATGCCAGCAATGGATATAATCcg   >  1:6449601/1‑36 (MQ=255)
                                  tACCGCCACTTAATGCCAGCAATGGATATAATCcgg  >  1:401179/1‑36 (MQ=255)
                                  tACCGCCACTTAATGCCAGCAATGGATATAATCcgg  >  1:1028815/1‑36 (MQ=255)
                                  cACCGCCACTTAATGCCAGCAATGGATATAATCcgg  >  1:5638193‑M1/2‑36 (MQ=255)
                                   |                                  
CCTTCCAGCGCAGTAGTTAATGAGAAAATCAGAATACCGCCACTTAATGCCAGCAATGGATATAATCCGG  >  REL606/1241494‑1241563

Base quality scores:  ATCG  < 3 ≤  ATCG  < 13 ≤  ATCG  < 24 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG