diff mergegffs.py @ 12:e188a09bf012 draft

planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/mergegffs commit 96fe2c546e9e6a58707cfaf6e6435ed4289c55a2
author eduardo
date Sat, 17 Jun 2017 07:52:11 -0400
parents 791fae70fa3c
children 742d726397b3
line wrap: on
line diff
--- a/mergegffs.py	Sun Apr 09 10:21:46 2017 -0400
+++ b/mergegffs.py	Sat Jun 17 07:52:11 2017 -0400
@@ -25,16 +25,16 @@
         db = gffutils.create_db(args.peptides, dbfn=args.database, force=True, keep_order=True,merge_strategy='merge', sort_attribute_values=True)
         if args.blastp is not None:
             os.system("convert2gff.py -i "+args.blastp+" -g  "+args.peptides+" -p blastp -t protein_match -T  -x > blastp.gff")
-            db.update("blastp.gff")
+            db=db.update("blastp.gff")
         if args.signalP is not None:
             os.system("convert2gff.py -i "+args.signalP+" -g "+args.peptides+" -T -p signalP -t signalpep > signalp.gff")
-            db.update("signalp.gff")
+            db=db.update("signalP.gff")
         if args.tmhmm is not None:
             os.system("convert2gff.py -i "+args.tmhmm+" -g "+args.peptides+" -T -p tmhmm -t trans_helix > tmhmm.gff")
-            db.update("tmhmm.gff")
+            db=db.update("tmhmm.gff")
         if args.pfam is not None:
             os.system("pfam2gff.py -i "+args.pfam+" -g "+args.peptides+" -T > PFAM.gff")
-            db.update("PFAM.gff")
+            db=db.update("PFAM.gff")
 
         with open(args.output, 'w') as fout:
             for f in db.all_features():