New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 3615477 =121 (1.180)11 (0.110) 4/48 NT 10.1% coding (3189/3474 nt) bcsC cellulose synthase subunit
?REL606 = 3947083 74 (0.720)coding (724/1656 nt) aslA acrylsulfatase‑like enzyme

CACGCACCAAAACC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  REL606/3615490‑3615477
‑‑‑‑‑‑‑‑‑‑‑‑‑‑CTGCCGTTCAGCAAAGATGACGTTCATGCCGTGC  <  REL606/3947083‑3947050
                                                
CACCCACCAAAACCCTGCCGTTCAGCAAAGATGAC               <  1:5238842/35‑1
 cacCACCAAAACCCTGCCGTTCAGCAAAGATGACG              <  1:4704627/32‑1
 accCACCAAAACCCTGCCGTTCAGCAAAGATGACG              <  1:5558885/32‑1
   ccgCCAAAACCCTGCCGTTCAGCAAAGATGACGTT            <  1:16076943/32‑1
   cCACCAAACCCCTGCCGTTCAGCAAAGATGACGTT            <  1:15384373/34‑1
   cCACCAAACCCCTGCCGTTCAGCAAAGATGACGTT            <  1:17303595/34‑1
   aaACCAAAACCCTGCCGTTCAGCAAAGATGACGTT            <  1:13041534/33‑1
      cccAAACCCTGCCGTTCAGCAAAGATGACGTTCAT         <  1:9874127/32‑1
      CCAAAACCCTGCCGTTCAGCAAAGATGACGTTCAT         <  1:6753747/35‑1
       CAAAACCCTGCCGTTCAGCAAAGATGACGTTCATG        <  1:17485750/35‑1
        AAACCCCTGCCGTTCAGCAAAGATGACGTTCATGC       <  1:12826600/35‑1
        AAAACCCTGCCGTTCAGCAAAGATGACGTTCATGC       <  1:741744/35‑1
         accCCCTGCCGTTCAGCAAAGATGACGTTCATGCC      <  1:13436441/32‑1
         aacCCCTGCCGTTCAGCAAAGATGACGTTCATGCC      <  1:230793/32‑1
         aacCCCTGCCGTTCAGCAAAGATGACGTTCATGCC      <  1:10025650/32‑1
         AAACCCTGCCGTTCAGCAAAGATGACGTTCATGCC      <  1:4124910/35‑1
          ccCCCTGCCGTTCAGCAAAGATGACGTTCATGCCG     <  1:10053253/33‑1
          acCCCTGCCGTTCAGCAAAGATGACGTTCATGCCG     <  1:16437497/33‑1
          acCCCTGCCGTTCAGCAAAGATGACGTTCATGCCG     <  1:17752810/33‑1
           cCCCTGCCGTTCAGCAAAGATGACGTTCATGCCGT    <  1:3412097/34‑1
           cCCCTGCCGTTCAGCAAAGATGACGTTCATGCCGT    <  1:16518166/34‑1
           ACCCTGCCGTTCAGCAAAGATGACGTTCATGCCGT    <  1:3757545/35‑1
            CCCTGCCGTTCAGCAAAGATGACGTTCATGCCGTG   <  1:4069072/35‑1
            CCCTGCCGTTCAGCAAAGATGACGTTCATGCCGTG   <  1:16168894/35‑1
            CCCTGCCGTTCAGCAAAGATGACGTTCATGCCGTG   <  1:12131706/35‑1
            CCCTGCCGTTCAGCAAAGATGACGTTCATGCCGTG   <  1:6411379/35‑1
            CCCTGCCGTTCAGCAAAGATGACGTTCATGCCGTG   <  1:6600737/35‑1
            CCCTGCCGTTCAGCAAAGATGACGTTCATGCCGTG   <  1:10309945/35‑1
            CCCTGCCGTTCAGCAAAGATGACGTTCATGCCGTG   <  1:757128/35‑1
             CCTGCCGTTCAGCAAAGATGACGTTCATGCCGTGC  <  1:3672730/35‑1
             CCTGCCGTTCAGCAAAGATGACGTTCATGCCGTGC  <  1:4638967/35‑1
             CCTGCCGTTCAGCAAAGATGACGTTCATGCCGTGC  <  1:20449305/35‑1
                                                
CACGCACCAAAACC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  REL606/3615490‑3615477
‑‑‑‑‑‑‑‑‑‑‑‑‑‑CTGCCGTTCAGCAAAGATGACGTTCATGCCGTGC  <  REL606/3947083‑3947050

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