bioperl 格式化genebank的输出
代码如下:
use Bio::SeqIO; use Bio::DB::GenBank; use Bio::DB::Query::GenBank;my $db_obj = Bio::DB::GenBank->new; my $seq_obj = $db_obj->get_Seq_by_acc('JN093905'); my $id = $seq_obj->display_id(); my $organ; foreach my $feat_object ($seq_obj->get_SeqFeatures) {if ($feat_object->primary_tag eq "source") {($organ) = $feat_object->get_tag_values('organism');}if ($feat_object->primary_tag eq "CDS") {my ($gene_name) = $feat_object->get_tag_values('gene');next if $gene_name ne 'nifH';my ($seq) = $feat_object->get_tag_values('translation');my ($pro) = $feat_object->get_tag_values('product');$seq = lc($seq);my $len = (length($seq) + 1) * 3;print qq{>$id coded_by=<1..>$len,organism=$organ,definition=$pro\n$seq\n};} }
运行结果如下:
>JN093905 coded_by=<1..>330,organism=uncultured Trichodesmium sp.,definition=dinitrogenase reductase rlilnakaqttvlhvaaergavedveldevlkpgfggikcvesggpepgvgcagrgiitainfleeegaytdldfvsydvlgdvvcggfampirenkaqeiyivcsgem
bioperl 格式化genebank的输出相关推荐
- html+注释格式化,使用xml注释来生成格式化的html输出
我试图从我在xml文件中的注释中生成一个格式良好的html文档.目前我有一个xml文件,用于生成xml表格的html列表.为了让我添加有关表格的评论,我手动将注释添加到输出html文件中.使用xml注 ...
- (03)格式化输入和输出
一.概念简介 基础学习阶段:键盘输入,屏幕输出. 从键盘输入数据,以及实现简单的人机交互. 通过格式化的方式输出 整数.浮点数.字符.字符串,以及对数据进行对齐. 1.标准输入输出 在C语言中,有三个 ...
- 文件格式化输入和输出
随时随地技术实战干货,获取项目源码.学习资料,请关注源代码社区公众号(ydmsq666) 在控制台操作时,使用的格式化输入和输出为scanf和printf,那么对文件的IO操作也可以使用fscanf和 ...
- GDB 格式化结构体输出
set print address set print address on 打开地址输出,当程序显示函数信息时,GDB会显出函数的参数地址.系统默认为打开的, show print address ...
- linux gdb 格式化结构体输出
set print address set print address on 打开地址输出,当程序显示函数信息时,GDB会显出函数的参数地址.系统默认为打开的, show print address ...
- PAT甲级1016 Phone Bills :[C++题解]字符串处理(复杂题)(C语言格式化读入、输出很便利!!!)
文章目录 题目分析 题目链接 题目分析 原题: 长途电话公司按以下规则向客户收费: 拨打长途电话每分钟要花费一定的费用,具体收费取决于拨打电话的时间. 客户开始拨打长途电话的时间将被记录,客户挂断电话 ...
- C++ Primer 5th笔记(chap 17 标准库特殊设施)未格式化的输入/输出操作
1. 格式化IO 输入和输出运算符(<< 和>>)根据读取或写入的数据类型来格式化它们. 输入运算符忽略空白符 输出应用补白 精度等规则操作 2. 未格式化 IO (unfor ...
- oracle 格式化报表输出,perl的格式化(Format)报表输出实现代码
perl有最好的文本数据处理能力.这是大家都知道的.在perl本身有一个别的软件没有的小功能,就是Perl格式.它相当于简单的命令行报表和图表输出.这个需要先声明格式.创建格式的文件句柄.然后由@,^ ...
- php 输出text格式化,php printf() 输出格式化的字符串,phpprintf
php printf() 输出格式化的字符串,phpprintf php printf() 函数用于输出格式化的字符串,本文章向码农介绍php printf()函数的使用方法和基本使用实例,感兴趣的码 ...
最新文章
- Qt实现 指针式时钟+动态时钟 (详细注释)
- WebMisDeveloper4.2.0面世
- 微信可能会颠覆快递的运营模式
- 黑客数字雨html单页,Hei客帝国数字雨.html
- 数据分析与挖掘实战-窃电漏电用户的发现
- BNUOJ 4215 最长公共连续子序列
- 三星为什么要卖芯片?手机干不过华为小米,半导体好挣钱!
- C++标准库函数begin和end函数
- java的equals方法_Java Duration类| 带示例的equals()方法
- c语言sort函数排序二维数组,js 二维数组排序sort()函数
- 作者:​冯景华(1984-),男,国家超级计算天津中心主任助理、系统管理部部长。...
- 说明多级指针与动态变量,静态变量关系的几个简单程序
- 八周二次课(5月14日)
- goaheadlinux移植_goahead
- 从谷歌浏览器复制不带样式_如何更新ETHERSCAN区块链浏览器上的代币信息
- python3、ipython3、setup-tools、pip等环境搭建详细总结
- 内部跳转(请求转发)和外部跳转(重定向)的区别?
- c++时间片轮转rr进程调度算法_进程,线程基础(—)
- 隐马尔可夫的知识要点
- Linux的 7个运行级别
热门文章
- directadmin php.ini 修改,如何快速修改DirectAdmin面板默认2222端口
- 提高开发效率之安卓模板(上面有四种模板的教程,我之前会两种,看完之后还是只会两种2333)
- QT5界面操作3: QPainter和paintEvent(self, event)的配合
- 宜出行热力图怎么抓取_滴滴听不到单怎么办?七个小技巧
- 3des java ecb_PHP 3DES加密 与JAVA通用 加密模式:ECB
- 目标检测:NMS和计算mAP时的置信度阈值和IoU阈值
- 天河二号 服务器 操作系统,天河二号云服务器
- java comet_用java实现comet,基于 HTTP长连接的实现,用于从服务端实时发送信息到客户端...
- php中怎么定义page,css @page的使用与定义详解
- 【c++】29.设计模式总结