-
Notifications
You must be signed in to change notification settings - Fork 2
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Confusion with tgls file #8
Comments
Hi Erick,
Well, I've been trying to remember what I was thinking back when I
implemented this one (quite a few years ago). I recall going off of a VCF
spec document, but you're of course right the PL entry given in VCF files
has 3 values. Well, the best advice I can give at the moment is to actually
use GQ instead of PL. If you don't have GQ values called in your data, it
should be the second highest number among the three PL. In practice, our
lab has only ever used GQ.
…-Zachary
On Tue, Jun 11, 2024 at 2:43 PM erickfigue ***@***.***> wrote:
Hi I'm Erick,
I've been using Garlic with only tped and tfam files, I'm trying to
include the tgls file, but something that is not clear for me is the
likelihood value format. I'm trying with the file in phred scale (PL),
since I'm getting it from the vcf using bcftools. In the example tgls file
there is only one value for each individual, I assume is the GL for the
called genotype. In my case I'm getting three values for the PL data from
my vcf file, for each individual (each possible genotype), and Garlic
accept it with no error. My doubt is if the program is identifying the
correct value from the three scores in the input or I need to extract only
the PL score for the called genotype.
Thanks in advance.
—
Reply to this email directly, view it on GitHub
<#8>, or unsubscribe
<https://github.com/notifications/unsubscribe-auth/ABAKRQVZCEKD7EL43ETCVZLZG5AM5AVCNFSM6AAAAABJE4STK6VHI2DSMVQWIX3LMV43ASLTON2WKOZSGM2DOMBWGQYDANI>
.
You are receiving this because you are subscribed to this thread.Message
ID: ***@***.***>
|
Understood. |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
Hi I'm Erick,
I've been using Garlic with only tped and tfam files, I'm trying to include the tgls file, but something that is not clear for me is the likelihood value format. I'm trying with the file in phred scale (PL), since I'm getting it from the vcf using bcftools. In the example tgls file there is only one value for each individual, I assume is the GL for the called genotype. In my case I'm getting three values for the PL data from my vcf file, for each individual (each possible genotype), and Garlic accept it with no error. My doubt is if the program is identifying the correct value from the three scores in the input or I need to extract only the PL score for the called genotype.
Thanks in advance.
The text was updated successfully, but these errors were encountered: