New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? AP009048 = 356873429 (0.520)3 (0.060) 3/90 NT 10.4% coding (1330/2037 nt) yihQ alpha‑glucosidase
?AP009048 = 3568756 26 (0.520)coding (1352/2037 nt) yihQ alpha‑glucosidase

GTCAGTGCCGAAATTATGCATAACGCCTGGCCTGCGCTGTGGGCGAAGTGTAACTACGAAGCCCTTGAAGAAA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AP009048/3568662‑3568734
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑aagaaaAGGATCTCGCCGAGCTTGCCCGTTTCTTCAAGGGCTTCGTAGTTACACTTCGCCCACAGCGCAGG  <  AP009048/3568756‑3568692
                                                                                                                                          
GTCAGTGCCGAAATTATGCATAACGCCTGGCCTGCGCTGTGGGCGAAGTGTAACTACGAAGCCCTTGAAGA                                                                     <  1:66026/71‑1
 TCAGTGCCGAAATTATGCATAACGCCTGGCCTGCGCTGTGGGCGAAGTGTAACTACGAAGCCCTTGAAGAA                                                                    <  1:2825783/71‑1
   AGTGCCGAAATTATGCATAACGCCTGGCCTGCGCTGTGGGCGAAGTGTAACTACGAAGCCCTTGAAGA                                                                     <  1:3641234/68‑1
     TGCCGAAATTATGCATAACGCCTGGCCTGCGCTGTGGGCGAAGTGTAACTACGAAGCCCTTGAAGAA                                                                    <  1:4583186/67‑1
                       CGCCTGGCCTGCGCTGTGGGCGAAGTGTAACTACGCAGCCCTTGAAGAAA                                                                   >  1:4168233/1‑50
                              CCTGCGCTGTGGGCGAAGTGTAACTACGAAGCCCTTGAAGAAAAGGATCTCGCCGAGCTTGCCCGTTTC                                         <  1:3156478/69‑1
                                TGCGCTGTGGGCGAAGTGTAACTACGAAGCCCTTGAAGAAAAGGATCTCGCCGAGCTTGCCCGTTTCTTC                                      <  1:3416444/70‑1
                                   GCTGTGGGCGAAGTGTAACTACGAAGCCCTTGAAGAAAAGGATCTCGCCGAGCTTGCCCGTTTCTTCA                                     >  1:4964538/1‑68
                                                                   AAGAAAAGGATCTCGCCGAGCTTGCCCGTTTCTTCAAGGGCTTCGTAGTTACACTTCGCCCACAGCGCAGG  >  1:3828661/1‑71
                                                                   AAGAAAAGGATCTCGCCGAGCTTGCCCGTTTCTTCAAGGGCTT                              <  1:3119054/43‑1
                                                                    AGAAAAGGATCTCGCCGAGCTTGCCCGTTTCTTCACGGGCTTCGTAGTTACACTTCGCCCACAGCGCAGG  >  1:2024156/1‑70
                                                                    AGAAAAGGATCTCGCCGAGCTTGCCCGTTTCTTCAAGGGCTTCGTAGTTACACTTCGCCCGCAGCGCAGG  >  1:21387/1‑70
                                                                    AGAAAAGGATCTCGCCGAGCTTGCCCGTTTCTTCAAGGGCTTCGTAGTTACACTTCGCCCACA         >  1:1041368/1‑63
                                                                    AGAAAAGGATCTCGCCGAGCTTGCCCGTTTCTTCAAGGGCTTCGTAGTTACACTTCGCCCACA         >  1:396153/1‑63
                                                                      AAAAGGATCTCCCCGAGCTTGCCCGTTTCTTCAAGGGCTTCGTAGTTACACT                  >  1:4994758/1‑52
                                                                                                                                          
GTCAGTGCCGAAATTATGCATAACGCCTGGCCTGCGCTGTGGGCGAAGTGTAACTACGAAGCCCTTGAAGAAA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AP009048/3568662‑3568734
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑aagaaaAGGATCTCGCCGAGCTTGCCCGTTTCTTCAAGGGCTTCGTAGTTACACTTCGCCCACAGCGCAGG  <  AP009048/3568756‑3568692

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