- #!/bin/sh
- #PBS -N seqERRor
- #PBS -j oe
- #PBS -l ncpus=20
- #PBS -l nodes=1
- #PBS -l mem=40G
- source /home/liuxiangqiong/miniconda3/bin/activate base
- target=/cgdata/liuxiangqiong/work62pancancer/pipeline/v0/refdata/NanOnco_Plus_Panel_v2.0_Covered_b37_cg.parY2X.sort.bed
- ###采用fgbio计算测序错误率
- outputdir2=${inputpath}/8Fastqc
- fgbio=/cgdata/liuxiangqiong/software/fgbio/fgbio-1.3.0.jar
- ref=/cgdata/Database/GATK/b37/human_g1k_v37_decoy.fasta
- dbsnp138=/cgdata/Database/GATK/b37/dbsnp_138.b37.vcf
- java -jar $fgbio ErrorRateByReadPosition -i ${bam_dir}/${sample}_clean.bam --collapse=false -r $ref -v $dbsnp138 -o $outputdir2/$sample
|