SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
How to detect Heterozygous (double) peaks from Sanger shuang Bioinformatics 10 01-25-2013 01:22 PM
samtools/mpileup heterozygous SNPs calling combiochem Bioinformatics 4 08-02-2011 08:05 AM
calling Heterozygous SNPs with samtools mpileup egatti Bioinformatics 1 07-21-2011 09:16 AM
Hi all! [SNP calls from Illumina w/o reference] gdilla Introductions 0 08-13-2010 01:29 PM
Manually correct heterozygous indels captainobvious Bioinformatics 2 03-03-2009 11:07 AM

Reply
 
Thread Tools
Old 01-03-2011, 05:34 AM   #1
Chiel
Junior Member
 
Location: Netherlands

Join Date: Sep 2008
Posts: 2
Default None reference heterozygous calls

Hello all,

In a single sample resequencing project we are trying to extract snps from a BAM file. We are also interested in heterozygous snps that differ from the reference allele. Say in a case where the reference shows a T, in theory we can sequence a C/G in a individual. (or S in IUPAC coding.)

I used mpileup as described on the samtools page, but only supplied a single sample. (Normal pileup deprecated) When I see this profile in the pilup: ggggggCccggg

This is what the vcf file tells:
chrY position . T G,C 29.3 . DP=12;AF1=1;CI95=0.5,1;DP4=0,0,1,7;MQ=41 PL 75,28,64,13,0,64

Yet, the vcf documentation states: "ALT comma separated list of alternate non-reference alleles called on at least one of the samples.."

So the first question: How can you still discriminate between samples if you would want to do so?

Second: The outputed format is not accepted by seattleseq snp annotation which we would like to do. How can we fix this, or are there other methods we can try?

I feel the vcf format may be a helpfull addition to a standard format. Yet I get the idea there are still some inconsistencies which makes pipelining these data a frustrating job. Hopefully this thread is welcome in the bioinformatics forum. Any comments are welcome, thanks!

Last edited by Chiel; 01-03-2011 at 05:35 AM. Reason: Title typo
Chiel 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 01:23 PM.


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