Mercurial > repos > eduardo > mergegffs
comparison mergegffs.py @ 14:5e27f44500d6 draft default tip
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/mergegffs commit 230f5fe14cc469e56626201a5c377686976d81fc-dirty
| author | eduardo |
|---|---|
| date | Mon, 26 Jun 2017 03:49:04 -0400 |
| parents | 742d726397b3 |
| children |
comparison
equal
deleted
inserted
replaced
| 13:742d726397b3 | 14:5e27f44500d6 |
|---|---|
| 23 | 23 |
| 24 | 24 |
| 25 db = gffutils.create_db(args.peptides, dbfn=args.database, force=True, keep_order=True,merge_strategy='merge', sort_attribute_values=True) | 25 db = gffutils.create_db(args.peptides, dbfn=args.database, force=True, keep_order=True,merge_strategy='merge', sort_attribute_values=True) |
| 26 if args.blastp is not None: | 26 if args.blastp is not None: |
| 27 os.system("convert2gff.py -i "+args.blastp+" -g "+args.peptides+" -p blastp -t protein_match -T -x > blastp.gff") | 27 os.system("convert2gff.py -i "+args.blastp+" -g "+args.peptides+" -p blastp -t protein_match -T -x > blastp.gff") |
| 28 db=db.update("blastp.gff") | 28 db.update("blastp.gff") |
| 29 if args.signalP is not None: | 29 if args.signalP is not None: |
| 30 os.system("convert2gff.py -i "+args.signalP+" -g "+args.peptides+" -T -p signalP -t signalpep > signalp.gff") | 30 os.system("convert2gff.py -i "+args.signalP+" -g "+args.peptides+" -T -p signalP -t signalpep > signalp.gff") |
| 31 db=db.update("signalp.gff") | 31 db.update("signalp.gff") |
| 32 if args.tmhmm is not None: | 32 if args.tmhmm is not None: |
| 33 os.system("convert2gff.py -i "+args.tmhmm+" -g "+args.peptides+" -T -p tmhmm -t trans_helix > tmhmm.gff") | 33 os.system("convert2gff.py -i "+args.tmhmm+" -g "+args.peptides+" -T -p tmhmm -t trans_helix > tmhmm.gff") |
| 34 db=db.update("tmhmm.gff") | 34 db.update("tmhmm.gff") |
| 35 if args.pfam is not None: | 35 if args.pfam is not None: |
| 36 os.system("pfam2gff.py -i "+args.pfam+" -g "+args.peptides+" -T > PFAM.gff") | 36 os.system("pfam2gff.py -i "+args.pfam+" -g "+args.peptides+" -T > PFAM.gff") |
| 37 db=db.update("PFAM.gff") | 37 db.update("PFAM.gff") |
| 38 | 38 |
| 39 with open(args.output, 'w') as fout: | 39 with open(args.output, 'w') as fout: |
| 40 for f in db.all_features(): | 40 for f in db.all_features(): |
| 41 fout.write(str(f) + '\n') | 41 fout.write(str(f) + '\n') |
| 42 if __name__ == "__main__": | 42 if __name__ == "__main__": |