synmuv/sfs_invariant_sites

Calculate Site Frequency Spectra With and Without Invariant Sites

SFS_Invariant_Sites

Calculate Site Frequency Spectra With and Without Invariant Sites

All scripts used to generate the spectra are currently located in the repo, but more work needs to be done to make the scripts more extensible.

Current Workflow

1. Run GENERATE_SFS_FILES.sh. This takes a VCF and generates a processed data set that is used by Generate_Spectra.R to generate spectra files. This script needs to be updated in the following ways (currently doesn't run as a script, but all the commands are there):

• Modify paths to work in the github repository file structure.
• Modify script to take VCF, sample names file, and ancestor name as an input.

2. Run Invariant_SFS.R. This script takes a spliced CDS fasta file and the output from GENERATE_SFS_FILES.sh to generate counts of 0- and 4-fold sites that do are invariant across the tested population. This script needs to be updated to work within the repository file structure.
3. Run Generate_Spectra.R. This generates .sfs files used for DFE analysis. The parameters for this script are:

• The output from GENERATE_SFS_FILES.sh: SFS_INPUT.tsv
• no_indel or indel to include indels in the spectra
• The output from Invariant_SFS.R: invariant_site_by_region.tsv