ANGSD: Analysis of next generation Sequencing Data

Latest tar.gz version is (0.934/0.935 on github), see Change_log for changes, and download it here.

Direct Ancestry

From angsd
Jump to: navigation, search

Based on the 2dsfs for two individuals it is possible to estimate the direct ancestry

angsd -i ind1.bam -gl 1 -anc ancestral.fasta -dosaf 1 -out ind1
angsd -i ind2.bam -gl 1 -anc ancestral.fasta -dosaf 1 -out ind2
realSFS ind1.saf.idx ind2.saf.idx -P 32 >2dsfs
Rscript direct_ancestry.R sfsfile='2dsfs'