use gatk phased vcf as input file of HapCut2
0
0
Entering edit mode
3.0 years ago
Peter Chung ▴ 140

I am new in sequencing analysis and would like to do haplotype blocks from the result of WGS. But I have some questions on the HAPCUT2.

I generated a vcf from gatk and then I phased it using gatk ReadBackedPhasing into phased vcf.

Commands below are from HAPCUT2:

./build/extractHAIRS [options] --bam reads.sorted.bam --VCF variants.VCF --out fragment_file
./build/HAPCUT2 --fragments fragment_file --vcf variantcalls.vcf --output haplotype_output_file

I saw there are two vcf files in different names, so which one should I input my phased vcf and what is another vcf file ? btw, is it better to put a phased vcf in HAPCUT2 ?

Thank you.

HapCut2 gatk phasing vcf • 1.2k views
ADD COMMENT
0
Entering edit mode

Please use the formatting bar (especially the code option) to present your post better. I've done it for you this time.
code_formatting

ADD REPLY
0
Entering edit mode

hello, did you figure it out which vcf file should be used. I only have one vcf file output from the WGS analysis

ADD REPLY

Login before adding your answer.

Traffic: 2507 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6