diff --git a/cross_species_comparison/scripts/lastz.sh b/cross_species_comparison/scripts/lastz.sh new file mode 100644 index 0000000000000000000000000000000000000000..6eed0813d7e26b4e0eaede38181161f3e14ca3ca --- /dev/null +++ b/cross_species_comparison/scripts/lastz.sh @@ -0,0 +1,17 @@ +##!/usr/bin/env bash + +## This script loops through all scaffolds in the dunnart genome and +## generates a separate lastz command for each scaffold +## These commands are then used in a slurm array script to run jobs in parallel + +TRA=($(for file in *.fa; do echo $file |cut -d "." -f 1;done)) + +echo ${TRA[@]} + +for tr in ${TRA[@]}; + +do + +echo 'lastz_32 /data/projects/punim0586/lecook/chipseq-pipeline/cross_species/data/genomes/mm10.fa[multiple] /data/projects/pun$ + +done diff --git a/cross_species_comparison/scripts/repeatMasker.sh b/cross_species_comparison/scripts/repeatMasker.sh new file mode 100644 index 0000000000000000000000000000000000000000..fc918a4e893c3bb793446ad89dc9eef72bfdf428 --- /dev/null +++ b/cross_species_comparison/scripts/repeatMasker.sh @@ -0,0 +1,18 @@ +##!/usr/bin/env bash + +## This script loops through all scaffolds in the dunnart genome and +## generates a separate RepeatMasker command for each scaffold +## These commands are then used in a slurm array script to run jobs in parallel + +TRA=($(for file in *.fa; do echo $file |cut -d "." -f 1;done)) + +echo ${TRA[@]} + +for tr in ${TRA[@]}; + +do + +echo 'RepeatMasker -q -xsmall smiCra1/'${tr}'.fa -default_search_engine hmmer -trf_prgm /home/lecook/.conda/envs/wga/bin/trf -hmmer_dir /home/lecook/.conda/envs/wga/bin/' + + +done