Questions tagged [snp]
Single Nucleotide Polymorphisms. How to find them? How to correlate them with other data? How to represent them? What is their distribution?
164
questions
2
votes
1
answer
27
views
Does using a microarray chip that matches the chip used in training lead to higher explained variance of a polygenic score?
I'm trying to replicate an existing polygenic score (i.e. test the accuracy in a new sample), and want to know if matching the original study's microarray chip will improve the accuracy (that is, ...
- 513
1
vote
0
answers
58
views
SNP How to detect regions of erratic pairwise alignment?
I have clusters of DNA sequences that contain some mutations. All the sequences in a cluster should contain the same mutations.
As the mutations aren't known, those clusters of sequences were pairwise-...
- 181
1
vote
2
answers
63
views
Why does my HWE QQ plot have extreme deviation and what does it mean?
This question was also asked on Reddit
I have recently completed my thesis and one of the comments was that I report on why this graph looks this way. I have tried to find a reason but the closest I ...
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
0
votes
0
answers
183
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
1
answer
139
views
Getting VCF file that contain common SNPs from 6 VCF file using isec
I have 6 VCF files, where I would like to obtain the SNPs that are common (by position) in all the 6 files. I have tried this command
...
2
votes
0
answers
20
views
How do you go about grabbing proper genotyping infomation from dbSNP using entrez API?
I made an R function to grab rsID genotypes from dbSNP using the entrez API. For example, here is the xml output that I am trying to parse:
https://eutils.ncbi.nlm.nih.gov/entrez/eutils/efetch.fcgi?db=...
- 31
2
votes
2
answers
116
views
How can I detect the Apo B gene mutations in my raw dna data?
Familial defective apolipoprotein B-100 (FDB) is an autosomal dominant disorder associated with increased plasma cholesterol concentration and may thus increase the risk of premature coronary heart ...
- 121
0
votes
0
answers
52
views
Filter VCF file for most differentiated SNPs (Fst) in VCFTools
I have a filtered VCF file that I want to subset for the 500 most differentiated SNPs based on estimated Fst values (calculated in vcftools). I am not completely sure how to approach this. I have ...
- 31
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 ...
0
votes
4
answers
179
views
Can you infer orientation from vcf file data?
I have a dante labs dataset that I would like to cross reference with SNPedia but I am unclear about orientation. I've read the info on SNPedia and other places about it but I just don't get it yet. I ...
- 131
-1
votes
1
answer
52
views
Variants list to aminoacids
I have the sequencing data of a Monkeypox virus sample. I have the list of variants extracted from Nextclade and I want to find the corresponding aminoacids substitution.
For example I have the ...
0
votes
1
answer
59
views
Help with analyzing SNP data for three patient groups
I am currently working on a project where I have genetic data for patients classified into three different groups. Specifically, I have information on 26 SNPs for each patient. I have recorded each ...
0
votes
0
answers
52
views
SKAT method in R
I have data that looks like this: 3 column SNPs their gene based on Annovar and a p-value for every SNP. What I would like is to aggregate the p values for every gene.
...
- 29
1
vote
0
answers
42
views
Project PLINK eigenvector from one group on top of another group
I have 1 .bcf file called with individuals from two groups, one which is my samples, another which is a reference panel. Due to my samples having a limited amount ...
- 85
1
vote
1
answer
152
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
2
votes
1
answer
259
views
Getting nearest gene from SNP data using SNP ID, CHR, and BP
I have data from a GWAS that provides me with the SNP, Chromosome, and base-pair. My data set has thousands of SNPs. What is the easiest way to find the nearest gene for each SNP using this ...
- 21
2
votes
1
answer
123
views
Arlequin - How to assign samples to populations?
I have converted a VCF file to the arlequin file format. It consists of 20 samples that should belong to three populations, along with a few thousands SNPs. However, when I load the arlequin file into ...
- 21
1
vote
1
answer
520
views
Retrieve rsIDs from chromosome positions using MRutils::get_rsid_from_position
I have 5609 chromosome positions with both ref and alt alleles. I want to retrieve rsIDs for the same. I tried get_rsid_from_position() of the ...
- 11
2
votes
2
answers
53
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
133
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....
0
votes
1
answer
44
views
SNPs location out side the gene
I downloaded exome data from gnomad to get the allele frequency of some SNPs that I'm analyzing. The SNPs are from external resource. Some SNPs (in my data) are not present in the exome data (in ...
2
votes
1
answer
687
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 ...
- 923
4
votes
1
answer
109
views
Find SNPs in yeast genomes
I'm a new Bioinformatic scientist working for a yeast genetics company.
Objective To create a database of yeast genomes from NCBI and identify SNP variants.
In my pipeline
FastQC,
Trimmomatic,
BWA
...
- 923
2
votes
0
answers
102
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
112
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 ...
2
votes
1
answer
55
views
What tool I can use to map short-reads sequences to reference genome and get specific mapped size
I have about 300 90-bp sequences which I would like to map to a reference genome to make it longer to a 300bp sequence, wherein the 90bp is at the middle.
Anyone knows what bioinformatics tool I can ...
- 21
1
vote
0
answers
48
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
2
answers
48
views
BAM files with no RNAME and POS, how to map contents to SNPs?
I have a set of 4 .bam files containing the exome of an individual, around 400 MB each. I used samtools to generate a 2.4 GB .sam file out of one of the .bam files, and I found it contains lines with ...
3
votes
1
answer
78
views
Imputing small region of the genome
If I'm looking for a specific SNP in my SNP-Chip data and it isn't there, are there any tools that let me quickly impute that SNP from surrounding SNPs rather than running a lengthy 'whole chromosome' ...
- 470
1
vote
1
answer
65
views
Reference variant detected as altered one in bam file
I received (from manufacturer) several .bam files and I used four callers (samtools, freebayes, haplotypecaller, deepvariant) to find some sequence variants. In obtained .vcf files, I took a closer ...
- 206
0
votes
1
answer
64
views
Regression plot of a continuous trait - is there a binary equivalent?
Within the GWAS wikipedia page, you can view the following plot:
"Illustration of a simulated genotype by phenotype regression for a
single SNP. Each dot represents an individual. A GWAS of a ...
- 155
1
vote
0
answers
37
views
RNAseq SNP discovery: deciding upon filters and dealing with allele expression bias
I am working with non-model plant RNA samples which we have been deep sequenced and analysed using STAR aligner under default parameters.
Aim We would like to conduct SNP discovery of these samples.
...
- 11
2
votes
2
answers
109
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 (...
- 470
2
votes
0
answers
173
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,783
3
votes
1
answer
38
views
SNPs in LD in which populations?
I have 3 SNPs that are in LD in 'several' populations here:
https://www.ebi.ac.uk/gwas/variants/rs13078854#ld_panel_label
It's a pain to look at all (30) populations one by one to manually check if ...
- 470
1
vote
1
answer
862
views
Filtering a vcf with a text file of SNP rsIDs
I have a vcf file containing an ID field with a snps and a list of human SNPs from dbSNP.I want to remove whatever in dbSNP from vcf file. The text file is like this:
...
- 1,981
-3
votes
1
answer
581
views
Extracting a column from a vcf file
I have a zipped vcf file of dbSNP hg38 version
No space left on my device to unzip that
I want to extract a column from that
I have tried this
...
- 1,981
1
vote
1
answer
502
views
imputed dosage values for vcf files
I have one more question related to genotype file. I submitted the job in michigan imputation server.
and my imputed dosage for chr 22 file is:
22 16050435 22:16050435:T:C T C . PASS
AF=0.00098;MAF=0....
- 15
0
votes
2
answers
130
views
VCF file inspector app with GUI
Does anyone know a good, free-to-use VCF file inspector app (Mac OS, with GUI) where I can browse the content of a VCF file similarly to bamseek*?
I would like to browse SNP annotation split by ...
- 123
1
vote
0
answers
41
views
How to split genotypes into two groups from a SNP txt file contains -1, 0 and 1 in R?
I have a SNP file(it is 96 obs. of 1178 variables with -1, 1 and 0 and in some places NA) and I need to calculate Jaccard genetic distance from them in R. Is there any way to convert these 3 values to ...
0
votes
2
answers
145
views
HG19 Position meaning
I am looking at genomic data (HG19). I have many SNPs, their chromosomes and positions.
I want to look at certain SNP (suppose its chromosome is 1 and position 77,226,919), and extract all SNPs in ...
- 231
1
vote
1
answer
64
views
Given SNP-Chip data for a population, what tool should I use to reconstruct haplotypes?
I think Eagle or haploSep will do the job of reconstructing haplotype data, but I'm not sure what the 'standard' or best practice tool to use is.
I have a VCF/PLINK format file of ~8,000 individuals ...
- 470
0
votes
1
answer
40
views
What to use to determine SNPs in a population? How to find significant SNP variation between generations?
I have a population in the hundreds and did an Illumina pool-seq (ended up with ~100bp reads for the whole population).
I have the reference genome for the species and now need to align all the ...
- 249
0
votes
0
answers
63
views
Combine all alternate alleles
I have a vcf with multi-allelic SNPs, e.g. REF: A and ALT: G, C.
I would like to combine all reference alleles to the most ...
- 101
1
vote
0
answers
255
views
Extract read names and the associated nucleotides on specific positions from a BAM file (in R)
Let's assume I have a BAM file and several positions that I would like to examine more closely in this alignment. My goal is to find out whether these positions are on the same reads and which ...
- 11
0
votes
1
answer
84
views
Snakemake rule issue with bcftools view filtering
I am creating a bioinformatics pipeline using Snakemake doing everying from adapter trimming of raw sequencing data in .fq to mapping, SNP calling and SNP filtering.
The issue i am having is when i ...
- 85
1
vote
1
answer
30
views
A unified database for CNV, SNP, Indel and MSI
I am looking for a database or different databases where I can find information on different gene variants for a gene. As an example if I enter PPARG, I could be able to see SNPs, CNV, InDels and MSI.
...
- 227
1
vote
1
answer
54
views
Output of allelic association doesn't write the rsID [closed]
I'm using plink (1.9b5) to do allelic association. My problem is that my output does not write the rsID (SNP), it just writes a dot.
Output:
...
- 11
1
vote
2
answers
69
views
Is there a way to calculate the Ka/Ks value in R/unix?
I am investigating the rate of synonymous non synonymous substitutions and would like to find the Ka/Ks value also known as dn and ds between multiple alleles? I want to investigate multiple alleles ...
- 95