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

AAATTCACTTTCCTTTAAAAAACAAAAAACCGCCA.AAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGC  >  REL606/2574477‑2574546
                                  |                                    
aaaTTCACTTTCATTAAAAAAAAAAAAAACCGCCa.a                                    <  1:119869/36‑1 (MQ=25)
  aTTCACTTTCCTTTAAAAAACAAAAAACAGCCa.aaa                                  >  1:6212548/1‑36 (MQ=255)
        tttCCTTTAAAAAACAAAAAACCGCCA.AAATCAGGc                            >  1:3532640/1‑36 (MQ=255)
            cTTTAAAAAACAAAAAACCGCCA.AAATCAGGCGGtt                        >  1:6674447/1‑36 (MQ=255)
             tttAAAAAACAAAAAACCGCCA.AAATCAGGCGGttt                       >  1:403742/1‑36 (MQ=255)
             tttAAAAAACAAAAAACCGCCA.AAATCAGGCGGttt                       >  1:3603499/1‑36 (MQ=255)
              ttAAAAAACAAAAAACCGCCA.AAATCAGGCGGtttt                      >  1:4244434/1‑36 (MQ=255)
              ttAAAAAACAAAAAACCGCCA.AAATCAGGCGGtttt                      >  1:3434321/1‑36 (MQ=255)
              attaaaaaaaaaaaacCGCCA.AATTC.GGCGGttttt                     <  1:7642955/33‑1 (MQ=17)
                aaaaaaaaaaaaacCGCCAgcagaaacttgtgctaa                     <  1:5743810‑M1/36‑18 (MQ=255)
                aaaaaaaaaaaaacCGCCA.AATTC.GGCGGtttttta                   <  1:1705625/36‑2 (MQ=17)
                aaaaaaaaaaaaacCGCCA.AAATCAGGCGGtttttt                    <  1:1193192/36‑1 (MQ=255)
                 aaaaaaaaaaaacCTCCA.AAATCAGGCGGTTTTttg                   <  1:6652018/36‑1 (MQ=25)
                 aaaaaaaaaaaacCGCCAgcagaaacttgtgctaag                    <  1:4104177‑M1/36‑19 (MQ=255)
                  aaaacaaaaaacCGACA.AAATCAGGCGGtttttttt                  >  1:2141702/1‑34 (MQ=37)
                  aaaaaaaaaaacCGCCAgcagaaacttgtgctaagg                   <  1:5339383‑M1/36‑20 (MQ=255)
                  aaaaaaaaaaacCGCCA.AAATCAGGCGGTTTTttgt                  <  1:2916755/36‑1 (MQ=255)
                   aaaaaaaaaacCGCCAtagaaaggaccaaagaatga                  <  1:6178424‑M1/36‑21 (MQ=255)
                    acaaaaaaaccgcCAtagaaaggaccaaagaatgac                 <  1:4472496‑M1/33‑22 (MQ=255)
                    aacaaaaaacCGCCA.AAATCAGGCGGTTTTttgttg                >  1:5296552/1‑36 (MQ=255)
                     ataaaaaaccgCCAtagaaaggaccaaagaatgaca                <  1:6503580‑M1/34‑23 (MQ=255)
                     acaaaaaacCGCCA.AAATCAGGCGGTGTTttttttc               >  1:899079/1‑34 (MQ=21)
                     aaaaaaaaccgCCAtagaaaggaccaaagaatgaca                <  1:5330331‑M1/34‑23 (MQ=255)
                     aaaaaaaaccgCCAgcagaaacttgtgctaaggatt                <  1:4115681‑M1/34‑23 (MQ=255)
                     aaaaaaaacCGCCAGTAAACCGGCGGTGTAttgttg                >  1:7640553/3‑36 (MQ=25)
                      aaaacaaccGCCAtagaaaggaccaaagaatgacac               <  1:4552154‑M1/35‑24 (MQ=255)
                       aaaaaccCGCCA.AAATCAGGCGGTTTTTTGTTGCTg             <  1:3215682/36‑1 (MQ=255)
                       aaaaaacCGCCAtagaaaggaccaaagaatgacact              <  1:6044528‑M1/36‑25 (MQ=255)
                       aaaaaacCGCCAgcagaaacttgtgctaaggatttc              <  1:6882549‑M1/36‑25 (MQ=255)
                        aaaaacCGCCAtagaaaggaccaaagaatgacactt             <  1:820405‑M1/36‑26 (MQ=255)
                        aaaaacCGCCA.AAATCAGGCGGTTTTTTGTTTCttg            >  1:5133876/1‑34 (MQ=25)
                         aaaccCGCCAtagaaaggaccaaagaatgacacttt            <  1:4291462‑M1/36‑27 (MQ=255)
                         aaaacCGCCAgcagaaacttgtgctaaggatttctt            <  1:1599010‑M1/36‑27 (MQ=255)
                         aaaacCGACAtagaaaggaccaaagaatgacacttt            <  1:5449419‑M1/36‑27 (MQ=255)
                          aacccgccAtagaaaggaccaaagaatgacactttc           <  1:1222132‑M1/33‑28 (MQ=255)
                          aacccGCCA.AAATCAGGCGGTTTTTTGTTGCTGGTc          <  1:2081787/33‑1 (MQ=255)
                          aaacCGCCAtagaaaggaccaaagaatgacactttc           <  1:4397841‑M1/36‑28 (MQ=255)
                          aaacCGCCAtagaaaggaccaaagaatgacactttc           <  1:3650358‑M1/36‑28 (MQ=255)
                          aaacCGCCA.AAATCAGGCGGTTTTTTGTTGCTGGTc          >  1:2464715/1‑36 (MQ=255)
                           aacCGCCAtagaaaggaccaaagaatgacactttca          >  1:217015‑M1/1‑8 (MQ=255)
                           aacCGCCAtagaaaggaccaaagaatgacactttca          >  1:5404714‑M1/1‑8 (MQ=255)
                           aacCGCCAtagaaaggaccaaagaatgacactttca          >  1:7642074‑M1/1‑8 (MQ=255)
                           aacCGCCAtagaaaggaacaaagaatgacactttca          >  1:580295‑M1/1‑8 (MQ=255)
                           aacCGCCA.AAATCAGGCGGTTTTTTTTTTCTTGTcc         >  1:4951115/1‑36 (MQ=21)
                           aacCGCCA.AAATCAGGCGGTTTTTTTTTGCTGTTca         >  1:758734/1‑35 (MQ=21)
                           aacCGCCA.AAATCAGGCGGTTTTTTGTTTCTGGTcc         >  1:2919651/1‑36 (MQ=255)
                            acCGCCAtagaaaggaccaaagaatgacactttcag         >  1:4624905‑M1/1‑7 (MQ=255)
                            acCGCCAgcagaaacttgtgctaaggatttcttcat         >  1:4375433‑M1/1‑7 (MQ=255)
                            acCGCCAgcagaaacttgtgctaaggatttcttcat         >  1:3710895‑M1/1‑7 (MQ=255)
                            acCGCCAgcagaaacttgtgctaaggatttcttcat         >  1:5610708‑M1/1‑7 (MQ=255)
                            acCGCCAgcagaaacttgcgctaaggatttcttcat         >  1:4794151‑M1/1‑7 (MQ=255)
                              cGCCAtagaaaggaccaaagaatgacactttcagcg       >  1:956339‑M1/1‑5 (MQ=255)
                              cGCCAtagaaaggaccaaagaatgacactttcagcg       >  1:6011815‑M1/1‑5 (MQ=255)
                              cGCCAtagaaaggaccaaagaatgacactttcagcg       >  1:6091859‑M1/1‑5 (MQ=255)
                              cGCCAtagaaaggaccaaagaatgacactttcagcg       >  1:6980178‑M1/1‑5 (MQ=255)
                              cGCCAtagaaaggaccaaagaatgacactttcagcg       <  1:6759461‑M1/36‑32 (MQ=255)
                              cGCCAgcagaaacttgtgctaaggatttcttcatgt       <  1:6190677‑M1/36‑32 (MQ=255)
                               gcaagcagaaacttgtgctaaggatttcttcatgtc      >  1:6635935‑M1/4‑4 (MQ=255)
                               gCCAtagaaaggaccaaagaatgacactttcagcgc      >  1:991862‑M1/1‑4 (MQ=255)
                               gCCAtagaaaggaccaaagaatgacactttcagcgc      >  1:1232995‑M1/1‑4 (MQ=255)
                               gCCAtagaaaggaccaaagaatgacactttcagcgc      >  1:2014258‑M1/1‑4 (MQ=255)
                               gCCAtagaaaggaccaaagaatgacactttcagcgc      >  1:2705456‑M1/1‑4 (MQ=255)
                               gCCAtagaaaggaccaaagaatgacactttcagcgc      <  1:5890457‑M1/36‑33 (MQ=255)
                               gCCAtagaaaggaccaaagaatgacactttcagcgc      >  1:6219532‑M1/1‑4 (MQ=255)
                               gCCAgcagaaacttgtgctaaggatttcttcatgtc      >  1:454975‑M1/1‑4 (MQ=255)
                               gCCAgcagaaacttgtgctaaggatttcttcatgtc      <  1:7668201‑M1/36‑33 (MQ=255)
                               gCCAgcagaaacttgtgctaaggatttcttcatgtc      <  1:2803911‑M1/36‑33 (MQ=255)
                                ccAtagaaaggaccaaagaatgacactttcagcgca     >  1:2278475‑M1/1‑3 (MQ=255)
                                ccAtagaaaggaccaaagaatgacactttcagcgca     <  1:3670860‑M1/36‑34 (MQ=255)
                                ccAtagaaaggaacaaagaatgacactttcagcgca     <  1:3594568‑M1/36‑34 (MQ=255)
                                ccAgcagaaacttgtgctaaggatttcttcatgtca     <  1:6983618‑M1/36‑34 (MQ=255)
                                ccAgcagaaacttgtgctaaggatttcttcatgtca     <  1:3595856‑M1/36‑34 (MQ=255)
                                ccAgcagaaacttgtgctaaggatttcttcatgtca     <  1:7919067‑M1/36‑34 (MQ=255)
                                ccA.AAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTc    <  1:4975372/36‑1 (MQ=255)
                                 cagcagaaacttgtgctaaggatttcttcatgtcat    <  1:2303784‑M1/36‑35 (MQ=255)
                                 cAtagaaaggaccaaagaatgacactttcagcgcag    >  1:3374399‑M1/1‑2 (MQ=255)
                                 cAtagaaaggaccaaagaatgacactttcagcgcag    <  1:350899‑M1/36‑35 (MQ=255)
                                 cAtagaaaggaccaaagaatgacactttcagcgcag    <  1:3522576‑M1/36‑35 (MQ=255)
                                 cAtagaaaggaccaaagaatgacactttcagcgcag    <  1:4838551‑M1/36‑35 (MQ=255)
                                  a.aaaTCAGGCGGTTTTTTGTTGCTGGTCCGGTTcgc  <  1:7127559/36‑1 (MQ=255)
                                  |                                    
AAATTCACTTTCCTTTAAAAAACAAAAAACCGCCA.AAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGC  >  REL606/2574477‑2574546

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