Skip to content

Scripts related to wild guppies expansion load estimation.

Notifications You must be signed in to change notification settings

0-Ioniel-0/Expantion_Load

Folders and files

NameName
Last commit message
Last commit date

Latest commit

 

History

12 Commits
 
 
 
 

Repository files navigation

Expansion Load

The repository name is misspelled because, this, unfortunately, is how it appears in my PhD thesis.

Scripts related to wild guppies genetic load estimation.

This repository contains scripts regarding:

  1. reads trimming
  2. reads alignment
  3. bam files manipulation
  4. variant calling
  5. filtering by quality, ancestral alleles and conservation score
  6. calculating genetic load
  7. estimating effective populations size and more...

!!! For details on the workflow of this project, please look into flowchart.png figure. !!!

Note that most of the scripts need to be edited to change paths to files and programs. DP thresholds, number of individuals, individuals' positions in the VCF files etc. should also be modified accordingly.

These scripts were run with:

Python (3.9)
FastQC (0.11.8)
Trimmomatic (0.39)
bwa (0.7.10)
samtools (1.6.0)
bcftools (1.9)
gatk (4.1.4.1)
vcftools (0.1.14)
R(4.3.3)
pixy(1.2.7.beta1)
admixture(1.3.0)
SNPeff(5.0)
plink(1.90)

For help or comments, contact: k.burda[a]amu.edu.pl

About

Scripts related to wild guppies expansion load estimation.

Resources

Stars

Watchers

Forks

Releases

No releases published

Packages

No packages published