Eagle phasing and minimac imputation
0
0
Entering edit mode
7.1 years ago
dp0b ▴ 80

Hi,

I wish to use eagle phasing and minimac imputation to impute to whole genome sequence.

I have successfully phased my genotypes using eagle but Im struggling to read them into minimac. Minimac will tell me there are no marker overlaps between the target and reference, when there definitely is. Can minimac read in eagle phased files directly?

My original files for phasing with eagle were plink binary files and when I specify vcf format for minimac they read in but then I get the no overlap error I was hoping someone could help me with the input line; I have tried the following;

./minimac --vcfReference --refHaps phased_seq.haps --haps impute_ani_phased.haps --snps target.snps --rs --rounds 5 --states 200 --prefix minimac_imputed

./minimac --vcfReference --refHaps phased_seq.haps --shape_haps impute_ani_phased.haps --chr 29 --sample impute_ani_phased.sample --rs --rounds 5 --states 200 --prefix minimac_imputed

the target snps was a file in the format

29:100
29:2000

and my phased files were in the format

29 29:16860 16860 A G 1 1 1 1
29 29:168604 168604 G A 1 1 0 1
29 29:168680 16860369 G A 1 1 1 1

Many thanks

eagle minimac imputation phasing genotypes • 3.0k views
ADD COMMENT
0
Entering edit mode

I am wondering how the issue is solved?

ADD REPLY

Login before adding your answer.

Traffic: 1830 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