vcflib GPAT++: How to select all individuals as targets?
1
0
Entering edit mode
6.1 years ago
spiral01 ▴ 110

I am attempting to run the iHS tool from the GPAT++ suite. It requires the target individuals to be listed as a zero based, comma separated list. I would like to apply this to all individuals in the 1000 genomes phase 3 vcf files. My command is:

iHS --target <(seq -s ',' 0 2503) --file ../chr1.vcf.gz --region 1:860694-861694 --type GT

Here I am using a substitution in bash to create a comma separated sequence of all individuals. However, GPAT rejects this:

INFO: there are 1 individuals in the target
INFO: target ids: /dev/fd/63
INFO: file: ../chr1.ancestral.vcf.gz
INFO: set seqid region to : 1:860694-861694
FATAL: target option is required -- or -- less than two individuals in target

Is there a way of selecting all individuals? The --target argument is required to run the program.

Just to confirm, I have tried running iHS with the argument --target 0,1,2,3,4,5,6 just to confirm there are no other issues with my input, and this runs absolutely fine.

SNP vcf vcflib GPAT++ • 1.2k views
ADD COMMENT
1
Entering edit mode
6.1 years ago
spiral01 ▴ 110

The argument requires a string as input. Therefore:

sequence="$(seq -s ',' 0 2503)"
echo $sequence

produces the sequence as one long string of numbers.

ADD COMMENT

Login before adding your answer.

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