What happens if you don't put that option, what is the default behaviour?
• 0 views
•
link
written
by
jxiang15 •
What happens to heterozygous sites when you go from reference sequence to sequence modified by variants?
In a previous post, it was recommended to use either vcftools or FastaAlternateReferenceMaker if you have a reference sequence and a variant file, and you to get a new FASTA file.
However, with the 1000 genomes data, the data is phased. So at heterozygous sites, should the ALT allele be substituted or should the REF allele be left in the sequence?
Thanks in advance
• 2,737 views
•
link
updated
by
Ram
•
written
by
jxiang15 •
•
written
by
jxiang15 •
1 answer
FastaAlternateReferenceMaker's --useIUPAC option will place the IUPAC code wherever there is heterozygous state.
• 0 views
•
link
updated
by
Ram
•
written
by
GouthamAtla
•
written
by
GouthamAtla
Log in to answer this question.
More posts like this
-
How to Resolve Reference Mismatch While Extracting Gene Variants from 1000 Genomes Project Phase 3 …
written by Rohan •
-
Non Variant Sites in GVCF File
written by saamhasan55 •
-
Building vcf from most common alleles
written by Eugene •
-
filter vcf to include only sites which have at least one homozygous alt call
written by grey •
-
Vcftools r2 Error: Insufficient sites remained after filtering or No data left for analysis!
written by hellbio
-
Substituting heterozygous sites by randomly sampling one allele
written by outlier95 •
-
BCFTools view equivelent of GATK -T SelectVariants -sn sample_x --excludeNonVariants
written by William
-
how to add QualityOfDepth(QD) annotation to freebayes vcf
written by William
-
How to identify Homozygous or Heterozygous variants from the annotation file ?
written by deepue
-
Ped File Individual Id Same As Bam Rg Sample Name For Gatk?
written by William
