SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
annovar output - the meaning of REJECT? dkrtndhkd Bioinformatics 7 03-25-2014 01:59 AM
BFAST and Variant Calling nexgengirl Bioinformatics 8 01-02-2013 06:03 AM
variant calling with SAMtools krueml Bioinformatics 2 09-28-2012 06:36 AM
annovar output - question dkrtndhkd Bioinformatics 3 04-29-2012 10:23 AM
variant calling kjaja Bioinformatics 1 11-04-2011 08:16 AM

Reply
 
Thread Tools
Old 01-14-2013, 02:53 AM   #1
MrFronk
Junior Member
 
Location: Mordor

Join Date: Jan 2013
Posts: 3
Question Weird annovar output after variant calling

Hello all,

I ran Annovar recently and I'm currently analyzing the results.
I gained an output file (tab delimited file) with 11 colums.
An example entry is given here:
avsift 0.53 chr1 10335315 10335315 G A hom 28.1 6 20

I followed the steps on this website:
http://www.openbioinformatics.org/an...ut.html#pileup
where i went to link 2 "Samtools genotype-calling pileup format"
and then looked at the 3rd code block.

The example stated on their website is pretty straight forward and easy to understand. The problem however is that in my case i have 11 instead of 9 columns, and the last column (used for the number of mutations found in the total amount of reads for that position) is 99% of the time '20'.

This is also the case for a read depth lower then 20.
(So there are cases of a nucleotide found 4 times total, with 20 of those being labeled as SNP...)


My question is: Did other people run into this or do people know how to fix this?

My current sollution would be a workaround to make a script which searches for the SNP in the annovar output format and then generate the last column by searching for the given position in the mPilup file.
MrFronk is offline   Reply With Quote
Old 01-17-2013, 04:30 AM   #2
MrFronk
Junior Member
 
Location: Mordor

Join Date: Jan 2013
Posts: 3
Default

Pease, anyone?
MrFronk is offline   Reply With Quote
Reply

Tags
annovar, column, mpileup, 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 08:50 PM.


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