SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
BWA - fail to infer insert size: too few good pairs id0 Bioinformatics 17 06-10-2014 07:27 AM
BWA 0.6.1 fail to insert infer isize: weird pairing + Segmentation Fault attilav Bioinformatics 5 02-20-2012 02:54 PM
How to infer Illumina paired-end strand specificity from SAM output? David Harmin Bioinformatics 0 02-16-2011 08:34 AM

Reply
 
Thread Tools
Old 10-05-2012, 12:46 AM   #1
DonJoe
Junior Member
 
Location: Sweden

Join Date: Oct 2012
Posts: 1
Default Using KING to infer kinship

Greetings!

I'm trying to use the software KING (http://people.virginia.edu/~wc9c/KING/) to check kinship between exome sequence samples. I'm not getting the expected results so I think I might be doing something wrong.

I have followed the procedure posted here: http://cphg.virginia.edu/quinlan/?p=300 to run the progran. My VCF file is created by running GATK UnifiedGenotyper with the options --dbsnp and --output_mode EMIT_VARIANTS_ONLY

The resulting output from KING show negative values for the kinship coefficient, although I know some of the individuals are siblings. Should I generate my VCF file in some other manner, like generate calls at all sites? This produce huge outputs though.

Has anyone else used KING or done similar analysis before? Any help or hints is greatly appreciated!

----------------------

These are the GATK options I've used:

-T UnifiedGenotyper \
-R ref.fa \
-I inputbams.list \
-o GATK_snpcall.vcf \
--output_mode EMIT_VARIANTS_ONLY \
--validation_strictness LENIENT \
-dcov 200 \
-l INFO \
-nt 8 \
--dbsnp dbsnp_135.b37.vcf
DonJoe is offline   Reply With Quote
Old 09-16-2013, 10:12 AM   #2
Zigster
(Jeremy Leipzig)
 
Location: Philadelphia, PA

Join Date: May 2009
Posts: 116
Default

Yeah, I tried Aaron's recipe once and got very low numbers (although siblings were higher than strangers). I think the no-calls is screwing it up.
__________________
--
Jeremy Leipzig
Bioinformatics Programmer
--
My blog
Twitter
Zigster is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 04:04 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO