24小时热门版块排行榜    

查看: 2544  |  回复: 15
本帖产生 1 个 BioEPI ,点击这里进行查看
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

我要学perl么

新虫 (初入文坛)

[求助] 【求助】如何得到某一物种的所有蛋白翻译基因的不同位置(exon、intron)的GC含量

如题:
如何得到某一物种的所有蛋白翻译基因的不同位置(exon、intron)的GC含量,如下图



是要从网上下染色体的全部序列信息然后再自己计算吗?
【求助】如何得到某一物种的所有蛋白翻译基因的不同位置(exon、intron)的GC含量
图片1.png



[ Last edited by 我要学perl么 on 2013-8-14 at 22:30 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

我要学perl么

新虫 (初入文坛)

引用回帖:
7楼: Originally posted by wizardfan at 2013-08-19 05:06:58
random的那个要看paper里是怎么写的了
exon以外的不就是intron吗?
在下拉菜单里有exons的选项吧。...

你好!谢谢,Galaxy是非常好用的工具。
蛋白编码基因的内含子、外显子和上下游的GC含量我会计算了。但是对于lncRNA仍有点问题。
现在我有lncRNA的gtf和gff两种文件,对于如何计算GC含量我有2个想法,不知道哪个对

1. 第一种就是讲gtf转化为bed,然后使用Extract Features的Gene BED To Exon/Intron/Codon BED expander工具得到包含内含子的bed文件,然后再转化为gff, 以便用Feature Sequence得到fasta序列。
但是fasta文件为空,提示是“susScr3
290 warnings, 1st is: Unable to fetch the sequence from '9162341' to '543' for chrom '1'. Skipped 290 invalid lines, 1st is #4, "1        bed2gff        CUFF.1911.1 9162341        9162884        0        -        . CUFF.1911.1;"”

2.第二种就是将gtf过滤得到transcrip和exon的两种文件,然后用Operate on Genomic Intervals里的Subtract工具可以得到内含子的文件,Get flanks 工具可以得到上下游序列。然后就可以计算了。

感觉第二种方法是对的。
By the way,在用第一种方法的时候如果transcription和exon的起始终止位置一样的话,Gene BED To Exon/Intron/Codon BED expander工具判断内含子的话会出现start 和end相反的情况(比如说正链的内含子应该为0的,可是出现start为1300 end 为300)。不知道这样算不算错误?

写的有点乱,可能表达不清楚,见谅!
谢谢你的回复!
8楼2013-08-20 01:07:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 16 个回答

wizardfan

至尊木虫 (著名写手)

优秀版主

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
silicare: 金币+2, 辛苦了 2013-08-21 17:51:30
自己计算是最方便最自由的
如果是人等的注释很好的基因组,可以考虑用galaxy
2楼2013-08-15 08:08:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

我要学perl么

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by wizardfan at 2013-08-15 08:08:24
自己计算是最方便最自由的
如果是人等的注释很好的基因组,可以考虑用galaxy

我做的 是猪的。
galaxy可以用吗?
如果要自己计算的话是要下载全基因组序列,然后自己写perl代码?

另外有个小问题:就是计算GC含量的时候假如遇到这样的情况,ACGN,GC含量算50%吗?
3楼2013-08-15 08:42:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wizardfan

至尊木虫 (著名写手)

优秀版主

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
我要学perl么: 金币+5, ★★★很有帮助 2013-08-16 16:55:05
silicare: 金币+5, BioEPI+1 2013-08-21 17:51:44
引用回帖:
3楼: Originally posted by 我要学perl么 at 2013-08-15 08:42:36
我做的 是猪的。
galaxy可以用吗?
如果要自己计算的话是要下载全基因组序列,然后自己写perl代码?

另外有个小问题:就是计算GC含量的时候假如遇到这样的情况,ACGN,GC含量算50%吗?...

I think so
Go to website https://main.g2.bx.psu.edu/
click "get data" in the left panel, choose UCSC main
Select pig in the genome
Search for geecee to calculate gc content

If you are not sure, read the tutorials, particularly this one https://main.g2.bx.psu.edu/u/aun1/p/galaxy101
4楼2013-08-15 17:33:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 可靠消息,8月21号(明天)放榜 +11 医学老男孩 2026-08-20 13/650 2026-08-20 10:50 by anjeeshine
[基金申请] 时间戳变了,能看出什么问题? +13 基诺咪客 2026-08-17 17/850 2026-08-20 10:43 by wangyue2015
[基金申请] 只有每年这种时候来逛逛小木虫 +21 yaoyewhu2008 2026-08-20 21/1050 2026-08-20 10:37 by Vivilian
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +15 Ldrop2023 2026-08-13 19/950 2026-08-20 10:30 by 陈晨晨陈啊
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +55 医学老男孩 2026-08-13 125/6250 2026-08-20 09:02 by 医学老男孩
[基金申请] 今天基金会出结果吗?20260819 +14 kkkl_v 2026-08-19 15/750 2026-08-20 08:21 by sunzitan
[基金申请] filecode,4个jtjc了 +11 ziyangfang 2026-08-19 13/650 2026-08-19 21:28 by aasahr
[教师之家] 为什么余额宝的年化利率越来越低?主要原因有哪些? +5 瞬息宇宙 2026-08-15 5/250 2026-08-19 20:23 by super2002521
[基金申请] 今天放榜没戏了吧 +4 yuleib84 2026-08-19 5/250 2026-08-19 19:52 by hhs666
[基金申请] 今天放榜吗? +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] filecode=后面第一个是大写字母 +10 wangze12014 2026-08-14 12/600 2026-08-19 16:56 by 苦难博士
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 什么时候开奖? +6 CrisMessi 2026-08-18 6/300 2026-08-19 08:06 by Equinoxhua
[基金申请] 重要消息,中午系统在维护 +10 yuleib84 2026-08-18 11/550 2026-08-19 07:30 by 羽毛枫f
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
信息提示
请填处理意见