Skip to content
This repository has been archived by the owner on Jan 16, 2019. It is now read-only.
Paul Hoffman edited this page Jan 24, 2016 · 26 revisions

Question: Who do I cite ANGSD-wrapper?

Answer: Currently the appopriate citation is to a preprint (see below). The manuscript is being submitted for peer review. Depending on the methods used in ANGSD-wrapper, you should also cite the ANGSD paper and the papers detailing methods. See our References for details. Durvasula A, Kent TV, Hoffman PJ, Liu C, Kono TJY, Morrell PL, Ross-Ibarra J (2015) ANGSD-wrapper. PeerJ PrePrints 3: e1837

Question: Is there an example dataset I can used with ANGSD-wrapper?

Answer: A small dataset consisting of BAM files for 11 maize and 11 teosinte inbred lines is available on figshare. A sample of Tripsacum dactyloides that can be used as an outgroup is also provided.

Question: Will ANGSD-wrapper be updated on a regular basis?

Answer: Yes! ANGSD-wrapper is being updated to work with the latest version of ANGSD.

Question: Can I use my own version of ANGSD?

Answer: No, you should use the version of ANGSD included with ANGSD-wrapper.

Question: When I try to compile ANGSD, the following error is reported: /usr/bin/ld cannot find -lz

Answer: Generally this error appears when your system doesn't have the required version of zlib. If you are working on computing cluster, contact your system administrator.

Question: How do I call SNPs using ANGSD and ANGSD-wrapper?

Answer: ANGSD was designed to perform most population genetic analyses without the need to make SNP calls. That is, variants are dealt with probabilistically and many common descriptors of polymorphism in a population can be calculated without SNP calls. SNP calling is possible in ANGSD and with the help of ANGSD-wrapper. See the tutorial for details.

Question: Where can I find the step-by-step process for running ANGSD-wrapper?

Answer: This set of tutorials will demonstrate the most common analyses conducted with ANGSD and ANGSD-wrapper.

Question: Do I need an outgroup sequence to use ANGSD and ANGSD-wrapper?

Answer: It is necessary to have an outgroup sequence to infer the ancestral state of the mutations and to estimate the derived site frequency spectrum. However, many of the descriptive statistics calculated by ANGSD are not dependent on the ancestral state of mutations. If not outgroup sequence is available, include the path to the reference sequence when asked for the outgroup.

Question: Do I need to include inbreeding coefficients in my analyses?

Answer: Yes, ANGSD-wrapper requires the input of inbreeding coefficients in most analyses. If you prefer not to incorporate inbreeding, create a file with 0 on a new line for each individual in your population(s), i.e., assume complete outcrossing, and link to this file in the required places.

Question: ANGSD complained about not being able to find a chromosome even though it's in my reference sequence.

Answer: This issue is caused by out of order chromosomes in your reference. You can fix this in the following way:

for i in `seq 1 10`; do samtools faidx ref.fa $i >> sorted_ref.fa; done

Where 10 is the number of chromosomes. Then gzip your fasta and run

samtools faidx sorted_ref.fa.gz

to get your .fai file.