parsea_salida_blast

#!/usr/bin/python
import os
data=open('Fh_Sm', 'r')
datai=open('Sm_Fh', 'r')
asd=open('Fh_Sm.temp','w')
asdi=open('Sm_Fh.temp','w')
columna=[]
dat=[]

def soloalgunos(dat,temp,i,f):
    for linea in dat:
        col=linea.split('\t')
        if int(col[i]) < int(col[f]) and (int(col[i])-1)%3 == 0:
            temp.writelines(linea)
    temp.close()

ini=8
fin=9
soloalgunos(data,asd,ini,fin)
ini=6
fin=7
soloalgunos(datai,asdi,ini,fin)
#os.system("")
data.close()
asd.close()
asdi.close()
datai.close()

#saco las unicas de los temp o sea las que tienen mejor score
os.system("perl solouniq1.pl > Fh_Sm_uniq.temp")
os.system("perl solouniq2.pl > Sm_Fh_uniq.temp")