之前用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 弄的还是昏头昏脑的，每次都要错了再改。
以上就是本文的全部内容，希望对大家的学习有所帮助，也希望大家多多支持。
免责声明：本站文章均来自网站采集或用户投稿，网站不提供任何软件下载或自行开发的软件！ 如有用户或公司发现本站内容信息存在侵权行为，请邮件告知！ 858582#qq.com
《魔兽世界》大逃杀！60人新游玩模式《强袭风暴》3月21日上线
暴雪近日发布了《魔兽世界》10.2.6 更新内容，新游玩模式《强袭风暴》即将于3月21 日在亚服上线，届时玩家将前往阿拉希高地展开一场 60 人大逃杀对战。
艾泽拉斯的冒险者已经征服了艾泽拉斯的大地及遥远的彼岸。他们在对抗世界上最致命的敌人时展现出过人的手腕，并且成功阻止终结宇宙等级的威胁。当他们在为即将于《魔兽世界》资料片《地心之战》中来袭的萨拉塔斯势力做战斗准备时，他们还需要在熟悉的阿拉希高地面对一个全新的敌人──那就是彼此。在《巨龙崛起》10.2.6 更新的《强袭风暴》中，玩家将会进入一个全新的海盗主题大逃杀式限时活动，其中包含极高的风险和史诗级的奖励。
《强袭风暴》不是普通的战场，作为一个独立于主游戏之外的活动，玩家可以用大逃杀的风格来体验《魔兽世界》，不分职业、不分装备（除了你在赛局中捡到的），光是技巧和战略的强弱之分就能决定出谁才是能坚持到最后的赢家。本次活动将会开放单人和双人模式，玩家在加入海盗主题的预赛大厅区域前，可以从强袭风暴角色画面新增好友。游玩游戏将可以累计名望轨迹，《巨龙崛起》和《魔兽世界：巫妖王之怒 经典版》的玩家都可以获得奖励。