import os,re,sys
infile=file("test.fq.txt")
def qual_stat(qstr):
q20 = 0
q30 = 0
for q in qstr:
qual = ord(q) - 33# ord 返回對應(yīng)的ASCII數(shù)值
if qual >= 30:
q30 += 1
q20 += 1
elif qual >= 20:
q20 += 1
return q20, q30
arr=[]
for line in infile:
line=line.strip()
arr.append(line)
j=0
a,g,c,t,n,dna=0,0,0,0,0,0
total_count = 0
q20_count = 0
q30_count = 0
for i in range(len(arr)):
if arr[i].startswith("@"):
j+=1
#dna+=arr[i+1].count("A")+arr[i+1].count("G")+arr[i+1].count("C")+arr[i+1].count("T")
dna+=len(arr[i+1])
a+=arr[i+1].count("A")
g+=arr[i+1].count("G")
c+=arr[i+1].count("C")
t+=arr[i+1].count("T")
n+=arr[i+1].count("N")
#total_count += len(arr[i+3])
q20, q30 = qual_stat(arr[i+3])
#q20_count += q20
q30_count += q30
ap=float(a)*100/float(dna)
gp=float(g)*100/float(dna)
cp=float(c)*100/float(dna)
tp=float(t)*100/float(dna)
np=float(n)*100/float(dna)
q30=float(q30_count)*100/float(dna)
print "Total Reads:",j
print "Total Bases:",dna
print "Q30 Bases:",q30_count
print "q30 percents:",("%.2f" % q30)+"%"
print "A:",a,("%.2f" % ap)+"%"
print "G:",g,("%.2f" % gp)+"%"
print "C:",c,("%.2f" % cp)+"%"
print "T:",t,("%.2f" % tp)+"%"
print "N:",n,("%.2f" % np)+"%"
統(tǒng)計(jì)fastaq信息逊彭，包括Q30瓦戚、堿基含量信息
最后編輯于 ：
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
- 文/潘曉璐 我一進(jìn)店門兵罢，熙熙樓的掌柜王于貴愁眉苦臉地迎上來族壳，“玉大人，你說我怎么就攤上這事趣些。” “怎么了贰您？”我有些...
- 文/不壞的土叔 我叫張陵坏平，是天一觀的道長拢操。 經(jīng)常有香客問我，道長舶替，這世上最難降的妖魔是什么令境？ 我笑而不...
- 正文 為了忘掉前任，我火速辦了婚禮顾瞪，結(jié)果婚禮上舔庶，老公的妹妹穿的比我還像新娘。我一直安慰自己陈醒，他們只是感情好惕橙，可當(dāng)我...
- 文/花漫 我一把揭開白布。 她就那樣靜靜地躺著钉跷，像睡著了一般弥鹦。 火紅的嫁衣襯著肌膚如雪。 梳的紋絲不亂的頭發(fā)上爷辙，一...
- 文/蒼蘭香墨 我猛地睜開眼，長吁一口氣：“原來是場噩夢啊……” “哼歹颓！你這毒婦竟也來了坯屿？” 一聲冷哼從身側(cè)響起，我...
- 序言：老撾萬榮一對情侶失蹤巍扛，失蹤者是張志新（化名）和其女友劉穎领跛，沒想到半個(gè)月后，有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體撤奸，經(jīng)...
- 正文 獨(dú)居荒郊野嶺守林人離奇死亡吠昭，尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
- 正文 我和宋清朗相戀三年，在試婚紗的時(shí)候發(fā)現(xiàn)自己被綠了胧瓜。 大學(xué)時(shí)的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片矢棚。...
- 正文 年R本政府宣布兜粘，位于F島的核電站申窘，受9級特大地震影響，放射性物質(zhì)發(fā)生泄漏孔轴。R本人自食惡果不足惜剃法，卻給世界環(huán)境...
- 文/蒙蒙 一、第九天 我趴在偏房一處隱蔽的房頂上張望路鹰。 院中可真熱鬧贷洲，春花似錦、人聲如沸晋柱。這莊子的主人今日做“春日...
- 文/蒼蘭香墨 我抬頭看了看天上的太陽趣斤。三九已至俩块，卻和暖如春，著一層夾襖步出監(jiān)牢的瞬間浓领，已是汗流浹背玉凯。 一陣腳步聲響...
- 正文 我出身青樓，卻偏偏與公主長得像泪幌，于是被迫代替她去往敵國和親盲厌。 傳聞我的和親對象是個(gè)殘疾皇子，可洞房花燭夜當(dāng)晚...
推薦閱讀更多精彩內(nèi)容
- 錨的就是你 英語難學(xué)嗎右蒲？ 如果一個(gè)人認(rèn)為英語難學(xué)阀湿，恐怕他的英語水平不會(huì)太好。 從小被植入腦海中這樣的錨：英語難瑰妄、過...
