New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 294204 =76 (0.740)6 (0.070) 3/38 NT 10.8% intergenic (+474/‑53) ykgD/ykgE predicted DNA‑binding transcriptional regulator/predicted oxidoreductase
?REL606 = 2574914 39 (0.480)intergenic (+756/+62) sseA/sseB 3‑mercaptopyruvate sulfurtransferase/rhodanase‑like enzyme, sulfur transfer from thiosulfate

ATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCC  >  REL606/2574881‑2574948
                                 |                                  
ataataACCGGAGTCTCTCCGGCGACAATTTACTg                                   <  1:16887462‑M2/35‑2 (MQ=255)
ataataACCGGAGTCTCTCCGGCGACAATTTACTg                                   >  1:8075755‑M2/1‑34 (MQ=255)
ataataACCGGAGTCTCTCCGGCGACAATTTACTg                                   >  1:9205534‑M2/1‑34 (MQ=255)
ataataACCGGAGTCTCTCCGGCGACAATTTACTg                                   >  1:20457041‑M2/1‑34 (MQ=255)
 taataaCCGGAGTCTCTCCGGCGACAATTTACTgg                                  <  1:5248395‑M2/35‑3 (MQ=255)
 taataaCCGGAGTCTCTCCGGCGACAATTTACTgg                                  <  1:14193213‑M2/35‑3 (MQ=255)
 atctaaccgGAGTCTCTCCGGCGACAATTTACTgg                                  >  1:10838360‑M2/4‑33 (MQ=255)
  aataaCCGGAGTCTCTCCGGCGACAATTTACttgt                                 >  1:18892306/1‑35 (MQ=255)
  aataaCCGGAGTCTCTCCGGCGACAATTTACTggt                                 <  1:4752067‑M2/35‑4 (MQ=255)
  aataaCCGGAGTCTCTCCGGCGACAATTTACTggt                                 <  1:14252286‑M2/35‑4 (MQ=255)
   ataaCCGGAGTCTCTCCGGCGACAATTTACTggtg                                <  1:19160850‑M2/35‑5 (MQ=255)
   ataaCCGGAGTCTCTCCGGCGACAATTTACTggtg                                >  1:3449838‑M2/1‑31 (MQ=255)
   ataaCCGGAGTCTCTCCGGCGACAATTTACTggtg                                <  1:8847384‑M2/35‑5 (MQ=255)
   ataaCCGGAGTCTCTCCGGCGACAATTTACTggtg                                >  1:9203392‑M2/1‑31 (MQ=255)
    taaCCGGAGTCTCTCCGGCGACAATTTACTggtgg                               >  1:10772502‑M2/1‑30 (MQ=255)
     aaCCGGAGTCTCTCCGGCGACAATTTACTggtggt                              <  1:5962818/35‑1 (MQ=255)
     aaCCGGAGTCTCTCCGGCGACAATTTACTggtggt                              <  1:5888106/35‑1 (MQ=255)
     aaCCGGAGTCTCTCCGGCGACAATTTACTggtggg                              <  1:15218677‑M2/35‑7 (MQ=255)
      aCCGGGGTCTCTCCGGCGACAATTTACTGGTGGtt                             <  1:7006298/35‑1 (MQ=255)
      aCCGGAGTCTCTCCGGCGACAATTTACTggtgggg                             <  1:12226152‑M2/35‑8 (MQ=255)
      aCCGGAGTCTCTCCGGCGACAATTTACTGGTGGtt                             <  1:1085865/35‑1 (MQ=255)
      aCCGGAGTATCTCCGGCGACAATTTACTGGTGGtt                             <  1:11772424/35‑1 (MQ=255)
       ccGGAGTCTCTCCGGCGACAATTTACTggtggggg                            <  1:652926‑M2/35‑9 (MQ=255)
       ccGGAGTCTCTCCGGCGACAATTTACTGGTGGTTa                            <  1:19878649/35‑1 (MQ=255)
        cGGAGTCTCTCCGGCGACAATTTACTggtgggggg                           <  1:9264673‑M2/35‑10 (MQ=255)
         ggAGTCTCTCCGGCGACAATTTACTggtggggggg                          <  1:11719300‑M2/35‑11 (MQ=255)
         ggAGTCTCTCCGGCGACAATTTACTggtgggggga                          <  1:16044502‑M2/35‑11 (MQ=255)
         ggAGTCTCTCCGGAGACAATTTACTGGTGGTTaac                          <  1:20339306/35‑1 (MQ=255)
          gAGTCTCTCCGGCGACAATTTACTggtgggggggg                         <  1:13285518‑M2/35‑12 (MQ=255)
          gAGTCTCTCCGGCGACAATTTACTggtgggggaga                         <  1:9563628‑M2/35‑12 (MQ=255)
          gAGTCTCTCCGGCGACAATTTACTGGTGGTTaaca                         <  1:7706109/35‑1 (MQ=255)
          gAGTCTCTCCGGCGACAATTTACTGGTGGTTaaca                         <  1:6938457/35‑1 (MQ=255)
          gAGTCTATCCGGCGACAATTTACTGGTGGTTaaca                         <  1:589064/35‑1 (MQ=255)
           nGTCTCTCCGGCGACAATTTACTggtggggggaga                        <  1:10611803‑M2/35‑13 (MQ=255)
           aGTCTCTCCGGCGACAATTTACTggtgggggggaa                        <  1:11868187‑M2/35‑13 (MQ=255)
           aGTCTCTCCGGCGACAATTTACTGGTGGTTaacaa                        <  1:8395279/35‑1 (MQ=255)
            gTCTCTCCGGCGACAATTTACTGGTGGTTaacaac                       <  1:15832047/35‑1 (MQ=255)
            gTCTCTCCGGCGACAATTTACTGGTGGTTaacaac                       <  1:11660846/35‑1 (MQ=255)
            gTCTCTCCGGCGACAATTTACTGGTGGTTaacaac                       <  1:13213176/35‑1 (MQ=255)
             tctctcCGGCTACAATTTACTggtgggggggaaac                      <  1:2817470‑M2/35‑15 (MQ=255)
             tctctcCGGCGACAATTTACTggtgggggggaaac                      <  1:4171774‑M2/35‑15 (MQ=255)
             tctctcCGGCGACAATTTACTggtggggggcaaaa                      <  1:9226844‑M2/35‑15 (MQ=255)
             tctctcCGGCGACAATTTACTGGTGGTTAACAAcc                      <  1:6918343/35‑1 (MQ=255)
             tctctcCGGCGACAATTTACTGGTGGTTAACAAcc                      <  1:18841804/35‑1 (MQ=255)
              ctctcCGGCGACAATTTACTggtgggggggaaaca                     <  1:2815950‑M2/35‑16 (MQ=255)
              ctctcCGGCGACAATTTACTggtgggggagaaacc                     <  1:16335979‑M2/35‑16 (MQ=255)
              ctctcCGGCGACAATTTACTGGTGGTTAACAACCt                     <  1:7966367/35‑1 (MQ=255)
              ctctcCGGCGACAATTTACTGGTGGTTAACAACCt                     <  1:18523205/35‑1 (MQ=255)
              ctctcCGGCGACAATTTACTGGTGGTTAACAACCt                     <  1:7698008/35‑1 (MQ=255)
               tctcCGGCGACAATTTACTGGTGGTTAACAACCtt                    <  1:3636150/35‑1 (MQ=255)
               tctcCGGCGACAATTTACTGGTGGTTAACAACCtt                    <  1:4930015/35‑1 (MQ=255)
                 tcCGGCGACAATTTACTGGTGGTTAACAACCTTCa                  <  1:20215234/35‑1 (MQ=255)
                 tcCGGCGACAATTTACTGGTGGTTAACAACCTTCa                  <  1:19179586/35‑1 (MQ=255)
                 tcCGGCGACAATTTACTGGTGGTTAACAACCTTCa                  <  1:10644450/35‑1 (MQ=255)
                  ccGGCGACAATTTACTGGTGGTTAACAACCTTCag                 <  1:5621908/35‑1 (MQ=255)
                   cGGCGACAATTTACTGGTGGTTAACAACCTTCaga                <  1:1110159/35‑1 (MQ=255)
                   cGGCGACAATTTACTGGTGGTTAACAACCTTCaga                <  1:4699435/35‑1 (MQ=255)
                   cGGCGACAATTTACTGGTGGTTAACAACCTTCaga                <  1:16935371/35‑1 (MQ=255)
                    ggCGACAATTTACTGGTGGTTAACAACCTTCagag               <  1:1030323/35‑1 (MQ=255)
                     gCGACAATTTACTGGTGGTTAACAACCTTCAGagc              <  1:6645710/35‑1 (MQ=255)
                      cGACAATTTACTGGTGGTTAACAACCTTCAGagca             <  1:2596413/35‑1 (MQ=255)
                      cGACAATTTACTGGTGGTTAACAACCTTCAGagca             <  1:14887066/35‑1 (MQ=255)
                       gNCAATTTACTGGTGGTTAACAACCTTCAGagcan            <  1:15622328/35‑1 (MQ=255)
                       gACAATTTACTGGTGGTTAACAACCTTCAGagcag            <  1:14498479/35‑1 (MQ=255)
                       gACAATTTACTGGTGGTTAACAACCTTCAGagcag            <  1:15487104/35‑1 (MQ=255)
                        aCAATTTACTGGTGGTTAACAACCTTCAGagcagc           <  1:5847838/35‑1 (MQ=255)
                         cAATTTACTGGTGGTTAACAACCTTCAGAgcagca          <  1:6264347/35‑1 (MQ=255)
                         cAATTTACTGGTGGTTAACAACCTTCAGAgcagca          <  1:9156647/35‑1 (MQ=255)
                         cAATTTACTGGTGGTTAACAACCTTAAGAgcagca          <  1:15696950/35‑1 (MQ=255)
                          aaTTTACTGGTGGTTAACAACCTTCAGAGCAGCaa         <  1:17842839/35‑1 (MQ=255)
                          aaTTTACTGGTGGTTAACAACCTTCAGAGCAGCaa         <  1:11023332/35‑1 (MQ=255)
                          aaTTTACTGGTGGTTAACAACCTTCAGAGCAGCaa         <  1:10570888/35‑1 (MQ=255)
                           aTTTACTGGTGGTTAACAACCTTCAGAGCAGCAAg        <  1:7916359/35‑1 (MQ=255)
                           aTTTACTGGTGGTTAACAACCTTCAGAGCAGCAAg        <  1:1976177/35‑1 (MQ=255)
                           aTTTACTGGTGGTTAACAACCTTCAGAGCAGCAAg        <  1:9811072/35‑1 (MQ=255)
                            tttACTGGTGGTTAACAACCTTCAGAGCAGCAAGt       <  1:6422197/35‑1 (MQ=255)
                             ttACTGGTGGTTAACAACCTTCAGAGCAGCAAGTa      <  1:8898270/35‑1 (MQ=255)
                                cTGGTGTTTAACAACCTTCAGAGCAGCAAGTAAGc   <  1:18137186/35‑1 (MQ=255)
                                cTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGc   >  1:838772/1‑35 (MQ=255)
                                cTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGc   >  1:20560270/1‑35 (MQ=255)
                                cTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGc   >  1:4516512/1‑35 (MQ=255)
                                 tggtggTTAACAACCTTCAGAGCAGCAAGTAAGcc  >  1:18293664/1‑35 (MQ=255)
                                 tggtggTTAACAACCTTCAGAGCAGCAAGTAAGcc  <  1:1091053/35‑1 (MQ=255)
                                 |                                  
ATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCC  >  REL606/2574881‑2574948

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.