Skip to content
Snippets Groups Projects
Commit 726c47cd authored by Remi  PLANEL's avatar Remi PLANEL
Browse files

features qualifiers in object

parent 7f69bdb7
No related branches found
No related tags found
No related merge requests found
......@@ -76,6 +76,7 @@ def from_seq(config, target, seq_format, output_file): # genome,
for record in SeqIO.parse(target, seq_format):
if config.verbose:
click.secho(" - search guide RNAs for %s " % record.id, fg=fg)
# g = SeqIO.parse(genome, "genbank")
guide_rnas = on_target_predict(str(record.seq))
write_guide_rnas(guide_rnas, output_file, record.id)
......
......@@ -72,12 +72,12 @@ def on_target_predict(seq, genome=None):
off_target_list = []
features = off_target_df.loc[0:, "features"]
for feat in features:
features_list = []
for x in feat:
features_dict = {}
for k, feat in x.qualifiers.items():
if k != "translation":
features_list.append((k, feat))
off_target_list.append(features_list)
features_dict[k] = " :: ".join(feat)
off_target_list.append(features_dict)
target.update({"off_targets": off_target_list})
else:
target.update({"off_targets": []})
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment