Comments (2)
Ok, that is a documentation bug in the header of the vcf. I’ve changed this in the latest commit on github.
It is loglikeratio to the most likely, and scaled as log10.
Let me know if you find other issues with the vcf. I dont think this is the most used part of angsd.
On 22 Sep 2015, at 01:47, Evan McCartney-Melstad [email protected] wrote:
Hello,
I've been trying to output a VCF file to look at IBD between samples, and I was curious if the GL's in the VCF file are scaled in any way. It says they have been scaled to add up to 1 in the header of the VCF, but they don't (the GP fields do sum to 1).Here's the command I used:
angsd0.902/angsd -bam bamlist.txt -out OUTPUT -uniqueOnly 1 -only_proper_pairs 1 -remove_bads 1 -doCounts 1 -nThreads 27 -doMaf 1 -doMajorMinor 2 -GL 2 -minMapQ 20 -minQ 30 -dovcf 1 -doPost 2And here's what my VCF looks like:
##fileformat=VCFv4.2(angsd version)
##FORMAT=
##FORMAT=
##FORMAT=
##FORMAT=
#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT ind0 ind1
scaffold_0 1616 . T A . PASS . GP:GL 0.666631,0.333333,0.000035:0.000000,-0.301007,-4.277052 0.666578,0.333333,0.000088:0.000000,-0.300972,-3.876948Thanks very much for any insight!
Evan—
Reply to this email directly or view it on GitHub #18.
from angsd.
Ok, I'm closing this issue. Feel free to reopen if needed
from angsd.
Related Issues (20)
- Calculating folded SFS to account for highly polymorphic regions
- ANGSD to calculate theta, is the result calculated log with 10 as the base or log with e as the base.
- Getting allele frequencies for different populations cutting a beagle file or using a -sites file?
- In version 0.934 how to convert haplo to tped? haploToPlink command no longer works
- Segmentation Fault HOT 1
- problems with estAvgError.R HOT 1
- Length of fastafile vs length of chr in BAM header HOT 1
- Contamination test, Error in dimnames(x)
- Population level theta with saf2theta
- [loadChr] Error loading fasta info from chr
- Missing windows in thetaStat
- Number of numbers output by realSFS HOT 1
- merging binary all 10 log genotype likelihood
- angsd: realloc() error
- 0 sites retained after filtering when using -sites HOT 1
- estAvgError.R script error
- -argument is unknown will exit HOT 1
- Killed HOT 1
- Problems with unsorted saf file chromoname HOT 1
- SAF file out of sync
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from angsd.