24小时热门版块排行榜    

查看: 2357  |  回复: 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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料调剂 +3 匹克i 2026-03-23 3/150 2026-03-23 23:18 by peike
[考研] 265求调剂 +10 梁梁校校 2026-03-17 10/500 2026-03-23 21:17 by 一切OK
[考研] 0703化学求调剂 +4 奶油草莓. 2026-03-22 5/250 2026-03-23 19:37 by pswait
[考研] 336化工调剂 +4 王大坦1 2026-03-23 5/250 2026-03-23 18:32 by allen-yin
[考研] 一志愿070300浙大化学358分,求调剂! +4 酥酥鱼.. 2026-03-21 4/200 2026-03-23 08:12 by Iveryant
[考研] 352求调剂 +3 大米饭! 2026-03-22 3/150 2026-03-22 23:28 by king123!
[考博] 招收博士1-2人 +3 QGZDSYS 2026-03-18 4/200 2026-03-22 10:25 by QGZDSYS
[考研] 0856材料专硕353求调剂 +4 NIFFFfff 2026-03-20 4/200 2026-03-22 09:49 by 2026paper
[考研] 资源与环境 调剂申请(333分) +5 holy J 2026-03-21 5/250 2026-03-21 22:42 by Catalysis25
[考研] 一志愿东华大学控制学硕320求调剂 +3 Grand777 2026-03-21 3/150 2026-03-21 19:23 by 简之-
[考研] 一志愿南大,0703化学,分数336,求调剂 +3 收到VS 2026-03-21 3/150 2026-03-21 18:42 by 学员8dgXkO
[考研] 【考研调剂】化学专业 281分,一志愿四川大学,诚心求调剂 +11 吃吃吃才有意义 2026-03-19 11/550 2026-03-21 18:23 by 学员8dgXkO
[考研] 一志愿天津大学化学工艺专业(081702)315分求调剂 +12 yangfz 2026-03-17 12/600 2026-03-21 03:30 by JourneyLucky
[考研] 307求调剂 +3 wyyyqx 2026-03-17 3/150 2026-03-21 03:20 by JourneyLucky
[考研] 二本跨考郑大材料306英一数二 +3 z1z2z3879 2026-03-17 3/150 2026-03-21 02:29 by JourneyLucky
[考研] 一志愿重庆大学085700资源与环境专硕,总分308求调剂 +3 墨墨漠 2026-03-18 3/150 2026-03-21 00:39 by JourneyLucky
[考研] 考研调剂求学校推荐 +3 伯乐29 2026-03-18 5/250 2026-03-20 22:59 by JourneyLucky
[考研] 261求B区调剂,科研经历丰富 +3 牛奶很忙 2026-03-20 4/200 2026-03-20 19:34 by JourneyLucky
[考研] 288求调剂,一志愿华南理工大学071005 +5 ioodiiij 2026-03-17 5/250 2026-03-19 18:22 by zcl123
[考研] 一志愿福大288有机化学,求调剂 +3 小木虫200408204 2026-03-18 3/150 2026-03-19 13:31 by houyaoxu
信息提示
请填处理意见