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.

Direct Ancestry

From angsd
Jump to navigation Jump to search
The printable version is no longer supported and may have rendering errors. Please update your browser bookmarks and please use the default browser print function instead.

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'