Revision b32b6ebd11b49349893ec69fc4788cf7ede26003 authored by mishashe on 07 June 2021, 16:47:27 UTC, committed by GitHub on 07 June 2021, 16:47:27 UTC
1 parent 023b4e1
ConnectNRToSEED.py
#Input file with connections of NR accession numbers to SEED ontology
#NR accession SSIDs Level1 level2 level3 level4
#Output present in the dataset accessions which can be connected to SEED
SS = {}
with open('SEED_annotation') as file1:
for line in file1:
dat = line.strip().split('\t')
SS[dat[0]] = '\t'.join(dat[2:])
out = open('matches_SEED','w')
with open('alignment.txt') as file1:
for line in file1:
dat = line.strip().split('\t')
if float(dat[-1]) > 50:
if dat[0] in SS:
out.write(SS[dat[0]] + '\n')
out.close()
Computing file changes ...