## ------------------------------------------------------------------------ library(alpineData) ## ------------------------------------------------------------------------ ERR188297() ERR188088() ERR188204() ERR188317() ## ------------------------------------------------------------------------ eh <- ExperimentHub() query(eh, "ERR188") eh[["EH166"]] ## ------------------------------------------------------------------------ library(GenomicAlignments) gap <- ERR188297() barplot(sort(table(seqnames(gap))[1:25], decreasing=TRUE), las=3, main="Distribution of reads") ## ------------------------------------------------------------------------ gap1 <- gap[seqnames(gap) == "1"] starts <- start(first(gap1)) par(mfrow=c(2,2)) hist(starts,col="grey")