| 查看: 2337 | 回復(fù): 15 | ||
| 本帖產(chǎn)生 1 個(gè) BioEPI ,點(diǎn)擊這里進(jìn)行查看 | ||
我要學(xué)perl么新蟲(chóng) (初入文壇)
|
[求助]
【求助】如何得到某一物種的所有蛋白翻譯基因的不同位置(exon、intron)的GC含量
|
|
|
如題: 如何得到某一物種的所有蛋白翻譯基因的不同位置(exon、intron)的GC含量,如下圖 是要從網(wǎng)上下染色體的全部序列信息然后再自己計(jì)算嗎? 圖片1.png [ Last edited by 我要學(xué)perl么 on 2013-8-14 at 22:30 ] |
至尊木蟲(chóng) (著名寫(xiě)手)
新蟲(chóng) (初入文壇)
至尊木蟲(chóng) (著名寫(xiě)手)
|
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 |
新蟲(chóng) (初入文壇)
新蟲(chóng) (初入文壇)
至尊木蟲(chóng) (著名寫(xiě)手)
新蟲(chóng) (初入文壇)
|
你好!謝謝,Galaxy是非常好用的工具。 蛋白編碼基因的內(nèi)含子、外顯子和上下游的GC含量我會(huì)計(jì)算了。但是對(duì)于lncRNA仍有點(diǎn)問(wèn)題。 現(xiàn)在我有l(wèi)ncRNA的gtf和gff兩種文件,對(duì)于如何計(jì)算GC含量我有2個(gè)想法,不知道哪個(gè)對(duì) 1. 第一種就是講gtf轉(zhuǎn)化為bed,然后使用Extract Features的Gene BED To Exon/Intron/Codon BED expander工具得到包含內(nèi)含子的bed文件,然后再轉(zhuǎn)化為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過(guò)濾得到transcrip和exon的兩種文件,然后用Operate on Genomic Intervals里的Subtract工具可以得到內(nèi)含子的文件,Get flanks 工具可以得到上下游序列。然后就可以計(jì)算了。 感覺(jué)第二種方法是對(duì)的。 By the way,在用第一種方法的時(shí)候如果transcription和exon的起始終止位置一樣的話,Gene BED To Exon/Intron/Codon BED expander工具判斷內(nèi)含子的話會(huì)出現(xiàn)start 和end相反的情況(比如說(shuō)正鏈的內(nèi)含子應(yīng)該為0的,可是出現(xiàn)start為1300 end 為300)。不知道這樣算不算錯(cuò)誤? 寫(xiě)的有點(diǎn)亂,可能表達(dá)不清楚,見(jiàn)諒! 謝謝你的回復(fù)! |
至尊木蟲(chóng) (著名寫(xiě)手)
新蟲(chóng) (初入文壇)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿南昌大學(xué),327分,材料與化工085600 +5 | Ncdx123456 2026-03-19 | 5/250 |
|
|---|---|---|---|---|
|
[考研]
|
.6lL 2026-03-18 | 7/350 |
|
|
[論文投稿]
申請(qǐng)回稿延期一個(gè)月,編輯同意了。但系統(tǒng)上的時(shí)間沒(méi)變,給編輯又寫(xiě)郵件了,沒(méi)回復(fù)
10+3
|
wangf9518 2026-03-17 | 4/200 |
|
|
[考研] 288求調(diào)劑 +15 | 于海海海海 2026-03-19 | 15/750 |
|
|
[考研] 復(fù)試調(diào)劑 +4 | z1z2z3879 2026-03-14 | 6/300 |
|
|
[考研] 085601材料工程專碩求調(diào)劑 +10 | 慕寒mio 2026-03-16 | 10/500 |
|
|
[考研] 一志愿福大288有機(jī)化學(xué),求調(diào)劑 +3 | 小木蟲(chóng)200408204 2026-03-18 | 3/150 |
|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 6/300 |
|
|
[考研] 297求調(diào)劑 +8 | 戲精丹丹丹 2026-03-17 | 8/400 |
|
|
[考研] 一志愿西南交大,求調(diào)劑 +4 | 材化逐夢(mèng)人 2026-03-18 | 4/200 |
|
|
[考研] 0703化學(xué)求調(diào)劑 總分331 +3 | ZY-05 2026-03-13 | 3/150 |
|
|
[考研] 環(huán)境工程調(diào)劑 +8 | 大可digkids 2026-03-16 | 8/400 |
|
|
[考研] 278求調(diào)劑 +5 | 煙火先于春 2026-03-17 | 5/250 |
|
|
[基金申請(qǐng)] 被我言中:新模板不強(qiáng)調(diào)格式了,假專家開(kāi)始管格式了 +4 | beefly 2026-03-14 | 4/200 |
|
|
[考研] 301求調(diào)劑 +4 | A_JiXing 2026-03-16 | 4/200 |
|
|
[考研] 332求調(diào)劑 +6 | Zz版 2026-03-13 | 6/300 |
|
|
[考研] 275求調(diào)劑 +4 | 太陽(yáng)花天天開(kāi)心 2026-03-16 | 4/200 |
|
|
[考研] 321求調(diào)劑 +5 | 大米飯! 2026-03-15 | 5/250 |
|
|
[考研] 求老師收留調(diào)劑 +4 | jiang姜66 2026-03-14 | 5/250 |
|
|
[考研] 289求調(diào)劑 +4 | 這么名字咋樣 2026-03-14 | 6/300 |
|