diff blastxml_to_gapped_gff3.py @ 40:d5e8d72c432a draft

planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/jbrowse commit 7be399c9ee10ecd4b076cc1df775f585c3eb0eeb
author iuc
date Wed, 01 Mar 2017 23:06:18 -0500
parents 154510b45303
children c1451868f8d4
line wrap: on
line diff
--- a/blastxml_to_gapped_gff3.py	Thu Feb 09 12:36:02 2017 -0500
+++ b/blastxml_to_gapped_gff3.py	Wed Mar 01 23:06:18 2017 -0500
@@ -1,12 +1,10 @@
 #!/usr/bin/perl
-import argparse
-import copy
-import logging
 import re
 import sys
-
+import copy
+import argparse
 from BCBio import GFF
-
+import logging
 logging.basicConfig(level=logging.INFO)
 log = logging.getLogger(name='blastxml2gff3')
 
@@ -220,6 +218,7 @@
         else:
             log.warn("Bad data: \n\t%s\n\t%s\n\t%s\n" % (query, match, subject))
 
+
         if strict_m:
             if ret == '=' or ret == 'X':
                 ret = 'M'
@@ -253,7 +252,7 @@
 
 if __name__ == '__main__':
     parser = argparse.ArgumentParser(description='Convert Blast XML to gapped GFF3', epilog='')
-    parser.add_argument('blastxml', type=open, help='Blast XML Output')
+    parser.add_argument('blastxml', type=file, help='Blast XML Output')
     parser.add_argument('--min_gap', type=int, help='Maximum gap size before generating a new match_part', default=3)
     parser.add_argument('--trim', action='store_true', help='Trim blast hits to be only as long as the parent feature')
     parser.add_argument('--trim_end', action='store_true', help='Cut blast results off at end of gene')