Genotypes in vcf files
Hi! I have a question my VCF files have weird genotypes in it like some with 0/2 or 0/3, some with 1/2, 1/9, 2/2, 2/1, 4/13, 6/6. Should I only keep 1/1 as homozygotes?
genotypes
exome
• 1,213 views
•
link
updated
by
jena •
•
written
by
mmghaffar29 •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Understanding Uncommon Genotype Formats in VCF Files Generated by nf-core RNAvar Pipeline
written by ASid •Hello BioStars Community, I've been working with VCF files generated by the nf-core RNAvar pipeline, which employs the GATK best practices for RNA-seq variant calling. …
-
Removing indels in VCF file
written by ColaptesHello, I am trying to do something very simple, but running into confusing behaviour. I have a VCF file of multiple samples and want to …
-
Count columns and delet specific one
written by salvatore_d •Hi All I've merged a lot of different files with different columns number (12 and 13) and now I'm looking for a way to delete …
-
4/6 genotype in VCF
written by ajuiwl •Hi all, I have been doing variant calling analysis in human cancer exomes and some variants in the VCF file have genotype like 4/6 or …
-
Compute correlation for depth of coverage between 2 datasets with different length in R
written by AnaHello everyone, I am trying to compute correlation between depth of coverage for individuals with homozygote genotypes vs. individuals with heterozygote genotypes. here is the …
-
Coverage of methylation txt files
written by Pin.BioinfHello, I have been provided with methylation txt files for target_regions that look like this: chr position meth unmeth total degree.of.meth chr1 1058346 1 1 …
-
Unlisting lists inside data frame and collapsing their values in R
written by arronarHello I have this data.frame and as you can see there are lists inside some cells. myList1 <- list() myList1[[1]] <- 0 myList1[[2]] <- list(3) …
-
Help with HapTree reads input format
written by phoenixice154 •Dear All, I am going to try HapTree to phase the SNP/MNP/InDel/Complex in my VCF, and want to give HapTree a try. But I really …
-
Can these fusion predictions by tophat be trusted
written by DataminerHi, I have got these results from tophat-fusion: ``` chr15-chr17 74317267 38506035 ff 1 13 1 7 43 33 10.000000 @ 13 26 38 51 …
-
how to compare sets using python (dealing with PDB file)
written by Jason Lin •Hi all, Sorry to bother you all again. so I have a text file which contains the PDBID and corresponding missing coordinates from PDB file. …
Please provide a snippet of a few lines from your file, and also give some more details on the source and contents of your VCF.
Hi and welcome. How many alternative alleles do you see in the ALT column at these positions? The encoding is as follows: 0=REF, 1=1st ALT, 2=2nd ALT and so on. So for example the genotype 0/3 means a heterozygous individual with one reference allele (REF) and one allele describer by the third entry in the ALT column.
Since there are only 4 possible bases, the other alleles are typically indels. You may wish to remove those (with packages like
bcftools) or keep them, depending on the questions of your study.