DU7:如何對(duì)小RNA測序數(shù)據(jù)進(jìn)行長度分布統(tǒng)計(jì)

小RNA測序數(shù)據(jù)長度分布規(guī)律:
1懒浮、24nt>21nt>22nt>23nt
2锚国、>60%數(shù)據(jù)在20-24nt之間

一会喝、數(shù)據(jù)為fasta格式列荔,首先寫一個(gè)腳本統(tǒng)計(jì)各種長度的sRNA的數(shù)量

#!/usr/bin/env python3
import sys
import collections
inFile = open(sys.argv[1],'r')
outFile = open ('sRNA_count.csv', 'w')
lenlist = []
while True:
    line = inFile.readline()
    if not line:break
    if ">" not in line:
        line = line.rstrip()
        lenlist.append(len(line))
lenlist.sort()
lencount = collections.Counter(lenlist)
for length in lencount:
       outFile.write(str(length) + "\t" + str(lencount[length]) + "\n")
inFile.close()

#運(yùn)行命令
python sRNA_count.py sample.sRNA.data.fa
sRNA_count.csv

二、手動(dòng)將sRNA_count.csv進(jìn)行分列加表頭

sRNA_count.csv

三肾请、用R語言ggplot2繪制直方圖

install.packages('gcookbook')
library(ggplot2)
library(gcookbook)
sRNACount <- read.csv("J:/myProject/sRNA_count/sRNA_count.csv", header = TRUE)
sRNACount
ggplot(sRNACount, aes(x=sRNA, y=NUM)) + scale_y_log10() + xlim(15,30) + geom_bar(stat="identity", fill="lightblue", colour="black")
ggplot(sRNACount, aes(x=sRNA, y=NUM)) + scale_y_log10() + scale_x_continuous(limits=c(14,31),breaks=seq(15,30,1)) + geom_bar(stat="identity", fill="lightblue", colour="black")
sRNA長度分布圖
最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
  • 序言:七十年代末留搔,一起剝皮案震驚了整個(gè)濱河市,隨后出現(xiàn)的幾起案子铛铁,更是在濱河造成了極大的恐慌隔显,老刑警劉巖,帶你破解...
    沈念sama閱讀 217,277評(píng)論 6 503
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件饵逐,死亡現(xiàn)場離奇詭異括眠,居然都是意外死亡,警方通過查閱死者的電腦和手機(jī)倍权,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 92,689評(píng)論 3 393
  • 文/潘曉璐 我一進(jìn)店門掷豺,熙熙樓的掌柜王于貴愁眉苦臉地迎上來捞烟,“玉大人,你說我怎么就攤上這事当船√饣” “怎么了?”我有些...
    開封第一講書人閱讀 163,624評(píng)論 0 353
  • 文/不壞的土叔 我叫張陵德频,是天一觀的道長苍息。 經(jīng)常有香客問我,道長壹置,這世上最難降的妖魔是什么竞思? 我笑而不...
    開封第一講書人閱讀 58,356評(píng)論 1 293
  • 正文 為了忘掉前任,我火速辦了婚禮钞护,結(jié)果婚禮上盖喷,老公的妹妹穿的比我還像新娘。我一直安慰自己患亿,他們只是感情好传蹈,可當(dāng)我...
    茶點(diǎn)故事閱讀 67,402評(píng)論 6 392
  • 文/花漫 我一把揭開白布。 她就那樣靜靜地躺著步藕,像睡著了一般惦界。 火紅的嫁衣襯著肌膚如雪。 梳的紋絲不亂的頭發(fā)上咙冗,一...
    開封第一講書人閱讀 51,292評(píng)論 1 301
  • 那天沾歪,我揣著相機(jī)與錄音,去河邊找鬼雾消。 笑死灾搏,一個(gè)胖子當(dāng)著我的面吹牛,可吹牛的內(nèi)容都是我干的立润。 我是一名探鬼主播狂窑,決...
    沈念sama閱讀 40,135評(píng)論 3 418
  • 文/蒼蘭香墨 我猛地睜開眼,長吁一口氣:“原來是場噩夢(mèng)啊……” “哼桑腮!你這毒婦竟也來了泉哈?” 一聲冷哼從身側(cè)響起,我...
    開封第一講書人閱讀 38,992評(píng)論 0 275
  • 序言:老撾萬榮一對(duì)情侶失蹤破讨,失蹤者是張志新(化名)和其女友劉穎丛晦,沒想到半個(gè)月后,有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體提陶,經(jīng)...
    沈念sama閱讀 45,429評(píng)論 1 314
  • 正文 獨(dú)居荒郊野嶺守林人離奇死亡烫沙,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點(diǎn)故事閱讀 37,636評(píng)論 3 334
  • 正文 我和宋清朗相戀三年,在試婚紗的時(shí)候發(fā)現(xiàn)自己被綠了隙笆。 大學(xué)時(shí)的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片锌蓄。...
    茶點(diǎn)故事閱讀 39,785評(píng)論 1 348
  • 序言:一個(gè)原本活蹦亂跳的男人離奇死亡升筏,死狀恐怖,靈堂內(nèi)的尸體忽然破棺而出瘸爽,到底是詐尸還是另有隱情仰冠,我是刑警寧澤,帶...
    沈念sama閱讀 35,492評(píng)論 5 345
  • 正文 年R本政府宣布蝶糯,位于F島的核電站,受9級(jí)特大地震影響辆沦,放射性物質(zhì)發(fā)生泄漏昼捍。R本人自食惡果不足惜,卻給世界環(huán)境...
    茶點(diǎn)故事閱讀 41,092評(píng)論 3 328
  • 文/蒙蒙 一肢扯、第九天 我趴在偏房一處隱蔽的房頂上張望妒茬。 院中可真熱鬧,春花似錦蔚晨、人聲如沸乍钻。這莊子的主人今日做“春日...
    開封第一講書人閱讀 31,723評(píng)論 0 22
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽银择。三九已至,卻和暖如春累舷,著一層夾襖步出監(jiān)牢的瞬間浩考,已是汗流浹背。 一陣腳步聲響...
    開封第一講書人閱讀 32,858評(píng)論 1 269
  • 我被黑心中介騙來泰國打工被盈, 沒想到剛下飛機(jī)就差點(diǎn)兒被人妖公主榨干…… 1. 我叫王不留析孽,地道東北人。 一個(gè)月前我還...
    沈念sama閱讀 47,891評(píng)論 2 370
  • 正文 我出身青樓只怎,卻偏偏與公主長得像袜瞬,于是被迫代替她去往敵國和親。 傳聞我的和親對(duì)象是個(gè)殘疾皇子身堡,可洞房花燭夜當(dāng)晚...
    茶點(diǎn)故事閱讀 44,713評(píng)論 2 354