ANGSD

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

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