How to extract allel, Genotype from vcf file using python or other language for 23GB files?
Well, I am able to right script to get allele but for large VCF file it's difficult? what should other possible way to get allele, Genotype information?
•
written
by
sukhjindermultani85 •
2 answers
See bcftools query.
EDIT: WIth bcftools query you can print any information you like. So in your case e.g.:
$ bcftools query -f '%CHROM %POS %REF %ALT [ %GT]\n' input.vcf
The output looks now like this:
chr1 10177 ACC ACCC 0/1
chr1 10327 T C 0/0
chr1 10352 TAC TAAC 1/1
chr1 12783 G A 1/1
fin swimmer
Hello Ram,
if an "answer" is just intended for full copy&paste solution then my post is indeed more a comment. But I thought that telling the tool with it's subcommand and linking to the good manual is an answer enough.
I extended my post now to an full answer :)
cpad was faster than me, right. I didn't saw his answer as I haven't reload the page.
fin swimmer
Life is not about fast pace. Those who take their time will metaphorically 'make it' first.
Hi finswimmer! Can I ask for the trick to convert the output to symbolic genotypes? for your example:
chr1 10177 ACC ACCC ACC/ACCC
chr1 10327 T C T/T
chr1 10352 TAC TAAC TAAC/TAAC
chr1 12783 G A A/A
Searched for a whole, but just did not have my luck.
Hello yifangt86 ,
that's also described in the manual I've linked to:
$ bcftools query -f '%CHROM %POS %REF %ALT [ %TGT]\n' input.vcf
fin swimmer
Extracting genotype information using R.
library(vcfR)
vcf <- read.vcfR(vcf_file, verbose = FALSE )
gt <- extract.gt(vcf, element = c('GT'), as.numeric = TRUE)
For python take a look at the following article.
http://alimanfoo.github.io/2017/06/14/read-vcf.html
Genotypes can also be extracted using SnpSift.jar in snpEff using the following command.
java -jar ../snpEff/SnpSift.jar extractFields annotated.vcf CHROM POS REF ALT "GEN[*].GT" > output.tsv
Doesn't look like vcfR does streaming read, so I would not recommend it as it's not a great idea to build an in-memory object of an entire VCF file. A better strategy would be to use closer-to-bare-metal tools such as bcftools to extract information, then use R or Python to compute on extracted information.
Log in to answer this question.
-
How can I view the information in an alternative allele?
written by iarmir •
-
using Haplotypearray in scikit-allel
written by peter •
-
How to get REF and ALT alleles from a genotype data?
written by biotechnology415 •
-
is it possible to know individual genotype(alleles) by SNP in plink?
written by nayeona.hi
-
SNP density per genes
written by rthapa •
-
Merging and converting multiple vcf files into a SNP array
written by Arko •
-
SnpMatrix from VCF file
written by Dhana
-
Adding allele frequency to the genotype field of a VCF file
written by Nonny27 •
-
Extract each annotation in 'info' column in snpEff output
written by jan
-
Extracting Genotype Information From Vcf
written by rickyflintoff •
how about VCFtools?
Why is this a
toolpost? A question about tools should be aquestion-type post, not atool-type post.What have you tried?
May help the user (AWK ideas):A: How to get sample names and genotype for SNP in multi-sample VCF fileActually, I have a Python script that can parse a VCF, in fact: Filtering VCF with python
Why have you replied to my comment, Kevin?
Did not want to create yet another 4th and independent comment
You can take a look at this two scripts wrote in python to split a vcf and select what you want : A: VCF file help and C: parsing vcf file
try
bcftools query.