bioconductor package for BWA-MEM ?
1
0
Entering edit mode
K ▴ 50
@k-8495
Last seen 3.1 years ago
United States

Hello,

I was wondering if there is a R/Bioconductor package that can run BWA-mem alignment (just like how RSamtools is the Bioconductor implementation of Samtools).

I am looking to consolidate all my analyses in R/Bioconductor, instead of having half my things in unix and other parts in R.

Thanks,

bwa ngs • 2.8k views
ADD COMMENT
0
Entering edit mode
@gordon-smyth
Last seen 47 minutes ago
WEHI, Melbourne, Australia

See the BiocViews page for "alignment":

 http://www.bioconductor.org/packages/release/BiocViews.html#___Alignment

As far as I know, there is no Bioconductor implementation of BWA.

The Rsubread package provides a native implementation of the subread aligner.

The QuasR and Rbowtie packages provide R wrappers for the bowtie and SpliceMap aligners.

ADD COMMENT
0
Entering edit mode

Thanks for the response. I will go through the packages you suggested.

ADD REPLY
0
Entering edit mode

There is also gmapR, too ;-)

ADD REPLY
0
Entering edit mode
Also, if the variant tolerance of BWA is important to you then the gsnap aligner from the gmapR package might be a great alternative. If you need BWA directly and you don't mind using a system command-based interface then the systemPipeR package might provide some solutions. The main and VAR-Seq vignettes give some examples how to run bwa mem. Thomas
ADD REPLY
0
Entering edit mode

Thanks ! The systemPipeR might be exactly what I am looking for...I will definitely try this out. Thanks for all the suggestions.

ADD REPLY

Login before adding your answer.

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