stretch 3/5p: anchor_end() shortens otherwise

plyranges=1.0.3 reported at tidyomics/plyranges#49
mschubert · Aug 7, 2018 · 2a1394e · 2a1394e
1 parent d1e2ba6
commit 2a1394e
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/cis_analysis/plot_ins.r b/cis_analysis/plot_ins.r
@@ -59,8 +59,8 @@ grtrack = GeneRegionTrack(txdb, name=args$gene)
 gtrack = GenomeAxisTrack(name="GRCm38.92")
 genes = seq$coords$gene(dset="mmusculus_gene_ensembl", granges=TRUE) # get this out of txdb obj?
 region = genes[genes$external_gene_name == args$gene] %>%
-    anchor_start() %>% stretch(as.integer(args$flank)) %>%
-    anchor_end() %>% stretch(as.integer(args$flank))
+    anchor_5p() %>% stretch(as.integer(args$flank)) %>%
+    anchor_3p() %>% stretch(as.integer(args$flank))
 aw = width(region) / 20
 
 dna_ins = io$load(args$dna_ins) %>%