Does any one here have ever used Bioscope with Cufflinks to do DE analysis? I am trying to use the BAM files generated by Bioscope to do it but all my output files have zero expression values (fpkm 0.0). I am using the recommended parameter "--library-type fr-secondstrand" and I also tried to convert the BAM to SAM and add the XS:A[+-] flag but it didn't work also. I don't know if it is the cause but I have two of these errors: "SAM error on line nnnnnnn: CIGAR op has zero length". Thank you
•
written
by
Ramon Vidal •
1 answer
Doesn't directly answer your question, but I use following workflow to do DE analysis with SOLiD data:
- Do the alignment using LifeScope's Whole Transcriptome Mapping Analysis (I only need the BAM file)
- Do the gene count using HTSeq-count (http://www-huber.embl.de/users/anders/HTSeq/doc/overview.html)
- Do the DE analysis using DESeq2 (http://www.bioconductor.org/packages/2.12/bioc/html/DESeq2.html)
I very much prefer DESeq over Cufflinks; see for example Bullard et al., 2010 why the RPKM approach is not adivsable.
Log in to answer this question.
-
How To Sort Sam File
written by Chen Sun
-
Sam To Bam - Loss Of Data Or Just Great Compression?
written by Clare
-
Can Cufflinks Handle Bam Files From Outside?
written by thejustpark
-
Solution To Find Out Read Orientation As Well As Strand To Which Reads Align From Sam File
written by k2bhide
-
Samtools View Of A Sam File Entire Chromosome
written by dfernan
-
Error In Converting Sam To Bam By Samtools
written by Raghav
-
Problem Running Pindel With Bam Files
written by xenophiliuslovegood •
-
Cuffdiff Error: Cuffdiff Requires At Least 2 Sam Files
written by Lisanne •
-
Soap Result Converted To Bam Not Working In Samtool
written by Sarwar •
-
Converting Bam Or Sam Files Into Bed Files
written by Nancy •
There is also a thread here on SeqAnswers with similar issues: http://seqanswers.com/forums/showthread.php?t=3551
I always used TopHat with SOLiD RNA-seq data and it worked to get alignments ready for CuffLinks.