New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? minE 2577253 =110 (1.860)6 (0.110) 3/94 NT 5.8% coding (846/1650 nt) treF cytoplasmic trehalase
?minE 2577286 = 95 (1.750)coding (813/1650 nt) treF cytoplasmic trehalase

CCTAATCAGGCCTATCGCCATGTTGTGCGGATGCCGGACGGATCGCTGCTCAAC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  minE/2577306‑2577253
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑caacATGGCGATAGGCCTGATTAGGGATTAACGATTCTGCACCGTCCATCCAGAAGGCATATTCCATTTT  >  minE/2577286‑2577351
                                                                                                                        
CCTAATCAGGCCTATCGCCATGTTGTGCGGATGCCGGACGGATCGCTGCTCAA                                                                     >  1:1103682/1‑53
CCTAATCAGGCCTATCGCCATGTTGTGCGGATGCCGGACGGATCGCTGCTCAA                                                                     >  1:453267/1‑53
             ATCGCCATGTTGTGCGGATGCCGGACGGATCGCTGCTCAA                                                                     >  1:2452967/1‑40
                 CCATGTTGTGCGGATGCCGGACGGATCGCTGCTCAACATGGCGATAGGCCTGATTAGGGATTAACGATT                                    <  1:100349/69‑1
                   ATGTTGTGCGGATGCCGGACGGTTCGCTGCTCAACATGGCGATAGGCCTGATTAGGGATTAACGATTCT                                  >  1:1248587/1‑69
                        GTGCGGATGCCGGACGGATCGCTGCTCAACATGGCGATA                                                           <  1:1303232/39‑1
                        GTGCGGATGCCGGACGGATCGCTGCTCAACATGGCGATA                                                           <  1:1778724/39‑1
                        GTGCGGATGCCGGACGGATCGCTGCTCAACATGGCGATA                                                           <  1:3457374/39‑1
                         TGCGGATGCCGGACGGATCGCTGCTCAACATGGCGATA                                                           <  1:1529023/38‑1
                                                   AACATGGCGATAGGCCTGATTAGGGATTAACGATTCTGCACCGTCCATCCAGAAGGCATATTCCATTTT  >  1:754071/1‑69
                                                   AACATGGCGATAGGCCTGATTAGGGATTAACGATTCTGCACCGTCCAT                       >  1:1970223/1‑48
                                                   AACATGGCGATAGGCCTGATTAGGGATTAACGATTCTGCACCGTCCAT                       >  1:3363242/1‑48
                                                   AACATGGCGATAGGCCTGATTAGGGATTAACGATTCTGCACCGTCCAT                       >  1:3584671/1‑48
                                                   AACATGGCGATAGGCCTGATTAGGGATTAACGATTCTGCACCGTCCAT                       >  1:371018/1‑48
                                                                                                                        
CCTAATCAGGCCTATCGCCATGTTGTGCGGATGCCGGACGGATCGCTGCTCAAC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  minE/2577306‑2577253
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑caacATGGCGATAGGCCTGATTAGGGATTAACGATTCTGCACCGTCCATCCAGAAGGCATATTCCATTTT  >  minE/2577286‑2577351

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