How to produce fasta sequence from haplotyping (illumina)

Hello,

I would like to produce a fasta file from haplotyping.

I have a reference sequence of 2kb and illumina PE reads, from a tetraploid organism. I see some SNPs in IGV and i would like to separate the variants in a fasta file, without parental information.

I produce VCF file with GATK haplotypecaller but failed to produce sequence with BCFtools and vcflib/vcf2fasta because my VCF file is not phased.

Could anyone help me for this ? Does anyone have a pipeline for that ?

Thanks, Luc

snp

• 1,524 views • link

updated 7.7 years ago by WouterDeCoster

• written 7.7 years ago by luc.cornet •

1 answer

For phasing I would suggest you take a look at WhatsHap or GATK ReadBackedPhasing.

• 1 views • link

written 7.7 years ago by WouterDeCoster

More posts like this

VCF file filter for breakpoints
written 3.4 years ago by EpiExplorer

Hi, I am trying to filter VCF file and extract breakpoints. Does anyone know how to do it using bcftools? I have extracted SV calls …
Phasing a partly-phased VCF
written 4.6 years ago by Sarah •

Hello, I am trying to phase a VCF file for a non-model organism that has no reference panel. I have a VCF file of genotypes …
Generating multiple species consensus alignments for tree topology tests
written 5.2 years ago by mglasena •

I have generated a VCF file with 9 samples, each being a separate species aligned to a common reference. I am hoping to generate multiple …
hapcut2 output and figures
written 6.9 years ago by rob234king

I have run Hapcut2 and have a VCF file. I can't see in the VCF where it has assigned the snp to one haplotype or …
call haplotypes in GATK HaplotypeCaller
written 7.3 years ago by Peter Chung

I am new in bioinformatics. I would like to call haplotypes in GATK HaplotypeCaller, however the output result is 0kb. Can anyone give me some …
Separate the Longest Strain in Bacteria
written 8.0 years ago by Shelle •

I have a FASTA file of bacteria which has multiple strains like below. I would like to go through the FASTA file and extract/separate the …
How to convert a raw vcf to a phased vcf
written 8.6 years ago by pinn

Hi, I generated a vcf file called variants using Freebayes. Further, I would like to perform normalisation for my vcf file I found a tool …
Sanger Imputation Service - The input file sanity check failed
written 11.0 years ago by alesssia

Dear all, I would like to impute a genotype dataset using the [Sanger Imputation Service][1]. I created the VCF file using PLINK (v1.90) and, as …
consensus sequence of complex estructural variants
written 11.5 years ago by user230613

Hi all , I would like to know if there is any tool to create a consensus FASTA sequence having a vcf file with complex …
Modifying A Gene Sequence With Personal Snps?
written 12.6 years ago by Chris Miller

<p>Given exome sequencing of a patient, along with a gene name, we'd like to create a fasta file of that gene's sequence, with the sequence …