Tweet
When I convert the bam file from the output of tophat 1.1.4 to a sam file, it loses the strand specificity ie,, XS:A+/- at the end. This renders it unusable for cufflinks. The conversion I used is samtools view -h -o out.sam in.bam. Any suggestions?
-
-
The latest version of cufflinks directly accepts the bam file produced by TopHat. you do not need to conver the bam file into sam file. I just did this a couple of days ago.
$./cufflinks -G refGene.gtf accepted_hits.bam
This works quite well -
I am having a similar problem with converted output .sam files. And when I input my .bam directly into cufflinks, I get the following within a second before the process terminates:
Do .bam files have a CIGAR string? Has anyone else encountered this problem?Counting hits in map
CIGAR op has zero length
CIGAR op has zero length
CIGAR op has zero length
Total map density: 0.000000
CIGAR op has zero length
CIGAR op has zero length
CIGAR op has zero lengthComment
-
The sequencing world is rapidly changing due to declining costs, enhanced accuracies, and the advent of newer, cutting-edge instruments. Equally important to these developments are improvements in sequencing analysis, a process that converts vast amounts of raw data into a comprehensible and meaningful form. This complex task requires expertise and the right analysis tools. In this article, we highlight the progress and innovation in sequencing analysis by reviewing several of the...05-06-2024, 07:48 AM -
The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist...04-22-2024, 07:01 AM
Comment