New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 2583271NA (NA)32 (0.530) 3/28 NT 49.9% intergenic (‑441/+11) iscR/yfhQ DNA‑binding transcriptional repressor/predicted methyltransferase
?REL606 = 3686161 55 (0.540)coding (401/741 nt) yiaT hypothetical protein

GCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATT  >  REL606/2583253‑2583305
                  |                                  
gCATTCCGCGTCGCATCAGGCAACGGCTGTTATTc                    <  1:4095101/35‑1 (MQ=255)
gCATTCCGCGTCGCATCAGGCAACGGCTGTTATTc                    <  1:13013891/35‑1 (MQ=255)
gCATTACGCGTCGCATCAGGCAACGGCTGTTATTc                    <  1:18623432/35‑1 (MQ=255)
 cATTACGCGTCGCATCAGGCAACGGCTGTTATTCg                   <  1:18734152/35‑1 (MQ=255)
 cATTACGCGTCGCATCAGGCAACAGCTGTTATTCg                   >  1:434759/1‑35 (MQ=255)
     ccgcgTCGCATCAGGCAACGGCTGTTATTCGGCCt               <  1:16579752/34‑1 (MQ=255)
      cgcgTCGCATCAGGCAACGGCTGTTATTCGGCCtt              <  1:5334701/35‑1 (MQ=255)
      cgcgTCGCAGCGGGCACCGGCTGTTATTCGGCCtt              <  1:3216377/35‑1 (MQ=25)
       gcgTCGCATCAGGCACCGGCTGTTATTCGGCCTTg             <  1:532136/35‑1 (MQ=255)
       gcgTCGCATCAGGCACCGGCTGTTATTCGGCCTTg             <  1:12017351/35‑1 (MQ=255)
       gcgTCGCATCAGGCACCGGCTGTTATTCGGCCTTg             <  1:5911383/35‑1 (MQ=255)
       gcgTCGCATCAGGCAACGGCTGTTATTCGGCCTTg             <  1:10287183/35‑1 (MQ=255)
        cgTCGCATCAGGCAACGGCTGTTATTCGGCCTTGt            <  1:1651965/35‑1 (MQ=255)
        cgTCGCATCAGGCAACGGCTGTTATTCGGCCTTGt            <  1:11207719/35‑1 (MQ=255)
        cgTCGCAACGGGCACCGGCTGTTATTCGGCCTTGt            <  1:20178248/35‑1 (MQ=25)
         gTCGCATCAGGCAACGGCTGTTATTCGGCCTTGtt           <  1:3037436/35‑1 (MQ=255)
          tCGCATCAGGCACCGGCTGTTATTCGGCCTTGTTa          <  1:6529200/35‑1 (MQ=255)
          tCGCATCAGGCACCGGCTGTTATTCGGCCTTGTTa          <  1:16822470/35‑1 (MQ=255)
           cGCATCGGGCAACGGCTGTTATTCGGCCTTGTTAc         <  1:4649211/35‑1 (MQ=255)
           cGCATCGGGCAACGGCTGTTATTCGGCCTTGTTAc         <  1:11606113/35‑1 (MQ=255)
           cGCATCAGGCACCGGCTGTTATTCGGCCTTGTTAc         <  1:17547933/35‑1 (MQ=255)
           cGCATCAGGCACCGGCTGTTATTCGGCCTTGTTAc         <  1:4335041/35‑1 (MQ=255)
           cGCATCAGGCAACGGCTGTTATTCGGCCTTGTTAc         <  1:16105718/35‑1 (MQ=255)
           cGCATCAGGCAACGGCTGTTATTCGGCCTTGTTAc         <  1:15850060/35‑1 (MQ=255)
            gCATCAGGCACCGGCTGTTATTCGGCCTTGTTAcc        <  1:2637202/35‑1 (MQ=255)
            gCATCAGGCACCGGCTGTTATTCGGCCTTGTTAcc        <  1:16035348/35‑1 (MQ=255)
            gCATCAGGCACCGGCTGTTATTCGGCCTTGTTAcc        <  1:4971462/35‑1 (MQ=255)
            gCATCAGGCAACGGCTGTTATTCGGCCTTGTTAcc        <  1:15117135/35‑1 (MQ=255)
            gCATCAGGCAACGGCTGTTATTCGGCCTTGTTAcc        <  1:13671508/35‑1 (MQ=255)
            gCACCAGGCACCGGCTGTTATTCGGCCTTGTTAcc        <  1:4324402/35‑1 (MQ=37)
             cATCAGGCAACGGCTGTTATTCGGCCTTGTTACCt       <  1:18173289/35‑1 (MQ=255)
             cATCAGGCAACGGCTGTTATTCGGCCTTGTTACCt       <  1:6436923/35‑1 (MQ=255)
             cATCAGCCACCGGCTGTTATTCGGCCTTGTTACCt       <  1:9195406/35‑1 (MQ=37)
              aTCAGGCAACGGCTGTTATTCGGCCTTGTTACCtt      <  1:10791380/35‑1 (MQ=255)
              aTCAGCCAACGGCTGTTATTCGGCCTTGTTACCtt      <  1:8408600/35‑1 (MQ=255)
               tcgggCAACGGCTGTTATTCGGCCTTGTTACCttt     <  1:7675702/32‑1 (MQ=255)
               tCAGGCACCGGCTGTTATTCGGCCTTGTTACCttt     <  1:4547686/35‑1 (MQ=255)
               tCAGGCACCGGCTGTTATTCGGCCTTGTTACCttt     <  1:2736447/35‑1 (MQ=255)
               tCAGGCACCGGCTGTTATTCGGCCTTGTTACCttt     <  1:10900424/35‑1 (MQ=255)
               tCAGGCAACGGCTGTTATTCGGCCTTGTTACCttt     <  1:3162351/35‑1 (MQ=255)
               tCAGGCAACGGCTGTTATTCGGCCTTGTTACCttt     <  1:7916231/35‑1 (MQ=255)
               tCAGGCAACGGCTGTTATTCGGCCTTGTTACCttt     <  1:15720314/35‑1 (MQ=255)
                cAGGCACCGGCTGTTATTCGGCCTTGTTACCTTTa    <  1:3107534/35‑1 (MQ=255)
                cAGGCACCGGCTGTTATTCGGCCTTGTTACCTTTa    <  1:6666522/35‑1 (MQ=255)
                cAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTa    <  1:3903120/35‑1 (MQ=255)
                cAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTa    <  1:8823789/35‑1 (MQ=255)
                cAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTa    <  1:9327446/35‑1 (MQ=255)
                 aGGCAACGGCTGTTATTCGGCCTTGTTACCTTTAt   <  1:2397006/35‑1 (MQ=255)
                  ggCAACGGCTGTTATTCGGCCTTGTTACCTTTAtt  <  1:17414465/35‑1 (MQ=255)
                  |                                  
GCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATT  >  REL606/2583253‑2583305

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