Question: VCF file phasing by SHAPEIT
gravatar for seta
6 months ago by
seta1.2k wrote:

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 sequencing, the vcf file obtained by GATK. As I searched, SHAPEIT was mostly used; based on its manual, it requires genetic map for phasing, however, the provided link for genetic map is based on hapmap, hg37, which the link didn’t work (actually, an error says “The requested URL /genetics_software/shapeit/shapeit.html/files/genetic_map_b37.tar.gz was not found on this server”. Now, my questions are:

1) Could you please tell me where is genetic map?, I also need this map based on hg38, is there the genetic map for hg38, or how we can convert this map from hg37 to the related map for hg38?

2) In SHAPEIT manual, “read aware phasing” also described that takes bam and vcf file as input to extract the phase informative read (PIR) that used for phasing vcf file in the next step. So, the genetic map is no longer required, here. I think I should use this method (not that is based on genetic map) since I have the sequencing data, yes, am I right?

3) Also, please kindly let me know if it is possible to use a subset of interest from the vcf file of a given chromosome and the related part bam file extracted from whole bam file (so, not use whole bam and vcf file) for phasing?

Any suggestion and help would be highly appreciated.

ADD COMMENTlink modified 6 months ago • written 6 months ago by seta1.2k

I think you mean SHAPEIT ( ) instead of Shapiet, right?

ADD REPLYlink written 6 months ago by h.mon26k

Sorry, yes, I corrected it. Could you please kindly help me? I don't know why I cannot post the question in OXSTATGEN mailing list.

ADD REPLYlink modified 6 months ago • written 6 months ago by seta1.2k

I would contact the authors directly. They have a duty to follow-up on these things if they are still advertising their program for current usage. Feel free to point them to my comment, here.

ADD REPLYlink written 6 months ago by Kevin Blighe44k

Thanks, Kevin for your nice point. At the moment, I would like to try Eagle 2 due to the more accuracy and speed compared to Shapeit as its authors said.

ADD REPLYlink written 6 months ago by seta1.2k

Hi, have you found a solution? I need this genetic map for hg38 too. Thanks.

ADD REPLYlink written 4 months ago by janka20120

Please contact the authors of SHAPEIT.

ADD REPLYlink written 4 months ago by Kevin Blighe44k

I found the genetic map from 1KG Phase3 files under the IMPUTE2 ftp here:

ADD REPLYlink written 10 weeks ago by miaowzai140
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1632 users visited in the last hour