New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? AP009048 1933441 =53 (0.940)3 (0.060) 3/94 NT 5.1% coding (847/1179 nt) purT phosphoribosylglycinamide formyltransferase 2
?AP009048 1933468 = 63 (1.210)coding (874/1179 nt) purT phosphoribosylglycinamide formyltransferase 2

AAGGCACGTACATGCAGGGCAAACTCTGAGAGATCTTGAGAAATTAACGTCACCATCCCGGTATCATGTGGACG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AP009048/1933514‑1933441
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gacgTTAATTTCTCAAGATCTCTCAGAGTTTGCCCTGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCGGGA  >  AP009048/1933468‑1933537
                                                                                                                                                
AAGGCACGTACATGCAGGGCAAACTCTGAGAGATCTTGAGAAATTAACGTCACCATCCCGGTATCATGTGG                                                                           <  1:435249/71‑1
  GGCACGTACATGCAGGGCAAACTCTGAGAGATCTTGAGAAATTAACGTCACCATCCCGGTATCATGTGGAC                                                                         <  1:2133070/71‑1
             GCAGGGCAAACTCTGAGAGATCTTGAGAAATTAACGTCACCATCCCGGTATCATGTGGA                                                                          >  1:3713292/1‑59
                      ACTCTGAGAGATCTTGAGAAATTAACGTCACCATCCCGGTATCATGTGGAC                                                                         <  1:2582039/51‑1
                           GAGAGATCTTGAGAAATTAACGTCACCATCCCGGTATCATGTGGAC                                                                         <  1:310607/46‑1
                                                GTCACCATCCCGGTATCATGTGGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCTGC                                      >  1:1723118/1‑60
                                                    CCATCCCGGTATCATGTGGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCTGC                                      >  1:277680/1‑56
                                                          CGGTATCATGTGGACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCTGCATGTACGTGCCTTCCTCGGAC                 <  1:1708949/71‑1
                                                                      GACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCTGCATGTACGTGCCTTCCTCGGACTT               >  1:3276517/1‑61
                                                                      GACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCTGCATGTACGTGCCT                          <  1:2057839/50‑1
                                                                       ACGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCTGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCGG    >  1:1467947/1‑71
                                                                        CGTTAATTTCTCAAGATCTCTCAGAGTTTGCCCTGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCGGG   <  1:1316728/71‑1
                                                                         GTTAATTTCTCAAGATCTCTCAGAGTTTGCCCTGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCGGGA  >  1:4466996/1‑71
                                                                         GTTAATTTCTCAAGATCTCTCAGAGTTTGCCCTGC                                      <  1:281232/35‑1
                                                                                                                                                
AAGGCACGTACATGCAGGGCAAACTCTGAGAGATCTTGAGAAATTAACGTCACCATCCCGGTATCATGTGGACG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AP009048/1933514‑1933441
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gacgTTAATTTCTCAAGATCTCTCAGAGTTTGCCCTGCATGTACGTGCCTTCCTCGGACTTCCGGTTGGCGGGA  >  AP009048/1933468‑1933537

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