New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 1539621 =78 (0.760)8 (0.100) 4/38 NT 13.6% coding (51/2400 nt) dos cAMP phosphodiesterase, heme‑regulated
?REL606 = 3363882 40 (0.490)coding (771/1290 nt) rrmB 16S rRNA m5C967 methyltransferase, S‑adenosyl‑L‑methionine‑dependent

CATTTTTTTCCCCGCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  REL606/1539636‑1539621
‑‑‑‑‑‑‑‑‑‑‑ccgccGGGGGCGGCACAAAGATCCAAAATGTGTTCACCG  <  REL606/3363882‑3363849
                                                  
ttTTTTTTTCCCCGCCGGGGGCGGCACAAAGATCC                 <  1:16714079/33‑1
gcTTTTTTTCCCCGCCGGGGGCGGCACAAAGATCC                 <  1:13586159/33‑1
  TTTTTTTCCCCGCCGGGGGCGGCACAAAGATCCAA               <  1:8707507/35‑1
   ttgTTTCCCCGCCGGGGGCGGCACAAAGATCCAAA              <  1:15364491/32‑1
   TTTTTTCCCCGCCGGGGGCGGCACAAAGATCCAAA              <  1:20202255/35‑1
   TTTTTTCCCCGCCGGGGGCGGCACAAAGATCCAAA              <  1:527614/35‑1
    TTTTTCCCCGCCGGGGGCGGCACAAAGATCCAAAA             <  1:18958961/35‑1
    TTTTTCCCCGCCGGGGGCGGCACAAAGATCCAAAA             <  1:16223364/35‑1
     TTTTCCCCGCCGGGGGCGGCACAAAGATCCAAAAT            <  1:2635236/35‑1
          CCCGCCGGGGGCGGCACAAAGATCCAAAATGTGTT       <  1:8532131/35‑1
          CCCGCCGGGGGCGGCACAAAGATCCAAAATGTGTT       <  1:11898487/35‑1
           CCGCCGGGGGCGGCACAAAGATCCAAAATGTGTTC      <  1:18136557/35‑1
           CCGCCGGGGGCGGCACAAAGATCCAAAATGTGTTC      <  1:19729640/35‑1
            CGCCGGGGGCGGCACAAAGATCCAAAATGTGTTCA     <  1:4148125/35‑1
            CGCCGGGGGCGGCACAAAGATCCAAAATGTGTTCA     <  1:11383578/35‑1
             GCCGGGGGCTGCACAAAGATCCAAAATGTGTTCAC    >  1:5518761/1‑35
             GCCGGGGGCGGCACAAAGATCCAAAATGTGTTCAC    >  1:9953363/1‑35
             GCCGGGGGCGGCACAAAGATCCAAAATGTGTTCAC    <  1:10184431/35‑1
              CCGGGGGCGGCACAAAGATCCAAAATGTGTTCACC   <  1:19526442/35‑1
              CCGGGGGCGGCACAAAGATCCAAAATGTGTTCACC   >  1:4610040/1‑35
              CCGGGGGCGGCACAAAGATCCAAAATGTGTTCACC   <  1:13307573/35‑1
               CGGGGGCGGCACAAAGATCCAAAATGTGTTCACCG  <  1:18748705/35‑1
               CGGGGGCGGCACAAAGATCCAAAATGTGTTCACCG  <  1:17123629/35‑1
               CGGGGGCGGCACAAAGATCCAAAATGTGTTCACCG  >  1:8033970/1‑35
               CGGGGGCGGCACAAAGATCCAAAATGTGTTCACCG  >  1:10368893/1‑35
                                                  
CATTTTTTTCCCCGCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  REL606/1539636‑1539621
‑‑‑‑‑‑‑‑‑‑‑ccgccGGGGGCGGCACAAAGATCCAAAATGTGTTCACCG  <  REL606/3363882‑3363849

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