SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
No output from samtool's mpileup mikheyev Bioinformatics 6 03-18-2014 05:58 AM
Unstable output from mpileup/samtools KNS Bioinformatics 14 10-25-2011 01:38 AM
pysam IteratorSNPCall output different than the samtools mpileup aggp11 Genomic Resequencing 0 07-13-2011 10:25 AM
Cannot understand Tophat output... Help! gen2prot Bioinformatics 4 05-20-2010 07:05 AM

Reply
 
Thread Tools
Old 12-05-2010, 10:43 AM   #1
skblazer
Member
 
Location: Massachusetts

Join Date: Feb 2009
Posts: 50
Default How to understand the output of mpileup like this

vcf was generated by samtools-0.1.12
In the following line of output:

1 8075 . A C 99 . DP=7;AF1=0.9999;CI95=0.5,1;DP4=0,0,2,5;MQ=37 PL:GT:GQ 0,0,0:1/1:39 143,21,0:1/1:60

It tells me that the genotype likelihoods of sample1 are P(D|AA)=1, P(D|AC)=1, P(D|CC)=1?
skblazer 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 07:18 AM.


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