M2 Lab4 - areed37/Andrew_Reed_EPP_531 GitHub Wiki

Scaffolding

mkdir Juicer
cd Juicer/

#Hap 1
mkdir hap1
cd hap1/

mkdir references
cd references/
ln -s /pickett_sphinx/projects/EPP531_AGA/reed_EPP531/Juicer/hap1/references/MT_CP_removal/Sassafras_V1.0_with_Hi-C_80X_045.hic.hap2.p_ctg.fasta.filtered.filtered .

spack load /be6uuqp
bwa index Sassafras_V1.0_with_Hi-C_80X_045.hic.hap2.p_ctg.fasta.filtered.filtered

cd ..

mkdir restriction_sites
cd restriction_sites/
python /pickett_shared/software/juicer-1.6/misc/generate_site_positions.py DpnII sassafrass_hap1.p_ctg ../references/Sassafras_V1.0_with_Hi-C_80X_045.hic.hap2.p_ctg.fasta.filtered.filtered
awk 'BEGIN{OFS="\t"}{print $1, $NF}' sassafrass_hap1.p_ctg_DpnII.txt > sassafrass_hap1.p_ctg_DpnII.chrom.sizes
cd ..

mkdir fastq
cd fastq/
ln -s /pwd/salbidum01_1334141_S3HiC_R1.fastq.gz .
ln -s /pwd/salbidum01_1334141_S3HiC_R2.fastq.gz .
cd ..

bash /pickett_shared/software/juicer-1.6/AWS/scripts/juicer.sh \
        -g sassafrass_hap1.p_ctg \
        -D /pickett_shared/software/juicer-1.6/AWS/ \
        -t 10 \
        -d /pwd/Juicer/hap1/ \
        -z ./references/sassafrass_hap1.p_ctg.fasta \
        -p ./restriction_sites/sassafrass_hap1.p_ctg_DpnII.chrom.sizes \
        -y ./restriction_sites/sassafrass_hap1.p_ctg_DpnII.txt