Qualimap: java.lang.RuntimeException: The raw data doesn't exist for sample ?
2
6
Entering edit mode
8 months ago
dare_devil ★ 1.4k

I am trying to run qualimap on multiple bam file via command.

qualimap multi-bamqc -d bam_file.txt -gff target.interval_list

But it gives an error as follows:

qualimap multi-bamqc -d bam_file.txt -gff target.interval_list 
Java memory size is set to 120000M
Launching application...

OpenJDK 64-Bit Server VM warning: ignoring option MaxPermSize=1024m; support was removed in 8.0
QualiMap v.2.2.1
Built on 2016-10-03 18:14

Selected tool: multi-bamqc
Running multi-sample BAM QC

Checking input paths
Failed to run multi-bamqc
java.lang.RuntimeException: The raw data doesn't exist for sample: 001X1.sorted.deduped.bqsr.bam
Folder path:/media/molmed/Data/PJRCT1/bam
Please check raw data directory is present.

bam_file.txt conatins

head bam_file.txt

001X1.sorted.deduped.bqsr.bam   /media/molmed/Data/PJRCT1/bam
001X2.sorted.deduped.bqsr.bam   /media/molmed/Data/PJRCT1/bam
001M1.sorted.deduped.bqsr.bam   /media/molmed/Data/PJRCT1/bam
001M2.sorted.deduped.bqsr.bam   /media/molmed/Data/PJRCT1/bam
003X3.sorted.deduped.bqsr.bam   /media/molmed/Data/PJRCT1/bam
003M1.sorted.deduped.bqsr.bam   /media/molmed/Data/PJRCT1/bam

I also tried by giving the complete path to the file as /media/molmed/Data/PJRCT1/bam/001X1.sorted.deduped.bqsr.bam

and also by adding / to the folder ends as /media/molmed/Data/PJRCT1/bam/. But all gives the same error

qualimap bamqc multi-bamqc • 611 views
ADD COMMENT
6
Entering edit mode
8 months ago
dare_devil ★ 1.4k

Summary:

To run multi-bamqc in qualimap, you should have a bed file containg six columns: You can see the comments here

cat bed3.bed | perl -lane 'print "$F[0]\t$F[1]\t$F[2]\t.\t0\t."' > target.bed #this will contain 6 columns

The file containing name of the sample and paths should be as:

sample_1    /media/User/Data//bam/sample_1.bam #file name with extension `.bam`

and run the following command

qualimap multi-bamqc -r -d bam_file.txt -gff target.bed
ADD COMMENT
4
Entering edit mode
8 months ago
antonioggsousa ★ 2.2k

Not sure if it is related with copy/paste, but the file target.interval_list does not have any file extension.

Is this file a GFF/GTF or BED file?

Regarding the file bam_file.txt, the first column, I think it is the sample name that you want to give, not necessarily the same as the file name. In the second column you need to specify the bam file (this is my understanding from reading the docs - I never run this specific command), something like this:

sample_1   /media/molmed/Data/PJRCT1/bam/001X1.sorted.deduped.bqsr.bam

If you do not have BAM QC result you might need to specify the -r option:

Raw BAM files are provided as input. If this option is activated BAM QC process first will be run for each sample, then multi-sample analysis will be performed.

Source of this information.

I hope this helps,

António

ADD COMMENT
0
Entering edit mode

Thank you for the reply. I tried the same way you mentioned, but still gives the same error.. Interval file is in bed format.

ADD REPLY
1
Entering edit mode

The separator between both columns is a tab?

Did you try with and without the -r option?

ADD REPLY
2
Entering edit mode

I tried with -r and changing the bed file to six columns. It works!!.

If you post this as answer, I can accept and upvote

Thanks

ADD REPLY

Login before adding your answer.

Traffic: 1797 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6