from csamtools import *
+from ctabix import *
+import csamtools
+import ctabix
+from cvcf import *
+import cvcf
import Pileup
import sys
import os
self.parsers = parsers
self.stderr = []
- def __call__(self,*args, **kwargs):
- '''execute the samtools command
+ def __call__(self, *args, **kwargs):
+ '''execute a samtools command
'''
retval, stderr, stdout = csamtools._samtools_dispatch( self.dispatch, args )
if retval: raise SamtoolsError( "\n".join( stderr ) )
# Note that there is sometimes output on stderr that is not an error,
# for example: [sam_header_read2] 2 sequences loaded.
# Ignore messages like these
- stderr = [ x for x in stderr if not x.startswith( "[sam_header_read2]" ) ]
+ stderr = [ x for x in stderr \
+ if not (x.startswith( "[sam_header_read2]" ) or \
+ x.startswith("[bam_index_load]") or \
+ x.startswith("[bam_sort_core]") or \
+ x.startswith("[samopen] SAM header is present") )
+ ]
if stderr: raise SamtoolsError( "\n".join( stderr ) )
# call parser for stdout:
globals()[key] = SamtoolsDispatcher(cmd, parser)
# hack to export all the symbols from csamtools
-__all__ = csamtools.__all__ + [ "SamtoolsError", "SamtoolsDispatcher" ] + list(SAMTOOLS_DISPATCH) +\
- ["Pileup",]
-
+__all__ = \
+ csamtools.__all__ + \
+ ctabix.__all__ + \
+ cvcf.__all__ +\
+ [ "SamtoolsError", "SamtoolsDispatcher" ] + list(SAMTOOLS_DISPATCH) +\
+ ["Pileup" ]
+
+from version import __version__, __samtools_version__