Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
GNU Parallel Losing Records salamay Bioinformatics 0 12-09-2014 10:05 AM
Picardtools MarkDuplicates Question lre1234 Bioinformatics 0 09-11-2014 10:01 AM
Why am I losing up to 5 bases at start of reads? tjs7 Illumina/Solexa 6 10-03-2012 06:31 AM
What is the header format for Interval List in PicardTools? accipiter Bioinformatics 5 04-17-2012 03:37 PM
PicardTools IlluminaBaseCallsToSam Protaeus Bioinformatics 0 05-10-2011 11:06 AM

Thread Tools
Old 06-23-2016, 11:19 AM   #1
Junior Member
Location: USA

Join Date: Aug 2015
Posts: 7
Default picardtools LiftoverVcf losing majority of information

I'm currently attempting to lift over variants from a .vcf file mapped to the hg18 assembly to the chimpanzee panTro4 assembly (UCSC nomenclature used here). I'm using the picardtools LiftoverVcf to accomplish this, as well as the .chain file available on UCSC.

I'm noticing that only about 10% of variants are being lifted over, the rest are ending up in the 'regected' variants file. I'm wondering if this should be expected when lifting over from one organism to another, or if this is a strange occurrence.

Note: the original .vcf (which used hg18 assembly as the reference) was produced with short reads from chimpanzees. So, my guess (if it this is not a bug in the picardtools program) is that the vast majority of "variants" picked up in the .vcf are, in fact, just fixed differences between human and chimpanzee.
lstbl is offline   Reply With Quote

liftovervariants, variant calling

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 09:33 AM.

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