Questions tagged [vcf]
VCF, acronym for Variant Calling Format, is a text file type used to store information regarding genetic variants.
208
questions
0
votes
0
answers
27
views
Extracting information from VCF file for many specific positions in specific chromosomes?
I have an excel file that I created from VCF file for common SNPs across 6 samples. This excel have the chromosomes and the position of the SNPs only (see example table1)
Now I would like to obtain ...
1
vote
0
answers
38
views
What is the current state-of-the-art for classifying the ancestry groups of VCF files?
I'm wondering which tool or program is:
Most updated; doesn't use a superseded algorithm
Runs very quickly (ideally, under 30s per VCF file)
Able to be automated in some way for a large number of ...
- 96
2
votes
2
answers
65
views
Where do I get a large reference VCF?
I would like to download a large .vcf file containing many (hundreds or thousands) of samples. Ideally, I would download different population-specific .vcf files, but the ability to sort/filter by ...
- 96
3
votes
3
answers
298
views
Strategy for merging many VCFs
Can anyone suggest a strategy for speeding up VCF merging?
I have ~44,000 single-sample VCFs that I am trying to merge into a multi-sample VCF with bcftools merge, ...
- 643
0
votes
0
answers
50
views
How to convert VCF files to input file format required by XP-CLR1.0 software?
I attempt to do a selective sweep analysis by using the software XP-CLR1.0. My genotypic data is VCF format, so I wonder if someone with experience on this package could recommend any tool to convert ...
- 1
1
vote
0
answers
48
views
Title: Error merging VCF files: "The REF prefixes differ: G vs A"
I'm trying to merge two VCF files, but I'm running into an error message that says ...
- 11
2
votes
2
answers
181
views
How to generate rsid values locally from VCF file
I recently used Dante to obtain a WGS and have downloaded several files produced during this process. There is a file that ends '.vcf' which I did think was a VCF file but I'm not sure now after a ...
- 131
1
vote
1
answer
87
views
Converting a Dante Lab VCF file
Is there an easy way to grab the rsids (hg38) of a VCF file? I know one of the tabs can contain RSIDs, but this file doesn't contain any.
I have a software I created for mapping 23AndMe DTC testing to ...
- 13
3
votes
2
answers
80
views
Non-ACGTN alternate alleles in VCF
I was trying to left normalize my vcf file using following command:
bcftools norm -f reference.fasta -m -any input.vcf > output.norm.vcf
However this gives the ...
- 158
1
vote
1
answer
82
views
variant visualisation of VCF files for 30+ samples
I am relatively new to genome analysis and would like to compare 30+ vcf files, 1 sample per file. After filtering there are about ~20 variants per vcf file.
I would like to evaluate to mutational ...
- 158
0
votes
0
answers
64
views
How to Split VCF File into male and Female groups
Cross-posted on biostars
I have trio_example.vcf file and I want to split them into male and female groups. I tried this:
$ bcftools query -l trio_example.vcf
And ...
- 27
1
vote
1
answer
14
views
what does ALT ID DUP refer to when ALT ID for copy number are already specified?
ALT ID CN2, CN3,CN4,CN5 are regions of elevated copy number relative to the reference, Then what does the ID =DUP separately cover? I am new to these analyses. please help
...
- 13
2
votes
1
answer
242
views
How to Recalculate Allele Frequency from a VCF File?
I have a VCF file with 200 samples (mitochondrial genome of Plasmodium falciparum). Here is a pic to take a look at:
And a few relevant lines from the actual file:
...
- 41
3
votes
2
answers
237
views
Extract lines from a file based on a range of coordinates listed on a second file
I would be highly appreciated to have a script able to perform the following operation:
Only extract scaffold_1 lines from file1.vcf in which POS number (column 2) falls within column 2 and column 3 ...
- 33
2
votes
3
answers
108
views
Interpreting short indel calls in 1000 Genomes Project VCFs
Consider the following short indel polymorphism rs59679400 on chr7.
...
- 5,010
2
votes
1
answer
50
views
Availability of information on genes in Gnomad VCF data
I was wondering does the gnomad genome data that is downloaded in the vcf format on variants contain information on what is the nearest gene and is the genomic location available exome/intron?
if yes ...
2
votes
2
answers
50
views
VCF: false insertions?
I'm working with gvcf files, containing also non-variant positions. I have noticed some insertions that I'm not able to correctly interpret or handle. An example (selected samples from a huge vcf):
<...
- 123
0
votes
1
answer
77
views
Maintaining SNP names when calculating allele frequencies from a VCF file on plink 1.9
I downloaded the following vcf file from https://www.internationalgenome.org/.
http://ftp.1000genomes.ebi.ac.uk/vol1/ftp/release/20130502/ALL.chr22.phase3_shapeit2_mvncall_integrated_v5b.20130502....
1
vote
1
answer
285
views
Need to convert .bed to .vcf. Can the reference build (and necessary .fasta) be determined for a .bed file?
I have some .bed files (and .bim and .fam files) containing data for a number of different samples, and I need to convert them to .vcf.
I've found bed2vcf which is ...
0
votes
0
answers
64
views
Help sought with GWAS and vcf files, lack phenotype labels
This question has also been asked on Biostars
Hi, I am very new to this area, and I am taking a class about bioinformatics. For an independent project assignment, I need to do a GWAS. I am using the ...
3
votes
1
answer
214
views
Filter with bcftools
I want to filter a SNP, specifically CHROM1:POS:630128. I tried to use bcftools for this.
I came up with the following variant ...
- 145
2
votes
1
answer
45
views
Can I discard missing alleles?
I am converting a biallelic VCF into a SQL table and one of my tables will be something like:
...
- 123
1
vote
1
answer
68
views
Is picking variants belonging to a specific group of genomes from VCF file possible?
GnomAD has publicly available VCF files that have variant data gathered from over 15k individuals. For a university project, I want to randomly select 1k of these individuals and get the variant data ...
- 13
3
votes
1
answer
138
views
PL and QUAL values on VCF file?
I want to filter my VCF file to include only the relevant information I need but I have some questions about the results I got.
This is what the first 10 lines of my VCF file looks like:
...
- 494
2
votes
0
answers
78
views
BAF and LRR calculation and usage
I'm working with a vcf file, in particular I'm trying to identify a gene duplication in a sequenced genome.
I studied bcftools cnv command, but for this tool is needed a vcf file containing a BAF and ...
- 21
1
vote
0
answers
44
views
Benchmarking for variant identification using RNA-seq data
I am in need to benchmark the variant identification pipeline which uses RNA seq data alone without any matched-normal. I would like to know the reference dataset (and the pipeline on which the ...
- 11
1
vote
1
answer
73
views
Is it common to get different number of SNVs+Indels across samples from vcf files generated using GATK and DRAGEN (counts are higher for GATK)?
We have 2 vcf (Whole Exome Sequencing (WES) data; germline samples) files (e.g., vcf_1 and vcf_2). vcf_1 was generated (Ref. genome: hg38) using the GATK pipeline for 250 children and their parents ...
- 11
2
votes
1
answer
338
views
Compare VCF files
I have two strains, X and Y, and I would like to compare the variants found in both. I have my own parent strain, Z, that was sequenced and assembled to be used as the reference genome.
I first ...
- 494
0
votes
1
answer
137
views
2
votes
0
answers
100
views
SNPs with high population differentiation from 1k Genome dataset
I am trying to reproduce the results from this paper "Human genomic regions with exceptionally high levels of population differentiation identified from 911 whole-genome sequences".
...
- 21
1
vote
0
answers
83
views
Production database system for storing vcf records
I’m standing up a production bioinformatics pipeline that’ll be producing standard VCF files. I’d like to store this data in an SQL database that ingest the VCFs and allow me and my team to easily ...
- 181
1
vote
0
answers
95
views
Converting from VCF to PLINK while maintaining the SNP names from the VCF file
Just as the title says, I am working in PLINK 1.9 trying to convert the following VCF into a PLINK file (binary or not either is fine), while maintaining the SNP names from the VCF file. I am on ...
1
vote
0
answers
98
views
Polygenic Risk score calculation using vcf files
I have about 750 vcf files and I am looking to generate a Polygenic Risk Score for each individual. I am planning to use either PLINK or LDPhred to generate a Risk Score. I have phenotype information ...
- 11
1
vote
2
answers
89
views
How can I get latest .vcf files with annotation data?
I'd like to perform the annotation (1000genomes, COSMIC etc.) of my variants using SnpSift and SnpEff, however so far all I get are vcf files for separate chromosomes:
http://ftp.1000genomes.ebi.ac.uk/...
- 206
3
votes
3
answers
192
views
Compare my VCF to gnomAD variants
I have a VCF with small variant calls against HG38 and I would like to determine which of those calls are present in the gnomAD database.
Is there an existing tool that can do this? Should I be ...
- 133
1
vote
1
answer
101
views
Are the GIAB reference genomes (HG00N) fully phased?
Specifically, the NISTv4.2.1 HG002 dataset here (more info on Genome in a Bottle). According to their README, the HG002 v4.2.1 VCF should be fully phased.
For HG001, HG002, and HG005, we have phased ...
- 282
1
vote
0
answers
40
views
Weird data in VCF file
I am using GATK HaploTypeCaller for SNP calling. I have 300 samples. As far as I understand that if a sample is missing a SNP it should be denoted by ./., however, I am getting a lot of values like 0/...
0
votes
1
answer
69
views
syndip dataset for benchmark variant
I have a question about syndip dataset : https://github.com/lh3/CHM-eval .
I'm struggling to find the syndip vcf.
In the release ( https://github.com/lh3/CHM-eval/releases ), we have a file named : <...
- 25
1
vote
1
answer
172
views
VCF or BAM file for raw data of gene test?
My friend has a VUS (Polymicrogyria with or without VEDS) mutation that was found in her whole-exome sequencing with respect to the phenotype given that time. At present, her doctors have more ...
4
votes
0
answers
115
views
Error occurrence after merging files with bcftools: wrong number of fields?
This question was also asked on Biostars
I have multiple vcf of CASES and CONTROLS variations annotated by VEP, SNPEff and SnpSift.
First pair vcf -> only variations| CASES and CONTROLS
Second ...
- 141
6
votes
1
answer
692
views
Why does GATK produce both 0/1 and 1/0 genotypes in the same file? Are the two not equivalent?
I have always thought that 1/0 and 0/1 in VCF genotype fields are equivalent. And yet, GATK uses both. For example, these are ...
- 9,352
0
votes
1
answer
36
views
Is it valid VCF not to 'squash' positions with more than one ALT allele?
I'm seeing output from PBWT that looks like this:
...
- 440
1
vote
1
answer
514
views
Problems when using vcf2maf
I am trying to use a commonly recommended package named vcf2mafto convert vcf file to maf file and is reliant on Ensembl's Variant Effect Predictor (VEP) here. However, if I follow the guidelines ...
- 135
2
votes
2
answers
84
views
Confusing result from Sanger Imputation Service (Eagle v2.4 for phasing, PBWT v3.1 for imputation)
I ran my 23andMe data through the Sanger Imputation Service that uses Eagle v2.4 for phasing, PBWT v3.1 for imputation.
However, some aspects of the results are very confusing to me. Perhaps see #3 (...
- 440
0
votes
1
answer
172
views
Range of values for VAF in heterozygous genotype
(VAF = variant allele frequency)
This question follows from this one.
What is the expected range of values for VAF of mutations with heterozygous genotype (i.e., 0/1 or 0|1)?
Is it possible to have 0% ...
- 1,773
2
votes
0
answers
121
views
Assess ploidy estimation through SNV and CNV results
I want to assess ploidy estimation (i.e., diploid vs tetraploid) in a set of samples for which I have both results from the variant caller (i.e., SNVs) and from the copy number caller.
What are the ...
- 1,773
2
votes
2
answers
3k
views
How to convert a BAM file into a VCF file
I know that there is a lot of methods to do so, but most of them need to have a reference file. (For example, the GATK force me to provide a reference.)
However, I think the BAM file I am working with ...
- 135
0
votes
1
answer
287
views
Error in installing vcf2maf
I tried to install vcf2maf following this tutorial from GitHub: https://github.com/mskcc/vcf2maf.
I get a repeat error when I try to test out vcf2maf using the code:...
- 135
0
votes
0
answers
36
views
When can we have duplicated variants in a VCF file, excluding artifacts?
I read these two questions on duplicated VCF variants: 1, 2.
However, it is not entirely clear to me when duplicated variants can arise (excluding merging artifacts). I want to understand if ...
- 1,773
0
votes
1
answer
202
views
Imputing 23andMe 3v data to v4 and v5?
The 23andMe version 3 data (v3) is now quite 'old' and new data is being collected on the updated v5 Illumina SNP-Chip (details here).
Is there an online service for imputing 23andMe v3 files to v4 ...
- 440