24小时热门版块排行榜    

Znn3bq.jpeg
汕头大学海洋科学接受调剂
查看: 2399  |  回复: 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的回帖

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的回帖

我要学perl么

新虫 (初入文坛)

引用回帖:
4楼: Originally posted by wizardfan at 2013-08-15 17:33:44
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 a ...


Thank you very much!
5楼2013-08-16 16:54:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

我要学perl么

新虫 (初入文坛)

引用回帖:
4楼: Originally posted by wizardfan at 2013-08-15 17:33:44
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 a ...

你好,我还想再问2个问题。
1. 我有包含lncRNA的exon的位置信息的GFF文件,怎么样才能得到他们的内含子及上下游的信息?另外在USCS main里position选项的define regions里它提示There is a limit of 1,000 defined regions. 是不是意思就是说不能超过1000条的lncRNA的exon?

2. 在1楼的那张图里,怎么样才能得到random的位置呢?

非常谢谢你的解答!
6楼2013-08-18 23:20:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wizardfan

至尊木虫 (著名写手)

优秀版主

引用回帖:
6楼: Originally posted by 我要学perl么 at 2013-08-18 23:20:39
你好,我还想再问2个问题。
1. 我有包含lncRNA的exon的位置信息的GFF文件,怎么样才能得到他们的内含子及上下游的信息?另外在USCS main里position选项的define regions里它提示There is a limit of 1,000 define ...

random的那个要看paper里是怎么写的了
exon以外的不就是intron吗?
在下拉菜单里有exons的选项吧。
7楼2013-08-19 05:06:58
已阅   回复此楼   关注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的回帖

wizardfan

至尊木虫 (著名写手)

优秀版主

引用回帖:
8楼: Originally posted by 我要学perl么 at 2013-08-20 01:07:22
你好!谢谢,Galaxy是非常好用的工具。
蛋白编码基因的内含子、外显子和上下游的GC含量我会计算了。但是对于lncRNA仍有点问题。
现在我有lncRNA的gtf和gff两种文件,对于如何计算GC含量我有2个想法,不知道哪个对 ...

我不是GALAXY专家,推荐你加入官方邮件列表参与讨论
9楼2013-08-20 07:24:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

我要学perl么

新虫 (初入文坛)

引用回帖:
9楼: Originally posted by wizardfan at 2013-08-20 07:24:54
我不是GALAXY专家,推荐你加入官方邮件列表参与讨论...

嗯。
谢谢你的帮助!
10楼2013-08-20 09:52:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 我要学perl么 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料相关专业344求调剂双非工科学校或课题组 +17 hualkop 2026-04-12 18/900 2026-04-13 22:42 by pies112
[考研] 271求调剂 +24 2261744733 2026-04-11 24/1200 2026-04-13 22:24 by pies112
[考研] 化工学硕294分,求导师收留 +25 yzyzx 2026-04-12 28/1400 2026-04-13 22:16 by zhq0425
[考研] 生物医药调剂|SCI中科院三区一作+多项科研成果 +9 likangxing 2026-04-07 12/600 2026-04-13 21:35 by 崆峒药童
[考研] 266调剂 +10 daya sun 2026-04-07 11/550 2026-04-13 10:12 by fenglj492
[考研] 一志愿0807 数一英一 313 有没有二轮调剂 +12 emokidd 2026-04-08 13/650 2026-04-13 08:32 by lhj2009
[考研] 322求调剂,08工科 +4 今天是个小号 2026-04-08 4/200 2026-04-13 00:20 by baobaoye
[考研] 291求调剂 +11 关忆北. 2026-04-09 12/600 2026-04-12 10:32 by 逆水乘风
[考研] 求调剂 +18 MAX怅惘 2026-04-09 20/1000 2026-04-11 23:31 by zhen~
[考研] 一志愿厦大0856,306求调剂 +15 Bblinging 2026-04-11 15/750 2026-04-11 22:53 by 314126402
[教师之家] 请问地理、遥感方面,可以做哪些横向项目啊,纵向完不成考核啊 +3 锦衣卫寒战 2026-04-07 5/250 2026-04-11 20:51 by 豫椒
[考研] 085410-273求调剂 +6 X1999 2026-04-10 6/300 2026-04-11 10:32 by Delta2012
[考研] 085402通信工程调剂,有4项学科竞赛国奖(电赛国二),硕士研究生调剂自荐信。 +5 m永o不v言o弃m 2026-04-09 5/250 2026-04-11 09:33 by zhq0425
[考研] 0854调剂 +4 长弓傲 2026-04-09 4/200 2026-04-11 09:18 by 猪会飞
[考研] 一志愿211,化学学硕,310分,本科重点双非,求调剂 +27 努力奋斗112 2026-04-07 30/1500 2026-04-10 15:06 by Kilig0317
[考研] 机械还有还有名额吗?太难了 +6 笑笑袁 2026-04-10 6/300 2026-04-10 11:54 by 高维春
[考研] 求调剂 +11 翩翩一书生 2026-04-09 13/650 2026-04-10 10:27 by liuhuiying09
[考研] 考研调剂 +13 冰冰,,, 2026-04-07 13/650 2026-04-09 17:01 by Lilly_Li
[考研] 求机械专硕297第二批调剂 +5 拾柒12。 2026-04-08 5/250 2026-04-09 16:43 by 允当适度
[考研] 086004 求调剂 309 +7 Yin DY 2026-04-08 7/350 2026-04-09 13:59 by Delta2012
信息提示
请填处理意见