Question: Cuffdiff doesn't accepted sorted bam files that worked in cufflinks
1
gravatar for lzswindon
3.0 years ago by
lzswindon10
lzswindon10 wrote:

I have successfully run through cufflinks and cuffmerge on my STAR aligned data, but am running into an error on cuffdiff:

cuffdiff merged.asm/merged.gtf \
H_400Pa_2079out/Aligned.out.sorted.bam.bam, H_400Pa_2085out/Aligned.out.sorted.bam.bam \
H_60kPa_2080out/Aligned.out.sorted.bam.bam, H_60kPa_2086out/Aligned.out.sorted.bam.bam

This error repeatedly appears on the screen:

BAM record error: found spliced alignment without XS attribute

I had encountered this readout when running cufflinks, but addressed it by using samtools sort on my original bam files. I'm not sure how to address this new occurrence, since my bam files are already sorted. Help would be appreciated! Thanks.

bam cuffdiff cufflinks • 1.2k views
ADD COMMENTlink modified 3.0 years ago • written 3.0 years ago by lzswindon10
0
gravatar for lzswindon
3.0 years ago by
lzswindon10
lzswindon10 wrote:

Ah, silly me, I left out -library-type fr-firststrand in my cuffdiff command, which I'd used in cufflinks. So full command should be:

cuffdiff -library-type fr-firststrand merged.asm/merged.gtf \
H_400Pa_2079out/Aligned.out.sorted.bam.bam, H_400Pa_2085out/Aligned.out.sorted.bam.bam \
H_60kPa_2080out/Aligned.out.sorted.bam.bam, H_60kPa_2086out/Aligned.out.sorted.bam.bam
ADD COMMENTlink written 3.0 years ago by lzswindon10
Please log in to add an answer.

Help
Access

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