Creating a coverage plot using BEDTools and R

One of my Top 10 posts is on creating a coverage plot using R. For that post I used CAGE data, which is a transcriptomic data set containing transcription start sites, and I used R exclusively for building a “coverage plot.” The main issue with that post was that the plots were density plots rather…

Continue Reading

BEDTools puzzles

Updated 2015 December 16th: the files for the puzzles can be downloaded from the Amazon cloud Roughly two weeks ago I came across this excellent BEDTools tutorial and saw some puzzles or homework questions at the end of the tutorial; naturally I tweeted about it: Oh this looks like fun; BEDTools puzzles http://t.co/UMowwdPLEf— Dave Tang…

Continue Reading

Defining genomic regions

Updated 2014 June 24th to use GENCODE version 19 RNA sequencing (RNA-Seq) reads are typically mapped back to the genome (or transcriptome in some cases) after sequencing. The next task is to annotate the reads, to see which regions the reads mapped to. Typically one creates an annotation file and compares the coordinates of the…

Continue Reading