Tips:

# 仅供 IEMB-1 用户参考。

# 进入登陆服务器页面后,所有命令可直接使用复制粘贴使用。

# 已通过基因组组装软件(如SPAdes或Megahit等)组装获得contigs水平的基因组。

# 不要修改脚本生成文件的文件名,可能会导致脚本无法识别文件的问题。

1. 使用MEME搜寻基因组端粒结构

1
2
3
4
5
6
7
8
cd /your genome path/

# 以`contigs.fasta`文件,搜寻该基因组的端粒结构
# 运行脚本
/apps/users/andrew/littletools/find_telomere.sh contigs.fasta

# 根据基因组大小不同,脚本运行时间从几分钟到几十分钟不等
# 结果保存在`5prime_telomere`和`3prime_telomere`文件夹中,两个文件夹中的motif应该互为反向互补序列,即为该种的端粒结构