diff blast_to_scaffold.py @ 3:e279e53a3fde draft

planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/blast_to_scaffold commit 030207144f0811822dbdda9a10e036ff8e794d7c
author drosofff
date Fri, 25 Mar 2016 19:14:20 -0400
parents 61731b911f4e
children 3288cc5a57e5
line wrap: on
line diff
--- a/blast_to_scaffold.py	Fri Mar 25 17:01:52 2016 -0400
+++ b/blast_to_scaffold.py	Fri Mar 25 19:14:20 2016 -0400
@@ -30,10 +30,6 @@
         '--blast-tab', dest="blast_tab", action="store", type=str, help="13-columns tabular blastn or tblastx output")
     the_parser.add_argument(
         '--output', action="store", type=str, help="output file path, fasta format")
-    the_parser.add_argument(
-        '--scaffold_prefix', action="store", type=str, help="the prefix that will be used for the header of the fasta scaffold")
-    the_parser.add_argument(
-        '--scaffold_suffix', action="store", type=str, help="the sufix that will be used for the header of the fasta scaffold")
     args = the_parser.parse_args()
     return args
     
@@ -102,13 +98,13 @@
         for i, nucleotide in enumerate(sequence):
             GuideDict[i+subjectStart] = nucleotide
             
-def finalAssembly (GuideDict, outputfile, prefix, suffix):
+def finalAssembly (GuideDict, outputfile):
     finalSeqList = []
     for keys in sorted(GuideDict):
         finalSeqList.append(GuideDict[keys])
     finalSequence = insert_newlines("".join(finalSeqList) )
     Out = open (outputfile, "w")
-    print >> Out, ">Scaffold_from_%s_guided_by_%s" % (prefix, suffix)
+    print >> Out, ">Scaffold"
     print >> Out, finalSequence
     Out.close()
     
@@ -118,7 +114,7 @@
     GuideDict = myGuide (args.guideSequence)
     blastlist = blatnInfo(args.blast_tab)
     updateGuide(blastlist, GuideDict, ContigsDict)
-    finalAssembly(GuideDict, args.output, args.scaffold_prefix, args.scaffold_suffix)
+    finalAssembly(GuideDict, args.output)
 
 if __name__ == "__main__":
     __main__()