官方github网址:https://github.com/AmpliconSuite/AmpliconSuite-pipeline
申请mosek个人免费学术许可:https://www.mosek.com/products/academic-licenses/
申请后放在~/mosek/下
安装流程:
conda create -n ampsuite && conda activate ampsuite
conda install -c bioconda -c conda-forge ampliconsuite
conda install -c mosek mosek
# then run the installer script to finalize the locations of the data repo and mosek license
wget https://raw.githubusercontent.com/AmpliconSuite/AmpliconSuite-pipeline/master/install.sh
chmod +x install.sh
./install.sh --finalize_only # -h to see options
额外自行下载数据库: https://refs.ampliconrepository.org/?prefix=data/module_support_files/AmpliconArchitect/
已经经过bwa比对,下载没有index的GRCh38.tar.gz即可
并行脚本:
AA_DATA_REPO=~/database/AmpliconArchitect_repo/
perl -ne 'chomp; next if /^$/; @a = split /\t/; print "AmpliconSuite-pipeline.py -s $a[0] -t 4 --normal_bam ../align/$a[1]_bqsr.bam --bam ../align/$a[0]_bqsr.bam --run_AA --run_AC > log/$a[0].log 2>&1 && echo $a[0] AmpliconSuite ok\n";' ../mutect2/sample_pair.txt > RunAmpliconSuite.sh
# 输出示例:AmpliconSuite-pipeline.py -s FETB06-BLPT-E -t 4 --normal_bam ../align/FETB06-N_bqsr.bam --bam ../align/FETB06-BLPT-E_bqsr.bam --run_AA --run_AC > log/FETB06-BLPT-E.log 2>&1 && echo FETB06-BLPT-E AmpliconSuite ok
nohup bash -c "cat RunAmpliconSuite.sh | parallel -j 4" > log/AmpliconSuite.log 2>&1 &