Cufflinks Or Alternative From Bwa Alignments?
3
0
Entering edit mode
10.6 years ago
Biostar User ▴ 360

Can I use Cufflinks starting from reads in sam/bam format aligned using bwa instead of bowtie? If not, what is the alternative to isoform assembly of RNAseq reads from bwa alignments?

cufflinks bwa • 4.4k views
ADD COMMENT
2
Entering edit mode
10.6 years ago

You probably mean Cufflinks rather than TopHat, right? TopHat is a mapper (e. g. FASTQ -> BAM), Cufflinks is a quantifier / assembler (BAM -> FPKM values and/or "transfrags").

You can use BWA-produced SAM/BAM files in Cufflinks, yes.

ADD COMMENT
0
Entering edit mode

Support this answer!

ADD REPLY
0
Entering edit mode

Yes, TopHat is a mapper... but uses bowtie to align the seeds to the genome. I think he is looking for a program that uses the same approach than TopHat, but using BWA instead of bowtie.

ADD REPLY
0
Entering edit mode

He specifically mentions isoform assembly of RNA-seq reads, and also writes about starting from aligned SAM/BAM reads.

ADD REPLY
0
Entering edit mode

Yes he did... But maybe he's a little confused about mapping/assembly process. Anyways, one of the two interpretations of the question is right, and either your answer or mine will help him :)

ADD REPLY
0
Entering edit mode

True :-) I agree that the original poster needs to get the process clarified.

ADD REPLY
0
Entering edit mode

I have edited the question now.

ADD REPLY
1
Entering edit mode
10.6 years ago
Geparada ★ 1.5k

You will have better result if you just use TopHat 2.0 because this version uses Bowtie 2.0 wich is better than BWA (see this resent published paper)

ADD COMMENT
0
Entering edit mode
10.6 years ago
Arun 2.4k

In theory you can: You could take your fastq and align it using bwa. And from the sam file you could filter out the aligned reads and re-construct a fastq file with the unaligned ones and supply that to tophat. However, as the first stage of mapping, tophat will use bowtie. But I'd guess that there shouldn't be many reads that are mapped using bowtie on unaligned reads.

But then, I also suggest (as Geparada) the approach of using tophat2 and bowtie2, although bowtie2 seems to be in beta version. I am testing this right now.

ADD COMMENT

Login before adding your answer.

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