Skip to content

failed to identify which build for 1000g reference build on chr17 #123

@celisungmail

Description

@celisungmail

Hello,

I tried to run SSIMP and took ALL.chr{CHRM}.phase3_shapeit2_mvncall_integrated_v5a.20130502.genotypes.vcf (hg19) as reference panel but I found the log file showed that SSIMP identified which_build=hg19 for 1000g reference build on chr17 as following. I am wondering whether there is a way to specify the build for both gwas and reference panel.
Thanks,

[ssimp-0.5.5 (2018-12-26)]
file_name: xxxxgwas_chr17.txt
... loading the 1.7GB database of positions under three builds. This will take about a minute. Loaded.
Estimating which build (hg18/hg19/hg37) of the reference panel and the GWAS file, in case it is necessary to modify the GWAS file to match the reference panel
some_records_from_each_chromosome.size():1000
(ref) :count_of_hg18_0based,count_of_hg19_0based,count_of_hg20_0based,count_of_hg18_1based,count_of_hg19_1based,count_of_hg20_1based: 8 , 8 , 0 , 935 , 935 , 0
gwas_all_chrpos.size():136940
(gwas):count_of_hg18_0based,count_of_hg19_0based,count_of_hg20_0based,count_of_hg18_1based,count_of_hg19_1based,count_of_hg20_1based: 4356 , 16654 , 4228 , 4608 , 130206 , 4202
which_build_gwas,which_build_ref: hg19_1 , hg18_1

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions