Error While Running Cuffdiff On Multiple Bam Files Of A Group
1
0
Entering edit mode
7.5 years ago
Jordan ★ 1.2k

Hi,

I'm trying to run cuffdiff on 6 samples, divided into two groups of 3 samples each - G1(Female), G2(Male). But, it does not run. I gave the following command:

cuffdiff /refs/mm9/Annotation/Genes/genes.gtf -L G1,G2 ~/Female/sampleA/accepted_hits.bam,~/Female/sampleC/accepted_hits.bam,~/Female/sampleE/accepted_hits.bam ~/Male/sampleG/accepted_hits.bam,~/Male/sampleI/accepted_hits.bam,~/Male/sampleK/accepted_hits.bam

I get this error:

You are using Cufflinks v2.1.1, which is the most recent release.
open: No such file or directory
File ~/Female/sampleC/accepted_hits.bam doesn't appear to be a valid BAM file, trying SAM...
Error: cannot open alignment file ~/Female/sampleC/accepted_hits.bam for reading

I have generated these BAM files using tophat2. They are mouse genome. But they seem to work when I give single bam file of the same group. Like the following command.

cuffdiff /refs/mm9/Annotation/Genes/genes.gtf -L G1,G2 ~/Female/sampleA/accepted_hits.bam ~/Male/sampleG/accepted_hits.bam

Can anyone tell me how to run multiple bam files of a group on cuffdiff?

cuffdiff tophat2 rna-seq • 5.6k views
ADD COMMENT
0
Entering edit mode
7.5 years ago

replace '~' by ${HOME}.

$ echo A,~/B
A,~/B

$ echo A,${HOME}/B
A,/home/lindenb/B
ADD COMMENT
0
Entering edit mode

That worked. Thanks! But out of curiosity, can you explain why that works?

ADD REPLY
1
Entering edit mode
ADD REPLY
0
Entering edit mode

That helped. Thanks again!

ADD REPLY

Login before adding your answer.

Traffic: 1320 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