之前用featureCount 處理得到結果，要提出第一列gene_id 和 readcount 列，首先軟件輸出的第一行默認是你使用的命令行，沒有用，用bash批量刪掉。
for i in `ls`;do sed -i '1d' $i;done
刪除當前文件夾下所有文件第一行。
其實提出兩列很簡單，不過我受夠了每次一個文件執行一次的煩。想搞成別的程序調用時命令行參數直接就行。第一次知道sys.argv這玩意，學到了。
我設置了 -i 輸入，-o 輸出 這兩個參數來判斷輸入輸出文件個數，不過對于錯誤命令行輸入的判斷還很欠缺，畢竟只有自己用。
上代碼：
import sys
leng=len(sys.argv)
for i in range(leng): # index the input file and output file location,maybe it's too complex
if sys.argv[i]=="-i":
it=i
if sys.argv[i]=="-o":
out=i
for i in range(it+1,out): # input file number(there should have same # of in and out file)
print(sys.argv[i])
file=open(sys.argv[i],'r')
f=open(sys.argv[i+out-it],'w')
count=[]
for lines in file.readlines():
two=[]
lines=lines.rstrip('\n')
b=lines.split('\t')
two.append(b[0])
if(b[6].startswith('.')): #chage path to count,some output problem
b[6]="count"
two.append(b[6])
s='\t'.join(two)
f.write(s+'\n')
f.close
其實python里面 str 和list 弄的還是昏頭昏腦的，每次都要錯了再改。
以上就是本文的全部內容，希望對大家的學習有所幫助，也希望大家多多支持腳本之家。
更多文章、技術交流、商務合作、聯系博主
微信掃碼或搜索：z360901061
微信掃一掃加我為好友
QQ號聯系： 360901061
您的支持是博主寫作最大的動力，如果您喜歡我的文章，感覺我的文章對您有幫助，請用微信掃描下面二維碼支持博主2元、5元、10元、20元等您想捐的金額吧，狠狠點擊下面給點支持吧，站長非常感激您！手機微信長按不能支付解決辦法：請將微信支付二維碼保存到相冊，切換到微信，然后點擊微信右上角掃一掃功能，選擇支付二維碼完成支付。
【本文對您有幫助就好】元
