Question: NGS DNA-seq analysis packages for mapping, calling cariants, annotation on R / bioconductor
gravatar for houkto
4.1 years ago by
houkto60 wrote:


I am looking for a complete DNA-seq analysis in R this including reads mapping, variants calling and annotation ?

So far I found RNA-seq and Chip-seq packages only?

I would be grateful if you can points me to DNA-seq packages to analyse whole exome or targeted genes data (fastq) in R and preferable work through tutorial if exists ?



next-gen dna-seq R • 2.0k views
ADD COMMENTlink modified 4.1 years ago by Sean Davis25k • written 4.1 years ago by houkto60

I really can't recommend using R for alignment and variant calling. R is best at data analysis, not batch processing of large datasets. That's why very few people try to do their whole pipeline in R and rather use R at the end.

BTW, for alignment you could use Rsubread.

ADD REPLYlink written 4.1 years ago by Devon Ryan89k

Thanks, for the info. I did realise the issue with big data and R. Nevertheless, using Linux OS (whether fresh install or as virtual image) is not practical when it comes to demonstrate to colleges how basic DNA-seq analysis works especially in old machines.

ADD REPLYlink written 4.1 years ago by houkto60
gravatar for Sean Davis
4.1 years ago by
Sean Davis25k
National Institutes of Health, Bethesda, MD
Sean Davis25k wrote:

Targeted sequencing or exome sequence alignment and variant calling is best done using non-R-based tools (at least right now).  One can use R to call system programs such as BWA and the GATK tools to do analysis up to VCF files.  At that point, the VariantAnnotation package becomes quite useful for QC, annotation, filtering, and visualization of results.  

ADD COMMENTlink written 4.1 years ago by Sean Davis25k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1495 users visited in the last hour