Calculation of the genotype likelyhood from raw data
Entering edit mode
12 months ago
vit.filippo ▴ 50

I am working on amplicon data from an aploid genome (same length, already aligned with a miltiple alignment). Since I cannot generate a .bam file I ve to call mutations from aligned fastq through base count generation. To validate a mutation I ve to calculate genotype likelihood, but I cannot find a simplified version of the formula, despite I found a 'simpler' formula that expoites counter for GL calculation ( Moreover, all these formulae are not well explained in the page, so I am wondering if someone has already tryied to calculate GL by hand ( having base counts, base qualities and mapping qualities). If possible, I would like to calculate GL in a 'samtools like' manner which is more conservative respect GATK and other tools.

genotype likelihood bam counts raw data samtools • 277 views

Login before adding your answer.

Traffic: 2324 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6