Question: Input for VariantsToBinaryPed
gravatar for archie
2.3 years ago by
archie90 wrote:

Hello everyone,

I predicted SNPs by GATK which are in VCF format. I want to create binary ped file which is intermediate input for another tool.

java -jar GenomeAnalysisTK.jar  -T VariantsToBinaryPed   -R reference.fasta   -V variants.vcf  -m metadata.fam   -bed output.bed  -bim output.bim  -fam output.fam

I have reference genome -R , -V vcf input file. But after doing exploration, I found .fam is need to be given as input to GATK which is metadata information . I do not have information of family_id, paternal id and maternal id. What can I provide here in order to run in GATK ??

My metadata format is

family_id   individual_id   paternal_id     maternal_id     sex         phenotyp    
                  F545                                                          1             1

I found link

The PED file has 6 fixed columns at the beginning followed by the SNP information. The columns should be separated by a whitespace or a tab. The first six columns hold the following information:
Family ID (if unknown use the same id as for the sample id in column two)
Sample ID
Paternal ID (if unknown use 0)
Maternal ID (if unknown use 0)
Sex (if unknown use 0)
Not used, set to 0
Rest of the columns: SNPs

Is it right If I provide 0 in my metadata file in Paternal, meteral ID. Also I am using the same Family ID and Sample ID based on above description. Please correct me If I am doing anything wrong. I am new in this area of data analysis.

Thank you in advance


plink gatk • 595 views
ADD COMMENTlink modified 12 months ago by chrchang5237.3k • written 2.3 years ago by archie90
gravatar for chrchang523
12 months ago by
United States
chrchang5237.3k wrote:

For future reference, the solution is to not use GATK for this; use PLINK --vcf instead. (The VariantsToBinaryPed function was removed from GATK 4 since it no longer has any practical value.)

ADD COMMENTlink written 12 months ago by chrchang5237.3k
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: 1888 users visited in the last hour