进行Blast比对,用参数-m 6 可以以列表的方式输出结果,结果中从左到右每一列的意义分别是: [00] Query id [01] Subject id [02] % identity [03] alignment length [04] mismatches [05] gap openings [06] q. start [07] q. end [08] s. start [09] s. end [10] e-value [11] bit score…
用bioperl 解析blast的默认输出结果, 整理成-m8格式的输出 #!/usr/bin/perl use Bio::SearchIO; my ($blast) = @ARGV; my $searchio = new Bio::SearchIO(-format => "blast", -file => "test.bls"); while (my $result = $searchio->next_result) { while (my $h…
D - Blast the Enemy! Time Limit:3000MS Memory Limit:0KB 64bit IO Format:%lld & %llu Submit Status Practice UVALive 4426 Description A new computer game has just arrived and as an active and always-in-the-scene player, you should finish it bef…
from: https://secure.clcbio.com/helpspot/index.php?pg=kb.page&id=113 Can I run a local BLAST search again multiple blast databases simultaneously? (blast, blastn, blastp, blastx, serveral databases)By default, you will be offered single BLAST databas…
以前有的是非完整时间写的博客,抽时间需要统一整理一下. 今天在重新装repeatmasker. 整个过程是这样的,有关联的事情有两个. 1. 装repeatmasker需要各种Prerequisites,其中就可能用到了blast,而之前一直找这个版本的blast,在ncbi硬是没有找到: For RMBlast ( NCBI Blast modified for use with RepeatMasker/RepeatModeler ) please go to our download pa…