Heterozygous individuals and FastaAlternateReferenceMarker · Biostar

This is a test version of Biostars. For the public version, visit https://www.biostars.org.

Heterozygous individuals and FastaAlternateReferenceMarker

Good afternoon,

I am having a doubt with the function FastaAlternateReferenceMarker from GATK. I already have a phased VCF file and when I want to obtain the fasta files for some specific genes (giving the target genes sequences and the VCF file to the function) I just obtain one of the sequences. This is okey for homozygous individuals, but I don't know how to obtain the two sequences for heterozygous individuals.

I've been checking the GATK website but I did not find anything, anybody else that had the same issue?

Thanks in advance!!

gatk fasta snp

• 659 views • link

written 6.2 years ago by gubrins

0 answers

No answers yet.

Log in to answer this question.

More posts like this

Multisample VCF to FASTA
written 6.2 years ago by gubrins

Good afternoon, I know there are already quite a lot of posts explaining how to convert a VCF file to Fasta (In fact, I was …
Fasta files from phased VCF
written 6.2 years ago by gubrins

Good afternoon, I have a phased VCF file with the variation of my target genes. I have the reference for each of the genes and …
Count total number of Homozygous/Heterozygous SNPs from VCF file
written 6.2 years ago by nataliagru1

Dear Community, I would like to count the total number of heterozygous and homozygous SNP's in my VCF file. I have read up on other …
quality check of bam files
written 6.2 years ago by gubrins

Good afternoon, I'm working with GATK to do the SNP calling of some target capture sequencing data. Right now I'm creating the bam files and …
Phasing with Beagle
written 6.3 years ago by gubrins

Good afternoon, I'm working with a lizard species and we did a target capture sequencing for the nuclear genes related in the OXPHOS process and …
Coverage per gene
written 6.3 years ago by gubrins

Hey, I'm working with bedtools trying to find out the coverage of several target regions. I have my bam files and my target regions in …
Fastq to VCF pipeline configuration for low coverage data
written 6.8 years ago by igor •

Hi, I have a low coverage (1.0x) fastq sample. I am following GATK best practices pipeline (bwa mem > Mark Duplicates > BQSR > Apply …
Extract all genotype counts from phased data in vcf files
written 8.4 years ago by spiral01

I have vcf files (1000 genomes) with phased data and would like to extract counts of the genotypes. Several methods were suggested here: https://www.biostars.org/p/218444/ However, …
vcf file filter
written 8.8 years ago by ATCG

I need to filter my vcf file to include variants with at least 30 individuals in each of the possible groups: major allele homozygous, heterozygous, …
difference in genotypes in a multisample vcf file
written 9.5 years ago by prasundutta87

Hi, I have two individuals of the same animal (one male and one female water buffalo). I have RNAseq data from both the animals and …