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.

ANGSD: Difference between revisions

From angsd
Jump to navigation Jump to search
 
(81 intermediate revisions by 2 users not shown)
Line 1: Line 1:




ANGSD is a software for analyzing next generation sequencing data. The software can handle a number of different input types from mapped reads to imputed genotype probabilities. Most methods take genotype uncertainty into account instead of basing the analysis on called genotypes. This is especially useful for low and medium depth data. The software is written in C++ and has been used on large sample sizes.


=About=
This program is not for manipulating BAM/CRAM files, but solely a tool to perform various kinds of analysis. We recommend the excellent program [http://samtools.sourceforge.net/ SAMtools] for outputting and modifying bamfiles.
 
ANGSD is a software for analyzing next generation sequencing data. The software can handle a number of different input types from mapped reads to imputed genotype probabilities. All methods talk genotype uncertainty into account instead of basis the analysis on called genotypes. This is especially useful for low and medium depth data. The software is written in C++ and can handle thousands of samples


ANGSD is also on github: https://github.com/ANGSD/angsd
<!-- **
==Overview of input and intermediary data==
==Overview of input and intermediary data==
The input and intermediary data structures of angsd.


<classdiagram type="dir:LR">
<classdiagram type="dir:LR">
Line 33: Line 35:
<classdiagram>
<classdiagram>
//[input data|beagle output{bg:orange}]->[genotype;probabilities]
//[input data|beagle output{bg:orange}]->[genotype;probabilities]
[genotype;probabilities]->[output|genotype calling;MAF estimates;associations{bg:blue}]
[genotype;probabilities]->[output|genotype calling;MAF estimates;associations;SFS{bg:blue}]
  </classdiagram>
  </classdiagram>
-->


 
=Synopsis=
===Commands===
gatk
 
<code lang=sh>
angsd.g++ -outfiles samtools  -GL 1 uppile -b ceu5.mapped.list
</code>
soapsnp
If no recalibration matrix exists these will be created first.
<code lang=sh>
angsd.g++ -outfiles soapsnp  -GL 0 uppile -b ceu5.mapped.list
</code>
gatk
 
<code lang=sh>
<code lang=sh>
angsd.g++ -outfiles gatk  -GL 2 uppile -b ceu5.mapped.list
./angsd [OPTIONS]
</code>
</code>


example of allele frequency estimated from genotype likelihoods with bam files as input using 10 threads


<code lang=sh>
<code lang=sh>
angsd.g++ -outfiles suyeon  -doCounts 1 -qs 20 -doLike 1 uppile -b ceu5.mapped.list
./angsd -out outFileName -bam bam.filelist -GL 1 -doMaf 1 -doMajorMinor 1 -nThreads 10
</code>
</code>


=Possible Bugs=
=Platform=
0.18 If we don't get any reliable genotypelikelihoods for a site, the site wont be included in the .mafs file. If this is want we want I can't say.
The program is developed on tested on a Linux system with gcc compiler. It compiles on OSX with clang, but OSX is not really that tested.
 
0.20 the major minor, doesn't work for the 3genotypes likelihood format (beagle)
 
0.21 When using soap files do we want to infer major/minor from generated likes or from counts?
 
0.24 Should we plugin the keepInd vector for the realSFS
 
angsd
 
0.02 modify the -doGeno to set -doMaf to avoid a segfault with pars->results->asso->freq
 
0.03 when dumping beagle we need to -doMaf 2 otherwise segfault
 
0.3 clean version with less features. The lost features will be reintroduced later.
 
=Wish list=
1. start stop to work with soap input (Not gonna happen I think)

Latest revision as of 14:07, 4 December 2015


ANGSD is a software for analyzing next generation sequencing data. The software can handle a number of different input types from mapped reads to imputed genotype probabilities. Most methods take genotype uncertainty into account instead of basing the analysis on called genotypes. This is especially useful for low and medium depth data. The software is written in C++ and has been used on large sample sizes.

This program is not for manipulating BAM/CRAM files, but solely a tool to perform various kinds of analysis. We recommend the excellent program SAMtools for outputting and modifying bamfiles.

ANGSD is also on github: https://github.com/ANGSD/angsd

Synopsis

./angsd [OPTIONS]

example of allele frequency estimated from genotype likelihoods with bam files as input using 10 threads

./angsd -out outFileName -bam bam.filelist -GL 1 -doMaf 1 -doMajorMinor 1 -nThreads 10

Platform

The program is developed on tested on a Linux system with gcc compiler. It compiles on OSX with clang, but OSX is not really that tested.

Retrieved from "https://www.popgen.dk/angsd/index.php?title=ANGSD&oldid=2637"