Question: Input for VariantsToBinaryPed
1
gravatar for archie
16 months ago by
archie90
India
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 https://easygwas.ethz.ch/faq/view/15/

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

Archana

plink gatk • 386 views
ADD COMMENTlink modified 4 weeks ago by chrchang5235.7k • written 16 months ago by archie90
0
gravatar for chrchang523
4 weeks ago by
chrchang5235.7k
United States
chrchang5235.7k 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 4 weeks ago by chrchang5235.7k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2266 users visited in the last hour