Stack Exchange Network

Stack Exchange network consists of 174 Q&A communities including Stack Overflow, the largest, most trusted online community for developers to learn, share their knowledge, and build their careers.

Visit Stack Exchange

Questions tagged [vcf]

VCF, acronym for Variant Calling Format, is a text file type used to store information regarding genetic variants.

1
vote
0answers
12 views

R package grur: missing_visualization() fails during redundancy analysis of .vcf

I am trying to use the R package grur to process GBS SNP data, which is installed from github as instructed. For starters, I am running the test code in this vignette. As suggested I am running it ...
3
votes
1answer
43 views

Script Giving error at paste command

I wrote script for pasting rsids on CADD output. Here is script. ...
4
votes
3answers
41 views

bcftools filtering all files in a directory

Probably a silly oversight on my part, but I'm trying to filter all the vcfs in a directory with bcftools using a simple loop. My basic command is working fine: ...
2
votes
4answers
30 views

how to filter a multi-VCF for sex-specific genotypes?

I have a single VCF file of variants for a population that includes males & females and I would like to pull out all variants that are sex-specific. For example, sites where all males are 0/1 and ...
5
votes
1answer
56 views

calling diploid SNVs from long reads

I'd like to call diploid SNV variants from long-reads data (~80SMRTcells PacBio). I have generated a draft reference genome for an indivudual from a heterozygous (~4%) species (Canu+Haplomerger2). ...
1
vote
1answer
23 views

How to paste RSIDs in CADD output

I want to paste RSIDs in CADD output as CADD does not give RSID column in its output. For this purpose I am using bedtools intersect to compare two files and have RSID column in my CADD file. This is ...
0
votes
0answers
14 views

How to use FastQTL with LeafCutter?

I am following along with the example data in the LeafCutter documentation and I am a bit stuck on the final step of using FastQTL to generate QTL maps. Documentation for FastQTL is scant and so are ...
1
vote
0answers
72 views

Combine VCF files

I have a problem with using rbind to combine VCF files using the library VariantAnnotation from Bioconductor. I am reading two VCF files, when I try to combine them in a certain order with 'rbind' I'...
4
votes
3answers
80 views

Sampling haplotypes

I am trying to simulate different genome of peoples, I have data (VCF files) of various genes from the 1000K Gene project. I want to simulate different whole genomes i.e generate a new population by ...
8
votes
1answer
55 views

Best practices for deciding if two structural variants are actually the same variant?

I know that I can use bcftools isec to compare two VCF files containing single nucleotide variants (SNVs); it will generate four possible output files: one with all ...
2
votes
1answer
62 views

Simplest way to work out structural variant type?

In VCF 4.2, a strucutral variant can be described with the BND keyword in SVTYPE. For example, the following example is an ...
4
votes
2answers
168 views

No variant found using GATK 4.0 HaplotypeCaller

I am doing variant calling on RNA-seq datasets from wheat which is hexaploid,the binary alignment (BAM) files were created using STAR version 2.6.0c and variant calling was done using GATK 4.0 ...
1
vote
1answer
37 views

High Phred Quality score VCF file

I am working on ION Torrent S5 XL and I am facing to an issue with VCF output generated by the variant caller of thermofisher. It stating on VCF 4.1 that Phred Quality Score is on a range of 0 to 255. ...
3
votes
1answer
45 views

How to convert files to ADAM format?

I would like to convert BAM and VCF files to ADAM format. How do I do that?
4
votes
1answer
51 views

Meaning of the FORMAT fields of the VCF file coming from GIAB project

After reading the GIAB paper in https://www.biorxiv.org/content/early/2018/05/25/281006 and its Figure 1, I am still having trouble understanding the data inside the GIAB VCF file for HG001 (...
1
vote
0answers
29 views

Normalize by counting CNV within a pool of VCF files

I’m (trying) using the GATK4 germline CNV calling pipeline. I successfully got 57 VCFs from my sample batch, called with segments (obtained by merging the contiguous intervals), like in a classic VCF :...
2
votes
2answers
107 views

