我使用的服务器,由于文件很多,在使用FOR循环的时候一直出错,老师可以帮忙看一下是哪里有语法错误吗?

/public/singularity exec /public/reseq.sif 

for i in FP180000114BR_L01_142 FP180000114BR_L01_143 FP180000114BR_L01_144 FP180000114BR_L01_146 FP180000114BR_L01_147 FP180000114BR_L01_148; do 

echo "RUN CMD: fastp --thread 30 --qualified_quality_phred 10 \

    --unqualified_percent_limit 50 \

    --n_base_limit 10 \

    -i $datadir/${i}_1.fq.gz \

    -I $datadir/${i}_2.fq.gz \

    -o ${i}_1.clean.fq.gz \

    -O ${i}_2.clean.fq.gz \

    --detect_adapter_for_pe \

    -h ${i}.html -j ${i}.json"


fastp --thread 30 --qualified_quality_phred 10 \

    --unqualified_percent_limit 50 \

    --n_base_limit 10 \

    -i $datadir/${i}_1.fq.gz \

    -I $datadir/${i}_2.fq.gz \

    -o ${i}_1.clean.fq.gz \

    -O ${i}_2.clean.fq.gz \

    --detect_adapter_for_pe \

    -h ${i}.html -j ${i}.json    

done

请先 登录 后评论

2 个回答

马杰宇

attachments-2023-12-skBokKaR65854ee718f39.png

请先 登录 后评论
omicsgene - 生物信息
擅长:重测序,遗传进化,转录组,GWAS

多行命令 建议写到,一个SH文件里面,运行

reseq.sif sh w.sh
请先 登录 后评论
  • 2 关注
  • 0 收藏,448 浏览
  • 马杰宇 提出于 2023-12-22 16:53

相似问题