syri软件的安装

目录

1. 下载syri软件

​2. 其他软件的安装和环境配置

3. plotsr软件的安装


1. 下载syri软件

网址:syri 该软件基于 3.5 版本 python,请提前使用 conda 创建 python 3.5 版本 环境并安装依赖模块

点击 view project on github 

点击 +9 release

我选择的是syri v1.4版本

 2. 其他软件的安装和环境配置

 创建3.5版本的python环境

## 创建3.5版本python环境
conda create -n py35 python=3.5
## 进入环境
conda activate py35
## 安装syri 依赖的python模块
conda install cython numpy scipy pandas=0.23.4 biopython \
psutil matplotlib=3.0.0
conda install -c conda-forge python-igraph
conda install -c bioconda pysam

安装syri

## 下载1.4版本syri
git clone https://github.com/schneebergerlab/syri.git
## 或者直接wget下载
# wget https://github.com/schneebergerlab/syri/archive/refs/\

# tar -xzvf syri.v1.4.zip
# mv v1.4 syri
cd syri
python setup.py install
## 添加可执行权限
chmod 755 ./syri/bin/*

3. plotsr软件的安装

plotsr: visualizing structural similarities and rearrangements between multiple genomes

plotsr 软件是可视化基因组结构重拍的工具。

其下载网址:plotsr

##解压
tar -xzvf plotsr1.1.1.tar.gz
cd plotsr1.1.1

然后输入命令测序一下:

##进入bin目录
plotsr

出现帮助文档则表明安装成功。

注意:由于本人当时是在创建的syri环境python=3.5版本中,无论利用 conda install plotsr 或 下载后手动安装均无法成功,原因是在当前syri-python=3.5环境中未能与最新版本的plotsr相兼容。当切换至base=python.3.10版本时,上述两种方法均可正常使用。

  • 1
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
data_dir='/public/work/Personal/wuxu/qiantao_17' for file1 in ${data_dir}/*.fasta; do for file2 in ${data_dir}/*.fasta; do if [ "$file1" != "$file2" ]; then touch snp_indel.end.sh && cat snp_indel.end.sh && \ export PATH=/public/work/Personal/pangshuai/software/conda/miniconda3/bin/:${PATH} && \ nucmer --mum -t 8 -g 1000 -p ${file1##*/}.${file2##*/}.ref_based.nucmer $file1 $file2 && \ delta-filter -1 -l 200 ${file1##*/}.${file2##*/}.ref_based.nucmer.delta > ${file1##*/}.${file2##*/}.ref_based.nucmer.delta.filter && \ dnadiff -d ${file1##*/}.${file2##*/}.ref_based.nucmer.delta.filter -p ${file1##*/}.${file2##*/}.ref_based.nucmer && \ show-coords -rcloT ${file1##*/}.${file2##*/}.ref_based.nucmer.delta.filter > ${file1##*/}.${file2##*/}.ref_based.nucmer.delta.filter.coords && \ show-coords -THrd ${file1##*/}.${file2##*/}.ref_based.nucmer.delta.filter > ${file1##*/}.${file2##*/}.ref_based.nucmer.delta.filter.syri.coords && \ show-snps -ClrTH ${file1##*/}.${file2##*/}.ref_based.nucmer.delta.filter > ${file1##*/}.${file2##*/}.ref_based.nucmer.delta.filter.snp && \ show-diff ${file1##*/}.${file2##*/}.ref_based.nucmer.delta.filter > ${file1##*/}.${file2##*/}.ref_based.nucmer.delta.filter.inv && \ perl /public/work/Pipline/Structural_Variation/pipeline/2.1.1/bin/filter_the_MUmmer_SNP_file.pl ${file1##*/}.${file2##*/}.ref_based.nucmer.delta.filter.snp ${file1##*/}.${file2##*/}.ref_based.nucmer.delta.filter.snp.SNPs ${file1##*/}.${file2##*/}.ref_based.nucmer.delta.filter.snp.Insertions ${file1##*/}.${file2##*/}.ref_based.nucmer.delta.filter.snp.Deletions 10000000 && \ touch snp_indel.end.tmp && \ mv snp_indel.end.tmp snp_indel.end && \ sleep 10 fi done done ,增加一个判断,使/public/work/Personal/wuxu/qiantao_17路径下以.fasta结尾的文件两两一组不分前后只组合一次,然后再执行touch 后面的代码
06-03
好的,您可以将原始代码修改为以下内容: ```bash data_dir='/public/work/Personal/wuxu/qiantao_17' files=( "${data_dir}"/*.fasta ) for (( i=0; i<${#files[@]}-1; i++ )); do for (( j=i+1; j<${#files[@]}; j++ )); do file1=${files[i]} file2=${files[j]} touch snp_indel.end.sh && cat snp_indel.end.sh && \ export PATH=/public/work/Personal/pangshuai/software/conda/miniconda3/bin/:${PATH} && \ if [ "$file1" != "$file2" ]; then nucmer --mum -t 8 -g 1000 -p ${file1##*/}.${file2##*/}.ref_based.nucmer $file1 $file2 && \ delta-filter -1 -l 200 ${file1##*/}.${file2##*/}.ref_based.nucmer.delta > ${file1##*/}.${file2##*/}.ref_based.nucmer.delta.filter && \ dnadiff -d ${file1##*/}.${file2##*/}.ref_based.nucmer.delta.filter -p ${file1##*/}.${file2##*/}.ref_based.nucmer && \ show-coords -rcloT ${file1##*/}.${file2##*/}.ref_based.nucmer.delta.filter > ${file1##*/}.${file2##*/}.ref_based.nucmer.delta.filter.coords && \ show-coords -THrd ${file1##*/}.${file2##*/}.ref_based.nucmer.delta.filter > ${file1##*/}.${file2##*/}.ref_based.nucmer.delta.filter.syri.coords && \ show-snps -ClrTH ${file1##*/}.${file2##*/}.ref_based.nucmer.delta.filter > ${file1##*/}.${file2##*/}.ref_based.nucmer.delta.filter.snp && \ show-diff ${file1##*/}.${file2##*/}.ref_based.nucmer.delta.filter > ${file1##*/}.${file2##*/}.ref_based.nucmer.delta.filter.inv && \ perl /public/work/Pipline/Structural_Variation/pipeline/2.1.1/bin/filter_the_MUmmer_SNP_file.pl ${file1##*/}.${file2##*/}.ref_based.nucmer.delta.filter.snp ${file1##*/}.${file2##*/}.ref_based.nucmer.delta.filter.snp.SNPs ${file1##*/}.${file2##*/}.ref_based.nucmer.delta.filter.snp.Insertions ${file1##*/}.${file2##*/}.ref_based.nucmer.delta.filter.snp.Deletions 10000000 && \ touch snp_indel.end.tmp && \ mv snp_indel.end.tmp snp_indel.end && \ sleep 10 fi done done ``` 这里使用了一个数组 `files` 存储了所有以 `.fasta` 结尾的文件路径,然后使用两个 `for` 循环对其进行两两组合,并且加上了判断条件,使得每个文件只会和其他文件组合一次。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值