NGS plot error
1
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
ChIP-Seq
next-gen-sequencing
• 5.0k views
•
link
updated 2.3 years ago by
Ram
44k
•
written 9.1 years ago by
Constantine
▴
290
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 clustering:
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.
•
link
updated 2.3 years ago by
Ram
44k
•
written 9.1 years ago by
Fidel
★
2.0k
Login before adding your answer.
Traffic: 2369 users visited in the last hour
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
"computeMatrix: error: argument : invalid choice: 'tss' (choose from 'scale-regions', 'reference-point')" now this doesn't work .
But when i using the something like this
It takes a lot of time but when i choose like chr21 this for my bed region it generates the matrix ...how to reduce the matrix generation time ?