New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 180751 =70 (0.680)8 (0.080) 5/46 NT 13.4% coding (554/801 nt) yadT vitamin B12‑transporter protein BtuF
?REL606 3966540 = 36 (0.370)intergenic (‑56/‑97) ECB_03694/ECB_03695 conserved hypothetical protein/Magnesium and cobalt transport protein corA

TTAAGCTCAATATCAATTCGAGAACAAATGCTAGGGTAAAAAAATTTTAAAAAAGAAGTACACTTCA  >  REL606/3966506‑3966572
                                  |                                
ttAAGCTCAATATCAATTCGAGAACAAATGCTAgg                                  >  1:3163084/1‑35 (MQ=255)
ttAAGCTCAATATCAATTCGAGAACAAATGCTAgg                                  >  1:13425670/1‑35 (MQ=255)
ctattcaaaaccaggtactcgaagtttgtggcggg                                  >  1:625100‑M2/35‑35 (MQ=255)
cgattcaaaaccaggtactcgaagtttgtggcggg                                  >  1:1610962‑M2/35‑35 (MQ=255)
cgattcaaaaccaggtactcgaagtttgtggcggg                                  >  1:8127048‑M2/35‑35 (MQ=255)
cgattcaaaaccaggtactcgaagtttgtggcggg                                  >  1:14235025‑M2/35‑35 (MQ=255)
cgattcaaaaccaggtactcgaagtttgtggcggg                                  >  1:12801506‑M2/35‑35 (MQ=255)
cgattcaaaaccaggtactcgaagtttgtggcggg                                  >  1:12748157‑M2/35‑35 (MQ=255)
cgattcaaaaccaggtactcgaagtttgtggcggg                                  >  1:10042545‑M2/35‑35 (MQ=255)
 tAAGCTCAATATCAATTCGAGAACAAATGCTAggg                                 <  1:3350878/35‑1 (MQ=255)
 tAAGCTCAATATCAATTCGAGAACAAATGCTAggg                                 >  1:9213594/1‑35 (MQ=255)
 gattcaaaaccaggtactcgaagtttgtggcgggg                                 >  1:7753169‑M2/34‑35 (MQ=255)
 gattcaaaaccaggtactcgaagtttgtggcgggg                                 >  1:935268‑M2/34‑35 (MQ=255)
 gattcaaaaccaggtactcgaagtttgtggcgggg                                 >  1:15895298‑M2/34‑35 (MQ=255)
  attcaaaaccatgtactcgaagtttgtggcggggg                                >  1:14300648‑M2/33‑34 (MQ=255)
  aaGCTCAATATCAATTCTAGAACAAATGCTAGGGt                                >  1:15809761/1‑35 (MQ=255)
  aaGCTCAATATCAATTCGAGAACAAATGCTAGGGt                                >  1:13731697/1‑35 (MQ=255)
  aaGCTCAATATCAATTCGAGAACAAATGCTAGGGt                                >  1:20366137/1‑35 (MQ=255)
  aaGCTCAATATCAATTCGAGAACAAATGCTAGGGt                                >  1:3226146/1‑35 (MQ=255)
  aaGCTCAATATCAATTCGAGAACAAATGCTAGGGt                                <  1:12479653/35‑1 (MQ=255)
  aaGCTCAATATCAATTCGAGAACAAATGCTAGGGt                                <  1:12106194/35‑1 (MQ=255)
   ttcaaaaccaggtactcgaagtttgtggcggggga                               >  1:4015929‑M2/32‑33 (MQ=255)
   aGCTCAATATCAATTCGAGAACAAATGCTAGGGTa                               >  1:19266818/1‑35 (MQ=255)
   aGCTCAATATCAATTCGAGAACAAATGCTAGGGTa                               <  1:12827075/35‑1 (MQ=255)
   aGCTCAATATCAATTCGAGAACAAATGCTAGGGTa                               >  1:16978888/1‑35 (MQ=255)
    tcaaaaccaggtactcgaagtttgtggcggGGAaa                              >  1:15257790‑M2/31‑35 (MQ=255)
    tcaaaaccaggtactcgaagtttgtggcggGGAaa                              >  1:5556191‑M2/31‑35 (MQ=255)
    tcaaaaccaggtactcgaagtttgtggcggGGAaa                              >  1:12288653‑M2/31‑35 (MQ=255)
    gCTCAATATCAATTCGAGAACAAATGCTAGGGTaa                              >  1:5910054/1‑35 (MQ=255)
    gCTCAATATCAATTCGAGAACAAATGCTAGGGTaa                              <  1:4809826/35‑1 (MQ=255)
    gCTCAATATCAATTCGAGAACAAATGCTAGGGTaa                              >  1:7152372/1‑35 (MQ=255)
    gCTCAATATCAATTCGAGAACAAATGCTAGGGTaa                              <  1:7615231/35‑1 (MQ=255)
    gCTCAATATCAATTCGAGAACAAATGCTAGGGTaa                              >  1:18196090/1‑35 (MQ=255)
     caaaaccaggtactcgaagtttgtggcggGGGaaa                             >  1:15545006‑M2/30‑35 (MQ=255)
     caaaaccaggtactcgaagtttgtggcggGGGaaa                             >  1:13617395‑M2/30‑35 (MQ=255)
      aaaaccaggtactcgaagtttgtggcggGGGaaaa                            >  1:13554901‑M2/29‑35 (MQ=255)
       cAATATCAATTCGAGAACAAATGCTAGGGTaaaaa                           >  1:3095775/1‑35 (MQ=255)
       cAATATCAATTCGAGAACAAATGCTAGGGTaaaaa                           <  1:1010554/35‑1 (MQ=255)
       cAATATCAATTCGAGAACAAATGCTAGGGTaaaaa                           >  1:7185305/1‑35 (MQ=255)
       cAATATCAATTCGAGAACAAATGCTAGGGTaaaaa                           >  1:3468897/1‑35 (MQ=255)
       cAATATCAATTCGAGAACAAATGCTAGGGTaaaaa                           <  1:14705699/35‑1 (MQ=255)
        aaccaggtactcgaagtttgtggcggGGGaaaaaa                          >  1:3230551‑M2/27‑35 (MQ=255)
        aaccaggtactcgaagtttgtggcggGGGaaaaaa                          >  1:17666541‑M2/27‑35 (MQ=255)
         accaggtactcgaagtttgtggcggGGgaaaaaac                         >  1:13777568‑M2/26‑34 (MQ=255)
          tatCAATTCGAGAACAAATGCTAGGGTAAAAAAAt                        >  1:17924276/1‑35 (MQ=255)
           atCAATTCGAGAACAAATGCTAGGGTAAAAAAAtt                       >  1:9072049/1‑35 (MQ=255)
            tCAATTCGAGAACAAATGCTAGGGTAAAAAAAttt                      <  1:5413452/35‑1 (MQ=255)
              aaTTCGAGAACAAATGCTAGGGTAAAAAAATTTTa                    <  1:14440523/35‑1 (MQ=255)
              aaTTCGAGAACAAATGCTAGGGTAAAAAAATTTTa                    >  1:5549352/1‑35 (MQ=255)
               tactcgaagtttgtggcggGGGAAAAAAATTTTaa                   >  1:8470961‑M2/20‑35 (MQ=255)
                ttCGAGAACAAATGCTAGGGTAAAAAAATTTTaaa                  <  1:4753988/35‑1 (MQ=255)
                actcgaagtttgtggcggGGGAAAAAAATTTTaaa                  >  1:9701890‑M2/19‑35 (MQ=255)
                 tCGAGAACAAATGCTAGGGTAAAAAAATTTTaaaa                 >  1:20614284/1‑35 (MQ=255)
                  cGAGAACAAATGCTAGGGTAAAAAAATTTTaaaaa                >  1:20188145/1‑35 (MQ=255)
                  cGAGAACAAATGCTAGGGTAAAAAAATTTTaaaaa                <  1:8448636/35‑1 (MQ=255)
                   gagaACAAATGCTAGGGTAAAAAAATTTTaaaaaa               >  1:8083545/1‑35 (MQ=255)
                   gagaACAAATGCTAGGGTAAAAAAATTTTaaaaaa               <  1:12928513/35‑1 (MQ=255)
                   gagaACAAATGCTAGGGTAAAAAAATTTTaaaaaa               >  1:4220690/1‑35 (MQ=255)
                   gagaACAAATGCTAGGGTAAAAAAATTTTaaaaaa               >  1:4533397/1‑35 (MQ=255)
                    agaACAAATGCTAGGGTAAAAAAATTTTAAAAaag              <  1:17624756/35‑1 (MQ=255)
                     gaACAAATGCTAGGGTAAAAAAATTTTAAAAaaga             >  1:12575806/1‑35 (MQ=255)
                      aaCAAATGCTAGGGTAAAAAAATTTTAAAAaagaa            >  1:7599844/1‑35 (MQ=255)
                      aaCAAATGCTAGGGTAAAAAAATTTTAAAAaagaa            >  1:19615919/1‑35 (MQ=255)
                      aaCAAATGCTAGGGTAAAAAAATTTTAAAAaagaa            >  1:1535954/1‑35 (MQ=255)
                       aCAAATGCTAGGGTAAAAAAATTTTAAAAaagaag           >  1:20081891/1‑35 (MQ=255)
                       aCAAATGCTAGGGTAAAAAAATTTTAAAAaagaag           >  1:14402625/1‑35 (MQ=255)
                        cAAATGCTAGTGTAAAAAAATTTTAAAAAAGAAGt          >  1:9018111/1‑35 (MQ=255)
                        cAAATGCTAGGGTAAAAAAATTTTAAAAAAGAAGt          >  1:4959230/1‑35 (MQ=255)
                        cAAATGCTAGGGTAAAAAAATTTTAAAAAAGAAGt          <  1:18842457/35‑1 (MQ=255)
                         aaaTGCTAGGGTAAAAAAATTTTAAAAAAGAAGTa         >  1:9810040/1‑35 (MQ=255)
                         aaaTGCTAGGGTAAAAAAATTTTAAAAAAGAAGTa         >  1:9179414/1‑35 (MQ=255)
                         aaaTGCTAGGGTAAAAAAATTTTAAAAAAGAAGTa         <  1:16432858/35‑1 (MQ=255)
                         aaaTGCTAGGGTAAAAAAATTTTAAAAAAGAAGTa         <  1:19902394/35‑1 (MQ=255)
                          aaTGCTAGGGTAAAAAAATTTTAAAAAAGAAGTac        >  1:19211580/1‑35 (MQ=255)
                           aTGCTAGGGTAAAAAAATTTTAAAAAAGAAGTaca       <  1:6047046/35‑1 (MQ=255)
                           aTGCTAGGGTAAAAAAATTTTAAAAAAGAAGTaca       >  1:10536293/1‑35 (MQ=255)
                           aTGCTAGGGTAAAAAAATTTTAAAAAAGAAGTaca       <  1:3374373/35‑1 (MQ=255)
                            tGCTAGGGTAAAAAAATTTTAAAAAAGAAGTacac      >  1:18170414/1‑35 (MQ=255)
                            tGCTAGGGTAAAAAAATTTTAAAAAAGAAGTacac      <  1:13299423/35‑1 (MQ=255)
                              cTAGGGTAAAAAAATTTTAAAAAAGAAGTACACtt    <  1:18274462/35‑1 (MQ=255)
                              cTAGGGTAAAAAAATTTTAAAAAAGAAGTACACtt    <  1:14337570/35‑1 (MQ=255)
                              cTAGGGTAAAAAAATTTTAAAAAAGAAGTACACtt    <  1:13659349/35‑1 (MQ=255)
                               tAGGGTAAAAAAATTTTAAAAAAGAAGTACACTTc   >  1:3942389/1‑35 (MQ=255)
                               tAGGGTAAAAAAATTTTAAAAAAGAAGTACACTTc   <  1:9243879/35‑1 (MQ=255)
                               tAGGGTAAAAAAATTTTAAAAAAGAAGTACACTTc   <  1:5532189/35‑1 (MQ=255)
                               tAGGGTAAAAAAATTTTAAAAAAGAAGTACACTTc   <  1:8908985/35‑1 (MQ=255)
                               tAGGGTAAAAAAATTTTAAAAAAGAAGTACACTTc   <  1:18454403/35‑1 (MQ=255)
                                aGGGTAAAAAAATTTTAAAAAAGAAGTACACTTca  <  1:13748002/35‑1 (MQ=255)
                                aGGGTAAAAAAATTTTAAAAAAGAAGTACACTTca  <  1:14500996/35‑1 (MQ=255)
                                  |                                
TTAAGCTCAATATCAATTCGAGAACAAATGCTAGGGTAAAAAAATTTTAAAAAAGAAGTACACTTCA  >  REL606/3966506‑3966572

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ 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.