//创建基础文件夹mkdir /work_dictionary/softwaremkdir /work_dictionary/scriptmkdir /work_dictionary/index//安装sratoolkit软件包wget -c https://ftp-trace.ncbi.nlm.nih.gov/sra/sdk/3.0.0/sratoolkit.3.0.0-centos_linux64.tar.gz//使用Xftp将circRNA分析软件传输到文件夹//下载并构建索引wget -c /pub/release-105/fasta/homo_sapiens/dna/Homo_sapiens.GRCh38.dna_sm.primary_assembly.fa.gzwget -c /pub/release-105/gtf/homo_sapiens/Homo_sapiens.GRCh38.105.gtf.gzbwa index -a bwtsw Homo_sapiens.GRCh38.dna_sm.primary_assembly.fa//创建工作项目文件夹mkdir /work_dictionarymkdir /work_dictionary/00.SRR_filesmkdir /work_dictionary/01.readsmkdir /work_dictionary/02.bwamkdir /work_dictionary/03.CIRI2mkdir /work_dictionary/04.CIRI_ASmkdir /work_dictionary/05.CIRI_pipelinemkdir /work_dictionary/06.CIRI_ro1mkdir /work_dictionary/07.bwa_mem_ro1mkdir /work_dictionary/08.CIRI_ro2mkdir /work_dictionary/09.CIRI_mergemkdir /work_dictionary/10.CIRI_vis//00.SRR_filescd mkdir /work_dictionary/00.SRR_filescp /work_dictionary/script/wget.sh /work_dictionary/00.SRR_files/wget.shqsub wget.sh//01.readscd input=/work_dictionary/00.SRR_filesls -lR |grep SRR|awk '{print $9}' > fastq_samples.listcp /work_dictionary/script/fastq_dump_sample.sh /work_dictionary/00.SRR_files/fastq_dump_sample.shsed -i 's/dictionary/xxx_circseq/g' fastq_dump_sample.shfor i in `cat fastq_samples.list`; do cp fastq_dump_sample.sh fastq_dump_${i}.sh; donefor i in `cat fastq_samples.list`; do sed -i 's/sample/'$i'/g' fastq_dump_${i}.sh; donefor i in `cat fastq_samples.list`; do qsub fastq_dump_${i}.sh; done//02.bwacd /work_dictionary/01.readsls -lR |grep _1.fastq.gz$|awk '{print $9}' > bwa_samples.list //$代表结尾sed -i 's/_1.fastq.gz//g' bwa_samples.list //用空白替换_1.fastq.gz这样的尾缀，等同于去除尾缀cp /work_dictionary/script/bwa_sample.sh /work_dictionary/01.reads/bwa_sample.sh //复制脚本到工作目录sed -i 's/dictionary/xxx_circseq/g' bwa_sample.sh //修改sample脚本中输入和输出的目录为当前工作目录for i in `cat bwa_samples.list`; do cp bwa_sample.sh bwa_${i}.sh; done //复制sample脚本并批量重命名脚本for i in `cat bwa_samples.list`; do sed -i 's/sample/'$i'/g' bwa_${i}.sh; done //遍历list中各行，替换sample脚本中的样本名称for i in `cat bwa_samples.list`; do qsub bwa_${i}.sh; done //批量提交作业//03.CIRI2cd /work_dictionary/02.bwals -lR |grep .sam$|awk '{print $9}' > CIRI2_samples.listsed -i 's/.sam//g' CIRI2_samples.listcp /work_dictionary/script/CIRI2_sample.sh /work_dictionary/02.bwa/CIRI2_sample.shsed -i 's/dictionary/xxx_circseq/g' CIRI2_sample.shfor i in `cat CIRI2_samples.list`; do cp CIRI2_sample.sh CIRI2_${i}.sh; donefor i in `cat CIRI2_samples.list`; do sed -i 's/sample/'$i'/g' CIRI2_${i}.sh; donefor i in `cat CIRI2_samples.list`; do qsub CIRI2_${i}.sh; done//04.CIRI_AScd /work_dictionary/03.CIRI2ls -lR |grep .ciri$|awk '{print $9}' > CIRI_AS_samples.listsed -i 's/.ciri//g' CIRI_AS_samples.listcp /work_dictionary/script/CIRI_AS_sample.sh /work_dictionary/03.CIRI2/CIRI_AS_sample.shsed -i 's/dictionary/xxx_circseq/g' CIRI_AS_sample.shfor i in `cat CIRI_AS_samples.list`; do cp CIRI_AS_sample.sh CIRI_AS_${i}.sh; donefor i in `cat CIRI_AS_samples.list`; do sed -i 's/sample/'$i'/g' CIRI_AS_${i}.sh; donefor i in `cat CIRI_AS_samples.list`; do qsub CIRI_AS_${i}.sh; done//05.CIRI_pipelinecd /work_dictionary/04.CIRI_ASls -lR |grep _library_length.list$|awk '{print $9}' > CIRI_pipeline_samples.listsed -i 's/_library_length.list//g' CIRI_pipeline_samples.listcp /work_dictionary/script/CIRI_pipeline_sample.sh /work_dictionary/04.CIRI_AS/CIRI_pipeline_sample.shsed -i 's/dictionary/xxx_circseq/g' CIRI_pipeline_sample.shfor i in `cat CIRI_pipeline_samples.list`; do cp CIRI_pipeline_sample.sh CIRI_pipeline_${i}.sh; donefor i in `cat CIRI_pipeline_samples.list`; do sed -i 's/sample/'$i'/g' CIRI_pipeline_${i}.sh; donefor i in `cat CIRI_pipeline_samples.list`; do qsub CIRI_pipeline_${i}.sh; done//06.CIRI_ro1cd /work_dictionary/05.CIRI_pipelinecp /work_dictionary/04.CIRI_AS/CIRI_pipeline_samples.list /work_dictionary/05.CIRI_pipeline/CIRI_ro1_samples.listcp /work_dictionary/script/CIRI_ro1_sample.sh /work_dictionary/05.CIRI_pipeline/CIRI_ro1_sample.shsed -i 's/dictionary/xxx_circseq/g' CIRI_ro1_samples.listfor i in `cat CIRI_ro1_samples.list`; do cp CIRI_ro1_sample.sh CIRI_ro1_${i}.sh; donefor i in `cat CIRI_ro1_samples.list`; do sed -i 's/sample/'$i'/g' CIRI_ro1_${i}.sh; donefor i in `cat CIRI_ro1_samples.list`; do qsub CIRI_ro1_${i}.sh; done//07.bwa_mem_ro1cd /work_dictionary/06.CIRI_ro1ls -lR |grep _ro1.fq$|awk '{print $9}' > bwa_mem_ro1_samples.listsed -i 's/_ro1.fq//g' bwa_mem_ro1_samples.listcp /work_dictionary/script/bwa_mem_ro1_sample.sh /work_dictionary/06.CIRI_ro1/bwa_mem_ro1_sample.shsed -i 's/dictionary/xxx_circseq/g' bwa_mem_ro1_sample.shfor i in `cat bwa_mem_ro1_samples.list`; do cp bwa_mem_ro1_sample.sh bwa_mem_ro1_${i}.sh; donefor i in `cat bwa_mem_ro1_samples.list`; do sed -i 's/sample/'$i'/g' bwa_mem_ro1_${i}.sh; donefor i in `cat bwa_mem_ro1_samples.list`; do qsub bwa_mem_ro1_${i}.sh; done//08.CIRI_ro2cd /work_dictionary/07.bwa_mem_ro1ls -lR |grep _ro1.sam$|awk '{print $9}' > CIRI_ro2_samples.listsed -i 's/_ro1.sam//g' CIRI_ro2_samples.listcp /work_dictionary/script/CIRI_ro2_sample.sh /work_dictionary/07.bwa_mem_ro1/CIRI_ro2_sample.shsed -i 's/dictionary/xxx_circseq/g' CIRI_ro2_sample.shfor i in `cat CIRI_ro2_samples.list`; do cp CIRI_ro2_sample.sh CIRI_ro2_${i}.sh; donefor i in `cat CIRI_ro2_samples.list`; do sed -i 's/sample/'$i'/g' CIRI_ro2_${i}.sh; donefor i in `cat CIRI_ro2_samples.list`; do qsub CIRI_ro2_${i}.sh; done//09.CIRI_mergecd /work_dictionary/08.CIRI_ro2ls -lR |grep _ro2_info.list$|awk '{print $9}' > CIRI_merge_samples.listsed -i 's/_ro2_info.list//g' CIRI_merge_samples.listcp /work_dictionary/script/CIRI_merge_sample.sh /work_dictionary/08.CIRI_ro2/CIRI_merge_sample.shsed -i 's/dictionary/xxx_circseq/g' CIRI_merge_sample.shfor i in `cat CIRI_merge_samples.list`; do cp CIRI_merge_sample.sh CIRI_merge_${i}.sh; donefor i in `cat CIRI_merge_samples.list`; do sed -i 's/sample/'$i'/g' CIRI_merge_${i}.sh; donefor i in `cat CIRI_merge_samples.list`; do qsub CIRI_merge_${i}.sh; done//10.CIRI_viscd /work_dictionary/09.CIRI_mergels -lR |grep _merge_circRNA_detail.anno$|awk '{print $9}' > CIRI_vis_samples.listsed -i 's/_merge_circRNA_detail.anno//g' CIRI_vis_samples.listcp /work_dictionary/script/CIRI_vis_sample.sh /work_dictionary/09.CIRI_merge/CIRI_vis_sample.shsed -i 's/dictionary/xxx_circseq/g' CIRI_vis_sample.shfor i in `cat CIRI_vis_samples.list`; do cp CIRI_vis_sample.sh CIRI_vis_${i}.sh; donefor i in `cat CIRI_vis_samples.list`; do sed -i 's/sample/'$i'/g' CIRI_vis_${i}.sh; donefor i in `cat CIRI_vis_samples.list`; do qsub CIRI_vis_${i}.sh; done
