Post date: Mar 25, 2016 9:55:10 PM
As a check, I want to see how similar our current genome from Dovetail is to what we had before. I am using mugsy to align the genomes (this is what I used for Lycaeides too). The submission script (below) and results are in (will be in) /uufs/chpc.utah.edu/common/home/u6000989/data/timema/tcrDovetail/MugsyComp.
#!/bin/sh
#SBATCH --time=72:00:00
#SBATCH --nodes=1
#SBATCH --ntasks=16
#SBATCH --account=gompert
#SBATCH --partition=kingspeak
#SBATCH --job-name=mugsy
#SBATCH --mail-type=FAIL
#SBATCH --mail-user=zach.gompert@usu.edu
echo ------------------------------------------------------
echo ------------------------------------------------------
echo SLURM: job identifier is $SLURM_JOBID
echo SLURM: job name is $SLURM_JOB_NAME
echo ------------------------------------------------------
bash
source /uufs/chpc.utah.edu/common/home/u6000989/source/mugsy_x86-64-v1r2.3/mugsyenv.sh
mkdir /scratch/general/lustre/mugsyTcr
cd /scratch/general/lustre/mugsyTcr/
/uufs/chpc.utah.edu/common/home/u6000989/source/mugsy_x86-64-v1r2.3/mugsy --directory /scratch/general/lustre/mugsyTcr --prefix mualnTcr /uufs/chpc.utah.edu/common/home/u6000989/data/timema/tcrDovetail/timema_26Feb2016_g9wYt.fasta /uufs/chpc.utah.edu/common/home/u6000989/data/timema/draft_genome/draft0.1/timemaGenome.fasta
rsync -avz /scratch/general/lustre/mugsyTcr /uufs/chpc.utah.edu/common/home/u6000989/data/timema/tcrDovetail/MugsyComp/