New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? AE015451 = 519589543 (1.050)4 (0.100) 4/320 NT 9.8% coding (620/690 nt) ybgJ allophanate hydrolase subunit
?AE015451 = 5195952 33 (0.850)coding (563/690 nt) ybgJ allophanate hydrolase subunit

CTCCACATCACCACCCATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AE015451/5195837‑5195895
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcgcacccCGGTGCGGTTGTTCGACCGTGAGCGTGAGGGTTACAGCCTCTTGCAA  <  AE015451/5195952‑5195906
                                                                                                          
CTCCACATCACCACCCATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCCGGTGCGGTTGTTCGACCGTGAGCGTGAGGGTTACAGCCTCT       >  1:606985/1‑101
CTCCACATCACCACCCATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCCGGTGCGGTTGTTCGACCGTGAGCGTGAGGGTTACAGCCTCT       <  2:606984/101‑1
                ATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCCGGTGCGGTTGTTCGACCGTGAGCGTGAGGGTTACAGCCTCTTGCA   >  1:590302/1‑89
                ATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCCGGTGCGGTTGTTCGACCGTGAGCGTGAGGGTTACAGCCTCTTGCA   <  2:590301/89‑1
                                                         gCCGGTGCGGTTGTTCGACCGTGAGCGTGAGGGTTACAGCCTCTTGCAA  <  2:14604/48‑1
                                                                                                          
CTCCACATCACCACCCATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AE015451/5195837‑5195895
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcgcacccCGGTGCGGTTGTTCGACCGTGAGCGTGAGGGTTACAGCCTCTTGCAA  <  AE015451/5195952‑5195906

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