Latest
Open
RNA-Seq
ChIP-Seq
SNP
Assembly
Tutorials
Tools
Jobs
Forum
Planet
All »
Home
Log In
Question: NGS plot error
0
gravatar for Constantine
3.4 years ago by
Constantine • 220
USA
Constantine • 220 wrote:

Hello 

I'm trying to plot the coverage of a bam file from a ChIP-seq experiment via NGSplot and I get the following error:

Configuring variables...Error in CheckRegionAllowed(reg2plot, default.tbl) :
  Unknown region specified. Must be one of: bed
Execution halted

This is the command I'm using:

ngs.plot.r -G mm10 -R tss -C myfile.bam -O mybamfile.ngs

If I change the -R from "tss" to "bed" and specify a bed file I get from the PeakCaller software, it works just fine. However, I want it to take into consideration only the bam file

Any ideas why I'm getting this error message?

Thanks

 
sequencing chip-seq next-gen • 2.6k views
ADD COMMENT • link •
modified 3.4 years ago by Fidel • 1.9k • written 3.4 years ago by Constantine • 220
1
gravatar for Fidel
3.4 years ago by
Fidel • 1.9k
Germany
Fidel • 1.9k wrote:

The problem could be that your region's file does not end in .bed?

Otherwise, I suggest you to use deepTools instead (note: I am part of the developing team). Simply:

pip install deepTools

bamCoverage -b myfile.bam -o myfile.bw

computeMatrix -R tss -S myfile.bw -o matrix

profile -m matrix -o profile.pdf

 

Also, once you have a matrix you can plot a heatmap and add some clustring:

heatmapper -m matrix --kmeans 2 -o heatmap.png

You can take a look at the documentation (https://github.com/fidelram/deepTools/wiki) for more information.
ADD COMMENT • link modified 3.4 years ago • written 3.4 years ago by Fidel • 1.9k
1

I'm mainly using deeptools (prefer it better) ;)

Just wanted to also use NGS to see whether I'm getting the same output...Guess I'll have to stick with deeptools then

ADD REPLY • link written 3.4 years ago by Constantine • 220
Please log in to add an answer.

Similar posts • Search »

Content
Help
Access
Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1849 users visited in the last hour