How to filter for “CAF=0.5747,0.4253”

In my gzipped VCF file I have many variants, and they report allele frequency such as: ...
2
votes
0answers
41 views

bcftools consensus not working when using -H flag

Sorry, I have been trying to work this out desparately all day. Biostars couldn't find any answer so I am cross posting here https://www.biostars.org/p/315458/#315686 . I am trying to get both the ...
2
votes
2answers
119 views

Protein sequence from patient data

Currently, I am working on NGS data and my aim is to get significance prediction of variants present in the vcf file. As we know about SIFT Score for significance score prediction, I am trying to ...
1
vote
1answer
126 views

update dbSNP ID

I have used dbSNP build 138 to tag 'rs' ID in my VCF file. But later I came to know that now dbSNP build 150 has been released. I want to update the rs IDs in my VCF file. I tried but I found a ...
1
vote
1answer
34 views

How to add OVERLAP_BP field to VEP output?

I'm using the Ensembl Variant Effect Predictor (VEP) to annotate VCF files, and I am trying to add some additional fields to the output, among these is OVERLAP_BP, ...
1
vote
0answers
72 views

Error filtering SNPs with Plink

Im trying to filter SNPs using PLINK with this command. Also, the .map and .ped files were converted to a plink format from a vcf file (...
3
votes
2answers
152 views

GATK CombineVariants complains the contig order in the VCF files

I have called variants on two strains of C. elegans separately. I now want to merge the VCF files into one using the following code: Create a sequence dictionary of the reference sequence Sort the ...
4
votes
1answer
74 views

Keep Format and Individual fields when annotating VCF with VEP

I'm currently updating my Variant Calling Pipeline by switching the VCF annotating software from Annovar to VEP for a variety of regions, not least how easy it is to annotate with HGVS notation and ...
4
votes
1answer
195 views

Merge 2 VCFs from different variant callers

I am working with WES data for detection of somatic variants and I have used two variant callers because no variant caller is complete in itself. I have used GATK Haplotypecaller for small variants ...
0
votes
0answers
42 views

Complete Genomics Junction file for SV Detection

At present, I have vcf file for each patient generated from tumor data. I am using virtual normal correction method for somatic mutation detection because I don't have normal data. For SNP or small ...
7
votes
2answers
613 views

Converting a BAM file into VCF

I have NGS illumina RNA-seq reads from M. musculus (mm10). I am trying to find variants along the strand portion of the reads in the refseq (mm10). I mapped a pair of sequence files and generated a ...
1
vote
1answer
70 views

How to get unique somatic mutations for each individual patients

I have exome data 32 patients. My aim is to identify unique somatic mutations in each patients. For this I have completed my pipeline for Whole Exome Data analysis using various tools like BWA mem, ...
4
votes
1answer
2k views

Converting VCF file to PLINK bed/bim/fam files

I am trying to find the best way to convert VCF files to PLINK binary bed/bim/fam files, but it seems like there are many varied ways to do this. (For example, using Plink 1.9 --vcf tag, bcftools, ...
4
votes
1answer
100 views

Can a data file in VCF format be converted into FASTA?

I'm considering purchasing the 'MyGenome' product by Veritas Genetics to analyze my genome for a project. I'd like the data to be in FASTA format, but Veritas only provides VCF data. Is it possible to ...
2
votes
0answers
43 views

Post ipyrad filtering of SNP loci from GBS, based on SNP quality

I want to use SNPs produced by ipyrad analysis of GBS data in RAxML to examine the monophyly of a highly-variable focal species and its phylogenetic relationships with about 40 congeners. An initial ...
2
votes
2answers
881 views

How to subset samples from a VCF file?

I have VCF files (SNPs & indels) for WGS on 100 samples, but I want to only use a specific subset of 10 of the samples. Is there a relatively easy way to pull out only the 10 samples, while still ...
3
votes
1answer
555 views

Tools to do VCF to MAF and MAF to VCF conversion?

Normally, I would use the vcf2maf scripts to convert a VCF to a MAF (or vice versa). This is great software, but on my system, perl scripts with dependencies are easy to break. (Here it uses VEP.) ...
3
votes
2answers
1k views

How to subset a VCF by chromosome and keep the header?

I would like to subset a VCF which only has chromosome 2. The problem with using various grep commands, e.g. ...
2
votes
3answers
1k views

How to concatenate “by chromosome”-VCFs?

I have a several VCFs which are VCF which only contain information by chromosome. That is, there's a chromosome 1 VCF (with only chr1), a chromosome 2 VCF (with only chr2), etc. I checked to make ...
6
votes
1answer
128 views

ethnicity check either from bam or vcf files

What tool could I use to check the ethnicity of a human bam or vcf file? I would like to use the results as a QC check to know whether a given sample or set of samples match the ethnicity information ...
1
vote
2answers
122 views

Approximate matches with bedtools intersect

I am using bedtools to test whether a set of variant calls produced by a new method is missing any variants from a set of known validated variants. ...
4
votes
1answer
46 views

Intersecting two different files with one “master” file based on different columns

I have the following sets of data: file1: 1 15776220 15776240 GTGACCAGCAGGTGTCTCTG 16855676 16855696 CTGTCCAGCAGAGGGCGGTG file2: ...
5
votes
2answers
102 views

Using column 2 of one file to match with two columns of another file, and append

I have file 1 like following: 1 15776220 15776240 GTGACCAGCAGGTGTCTCTG 16855676 16855696 CTGTCCAGCAGAGGGCGGTG And file 2 as following ...
5
votes
3answers
167 views

How to indicate the END of a haplotype block in VCF?

In VCF I know how to indicate that two genotypes are in the same phase by using consecutive "0|1" and "1|0" genotype fields, for example. However, how do I deal with the case that the first two ...
2
votes
1answer
466 views

Using Python, how to convert a pandas DataFrame into a VCF?

Let's say I have a pandas dataframe with fields CHROM, POS, ALT, ...
3
votes
1answer
102 views

Where can I get the population allele frequency vcf file?

I want to use GATK to Estimate cross-sample contamination for Whole Genome Sequencing data. The specific tool is ContEst and it is run with: ...
4
votes
1answer
59 views

Is it possible to check if a patient has the HLA-B27 antigen from his exome stored in a VCF file

Is it possible to check if patient has the HLAB27 antigen by comparing his exome stored in VCF with all known HLAB27 variants? For instance, rs4349859 is a known HLAB27 SNP. If I can find rs4349859 ...
2
votes
2answers
225 views

bcftools output vs. bgzip

I've got what seems like it should be a simple question, but I can't seem to figure it out from google. bcftools has 4 output options: ...
6
votes
2answers
2k views

Converting a VCF into a FASTA given a reference with Python, R

I am interested in converting a VCF file into a FASTA file given a reference sequence with Python or R. Samtools/BCFtools (Heng Li) provides a Perl script ...
5
votes
3answers
144 views

Any fast options to query large VCF bed intervals?

I'm doing some analysis and I need to subset a large VCF file (~8GB gziped) given a bed interval and identify within a subset of rsid. Unfortunately, both my normal choices to do this analysis (<...
6
votes
1answer
162 views

VCF merge containing CNV

How do I merge VCFs files containing CNVs? I use vcf-merge, a VCFtools function, and after bgzip and tabix, SAMtools, to index and tab separate variants, but I don't know if it is the right way.
5
votes
1answer
264 views

How to represent a deletion at position 1 in a VCF file?

I am writing a small script to write a VCF file from short multiple sequence alignments, and stumbled into an interesting corner case. This is the way a deletion should be represented in VCF format (...
4
votes
2answers
107 views

Is there a way to quickly verify the presence of some SNPs in Fastq files?

Is there a tool that can scan fastq files without assembling them for a custom list of user defined snps?
7
votes
3answers
360 views

Tools to create annotated table of variants from VCF

The problem: I have a VCF file, a reference genome, and a bunch of annotations for the reference (genes, repeat regions, etc.) as GFF or BED files. What I would like is a tool that takes all of this ...