发文章
发文工具
撰写
网文摘手
文档
视频
思维导图
随笔
相册
原创同步助手
其他工具
图片转文字
文件清理
留言交流
给定一个FASTA文件,他以多行形式存放数据,如下所示
请根据每个序列的ID中的染色体编号拆分数据,例如上面两个数据,需要新建两个文件,分别是chr7和chr8,将这两条序列分别输出到这两个文件中。
不限制语言,测试数据可以从网盘中下载。
链接:https://share./5ecURT0 密码:pzsiab
我调用kseq.h文件,编写了一个程序,代码如下
运行速度如下
来自: 生物_医药_科研 > 《待分类》
0条评论
发表
请遵守用户 评论公约
使用SOAPfuse进行融合基因的分析
使用SOAPfuse进行融合基因的分析。和其他软件类似,SOAPfuse也需要建立物种对应的数据库,在软件安装的source目录下,提供了建立数据库...
我的Python笔记·生信编程直播第二题
= ''''''''>'''''''': seq = seq + line.strip() elif line[0] == ''''''''>'...
HLA-VBSeq:对全基因组数据进行HLA分型
HLA-VBSeq:对全基因组数据进行HLA分型。将比对到HLA区域的reads和没比对上参考基因组的reads合并,命令如下cat R1.fq unmapped_R1.fq &g...
blast简介及格式解读及练习题
blast简介及格式解读及练习题。仅仅对对延伸匹配进行连接的区域(限制性区域),而不是整个矩阵,是blast 相对于其他算法速度提高的关键,是以牺牲对角线带以外的任何匹配信息为代价,因此并不能确保que...
序列提取 - 《-零命令行-生信下游数据分析》- 第一弹
序列提取 - 《-零命令行-生信下游数据分析》- 第一弹写在前面。序列文件,Fasta格式(任何Fasta格式的序列文件,如genome.fa, unigenes....
在文件后添加数据
在文件后添加数据Private Sub Command1_Click()Open "C:\2.txt" For Append As #1Write #1, Chr(13) &Chr(10) &"Hello World"Close #1End Sub.
基于gff3/gtf文件-批量提取启动子或CDS序列-任何人都可以
基于gff3/gtf文件-批量提取启动子或CDS序列-任何人都可以写在前面。这篇推文应是 《-零命令行-生信下游数据分析》的第五贴,主题定为基...
GEO的数据注释文件没有基因名肿么破?
GEO的数据注释文件没有基因名肿么破?我们在处理GEO芯片数据的时候,经常会碰到芯片的数据的注释文件没有提供基因名,只有基因的序列。...
生信分析人员数据处理脚本实战 | 生信菜鸟团
微信扫码,在手机上查看选中内容