Rsamtools to separate vcf file by chromosome?
I'm using Rsamtools to try and separate a vcf file by chromosome and while there are plenty of questions answered using samtools/vcftools/etc, I can't seem to adapt them to work in R. (I'm using Windows so I'm stuck to using R). I also have Rbamtools but that doesn't seem to have a split function.
Ideally I need this to loop as well.
Thanks in advance!
• 1,869 views
•
link
written
by
nadrianam •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Filter BAM by chromosome using PySam
written by pablosolar.r •
-
pca in r using plink files
written by stalberg •
-
How to align multiple fastq files with Hisat2 ?
written by sarai •
-
How to automate the conversion of a large number of .sra files to .fastq?
written by Jeff M •
-
Input own BAM files
written by siangoldie •
-
Does bamtools split preserve sort and rmdup?
written by ccnn •
-
Split bwa mem for MergeBamAlignment
written by dr.dimo •
-
Performing PGLS for every tree in a .tre file
written by davidstupski •
-
Reading a .ped file in R?
written by nadrianam •
-
Iterating Vcftools Commands With Shell Script
written by vitti •
It's good practice to show us what you tried, what didn't work and how it didn't work :)