# HG changeset patch # User drosofff # Date 1497174989 14400 # Node ID 74394e39ad22b9e7054af8d1e6c7f1415d5ea51d # Parent 933475df0fd2912103ec97a5672cc290ad27b5fa planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/MirParser commit 76b92fc0dd106f009b1df4221347fd892cfdd31f diff -r 933475df0fd2 -r 74394e39ad22 MirParser.py --- a/MirParser.py Wed Oct 21 11:25:01 2015 -0400 +++ b/MirParser.py Sun Jun 11 05:56:29 2017 -0400 @@ -8,7 +8,7 @@ import sys import subprocess -from smRtools import * +from smRtools import HandleSmRNAwindows IndexSource = sys.argv[1] ExtractionDirective = sys.argv[2] diff -r 933475df0fd2 -r 74394e39ad22 MirParser.xml --- a/MirParser.xml Wed Oct 21 11:25:01 2015 -0400 +++ b/MirParser.xml Sun Jun 11 05:56:29 2017 -0400 @@ -1,40 +1,38 @@ - + from aligment files - bowtie - pysam - numpy - R - biocbasics - samtools + bowtie + numpy + pysam + r-lattice - - MirParser.py - #if $refGenomeSource.genomeSource == "history": - $refGenomeSource.ownFile ## index source sys.arg[1] - --do_not_extract_index ## sys.argv[2] - #else: - #silent reference= filter( lambda x: str( x[0] ) == str( $input_list.dbkey ), $__app__.tool_data_tables[ 'bowtie_indexes' ].get_fields() )[0][-1] - $reference ## sys.argv[1] + + #end if + #for $i in $refGenomeSource.input_list + '$i' '$i.ext' "$i.name" ## sys.argv[9,10,11] modulo 3 + #end for + #silent plottingoption = $plotting.plottingOption + ]]> diff -r 933475df0fd2 -r 74394e39ad22 smRtools.py --- a/smRtools.py Wed Oct 21 11:25:01 2015 -0400 +++ b/smRtools.py Sun Jun 11 05:56:29 2017 -0400 @@ -142,26 +142,6 @@ self.alignedReads += 1 F.close() return self.instanceDict -# elif self.alignmentFileFormat == "sam": -# F = open (self.alignmentFile, "r") -# dict = {"0":"+", "16":"-"} -# for line in F: -# if line[0]=='@': -# continue -# fields = line.split() -# if fields[2] == "*": continue -# polarity = dict[fields[1]] -# gene = fields[2] -# offset = int(fields[3]) -# size = len (fields[9]) -# if self.size_inf: -# if (size>=self.size_inf and size<= self.size_sup): -# self.instanceDict[gene].addread (polarity, offset, size) -# self.alignedReads += 1 -# else: -# self.instanceDict[gene].addread (polarity, offset, size) -# self.alignedReads += 1 -# F.close() elif self.alignmentFileFormat == "bam" or self.alignmentFileFormat == "sam": import pysam samfile = pysam.Samfile(self.alignmentFile) @@ -184,22 +164,6 @@ self.alignedReads += 1 return self.instanceDict -# def size_histogram (self): -# size_dict={} -# size_dict['F']= defaultdict (int) -# size_dict['R']= defaultdict (int) -# size_dict['both'] = defaultdict (int) -# for item in self.instanceDict: -# buffer_dict_F = self.instanceDict[item].size_histogram()['F'] -# buffer_dict_R = self.instanceDict[item].size_histogram()['R'] -# for size in buffer_dict_F: -# size_dict['F'][size] += buffer_dict_F[size] -# for size in buffer_dict_R: -# size_dict['R'][size] -= buffer_dict_R[size] -# allSizeKeys = list (set (size_dict['F'].keys() + size_dict['R'].keys() ) ) -# for size in allSizeKeys: -# size_dict['both'][size] = size_dict['F'][size] + size_dict['R'][size] -# return size_dict def size_histogram (self): # in HandleSmRNAwindows '''refactored on 7-9-2014 to debug size_histogram tool''' size_dict={} @@ -361,24 +325,7 @@ for offset in range (min(dicsize.keys()), max(dicsize.keys())+1): dicsize[size] = dicsize.get(size, 0) # to fill offsets with null values return dicsize - -# def size_histogram(self): -# norm=self.norm -# hist_dict={} -# hist_dict['F']={} -# hist_dict['R']={} -# for offset in self.readDict: -# for size in self.readDict[offset]: -# if offset < 0: -# hist_dict['R'][size] = hist_dict['R'].get(size, 0) - 1*norm -# else: -# hist_dict['F'][size] = hist_dict['F'].get(size, 0) + 1*norm -# ## patch to avoid missing graphs when parsed by R-lattice. 27-08-2014. Test and validate ! -# if not (hist_dict['F']) and (not hist_dict['R']): -# hist_dict['F'][21] = 0 -# hist_dict['R'][21] = 0 -# ## -# return hist_dict + def size_histogram(self, minquery=None, maxquery=None): # in SmRNAwindow '''refactored on 7-9-2014 to debug size_histogram tool''' @@ -480,7 +427,6 @@ return ". | %s" % (freqDic["Trev"] / reverse_sum * 100) else: return "%s | %s" % (freqDic["Tfor"] / forward_sum * 100, freqDic["Trev"] / reverse_sum * 100) - def readplot (self): norm=self.norm diff -r 933475df0fd2 -r 74394e39ad22 smRtools.pyc Binary file smRtools.pyc has changed diff -r 933475df0fd2 -r 74394e39ad22 tool_dependencies.xml --- a/tool_dependencies.xml Wed Oct 21 11:25:01 2015 -0400 +++ b/tool_dependencies.xml Sun Jun 11 05:56:29 2017 -0400 @@ -1,21 +1,18 @@ - - + + - - + + - + - + - - -