Beagle log not matching SNP data
0
0
Entering edit mode
15 months ago
Gerard • 0

Hello,

I am trying to impute the genotypes of 32 individuals to WGS with Beagle 5.4 and I found inconsistences between beagle logs and my SNP data. I run Beagle with this:

java -Xss51m -Xmx64g -jar beagle.22Jul22.46e.jar gt=genotype_Chr1.vcf ref=training_ref_Chr1.vcf.gz out=genotypeChr1_imputed

These are the first lines of the log:

No genetic map is specified: using 1 cM = 1 Mb

Reference samples:                1,842
Study     samples:                   32

Window 1 [1:695-40000660]
Reference markers:            2,304,542
Study     markers:                    4

Which makes no sense since in my SNP data (genotype_Chr1.vcf) there are more than 4 SNPs in the 1:695-40000660 window. Here are the first lines of the corresponding .bim file :

1   Hapmap43437-BTA-101873  0   135098  G   A
1   ARS-BFGL-NGS-16466  0   267940  A   G
1   EuroGMD_SNP_BTA1_rs110119084    0   349826  A   G
1   EuroGMD_SNP_BTA1_rs41567297 0   351236  A   G
1   EuroGMD_SNP_BTA1_rs109566607    0   351516  A   C
1   EuroG10K_ARS-BFGL-NGS-105096    0   353745  A   G
1   Hapmap34944-BES1_Contig627_1906 0   393248  A   C

Is this supposed to happen? Why is Beagle only using 4 of the hundreds/thousands of markers present in this window? I'm fairly new to beagle, so I do not know if I'm doing something wrong.

Thanks

Beagle Imputation • 380 views
ADD COMMENT

Login before adding your answer.

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