Mercurial > repos > eduardo > mergegffs
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():
