This is a test version of Biostars. For the public version, visit https://www.biostars.org.
DeepTools Error: Be sure that the bam files you are using is indexed.

I am attempting to use the bamCorrelate function of deeptools to find the reproducibility between two chip seq replicates.

I have used samtools index BAMFILE to index my files, but I am still running into this error so I must not be doing this correctly.

Here is the command I used for deeptools

bamCorrelate bins \
  --bamfiles /Volumes/mydata/2015108/IDRM/2015_10_21_0/analysis_results/Sample_Pol_II_9_29/Alignments/Pol_II_9_29.mapped.rgid.sorted.filtered.dups_remove.bam.bai \
  /Users/jonathanreeder/carlos-projects/Dropbox/ChIP-Seq/hg19/old_data/Sample_Pol-II-Chip-N20/Alignments/Pol-II-Chip.mapped.rgid.sorted.filtered.dups_remove.bam.bai \
  -out pol-correlation

Here's the error:

Be sure that the bam file you areusing is indexed.
deeptools
• 4,006 views • link
updated by Ram Ram • written by dally dally

1 answer

You are giving it the index files (*.bam.bai), when you should probably be giving it the bam files (*.bam).

• 0 views • link
updated by Ram Ram • written by matted matted

Ahh I see. I fixed that and was able to run the command line just fine. However, I can't seem to get it to work when running the --corMethod pearson parameters. It keeps telling me that it is an unrecognized argument. But everything on the github has the parameters written exactly like this. Any ideas?

The --plotFile command appears to be unrecognized as well.

• 0 views • link
updated by Ram Ram • written by dally dally

I'm not sure. Maybe you have an old version installed?

• 0 views • link
updated by Ram Ram • written by matted matted

Log in to answer this question.

More posts like this
  • Same data, very different heatmaps
    written by Sam Sam
    Me and a ex lab member generated two heatmaps using the same Pol II dataset using HOMER and deeptools. I used deeptools, and the ex …