Questions tagged [variation]
Refers to genetic variation such as SNPs, indels, and structural variants within a population or between individuals.
46
questions
1
vote
0
answers
80
views
Eigen phred score interpretation
I need help. I am not able to find any reference for Eigen_phred_coding prediction tools score from dbNSFP.
I know that it may be scaled in a similar manner to CADD_phred - ranking raw values from ...
2
votes
0
answers
24
views
Feasible to find genetic variations of two samples using RNAseq data?
I have bulk RNAseq data from two strains of mice from Jackson Lab: C57BL/6 and B6.SJL.
The former expresses a Ptprc-b allele and the latter expresses a ...
- 237
3
votes
3
answers
129
views
When it is stated in your consumer genetics testing report that you have x% (lets say 1.3 percent...) Neanderthal DNA, what is actually meant?
This question has also been asked on Biology SE
When a raw number like 1.3 percent is stated, it is confusing without context,and can be misinterpreted. So does 1.3 percent mean that 1.3 percent of ...
1
vote
0
answers
57
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
2
answers
142
views
Multiple Sequence Alignment
I applied three multiple sequence alignment methods of muscle, mafft and clistal omega on 30 sequences and 400 sequences . In the case of 30 seq, the output of alignment from the 3 methods was in the ...
0
votes
1
answer
39
views
Getting all variants in a region from ensembl
I want to get all the variants (SNPs) in a given region of the human genome. Is there a way on how to do this with the ensembl-API? It should be a representation similar to this:
There I added dbSNP ...
3
votes
1
answer
491
views
How to estimate the phenotypic variation explained by top SNPs from a GWAS study?
I have conducted a large-scale GWAS study and got a few significantly associated SNPs. I used GEMMA with -lmm 1 options to run ...
- 75
5
votes
0
answers
51
views
Parse RNA variant effect annotations ("r." format)
I've got annotations for splicing variants in a format like this (this is one variant):
Variant: NM_004092.3:c.88+5G>A
Effect: Retention; r.87_88ins1_88+10:p.(Ala31Glufs*23)
I want to extract ...
- 51
2
votes
2
answers
447
views
Produce .vcf file of ALL mutations in .bam file
I got tasked with a problem, that i thought would be quite simple to solve, but turned out to be quite tricky. Our lab is running targeted mutagenesis experiments in yeast using crispr base editors. ...
1
vote
1
answer
64
views
BWA: Detecting Variation between Reference Genome and Short-Read Sequences
I need to identify all loci in the short-read sequence at which the number of microsatellite repeats (i.e. number of copies of "AA," "GTC," etc.) differ from the reference genome, ...
- 189
3
votes
2
answers
2k
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 ...
- 257
2
votes
1
answer
387
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?
...
- 29
3
votes
3
answers
5k
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 ...
- 465
4
votes
1
answer
193
views
Remove variants that do not map to human genome
[This question was also asked on Biostars]
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 ...
- 151
3
votes
1
answer
60
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 ...
- 31
1
vote
4
answers
1k
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 ...
- 5,060
0
votes
2
answers
43
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 ...
- 101
0
votes
2
answers
1k
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 ...
- 101
2
votes
1
answer
40
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)?
- 199
5
votes
2
answers
3k
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 ...
- 9,662
4
votes
1
answer
604
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:
...
- 9,662
3
votes
2
answers
275
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 ...
- 2,324
3
votes
1
answer
482
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.
...
- 255
6
votes
1
answer
7k
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 ...
- 199
3
votes
1
answer
105
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 ...
- 333
3
votes
1
answer
145
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 ...
- 105
2
votes
0
answers
48
views
Mutation detection using Varscan2 on RNA sequencing for estimating tumour clones with pyclone or other package
I would like to analyze my RNAseq profiles from bulk tissue samples (Paired-End, 50M reads/sample, tumour-normal pairs) with varscan2 to detect mutations.
Then I ...
- 21
5
votes
1
answer
858
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 ...
- 2,739
4
votes
2
answers
428
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 ...
- 560
4
votes
2
answers
840
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 ...
- 2,459
7
votes
1
answer
2k
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 ...
- 560
9
votes
2
answers
12k
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 ...
- 361
2
votes
1
answer
241
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, ...
- 560
2
votes
1
answer
524
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:
...
- 2,459
1
vote
1
answer
220
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:
...
- 967
2
votes
1
answer
43
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 ...
- 2,324
7
votes
1
answer
2k
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 (...
- 967
8
votes
4
answers
1k
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 ...
- 1,681
14
votes
7
answers
3k
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 ...
- 243
7
votes
1
answer
143
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....
- 71
3
votes
2
answers
159
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 ...
- 1,022
6
votes
2
answers
432
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 ...
- 1,681
13
votes
2
answers
383
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 ...
- 232
13
votes
2
answers
2k
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):
...
- 2,739
20
votes
1
answer
465
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 ...
- 788
13
votes
3
answers
366
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 ...
- 1,681