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?

Filter by
Sorted by
Tagged with
1
vote
0answers
24 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 ...
0
votes
1answer
23 views

Biological datasets for polythetic classification

Question: What are some instances of polythetic datasets in biology? In particular, I am looking for a dataset to benchmark a machine learning algorithm optimized via episodic training. On polythetic ...
0
votes
0answers
24 views

How to make CIRCOS plot of VCF file?

I have 3 VCF files, I want to make CIRCOS plot just like the following: https://www.nature.com/articles/s41598-020-70527-8/figures/4 I try to follow the tutorial but could not understand. The genome ...
-1
votes
1answer
60 views

Find common SNPs in multiple VCF files

I have 3 VCF files. A1.vcf A2.vcf A3.vcf I want to get the common SNPs that are present in all these three files. And output must be in vcf format. output: ...
0
votes
0answers
26 views

How to plot a graph of Linkage disequilibrium with limited information

I have a snp information in an excel file and I used SNPstat to do some analysis. But I also need a plot of LD for these snps. But my file is not formatted as the most programs required (.ped, hapmap, ...
-1
votes
1answer
24 views

How to know the number of SNP in my samples after running the population module in stacks

I am new to stacks for analyzing RAD-Seq data, I have run all the modules. I will like to know to get the total SNP from my sample after running the population module. Is there a command line to use ...
0
votes
1answer
28 views

Annotating gene names or gene IDs to a dataframe containing SNPs?

I have a large data-frame (excel file) of SNPs with genotyping data. I need to filter the data to get SNP information of a specific gene alone. The list is too long to do it manually. I was wondering ...
0
votes
1answer
39 views

Converting FASTA Sequences to SNP FASTA Format

I have two FASTA sequences. Each is from a different species' mitochondrial genome. The sequences are pasted below. ...
1
vote
1answer
31 views

Question about public availability of human SNP dataset with country of origin

Is there a publicly available data set for humans SNPs, preferably together with the country of origin? If so, could someone please point me towards it. I've looked through other questions on the site ...
0
votes
1answer
16 views

Denovo, Stacks: Getting an “ambiguous redirect” error

I tried running the following command ...
0
votes
1answer
64 views

Supplementary aligments in VAF

This question has also been asked on Biostars I have a doubt, are supplementary alignment usually considered when the variant allele frequency is calculated? Thanks a lot. In my case I have some ...
0
votes
0answers
8 views

Can we use reference strain as an out-group strain to root SNP based species tree?

I have generated a species tree based on whole-genome SNP using CSI_SNP_Phylogeny, wherein I used Reference strain as an outgroup to root my species tree (SNP fasta were downloaded and RaxML used to ...
0
votes
0answers
21 views

FMT/DP of missing genotypes in bcftools

I have a VCF file where the FMT/DP of missing genotypes has been recorded itself as missing data (".") When I filter ...
0
votes
0answers
75 views

Filtering samples by FORMAT attributes using bcftools

The -e and -i options of the bcftools filter command appear, by default, to only allow for ...
0
votes
0answers
20 views

How to assess how prone genes may be to acquire structural polymorphisms?

I have 5 strains of P.falciparum. Each FASTA file has all its annotated CDSs. After a first pre-processing phase, where I eliminated the strangest sequences (perhaps the longest or shorter ones, which ...
0
votes
1answer
16 views

How do I infer someone's haplotype of a gene from data on several of the gene's SNPs?

For a a few SNPs from given gene, I know there happen to be several haplotypes that most of the population has. But from a person's data on those SNPs, we have two alleles at each SNP, so it's not ...
1
vote
1answer
27 views

freebayes segmentation fault

I have a draft, Illumina-based, genome and bunch of reseq samples. I generated bam files using bowtie2 and marked duplicates using ...
1
vote
1answer
52 views

Biology behind PCA analysis based on SNP

What is the biology behind SNP based PCA analysis to study population structure? I am reading some articles where they used PCA analysis to compare isolates of drosophila that is collected from a ...
0
votes
0answers
20 views

How to control/normalize for number of reads when calling SNPs using RNA-Seq?

I used the GATK pipeline to call SNPs on males and females using RNA-Seq data. But the males have a higher read count (~43-46M reads) than the females (~40-42M reads). This causes SNP counts to be ...
0
votes
0answers
14 views

What number to give for `-ploidy` for durum genome in GATK?

I have some genome sequencing data from durum wheat. I'm using GATK to call variants. I got the genome sequence here: http://plants.ensembl.org/Triticum_turgidum/Info/Index Here is the reference ...
1
vote
1answer
39 views

Random GWAS data generator

I was wondering if there is a tool/script/program that randomly generates GWAS data. The purpose of such a tool would be to use it for educational purposes. So you generate some random .ped, .map and ...
1
vote
2answers
40 views

RAD Seq Data Analysis without barcode

I received a ddRAD Seq data from my supervisor without barcodes and restriction enzymes. I asked him for both and he said I don't need it since the data has been cleaned by the company. Now, I am ...
0
votes
1answer
26 views

Loss of predictive power of polygenic risk score when dataset contains missing variants

I am trying to calculate polygenic risk scores (PRS) scores for a new dataset. This dataset does not have all the variants that the PRS score needs. The PRS score I am interested in has 40 variants, ...
1
vote
1answer
42 views

Apply trained PRS on another dataset

I am using PRSice to compute the PRS over a train set and want to use the coefficient used on the train set to apply it on another set which I will call the test set. Once I compute the PRS I get a ...
1
vote
1answer
29 views

dbnSNP frequency anomalies

Sometimes dbSNP reports very different allele frequencies for different large-scale genome projects e.g. between 1000 Genomes and GnomAD rs11822440 1000Genomes A=0.4629 C=0.5371 GnomAD A=0.99997 ...
1
vote
0answers
27 views

GCF VS GCA or combination of both for pangenomic studies

I have just started to learn bioinformatics and pangenomics. So if this question seems to you pretty basic then I apologize in advance. As we know, NCBI for the genome database, there are two kinds of ...
0
votes
1answer
35 views

Inheritance annotation for nsSNPs

I am relatively new to SNP analysis. Is there a database to find annotation about the Mendelian inheritance of SNPs? I have a small list of nsSNPs (non-synonymous SNPs) and I need to find how they are ...
0
votes
1answer
50 views

Odds ratio and enrichment of SNPs in gene regions?

I did a QTL analysis with a panel of 7M SNPs, and want to analyze the enrichment of the significant qtl-SNPs in different genic regions (promoters, gene bodies, TFBS, etc.). A straightforward way to ...
2
votes
1answer
30 views

Core genome too short

I'm trying to create a snp distance matrix from M tuberculosis isolates to try and infer transmission networks. My plan was to use snippy to make a core genome ...
2
votes
2answers
56 views

Where can I find the genotype frequencies for the SNP's that are tested at 23andMe, Ancestry and FTDNA/MyHeritage?

I want to add the frequencies of the genotypes to a baseline of raw DNA data I've got from the DNA testing companies 23andMe, Ancestry and MyHeritage (tested at all). All in build 37 human genome. ...
1
vote
1answer
162 views

Convert VCF to genotype table

How can I convert a VCF file into a genotype table (SNP matrix)? I have this format: ...
0
votes
1answer
20 views

e-utils api: get chromosome position for each snp on a gene

So I have the following e-utils api-link: https://eutils.ncbi.nlm.nih.gov/entrez/eutils/elink.fcgi?dbfrom=gene&db=snp&id=5726&retmode=json it returns all the SNP-ids for the gene with ...
0
votes
1answer
29 views

How can I extract all known mutations of my BAM (or SNP/INDL files)?

I am using a Genome Explorer tool to see all the mutations on my own DNA. My particular interest is on listing the variants and get their names/ids. Here are a few screenshots: You can see that ...
2
votes
0answers
64 views

Global SNP frequency database?

Is NCBI's dbSNP the most comprehensive database of SNP frequency? I'd like to get the largest possible database of allele frequencies per 'population' or geographic region. For example, on average, ...
0
votes
0answers
21 views

General question: How to assign my SNP (Wheat, Illumina 15K) to specific chromosomes?

I have 2000 lines of wheat genotyped with a 15K illumina chip. I also have the name of the SNPs. I would like to know how could I find out which SNPs are in each chromosomes.
1
vote
1answer
39 views

Database of RNA-Seq and SNP genotyping data

I need to obtain publicly available RNA-seq + SNP genotyping data for a cohort - ideally one where the RNA-seq data has been sequenced from brain tissue. I have found CMC and GTEx to be sufficient ...
0
votes
1answer
47 views

How to analyze co-occurrence of multiple SNPs?

I am interested in 20 different SNPs that all are either As or Gs, and they all occur on the same chromosome. How can I assess the co-occurrence of these SNPs? In other words, I want to know, if SNP1 ...
2
votes
0answers
39 views

Subsetting SNPs with specific exclusion criteria using bcftools

I am trying to subset SNPs from 32 cultivars. I have the following exclusion criteria: >20% missing data SNPs present in less than 50% of genotypes Each genotype must have at least 50% of included ...
2
votes
1answer
119 views

Allele count is zero

I ran freebayes on bacterial data and got a vcf file. Some of the SNPs found have AC=0, but in the bam file I can see that there are alternative alleles, so in my ...
2
votes
1answer
91 views

How can I download the SNP genotyping information from GTEX?

I am new to the field of Bioinformatics. GTEx database is great; I know how to check whether a SNP is eQtl with the online tools of GTEx. But my question is how to download the right data to draw a ...
1
vote
1answer
52 views

REMC SNP annotations for transcriptome prediction - epiXcan

I have been trying to use a pipeline for predicting gene expression for a target gene (Although you can do it for many) - the name of this method is epiXcan. Here is a link to the paper: https://www....
1
vote
1answer
50 views

How to calculate Centimorgan lengths from run lengths in one-to-one dna profile comparison

For easier management of a database of dna profiles, I'm creating a tool to perform one-to-one matching on two profiles with snp genotype information. Basically, something that works similar to ...
0
votes
2answers
51 views

Q: What analyses can I perform with a completely phased genome assembly?

Phasing a genome is the process of determining which variants lie on which copy of each chromosome. For example, position 5,430,500 might have an A on the paternal chromosome, but a T base on the ...
0
votes
2answers
44 views

Variant calling for a subset of genes using whole gneome sequencing data

I have few 100 raw fastq files from whole-genome sequencing data and I would like to map these files to a set of genes only (and not whole genome) so as to find SNP's associated with them. Can anyone ...
0
votes
1answer
27 views

How to Convert From Gene Location to SNP Location

I am new to this. I have an SNP array composed of digits [0,1,2] for chromosome 21. How do I find the start:end locations for the region that is "21q22.11" step by step?
1
vote
0answers
170 views

Running IBD using PLINK --parallel (v2)

I have merged two huge .gen files with roughly 500k samples in total - around 3 TB. I wish to run an IBD (identity by descent) check to identify if there are ...
0
votes
2answers
55 views

Questions regarding the detection of variants (SNPs and Indels) in mammal

Greetings, I’m a newbie in bioinformatics working currently on my first project regarding detection of variants (SNPs and Indels) in mammal assembly genome. I used the following method to detect ...
0
votes
1answer
113 views

Simulating phenotype with the 1000 Genomes Project

I'm looking for a way to simulate phenotypes against a real SNP data source, such as the 1000 Genomes. It must be free for commercial purpose (Eg.: MIT license). Any recommendation? I'm trying to use ...
1
vote
0answers
31 views

How to segment genome into homozygote/heterozygote blocks based on RNA-Seq data

Given a RNA-Seq data of a progeny between two inbred lines (for example an F4 plant which was derived from two inbred lines (parents) and the F1 selfed 3 times; however this is probably not crucial ...
1
vote
2answers
70 views

1000 genomes missing SNPs on chr X in GRCh38

I am trying to get a reference together for imputation using the 1000 genomes data set. The GRCh38 version of the data set was made available earlier this year, but seems to be missing a significant ...