基因组与基因组之间的比较（水稻）

一、准备工具

两个基因组文件： fs32.fa r498.fa (作为参考基因组)

minimap2的安装：

1.conda 直接安装

conda install -c bioconda minimap2

2. 下载安装

wget https://github.com/lh3/minimap2/releases/download/v2.17/minimap2-2.17_x64-linux.tar.bz2

tar -jxvf minimap2-2.17_x64-linux.tar.bz2

cd minimap2-2.17_x64-linux/

二、基因组比较

minimap2 -ax asm5 --eqx r498.fa fs32.fa >fs32.sam

最终得到sam格式文件

1.安装软件

conda install -c bioconda syri

或者：

cd /home/lsh/biosoft/

git clone https://github.com/schneebergerlab/syri.git

python setup.py install

chmod+x syri/bin/syri syri/bin/chroder syri/bin/plotsr

2.python3 /home/lsh/biosoft/syri-master/syri/bin/syri -c fs32.sam -r r498.fa -q fs32.fa -k -F S

最终得到结果文件,其中

syri.out 可以用来画共线性图

syri.vcf 包含了SNP

3.画图

python /home/lsh/biosoft/syri-master/syri/bin/plotsr syri.out r498.fa fs32.fa

©著作权归作者所有,转载或内容合作请联系作者
【社区内容提示】社区部分内容疑似由AI辅助生成，浏览时请结合常识与多方信息审慎甄别。
平台声明：文章内容（如有图片或视频亦包括在内）由作者上传并发布，文章内容仅代表作者本人观点，简书系信息发布平台，仅提供信息存储服务。