New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 257451125 (0.520)6 (0.120) 5/56 NT 15.6% intergenic (+353/+465) sseA/sseB 3‑mercaptopyruvate sulfurtransferase/rhodanase‑like enzyme, sulfur transfer from thiosulfate
?REL606 4395051 = 40 (0.670)intergenic (‑126/‑56) yjfO/yjfP hypothetical protein/predicted hydrolase

AACCCACATCTCAACATTATGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATAC  >  REL606/4395016‑4395086
                                   |                                   
aaCCCACATCTCAACATTATGAACTAATACCGCcag                                     >  1:6275284/1‑36 (MQ=255)
 aCCCACATCTCAACATTATGAACTAATACCGCcagc                                    >  1:3014478/1‑36 (MQ=255)
 aCCCACATCTCAACATTATGAACTAATACCGCcagc                                    >  1:3002061/1‑36 (MQ=255)
  cccACATCTCAACATTATGAACTAATACCGCcagca                                   >  1:6674417/1‑36 (MQ=255)
  cccACATCTCAACATTATGAACTAATACCGCcagca                                   >  1:7478406/1‑36 (MQ=255)
    cacaTCTCAACATTATGAACTAATACCGCCAGCAGa                                 >  1:6430959/1‑36 (MQ=255)
    cacaTCTCAACATTATGAACTAATACCGCCAGCAGa                                 >  1:6331275/1‑36 (MQ=255)
    aacaTCTAAAAATTATAAACTAATACCGCCAGCAGa                                 <  1:7299691/35‑1 (MQ=21)
     acaTCTCAACATTATGAACTAATACCGCCAGCAGaa                                >  1:5447049/1‑36 (MQ=255)
     acaTCTCAACATTATGAACTAATACCGCCAGCAGaa                                >  1:5359771/1‑36 (MQ=255)
     acaTCTCAACATTATGAACTAATACCGCCAGCAGaa                                >  1:6944628/1‑36 (MQ=255)
         ctcAACATTATGAACTAATACCGCCAGCAGAAACtt                            >  1:2742735/1‑36 (MQ=255)
         ctcAACATTATGAACTAATACCGCCAGCAGAAACtt                            >  1:1793424/1‑36 (MQ=255)
          tcAACATTATGAACTAATACCGCCAGCAGAAACTtg                           >  1:3668173/1‑36 (MQ=255)
           cAACATTATGAACTAATACCGCCAGCAGAAACTtgt                          >  1:4425997/1‑36 (MQ=255)
           cAACATTATGAACTAATACCGCCAGCAGAAACTtgt                          >  1:2900592/1‑36 (MQ=255)
           cAACATTATGAACTAATACCGCCAGCAGAAACTtgt                          >  1:5500706/1‑36 (MQ=255)
           cAACATTATGAACTAATACCGCCAGCAGAAACTtgt                          >  1:4595514/1‑36 (MQ=255)
           cAACATTATGAACTAATACCGCCAGCAGAAACTtgt                          >  1:4766452/1‑36 (MQ=255)
            aaCATTATGAACTAATACCGCCAGCAGAAACTtgtg                         >  1:6282703/1‑36 (MQ=255)
             aCATTATGAACTAATACCGCCAGCAGAAACTTGTGc                        >  1:3637952/1‑36 (MQ=255)
             aCATTATGAACTAATACCGCCAGCAGAAACTTGTGc                        >  1:1909052/1‑36 (MQ=255)
             aCATTATGAACTAATACCGCCAGCAGAAACTTGTGc                        >  1:6440852/1‑36 (MQ=255)
              cATTATGAACTAATACCGCCAGCAGAAACTTGTGCt                       >  1:7469196/1‑36 (MQ=255)
               aTTATGAACTAATACCGCCAGCAGAAACTTGTGCTa                      >  1:7716825/1‑36 (MQ=255)
               aTTATGAACTAATACCGCCAGCAGAAACTTGTGCTa                      >  1:1987238/1‑36 (MQ=255)
                aaaaaaaaaaaaaccgccaGCAGAAACTTGTGCTaa                     <  1:5743810‑M2/17‑1 (MQ=255)
                 tATGAACTAATACCGCCAGCAGAAACTTGTGCTAAg                    >  1:2960273/1‑36 (MQ=255)
                 tATGAACTAATACCGCCAGCAGAAACTTGTGCTAAg                    >  1:7844912/1‑36 (MQ=255)
                 aaaaaaaaaaaaccgccaGCAGAAACTTGTGCTAAg                    <  1:4104177‑M2/18‑1 (MQ=255)
                  aaaaaaaaaaaccgccaGCAGAAACTTGTGCTAAgg                   <  1:5339383‑M2/19‑1 (MQ=255)
                  aTGAACTAATACCGCCAGCAGAAACTTGTGCTAAgg                   >  1:357201/1‑36 (MQ=255)
                  aTGAACTAATACCGCCAGCAGAAACTTGTGCTAAgg                   >  1:7645703/1‑36 (MQ=255)
                   tGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGa                  >  1:4752146/1‑36 (MQ=255)
                   tGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGa                  >  1:5091006/1‑36 (MQ=255)
                   tGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGa                  >  1:6249667/1‑36 (MQ=255)
                    gAACTAATACCGCCAGCAGAAACTTGTGCTAAGGAt                 >  1:3566216/1‑36 (MQ=255)
                     aaaaaaaaccgccaGCAGAAACTTGTGCTAAGGAtt                <  1:4115681‑M2/22‑1 (MQ=255)
                     aaCTAATACCGCCAGCAGAAACTTGTGCTAAGGAtt                >  1:1674514/1‑36 (MQ=255)
                      aCTAATACCGCCAGCAGAAACTTGTGCTAAGGAttt               >  1:4117700/1‑36 (MQ=255)
                       cTAATACCGCCAGCAGAAACTTGTGCTAAGGATttc              >  1:7007051/1‑36 (MQ=255)
                       aaaaaaccgccaGCAGAAACTTGTGCTAAGGATttc              <  1:6882549‑M2/24‑1 (MQ=255)
                        tAATACCGCCAGCAGAAACTTGTGCTAAGGATttct             >  1:5417405/1‑36 (MQ=255)
                         aaaaccgccaGCAGAAACTTGTGCTAAGGATttctt            <  1:1599010‑M2/26‑1 (MQ=255)
                          aTACCGCCAGCAGAAACTTGTGCTAAGGATttcttc           <  1:6296937/36‑1 (MQ=255)
                          aTACCGCCAGCAGAAACTTGTGCTAAGGATttctta           >  1:2387223/1‑35 (MQ=255)
                            cccGCCATCAGAAACTTGTGCTAAGTATTTCTTCAt         <  1:5278958/35‑1 (MQ=37)
                            accgccaGCAGAAACTTGTGCTAAGGATTTCTTCAt         >  1:5610708‑M2/8‑36 (MQ=255)
                            accgccaGCAGAAACTTGTGCTAAGGATTTCTTCAt         >  1:4375433‑M2/8‑36 (MQ=255)
                            accgccaGCAGAAACTTGTGCTAAGGATTTCTTCAt         >  1:3710895‑M2/8‑36 (MQ=255)
                            accgccaGCAGAAACTTGCGCTAAGGATTTCTTCAt         >  1:4794151‑M2/8‑36 (MQ=255)
                              cgccaGCAGAAACTTGTGCTAAGGATTTCTTCATGt       <  1:6190677‑M2/31‑1 (MQ=255)
                               gccaGCAGAAACTTGTGCTAAGGATTTCTTCATGTc      >  1:454975‑M2/5‑36 (MQ=255)
                               gccaGCAGAAACTTGTGCTAAGGATTTCTTCATGTc      <  1:7668201‑M2/32‑1 (MQ=255)
                               gccaGCAGAAACTTGTGCTAAGGATTTCTTCATGTc      <  1:2803911‑M2/32‑1 (MQ=255)
                               gcaaGCAGAAACTTGTGCTAAGGATTTCTTCATGTc      >  1:6635935‑M2/5‑36 (MQ=255)
                                ccaGCAGAAACTTGTGCTAAGGATTTCTTCATGTCa     <  1:7919067‑M2/33‑1 (MQ=255)
                                ccaGCAGAAACTTGTGCTAAGGATTTCTTCATGTCa     <  1:6983618‑M2/33‑1 (MQ=255)
                                ccaGCAGAAACTTGTGCTAAGGATTTCTTCATGTCa     <  1:3595856‑M2/33‑1 (MQ=255)
                                 caGCAGAAACTTGTGCTAAGGATTTCTTCATGTCAt    <  1:2303784‑M2/34‑1 (MQ=255)
                                   gcagAAACTTGTGCTAAGGATTTCTTCATGTCATAc  >  1:2876742/1‑36 (MQ=255)
                                   gcagAAACTTGTGCTAAGGATTTCTTCATGTCATAc  >  1:3028720/1‑36 (MQ=255)
                                   gcagAAACTTGTGCTAAGGATTTCTTCATGTCATAc  >  1:3173358/1‑36 (MQ=255)
                                   gcagAAACTTGTGCTAAGGATTTCTTCATGTCATAc  <  1:1856134/36‑1 (MQ=255)
                                   gcagAAACTTGTGCTAAGGATTTCTTCATGTCATAc  <  1:377857/36‑1 (MQ=255)
                                   |                                   
AACCCACATCTCAACATTATGAACTAATACCGCCAGCAGAAACTTGTGCTAAGGATTTCTTCATGTCATAC  >  REL606/4395016‑4395086

Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 27 ≤  ATCG  < 33 ≤  ATCG  < 37 ≤  ATCG