New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? AP009048 113675 =45 (0.770)3 (0.060) 3/94 NT 6.5% coding (232/1044 nt) guaC GMP reductase
?AP009048 113695 = 45 (0.830)coding (252/1044 nt) guaC GMP reductase

ACCATCACATGTTTCAGCACATCAGCGGAAGAATTGTTGATAAACGCTTGCCACTCTTCGACAGAATAGTG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AP009048/113745‑113675
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑agtgGCAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTG  >  AP009048/113695‑113762
                                                                                                                                           
ACCATCACATGTTTCAGCACATCAGCGGAAGAATTGTTGATAAACGCTTGCCACTCTTCGACAGAATAGTG                                                                      >  1:4393635/1‑71
                    ATCAGCGGAAGAATTGTTGATAAACGCTTGCCACTCTTCGACAGAATAGTG                                                                      <  1:1172200/51‑1
                                                 GCCACTCTTCGACAGAATAGTGGCAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATG                     >  1:2968983/1‑71
                                                    ACTCTTCGACAGAATAGTGGCAAGCGTTTATCAACAATT                                                  <  1:1433861/39‑1
                                                     CTCTTCGACAGAATAGTGGCAAGCGTTTATCAACAATTCTTCC                                             >  1:930399/1‑43
                                                                    GTGGCAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTG  <  1:4747314/71‑1
                                                                     TGGCAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGT         <  1:3972046/63‑1
                                                                                                                                           
ACCATCACATGTTTCAGCACATCAGCGGAAGAATTGTTGATAAACGCTTGCCACTCTTCGACAGAATAGTG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AP009048/113745‑113675
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑agtgGCAAGCGTTTATCAACAATTCTTCCGCTGATGTGCTGAAACATGTGATGGTTTCTACCGGTACGTCTG  >  AP009048/113695‑113762

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