![240](https://upload.jianshu.io/users/upload_avatars/23256478/1d7488b7-b781-4580-9ac5-80a99414cf66.jpg?imageMogr2/auto-orient/strip|imageView2/1/w/240/h/240)
將GeneID與GO號一一對應(yīng)的兩列表格,轉(zhuǎn)換為同一GeneID對應(yīng)GO號排一行(如圖),主要使用字典歉铝、列表完成轉(zhuǎn)換。 代碼如下
最常用的Linux循環(huán)應(yīng)該是for和while绍刮,記錄兩個直接輸入command的實例(輸出Ta開頭的所有文件名字): 開始正片:使用同一個que...
100天生信-Day11 最近一直在做濕實驗览爵,需要大量設(shè)計引物烈拒,發(fā)現(xiàn)primer3有python版本嚷量,可以批量設(shè)計陋桂,簡直神器。 參考教程:htt...
100天生信-Day10 記錄一下workflow
100天生信-Day9 所有高通量測序的下機(jī)數(shù)據(jù)一般都要進(jìn)行質(zhì)控津肛,檢測數(shù)據(jù)質(zhì)量章喉,并去除低質(zhì)量reads。之前比較常用的是fastQC質(zhì)控 + c...
100天生信-Day8 箱線圖繪制代碼記錄:
100天生信-Day7 ChIP-seq拉下來的reads經(jīng)過mapping生成.sam文件身坐,sam文件使用samtools轉(zhuǎn)bam,再sort...
100天生信-Day6 長數(shù)據(jù)和寬數(shù)據(jù)是滿足不同分析的兩種數(shù)據(jù)格式落包,可以很方便用R的reshape2包轉(zhuǎn)換部蛇,代碼如下: 參考教程:https:/...
100天生信-Day5 ChIP-seq、ATAC-seq等需要在全基因組上找到特定位置的結(jié)合/切割/富集的高通量測序技術(shù)咐蝇,都需要進(jìn)行reads...