Mercurial > repos > drosofff > blast_to_scaffold
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__()