Hi,
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?
Thanks
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.
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.