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.
Usage
Version | module name |
---|
0.933 | angsd/0.933 |
#!/bin/bash
#SBATCH -p lts
#SBATCH -t 60
#SBATCH -n 1
#SBATCH -N 1
echo "This examples downloads sample data if not present"
if [[ ! -d bams ]]; then
if [[ ! -f bams.tar.gz ]]; then
wget http://popgen.dk/software/download/angsd/bams.tar.gz
tar -xvzf bams.tar.gz
fi
module load samtools/1.10
for i in bams/*.bam
do
samtools index $i
done
ls bams/*.bam > bam.filelist
module unload samtools/1.10
fi
module load angsd
angsd -b bam.filelist -GL 1 -doMajorMinor 1 -doMaf 2 -P 5
angsd -b bam.filelist -GL 1 -doMajorMinor 1 -doMaf 2 -P 5 -minMapQ 30 -minQ 20 -minMaf 0.05
Links
For more information visit http://www.popgen.dk/angsd/index.php/ANGSD