Questions tagged [genotyping]
Determination of the genotype (for example SNP, allele, heterozygote mutant, etc.) of a sample.
28
questions
1
vote
1
answer
129
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
0
votes
0
answers
71
views
Combine gvcf for whole population or for different phenotypic groups?
I created 400 or so gvcf files (g.vcf.gz) using the gatk HaplotypeCaller. These files show populations for three distinct traits (HR, MR, HS).
Now, I want to use gatk GenomicDBImport to combine gvcfs ...
- 51
2
votes
2
answers
59
views
Reference for the statistical methodology used in joint genotype calling
I am searching for a good reference that explains the statistical methodology used in what is variously known as Joint Genotyping, Cohort Genotyping or Multi-sample Genotyping. One with a pseudo code ...
- 141
1
vote
1
answer
20
views
Definition of genotype in demuxlet
I am reading the Online methods of the demuxlet paper.
The genotype $g$ is taken from the set $\{0, 1, 2\}$, defined as "the true genotype of the sample corresponding to $c$-th droplet at $v$-th ...
- 1,783
2
votes
1
answer
54
views
What is the best way to process yeast genomes?
I have obtained several hundred raw, unassembled yeast genomes from NCBI and I am looking for advice on how to process the genomes for downstream analysis.
I have a reference genome (S288C) to use for ...
- 804
3
votes
1
answer
57
views
Is there an alternative to bulked segregant analysis for insects?
This strategy seems to be most commonly used for plants. When crossing animals isn't possible, how can I do a similar study to identify a particular locus responsible for a polymorphic trait (...
- 307
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
6
votes
1
answer
878
views
Why does GATK produce both 0/1 and 1/0 genotypes in the same file? Are the two not equivalent?
I have always thought that 1/0 and 0/1 in VCF genotype fields are equivalent. And yet, GATK uses both. For example, these are ...
- 9,422
1
vote
1
answer
446
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
1
vote
0
answers
94
views
Genotype quality in VCF for haploid calls
I cannot find in the VCF specs what is the default value of Genotype Quality GQ for haploid calls. In such case, is the default the missing value ...
- 1,783
1
vote
1
answer
152
views
Meaning of order of alleles in VCF GT field for unphased genotypes
Order of alleles in genotype field (on VCFs) can be used to infer haplotypes
However, what about unphased genotypes? I have a VCF file where I only have unphased genotypes:
...
- 1,783
2
votes
0
answers
496
views
What is the correct way of setting the genotype after splitting multi-allelic sites in a VCF file?
This question was also asked on Biostars
Consider this minimal VCF file which contains a multi-allelic call made by a somatic variant caller on a human tumor sample:
...
- 9,422
1
vote
1
answer
266
views
genotyping or variant calling in R
Are there any good options for calling variants or even just a pileup from R? R is not ideal for this, but I'd like to integrate with other functions.
I found deepSNV::bam2R which roughly does what I ...
- 2,159
2
votes
1
answer
52
views
Excess average estimated identical-by-descent in genotype data
(Cross-post with Biostars: https://www.biostars.org/p/470788/)
We have some genotype data that we are putting through quality control in PLINK 1.9. As part of this QC, we have limited the data to ...
- 133
1
vote
0
answers
88
views
Male DNA genotyped with Illumina GSAMD gives two alleles for certain loci within sex chromosomes
I have a tab-separated file (the columns are rsID CHROM POS GENOTYPE) generated by Diganomics with Illumina GSAMD (v2.2) from male DNA, and there are a few hundred loci within X or Y chromosome that ...
- 151
1
vote
1
answer
47
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 ...
- 836
1
vote
1
answer
74
views
Phasing partially phased genomic data
I am trying to create an integrated callset with coding regions from phased exomes and noncoding regions from a genotyping array. I would like to merge the data together in the correct phase. Since ...
- 841
2
votes
1
answer
62
views
Can someone recommend software for reference-based genotype-estimation?
Problem: I have low-coverage sequencing data from humans in vcf format. For each sample, I have PL scores / preliminary genotype calls from the software ATLAS (it's a bit like GATK for those who don't ...
- 1,621
1
vote
0
answers
558
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 ...
- 836
0
votes
2
answers
487
views
Proper use of BWA MEM on multiplexed GBS sample
I have a multiplexed lane of GBS sequencing reads as a fastq file. I understand the first step is to demultiplex and trim the adapter sequences from the reads.
This yields many individual fastq files ...
2
votes
1
answer
266
views
PLINK - Transposed BED file?
I'm currently working with the PLINK file system to store large amounts of genotype data.
The plain format consists of three files, two files for storing phenotype and marker information as well as ...
- 123
1
vote
3
answers
49
views
Genotyping technologies do not maintain phase?
Why exactly does genotyping not maintain phase? I guess I'm more confused on what's the difference between sequencing and genotyping, and as a result how phase is maintained in one, but not the other. ...
- 341
1
vote
0
answers
201
views
GBS: clustering of PCoA axes with R package mclust to describe identity by missingness
I have a SNP dataset with 56 populations and 430 samples that was produced by GBS sequencing. I need to identify groups of samples that have similar patterns of genotype missingness, since there is ...
- 183
3
votes
1
answer
84
views
using SNPs to identify mixed samples
Is there a way to identify mixed samples based on SNPs?
Example input (table of genotypes for multiple samples):
...
- 2,159
4
votes
4
answers
1k
views
Two variants associated with same chromosomal position in genotype data for one individual
I'm trying to play around with my family's raw genetic data from Ancestry.com so that I can do some genetic analyses, and I noticed in the AncestryDNA raw data files that some regions of the genome ...
- 465
4
votes
1
answer
88
views
Make ipyrad use cuda-enabled NVidia card on Ubuntu
I want to use ipyrad on a new Ubuntu machine that has an NVidia Quadro K2000 card with 384 cores. One can configure ipyrad to run on a linux cluster. Do I have any options to get ipyrad to access ...
- 183
3
votes
1
answer
391
views
Techniques for analyzing and quantifying sample bleed through in genotyping with Illumina
I am looking at the presence of viral genotypes within individual samples within an assay. Often times there is a sample whose read counts are firing off the charts and this sample tends to "bleed ...
- 255
21
votes
2
answers
4k
views
Single-sample vs. joint genotyping
I am trying to understand the benefits of joint genotyping and would be grateful if someone could provide an argument (ideally mathematically) that would clearly demonstrate the benefit of joint vs. ...
- 313