Mercurial > repos > drosofff > msp_blastparser_and_hits

comparison BlastParser_and_hits.py @ 5:a0dec1a0f2ef draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_blastparser_and_hits commit 59dcca02907c91f852e63db16acd72c9d79d6eb2
author drosofff
date Fri, 15 Jan 2016 12:51:19 -0500
parents 60b6bd959929
children 78c34df2dd8d
comparison
equal deleted inserted replaced
4:60b6bd959929 5:a0dec1a0f2ef
129 leftCoordinate = 1 129 leftCoordinate = 1
130 return getseq (fastadict, FastaHeader, leftCoordinate, rightCoordinate, polarity) 130 return getseq (fastadict, FastaHeader, leftCoordinate, rightCoordinate, polarity)
131 131
132 def outputParsing (F, Fasta, results, Xblastdict, fastadict, filter_relativeCov=0, filter_maxScore=0, filter_meanScore=0, filter_term_in="", filter_term_out="", mode="verbose"): 132 def outputParsing (F, Fasta, results, Xblastdict, fastadict, filter_relativeCov=0, filter_maxScore=0, filter_meanScore=0, filter_term_in="", filter_term_out="", mode="verbose"):
133 def filter_results (results, filter_relativeCov=0, filter_maxScore=0, filter_meanScore=0, filter_term_in="", filter_term_out=""): 133 def filter_results (results, filter_relativeCov=0, filter_maxScore=0, filter_meanScore=0, filter_term_in="", filter_term_out=""):
134 print "###", filter_term_in
135 for subject in results.keys(): 134 for subject in results.keys():
136 if results[subject]["RelativeSubjectCoverage"]<filter_relativeCov: 135 if results[subject]["RelativeSubjectCoverage"]<filter_relativeCov:
137 del results[subject] 136 del results[subject]
138 continue 137 continue
139 if results[subject]["maxBitScores"]<filter_maxScore: 138 if results[subject]["maxBitScores"]<filter_maxScore: