Episode #125 of the Stack Overflow podcast is here. We talk Tilde Club and mechanical keyboards. Listen now

Questions tagged [variation]

Refers to genetic variation such as SNPs, indels, and structural variants within a population or between individuals.

Filter by
Sorted by
Tagged with
1
vote
1answer
20 views

Find rsIDs for GRCh37 SNPs and rsIDs for GRCh38 SNPs and compare for overlap

I have a long list of variant IDs that I generated as a result of research done one GRCh37 genomes (e.g. 13_28025615_G_C_b37). I want to get their rsIDs and compare ...
1
vote
0answers
48 views

How to use -csvStats option?

I'm trying to get CSV files from snpEff -csvStats option and got stuck. I tried the code below and only got the vcf file with no contents. How can I use -csvStats option properly? ...
2
votes
3answers
224 views

Find overlap between VCF files

I have two VCF files and I want to compare the missing rate in each of them, but I want to only look at sites that are present in both files. How could I go about getting a list of positions that are ...
3
votes
1answer
78 views

Remove variants that do not map to human genome

I received an hg38 VCF file that's had variants imputed with 1000 genomes. I've encountered some issues with the VCF; REF alleles that do not align to a reference genome, ALT alleles that do not ...
3
votes
1answer
24 views

Database for germline copy number variations in cancer

I am interested in looking at germline copy number variations in individuals that are at high risk of developing cancer. Are there any databases where I can look, if a CNV seen in our test case has ...
1
vote
3answers
136 views

Index VCF file for rapid access by rsID?

With tabix one can index a VCF file for rapid variant retrieval based on genomic position. I'm wondering if there are any tools that will index a VCF file to allow rapid retrieval using rsIDs and/or ...
0
votes
2answers
40 views

Evidence of emergence of “genuine” novel protein interactions by aa mutation [closed]

Anyone knows of an example of a protein that, without coming from a recent duplication event, underwent a mutation(s) that caused it to have a novel interaction with a new ligand, substrate, other ...
0
votes
2answers
66 views

How to compute LD for pairs of variants with Plink

I have a set of SNPs (single nucleotide polymorphisms) { S1, S2, ..., SN } from approximately 200 humans. I wish to determine the linkage disequilibrium (LD) for a defined subset of these SNPs. By ...
2
votes
1answer
28 views

How to interpret amino acid representation

What is the difference between C383S and C383A? For C383S, does it mean Cysteine(C) position 383 in the sequence and its next amino acid connect to it is Serine(S)?
5
votes
2answers
389 views

How can I remove (non-trivial) duplicates from a VCF file?

This is related to the question I asked here. Consider a vcf file that contains duplicate variants, but where the duplicates aren't simply the same thing in the same notation but instead one is a ...
4
votes
1answer
109 views

Are duplicate variants against the VCF standard?

I have been trying to understand an error that the EBI's vcf_validator gives when run on my vcf file. Consider this minimal file: ...
3
votes
2answers
206 views

Selecting 65000 SNPs where AF is close to 0.5 in all or most populations

I am evaluating the tool somalier (https://github.com/brentp/somalier) and I need to create a list of about 65,000 SNPs where the allele frequency (AF) is as close to 0.5 as possible across the most ...
3
votes
1answer
100 views

Output from vcftools missingness

I'm new to data filtering on vcf data and vcftools. I performed variant calling on my dataset, CHR22, homo sapiens. I'd like to remove sites that are missing in more than 5% of individuals. ...
5
votes
1answer
141 views

SNP vs Point Mutation

What is the difference between a Single Nucleotide Polymorphism (SNP) and a point mutation? I am quite confused in understanding these term as both of them refer to one base difference from the ...
3
votes
1answer
47 views

multiple COSMIC id for the same mutation

I would like to know why there are different entries for the same mutation in COSMIC, for example COSM6954941 and COSM12833 both refer to ERBB4 c.908C>A. In this specific case the field Gene name in ...
3
votes
1answer
63 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 ...
2
votes
0answers
27 views

Mutation detection using Varscan2 on RNA sequencing for estimating tumour clones with pyclone or other package

I would like to analyze my RNA seq profiles from bulk tissue samples (Paired-End, 50M reads/sample, tumour-normal pairs) with varscan2 to detect mutations. Then I would like to use those detected ...
4
votes
1answer
273 views

Simplest way to work out structural variant type?

In VCF 4.2, a structural variant (SV) can be described with the BND keyword in SVTYPE. For example, the following example is an ...
3
votes
2answers
193 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 ...
4
votes
2answers
322 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 ...
6
votes
1answer
575 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 ...
7
votes
2answers
2k 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 ...
2
votes
1answer
105 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, ...
2
votes
1answer
162 views

Can I export the SNPs between whole genome alignments using the command line?

I am using Mauve to align two whole genomes. I use the following command to get the alignment in xmfa format: ...
1
vote
1answer
105 views

Write the output of the VariantAnnotation package to file

I'm having some problems in exporting the annotations of some variants from a VCF file, after using the "VariantAnnotation" package from bioconductor. Below the script I've used: ...
2
votes
1answer
38 views

remaining human genome variation that hasn't been sequenced?

Given all the genome variation information we have accumulated from resequencing human genomes up until now (early 2018), how much variation is left that hasn't been sequenced? EDIT: given one of the ...
5
votes
1answer
519 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 (...
8
votes
4answers
729 views

How to manipulate a reference FASTA or bam to include variants from a VCF?

I have some software which takes fastas as the input. I need to include SNVs and InDels from a VCF into the reference hg38 and then use this. The problem is, I don't know of an algorithmically sound ...
12
votes
3answers
976 views

Is there public RESTful api for Gnomad?

I currently find Harvard's RESTful API for ExAC extremely useful and I was hoping that a similar resource is available for Gnomad? Does anyone know of a public access API for Gnomad or possibly any ...
7
votes
1answer
135 views

What is a good pipeline for using public domain exomes as controls?

I'm currently attempting association analysis with an extremely small set of patient exomes (n=10), with no control or parental exomes available. Downloading the ExAC VCF of variant sites (http://exac....
3
votes
2answers
130 views

How can I find mutations associated with disease in human histone residues?

I would like to look if there are mutations in residues of human histones associated with any disease. For instance, if a mutation in residue K6 (lysine 6) of histone H2A1A is associated with any ...
6
votes
2answers
275 views

What is the distribution of indel sizes in a healthy human genome? of insertion:deletion ratios?

My understanding is that indels are from 1bp to 10Kb, and a healthy genome has ~400K-500K Indels. Surely most of these are small. What is the distribution of insertion sizes in a healthy human ...
12
votes
2answers
148 views

Remapping genomic coordinates to account for indels

I'm interested in obtaining coding sequences of my favourite gene in all individuals from the 1000Genomes (and similar projects). I use GATK to get the right subset of variants, vcf-consensus to map ...
12
votes
2answers
733 views

How to read structural variant VCF?

The IGSR has a sample for encoding structural variants in the VCF 4.0 format. An example from the site (the first record): ...
16
votes
1answer
323 views

The state, limitations and comparisons of large variant stores

Background: We're increasingly needing some way of storing lots of variant data associated with lots of subjects: think clinical trials and hospital patients, looking for disease-causing or relevant ...
12
votes
3answers
277 views

Given a VCF of a human genome, how do I assess the quality against known SNVs?

I'm looking for tools to check the quality of a VCF I have of a human genome. I would like to check the VCF against publicly known variants across other human genomes, e.g. how many SNPs are already ...