如何利用提供的位置信息提取CDS？求思路

justinlife 2016-11-11 03:27:07

RT
genebank文件中有CDS的具体位置
如：
CDS join(5059..5101,34957..35078,43500..43702,48107..48222,
52895..53024,55347..55474,55795..55934,56032..56185,
57126..57294,68281..68384,69526..69634,69896..70074,
71162..71268,71355..71487,71578..71682,77281..77391,
78328..78481,80626..80708,82092..82219,86980..87102,
89868..89979,90313..90412,91012..91117,101395..101517)
下面有整条DNA的序列1-104000+
请问有什么方法可以提取这个CDS的位置信息，然后获取对应的CDS序列？求大牛解答

...全文

425 1 打赏收藏转发到动态举报

写回复

1 条回复

切换为时间正序

请发表友善的回复…

发表回复

屎克螂 2016-11-15

打赏
举报

回复

你是不是发错了地方，我看着是一本无字天书

Fiori Elements是SAP Fiori应用程序开发的一种开发模式，它旨在提供一种快速、简单和一致的方式来创建SAP Fiori应用程序。它基于元数据驱动的开发模式，通过使用预定义的UI元素和模板，可以快速构建具有一致性和可...

最近在上生物信息学原理，打算记录一些课上的作业。第一次作业：如题。基本思路：　1.从GFF中读取CDS的起始终止位置以及正...利用起始/终止位置等信息从FNA文件中提取CDS序列。FNA格式见 http://boyun.sh.cn/bi...

gff文件当中存储了基因组当中所有基因的注释信息，如果想得到基因组当中所有基因的位置信息可以利用awk命令批量的提取，命令如下：$ grep -v '#' Arabidopsis_thaliana.TAIR10.41.gff3|awk -F "[\t=:;]" 'BEGIN{OFS=...

今天分享一个实用的数据提取思路：从VCF文件中提取指定列，比如有1000个样本，需要从中提取出200个制定样本，形成新的vcf文件。需要在Linux环境下进行操作，主要用到perl脚本和shell脚本。原始数据原始文件vcf ...

使用方法参考: 基于全基因组的基因...思路——首先根据已经获得的ID号从gff文件中获取染色体位置信息，然后再用bedtools工具根据得到的染色体位置信息来获取基因的序列，最终得到基因集。代码如下。 #!/bin/bash c=

37,719

社区成员

34,238

社区内容

发帖

与我相关

我的任务

社区管理员

加入社区

近7日
近30日
至今

加载中

查看更多榜单

试试用AI创作助手写篇文章吧

+ 用AI写文章