Question: Find SNPs among samples
0
gravatar for asgara
7 months ago by
asgara0
asgara0 wrote:

Hello!

I would like to calculate the number of SNPs in a sample (.vcf file) only located in specific genomic regions (which I have as .bed file). Also, having two samples (as 2 different .vcf files), I would like to do the same as above but reporting only the SNPs which are present in sample 2 and not in sample 1, and still only located in the genomic regions provided as .bed file.

They are small files and I was thinking if that can be solved with some Python code.

Any suggestions?

I hope the description of the problem was understandable and not too confusing.

Thanks!

python snp genome gene • 143 views
ADD COMMENTlink modified 7 months ago by lager0 • written 7 months ago by asgara0

have a look at bcftools view read the manual about --regions-file and --include . http://samtools.github.io/bcftools/bcftools.html

ADD REPLYlink written 7 months ago by Pierre Lindenbaum131k

Thanks for the answer! I know there are some existing tools which can probably do it smoothly, but I was wondering if the same results can be obtained also with some Python code.

ADD REPLYlink written 7 months ago by asgara0
0
gravatar for lager
7 months ago by
lager0
lager0 wrote:

try VCFtools

ADD COMMENTlink modified 7 months ago • written 7 months ago by lager0
Please log in to add an answer.

Help
Access

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