Compare commits
7 Commits
| Author | SHA1 | Date | |
|---|---|---|---|
| 343d9ec6df | |||
| 7a993ec841 | |||
| 6170fac081 | |||
| 764ee27c0d | |||
| 18d87aa1b1 | |||
| 1c91c27ee3 | |||
| 22d343b46a |
@@ -19,7 +19,7 @@ from os import path
|
|||||||
|
|
||||||
|
|
||||||
PACKAGE = "OBITools"
|
PACKAGE = "OBITools"
|
||||||
VERSION = "1.1.19"
|
VERSION = "1.1.22"
|
||||||
AUTHOR = 'Eric Coissac'
|
AUTHOR = 'Eric Coissac'
|
||||||
EMAIL = 'eric@coissac.eu'
|
EMAIL = 'eric@coissac.eu'
|
||||||
URL = 'metabarcoding.org/obitools'
|
URL = 'metabarcoding.org/obitools'
|
||||||
|
|||||||
+2
-2
@@ -352,7 +352,7 @@ if __name__=='__main__':
|
|||||||
|
|
||||||
search = lcsIteratorSelf(entries,db,options)
|
search = lcsIteratorSelf(entries,db,options)
|
||||||
|
|
||||||
print >>sys.stderr,'\nCache size : %d\n'
|
print >>sys.stderr,'\nCache size : %d\n' % __CACHE_SIZE__
|
||||||
|
|
||||||
|
|
||||||
for seq,best,match in search:
|
for seq,best,match in search:
|
||||||
@@ -442,7 +442,7 @@ if __name__=='__main__':
|
|||||||
|
|
||||||
|
|
||||||
writer(seq)
|
writer(seq)
|
||||||
print >>sys.stderr,'\n%5.3f% of the alignments was cached' % (__INCache__/(__INCache__+__OUTCache__)*100)
|
print >>sys.stderr,'\n%5.3f%% of the alignments was cached' % (__INCache__/(__INCache__+__OUTCache__)*100)
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|||||||
@@ -48,13 +48,14 @@ from obitools.options import getOptionManager, allEntryIterator
|
|||||||
from obitools.align import QSolexaReverseAssemble
|
from obitools.align import QSolexaReverseAssemble
|
||||||
from obitools.align import QSolexaRightReverseAssemble
|
from obitools.align import QSolexaRightReverseAssemble
|
||||||
from obitools.tools._solexapairend import buildConsensus
|
from obitools.tools._solexapairend import buildConsensus
|
||||||
from obitools.format.options import addInputFormatOption,addOutputFormatOption,\
|
from obitools.format.options import addOutputFormatOption,\
|
||||||
sequenceWriterGenerator
|
sequenceWriterGenerator
|
||||||
|
|
||||||
from itertools import chain
|
from itertools import chain
|
||||||
import cPickle
|
import cPickle
|
||||||
import math
|
import math
|
||||||
from obitools.fastq._fastq import fastqIterator
|
from obitools.fastq._fastq import fastqIterator # @UnresolvedImport
|
||||||
|
|
||||||
|
|
||||||
def addSolexaPairEndOptions(optionManager):
|
def addSolexaPairEndOptions(optionManager):
|
||||||
optionManager.add_option('-r','--reverse-reads',
|
optionManager.add_option('-r','--reverse-reads',
|
||||||
@@ -90,14 +91,6 @@ def addSolexaPairEndOptions(optionManager):
|
|||||||
const='illumina',
|
const='illumina',
|
||||||
help="input file is in fastq nucleic format produced by old solexa sequencer")
|
help="input file is in fastq nucleic format produced by old solexa sequencer")
|
||||||
|
|
||||||
# optionManager.add_option('--proba',
|
|
||||||
# action="store", dest="proba",
|
|
||||||
# metavar="<FILENAME>",
|
|
||||||
# type="str",
|
|
||||||
# default=None,
|
|
||||||
# help="null ditribution data file")
|
|
||||||
|
|
||||||
|
|
||||||
optionManager.add_option('--score-min',
|
optionManager.add_option('--score-min',
|
||||||
action="store", dest="smin",
|
action="store", dest="smin",
|
||||||
metavar="#.###",
|
metavar="#.###",
|
||||||
@@ -105,13 +98,6 @@ def addSolexaPairEndOptions(optionManager):
|
|||||||
default=None,
|
default=None,
|
||||||
help="minimum score for keeping aligment")
|
help="minimum score for keeping aligment")
|
||||||
|
|
||||||
# optionManager.add_option('--pvalue',
|
|
||||||
# action="store", dest="pvalue",
|
|
||||||
# metavar="#.###",
|
|
||||||
# type="float",
|
|
||||||
# default=None,
|
|
||||||
# help="maximum pvalue for keeping aligment")
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
def cutDirectReverse(entries):
|
def cutDirectReverse(entries):
|
||||||
@@ -222,6 +208,8 @@ if __name__ == '__main__':
|
|||||||
|
|
||||||
#WARNING TO REMOVE : DIRTY PATCH !
|
#WARNING TO REMOVE : DIRTY PATCH !
|
||||||
options.proba = None
|
options.proba = None
|
||||||
|
options.skip = None
|
||||||
|
options.only = None
|
||||||
|
|
||||||
|
|
||||||
options.sminL = None
|
options.sminL = None
|
||||||
|
|||||||
+1
-1
@@ -179,7 +179,7 @@ if __name__ == '__main__':
|
|||||||
i=0
|
i=0
|
||||||
for c in classes:
|
for c in classes:
|
||||||
i+=1
|
i+=1
|
||||||
progressBar(i,lclasses,False,"%15s" % c)
|
progressBar(i,lclasses,False,"%15s" % ("/".join(map(str,c)),))
|
||||||
seqs = classes[c]
|
seqs = classes[c]
|
||||||
sortclass(seqs, options)
|
sortclass(seqs, options)
|
||||||
if len(c)==1:
|
if len(c)==1:
|
||||||
|
|||||||
@@ -33,10 +33,6 @@ from obitools.format.sequence import skipOnErrorIterator, skipfirst, only
|
|||||||
from obitools import BioSequence
|
from obitools import BioSequence
|
||||||
from obitools.utils import FakeFile
|
from obitools.utils import FakeFile
|
||||||
|
|
||||||
from glob import glob
|
|
||||||
from test.test_compiler import Toto
|
|
||||||
|
|
||||||
|
|
||||||
def binarySequenceIterator(lineiterator):
|
def binarySequenceIterator(lineiterator):
|
||||||
|
|
||||||
f = FakeFile(lineiterator)
|
f = FakeFile(lineiterator)
|
||||||
@@ -263,9 +259,9 @@ def autoEntriesIterator(options):
|
|||||||
raise AssertionError,'file is not in fasta, fasta, embl, genbank or ecoPCR format'
|
raise AssertionError,'file is not in fasta, fasta, embl, genbank or ecoPCR format'
|
||||||
|
|
||||||
if reader==binarySequenceIterator:
|
if reader==binarySequenceIterator:
|
||||||
input = binarySequenceIterator(lineiterator)
|
input = binarySequenceIterator(lineiterator) # @ReservedAssignment
|
||||||
else:
|
else:
|
||||||
input = reader(chain([first],lineiterator))
|
input = reader(chain([first],lineiterator)) # @ReservedAssignment
|
||||||
|
|
||||||
return input
|
return input
|
||||||
|
|
||||||
@@ -313,11 +309,11 @@ def autoEntriesIterator(options):
|
|||||||
if options.skiperror:
|
if options.skiperror:
|
||||||
reader = skipOnErrorIterator(reader)
|
reader = skipOnErrorIterator(reader)
|
||||||
|
|
||||||
if options.skip is not None:
|
if hasattr(options, 'skip') and options.skip is not None:
|
||||||
print >>sys.stderr,"Skipping %d sequences" % options.skip
|
print >>sys.stderr,"Skipping %d sequences" % options.skip
|
||||||
reader = skipfirst(reader,options.skip)
|
reader = skipfirst(reader,options.skip)
|
||||||
|
|
||||||
if options.only is not None:
|
if hasattr(options, 'only') and options.only is not None:
|
||||||
print >>sys.stderr,"Analysing only %d sequences" % options.only
|
print >>sys.stderr,"Analysing only %d sequences" % options.only
|
||||||
reader = only(reader,options.only)
|
reader = only(reader,options.only)
|
||||||
|
|
||||||
@@ -325,7 +321,7 @@ def autoEntriesIterator(options):
|
|||||||
|
|
||||||
def sequenceWriterGenerator(options,output=sys.stdout):
|
def sequenceWriterGenerator(options,output=sys.stdout):
|
||||||
class SequenceWriter:
|
class SequenceWriter:
|
||||||
def __init__(self,options,file=sys.stdout):
|
def __init__(self,options,file=sys.stdout): # @ReservedAssignment
|
||||||
self._format=None
|
self._format=None
|
||||||
self._file=file
|
self._file=file
|
||||||
self._upper=options.uppercase
|
self._upper=options.uppercase
|
||||||
@@ -350,7 +346,7 @@ def sequenceWriterGenerator(options,output=sys.stdout):
|
|||||||
sys.exit(0)
|
sys.exit(0)
|
||||||
|
|
||||||
class BinaryWriter:
|
class BinaryWriter:
|
||||||
def __init__(self,options,file=sys.stdout):
|
def __init__(self,options,file=sys.stdout): # @ReservedAssignment
|
||||||
self._file=file
|
self._file=file
|
||||||
self._file.write("#!Pickle\n")
|
self._file.write("#!Pickle\n")
|
||||||
def put(self,seq):
|
def put(self,seq):
|
||||||
|
|||||||
@@ -1,5 +1,5 @@
|
|||||||
major = 1
|
major = 1
|
||||||
minor = 1
|
minor = 1
|
||||||
serial= '19'
|
serial= '22'
|
||||||
|
|
||||||
version = "%2d.%02d %s" % (major,minor,serial)
|
version = "%2d.%02d %s" % (major,minor,serial)
|
||||||
|
|||||||
Reference in New Issue
Block a user