Biostar Beta. Not for public use.
Rsubread versus BWA ( Burrows-Wheeler Alignmener) for whole genome alignment
Entering edit mode
12 months ago
tarek.mohamed • 250

Hi All,

I have some whole genome sequencing data that were generated from illumina hiseq platform.

as a first step, I will align the fastq file to a reference genome. Can I use Rsubread for that? I have used it for RNAseq mapping before and it was very good, I got 90% of the read were mapped to the ref genome. Does BWA have any advantage as compared to Rsubread.



Entering edit mode
16 months ago
Gordon Smyth • 780

This question was answered on Bioconductor:

In summary, Rsubread is a general purpose aligner and is fine for DNA as well as RNA. For DNA, use the Rsubread align() function with `type="dna".

More details about Rsubread are in: The article focuses on RNA-seq but DNA-seq is mentioned as well.

We compared Subread with BWA in the original Subread publication:


Login before adding your answer.

Similar Posts
Loading Similar Posts
Powered by the version 2.1