Questions tagged [genotyping]

Determination of the genotype (for example SNP, allele, heterozygote mutant, etc.) of a sample.

Filter by
Sorted by
Tagged with
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 ...
Shay Maor's user avatar
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 ...
Learner's user avatar
  • 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 ...
Mark's user avatar
  • 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 ...
gc5's user avatar
  • 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 ...
rimo's user avatar
  • 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 (...
Caterina's user avatar
  • 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. ...
Tom's user avatar
  • 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 ...
terdon's user avatar
  • 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....
Rhea Bedi's user avatar
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 ...
gc5's user avatar
  • 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: ...
gc5's user avatar
  • 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: ...
terdon's user avatar
  • 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 ...
burger's user avatar
  • 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 ...
Joni's user avatar
  • 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 ...
Maxim Masiutin's user avatar
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 ...
h3ab74's user avatar
  • 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 ...
Greg's user avatar
  • 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 ...
user438383's user avatar
  • 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 ...
h3ab74's user avatar
  • 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 ...
PlantGeek519's user avatar
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 ...
Scholar's user avatar
  • 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. ...
Jonathan's user avatar
  • 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 ...
Peter Pearman's user avatar
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): ...
burger's user avatar
  • 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 ...
Sarah's user avatar
  • 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 ...
Peter Pearman's user avatar
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 ...
quantik's user avatar
  • 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. ...
llevar's user avatar
  • 313