error when using shapeit2 phasing with vcf from GATK

my command line: shapeit.v2.r837.GLIBCv2.12.Linux.static/bin/shapeit --input-map /genetic_map_b37/genetic_map_chr3_combined_b37.txt --input-vcf filtered.vcf --output-max shapeit_phasing.chr3.phased

and I got the error: ERROR: Site pos=15274 alt=G,T is not bi-allelic!

I can't find any answer anywhere. 1. shapeit could accept vcf file, but it doesn't say what kind of vcf it can handle. 2. I tried to use GATK VariantsToBinaryPed to transform the vcf file, but I don't have the metaData(.fam) file since I only have one sample.

I'm sorry that I'm new in phasing. Can anybody give me hand? THANKS!

vcf shapeit phasing

• 2,480 views • link

updated 9.0 years ago by h.mon

• written 9.0 years ago by Loveinheaven •

0 answers

No answers yet.

More posts like this

Problem with SHAPEIT
written 4.0 years ago by dongxinyu510 •

My code is as follows,At the end it has an error and I don't know how to fix it: ``` shapeit \ -M genetic_map_hg38_withX.txt \ …
Convert SHAPEIT haps file to vcf
written 4.7 years ago by michael.flower.14 •

I've phased some sequencing using SHAPEIT, producing two files: 125QiPSC.haps 125QiPSC.log I would now like to subset the file for my SNPs of interest. To …
Error when Phasing with Beagle 5.2
written 5.0 years ago by mglasena •

I'm having trouble phasing a multi-sample (9-samples) vcf file produced by gatk HaplotypeCaller with Beagle 5.2. I do not have a genetic map or reference …
SHAPEIT v4 modify keep QUAL FORMAT and INFO field
written 5.2 years ago by nhaus •

Hello, I am using [shapeit v4][1] to phase my germline mutation calls which I got using GATKs HaplotypeCaller (WGS data). As a reference, I am …
Quality checking phasing
written 5.4 years ago by truebeliever24 •

Hi everyone, I recently phased my data using Shapeit2, and I am wondering how I would check the quality of the phasing process on my …
VCF file phasing by SHAPEIT
written 7.7 years ago by seta

Hi everybody, I would like to phase (just phasing, not imputation) vcf file containing about 1100 individuals (a given human population) derived from whole genome …
VariantsToBinaryPed tool missing from GATK4?
written 8.2 years ago by gaelgarcia05

Hi all, I've been prepping my data to input to [GATK's Variant Manipulation Tool `VariantsToBinaryPed`][1] , but I just realized that GATK4 doesn't list this …
Error while running BEAGLE for genotype imputation
written 10.3 years ago by Shab86

I am trying to run BEAGLE 4.1 for an imputation run. I have core exome chip data on variants of 20th chromosome in BED/BIM/FAM format, …
Phasing using SHAPEIT
written 10.9 years ago by genetic •

Hello, I need to use SHAPEIT for **phasing only** since I will conduct **CH (compound heterozygous) analysis for recessive rare variant**. I will not perform …
How to get unphased VCF file phased by shapeit?
written 11.7 years ago by Jie Ping •

I have a vcf file at hand and tried to get it phased by shapeit. The command line was referred [here][1]. [1]: https://mathgen.stats.ox.ac.uk/genetics_software/shapeit/shapeit.html#input like this: …