3933 个回答
利用Perl遍历哈希
处理fasta序列最好用bioperl包，非常好用，下面是代码部分，其实你不用遍历hash，直接边读序列边输出就好了： 也不清楚你要的意图，下面的代码是把fasta序列存储到%keep中，并输出到普通文本；遍历hash可参考：https://www.omicsclass.com/question/99 use Bio::SeqIO;use Bio::Seq;$in = Bio::SeqIO->new(-file =...
回答于 2018-07-25 13:01
在perl中怎么遍历哈希
不明白你要干啥，你是有一些基因ID，想把对应的基因序列提取出来到一个新文件，还是只是把fasta序列中的ID提取出来？问问题要明确，不要让人歧义： perl遍历hash： while(($key, $value) = each(%HASH)) { # do something with $key and $value} #----------------------------- foreach $key (keys %HASH) { $...
回答于 2018-07-24 23:35
pip install RSeQC安装时候报错，是不是底层的库不全？怎么解决...
最后一个错误有问题，GCC编译出错，你的系统环境不清楚，建议Google一下错误原因。 有时候也贴一下错误文字，方便别人帮你搜索一下。
回答于 2018-07-24 12:08
请问eclipse中的remote system 怎么安装
安装插件可参照下图，然后搜索remote system： 更详细的安装可看视频《基于eclipse搭建编程环境 perl+python+R》
回答于 2018-07-24 11:07
get-fa-by-id.pl脚本文件不能被打开
你的脚本名字输入错了，当然没有那个文件了，no such file。 检查下文件名字对不对？我记得是get_fa_by_id.pl 是下划线； 问题很简单，仔细些，linux给你提示你要看一下，这种问题很容易解决；
回答于 2018-07-23 10:50
McScan，跑出无结果是什么问题？用McScanX还是同样没有结果 附图
从截图看不出问题，我这只能推测：你的gff文件有可能有问题，检查一下： 1.gff文件是否包含两个物种的所有基因位置信息； 2.看你的染色体是数字，如果两个物种的染色体都是数字，应加个前缀对染色体进行区分； H1，H2...... T1,T2...... 3.用下MCScanX吧，更方便些：《MCScanX安装和使用》
回答于 2018-07-20 13:33
怎么查找人GRCh36/即hg18的基因版本注释
可到UCSC这里去下载对应的基因组和注释信息： http://hgdownload.soe.ucsc.edu/downloads.html 要浏览不同版本的基因组注释信息，可以在UCSC的这个地方： 地址：http://genome-asia.ucsc.edu/cgi-bin/hgGateway?hgsid=682957885_HpOgAL9oAwgRsbQzV0aH1BTKPNO1&redirect=manual&source=genome.ucsc.edu 操作...
回答于 2018-07-19 19:11
如何评价系统进化树的构建的结果
可以适当删除一些比对较差的GAP区域，保留比对结果中比较保守的区域，再进行系统进化树构建。至于系统进化树的好坏，可以变换下参数，再做几次，看看有没有符合预期的；
回答于 2018-07-19 10:51