将 paired count 和unpaired count 相加

尝试 count-based 算表达量过后，因为在map这一步一直是把过滤过后的paired 和unpaired reads 一起做的，后面count的时候就要出问题。发现可以用 samtools 从 bam 里分出来，于是就分开count，最后把reads 数相加。在用DESeq什么的。

首先，记下 samtools 命令

samtools view -bf  -1 all.bam > paired.bamsamtools view -bF -1 all.bam > unpaired.bam

然后是用上次那个提取1，6列的程序，现在就相加：

def readf(filename):lines = open(filename,'r').readlines()s=[]for i in lines:i=i.rstrip('\n')i_ = i.split('\t')s.append(i_[1])return ss_unpaired=readf("unpaired.txt")print(s_unpaired)f=open('paired.txt','r')file=open('out.txt','w')i=0head='Geneid\tcount\n'file.write(head)for lines in f.readlines():lines=lines.rstrip('\n')line=lines.split('\t')#print(line[0])if(i>0):line[1]=int(line[1])+int(s_unpaired[i])print(line)final=line[0]+'\t'+str(line[1])file.write(final+'\n')i+=1file.close