MUMmer 项目教程

MUMmer 项目教程

mummerMummer alignment tool项目地址:https://gitcode.com/gh_mirrors/mu/mummer

1. 项目的目录结构及介绍

MUMmer 项目的目录结构如下:

mummer/
├── bin/
├── docs/
├── src/
├── test/
├── README.md
├── LICENSE
└── INSTALL
  • bin/:包含编译后的可执行文件。
  • docs/: 包含项目的文档文件。
  • src/: 包含项目的源代码文件。
  • test/: 包含测试脚本和测试数据。
  • README.md: 项目的基本介绍和使用说明。
  • LICENSE: 项目的许可证文件。
  • INSTALL: 项目的安装指南。

2. 项目的启动文件介绍

MUMmer 项目的主要启动文件位于 bin/ 目录下,包括以下几个关键的可执行文件:

  • mummer: 用于快速比对 DNA 序列的主程序。
  • nucmer: 用于比对核苷酸序列的工具。
  • promer: 用于比对蛋白质编码序列的工具。

这些文件在编译安装后生成,用户可以通过命令行直接调用这些可执行文件来启动相应的比对任务。

3. 项目的配置文件介绍

MUMmer 项目没有传统的配置文件,其运行参数主要通过命令行参数进行设置。用户在启动 mummernucmerpromer 时,可以通过不同的命令行选项来调整比对的具体参数,例如比对的长度、线程数等。

例如,使用 nucmer 进行比对时,可以这样设置参数:

nucmer --maxmatch -p output reference_genome query_genome

其中,--maxmatch 是一个参数选项,-p output 指定输出文件的前缀,reference_genomequery_genome 分别是参考基因组和查询基因组的路径。

通过这些命令行参数,用户可以灵活地配置和调整 MUMmer 的比对行为。

mummerMummer alignment tool项目地址:https://gitcode.com/gh_mirrors/mu/mummer

翻译代码############# new the scripts for alignments ,change format and show the alignments ############ if (($MappingSoft eq "mummer") or ($MappingSoft eq "nucmer")) { #mummer-4.0.0/bin/nucmer --mum --mincluster 500 -t 30 Ref.AAfa RefBB.fa -p OUT #mummer-4.0.0/bin/delta-filter -1 -i 90 -l 2000 OUT.delta > OUT.filter1.delta #mummer-4.0.0/bin/show-coords -c -r OUT.filter1.delta > OUT.filter1.coords if ($MappingPara eq "") {$MappingPara = "--mum --mincluster 500 ";} open (OUTSH,">$OutPrefix.mapping.sh") || die "input file can't open $!"; print OUTSH "$nucmer $MappingPara -t $NumThreads $OutPrefix.A.fa $OutPrefix.B.fa -p $OutPrefix \n"; print OUTSH "$deltaFilter -1 -i 90 -l $MinAlnLen $OutPrefix.delta > $OutPrefix.filter.delta \n"; print OUTSH "$showcoords -c -r $OutPrefix.filter.delta > $OutPrefix.filter.coords\n"; print OUTSH "perl $0 Coords2Link $OutPrefix.filter.coords $MinAlnLen $OutPrefix.link \n"; print OUTSH "$NGenomeSyn -InConf $OutPrefix.conf -OutPut $OutPrefix.svg \n"; close OUTSH; system ("sh $OutPrefix.mapping.sh "); } else { if ($MappingPara eq "") {$MappingPara = " -x asm5 "; } open (OUTSH,">$OutPrefix.mapping.sh") || die "input file can't open $!"; print OUTSH "$minimap2 $MappingPara -t $NumThreads $OutPrefix.B.fa $OutPrefix.A.fa > $OutPrefix.paf \n"; print OUTSH "perl $0 Paf2Link $OutPrefix.paf $MinAlnLen $OutPrefix.link \n"; print OUTSH "$NGenomeSyn -InConf $OutPrefix.conf -OutPut $OutPrefix.svg \n"; close OUTSH ; system ("sh $OutPrefix.mapping.sh "); } print "\tALL done, see the xxx.png . you can optimized drawing by [NGenomeSyn] software\n"; print "\t optimized: [Filter] and [Merge] small syn blocks to big syn block\n\n";
06-09
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

罗琰锴

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值