ANGSD: Analysis of next generation Sequencing Data
Latest tar.gz version is (0.938/0.939 on github), see Change_log for changes, and download it here.
Introduction: Difference between revisions
Jump to navigation
Jump to search
No edit summary |
No edit summary |
||
| Line 4: | Line 4: | ||
./angsd -bam list.txt -GL 1 -r chr8:12500000 -out results -doGlf 1 | ./angsd -bam list.txt -GL 1 -r chr8:12500000 -out results -doGlf 1 | ||
</pre> | </pre> | ||
This is using the SAMtools model for calculating genotype likelihoods see [[Genotype likelihoods from alignments | genotype likelihoods]] for other models. Without the '-doGlf' parameter the genotype likelihoods will solely be estimated but will not be dumped. Therefore we also supply the | This is using the SAMtools model for calculating genotype likelihoods (see [[Genotype likelihoods from alignments | genotype likelihoods]] for other models). Without the '-doGlf' parameter the genotype likelihoods will solely be estimated but will not be dumped. Therefore we also supply the '-doGlf 1' which will dump the genotype likelihoods as binary see [[Genotype_likelihoods_from_alignments#output_genotype_likelihoods |Gl outputs]]. | ||
Revision as of 08:15, 23 October 2012
Assuming you have a list of bamfiles in in file: 'list.txt' and you want the genotype likelihoods for position 12500000 on chromosome 8.
./angsd -bam list.txt -GL 1 -r chr8:12500000 -out results -doGlf 1
This is using the SAMtools model for calculating genotype likelihoods (see genotype likelihoods for other models). Without the '-doGlf' parameter the genotype likelihoods will solely be estimated but will not be dumped. Therefore we also supply the '-doGlf 1' which will dump the genotype likelihoods as binary see Gl outputs.