123456789101112131415161718192021222324252627 |
- #!/bin/sh
- #PBS -N align_svcalling_noUMI
- #PBS -j oe
- #PBS -l ncpus=20
- #PBS -l nodes=1
- #PBS -l mem=20G
- #inputpath=$1
- #sample=$2
- source /home/liuxiangqiong/miniconda3/bin/activate base
- hg19=/cgdata/Database/GATK/b37/human_g1k_v37_decoy.fasta
- target=/cgdata/liuxiangqiong/work62pancancer/Client/v0/script/refdata/NanOnco_Plus_Panel_v2.0_Covered_b37_cg.parY2X.sort.bed
- #################################################################################1.trim#########################################################
- echo "#--------------------------------------STEP1.clean the raw data--------------------------------------#"
- fastp -i ${inputpath}/${sample}_R1.fastq.gz -I ${inputpath}/${sample}_R2.fastq.gz -o ${inputpath}/${sample}_clean_R1.fq.gz -O ${inputpath}/${sample}_clean_R2.fq.gz -c -l 50 -q 20 --umi --umi_loc read1 --umi_len 8 -j ${inputpath}/${sample}.fastp.json -h ${inputpath}/${sample}.fastp.html
- echo "#--------------------------------------STEP1.clean the raw data finished!---------------#"
- #################################################################################2.mapping#########################################################
- #2.1.speedseq align
- echo "#--------------------------------------STEP2.align to refdata--------------------------------------#"
- speedseq align -t 40 -o ${bam_dir}/${sample}_clean -R "@RG\tID:${sample}\tSM:${sample}\tLB:${sample}" ${hg19} ${inputpath}/${sample}_clean_R1.fq.gz ${inputpath}/${sample}_clean_R2.fq.gz
- echo "#--------------------------------------STEP2.align to refdata finished!--------------------------------------#"
- rm -rf ${bam_dir}/${sample}_clean.splitters.bam
- rm -rf ${bam_dir}/${sample}_clean.splitters.bam.bai
- rm -rf ${bam_dir}/${sample}_clean.discordants.bam
- rm -rf ${bam_dir}/${sample}_clean.discordants.bam.bai
|