Skip to content

When running ANGSD/realSFS, how can I exclude data from a specific chromosome, such as the X chromosome? #655

@kuangzhuoran

Description

@kuangzhuoran

My command is:
angsd -bam p1.bam.list -doSaf 1 -out p1 -anc anc.fa -GL 2 -P 4 -only_proper_pairs 1 -uniqueOnly 1 -remove_bads 1 -minQ 20 -minMapQ 3;
angsd -bam p2.bam.list -doSaf 1 -out p2 -anc anc.fa -GL 2 -P 4 -only_proper_pairs 1 -uniqueOnly 1 -remove_bads 1 -minQ 20 -minMapQ 3;
realSFS -P 30 p1.saf.idx p2.saf.idx > sfs;

It is certainly possible to remove the X chromosome data from the BAM files, but this would require additional time and storage. I would like to know whether ANGSD or realSFS provides an option to do this.

I have seen that realSFS allows using a specific chromosome with the -r chr1 option.

I tried removing the X chromosome sequence from the anc.fa file, but the software halts when it reaches the X chromosome. Unfortunately, the X chromosome is listed near the top in the chromosome order.

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions