SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



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

Reply
 
Thread Tools
Old 06-23-2016, 10:19 AM   #1
lstbl
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
Reply

Tags
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 06:38 AM.


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