X-Git-Url: http://woldlab.caltech.edu/gitweb/?a=blobdiff_plain;f=pysam%2Fpysam_util.h;h=c6184eb756cca0241dc06449b829f4198ebef02a;hb=8ab946e230447a6fe686edabc5e3cfa445da1721;hp=bfbd6dd0016b73ad94a9ad2d34c5811f969a0d82;hpb=6a7c3f175b210cc16d09a5e8e4c1d47333dbe1c6;p=pysam.git diff --git a/pysam/pysam_util.h b/pysam/pysam_util.h index bfbd6dd..c6184eb 100644 --- a/pysam/pysam_util.h +++ b/pysam/pysam_util.h @@ -1,6 +1,10 @@ #ifndef PYSAM_UTIL_H #define PYSAM_UTIL_H +////////////////////////////////////////////////////////////////// +// set pysam standard error to point to stream +FILE * pysam_set_stderr( FILE * f ); + ////////////////////////////////////////////////////////////////// ////////////////////////////////////////////////////////////////// ////////////////////////////////////////////////////////////////// @@ -38,5 +42,24 @@ bam1_t * pysam_bam_update( bam1_t * b, // translate a nucleotide character to binary code unsigned char pysam_translate_sequence( const unsigned char s ); +// defined in bam_import.c +extern unsigned char bam_nt16_table[256]; + +// translate a reference string *s* to a tid +int pysam_reference2tid( bam_header_t *header, const char * s ); + +// return number of mapped reads for tid +uint32_t pysam_get_mapped( const bam_index_t *idx, const int tid ); + +// return number of unmapped reads for tid +uint32_t pysam_get_unmapped( const bam_index_t *idx, const int tid ); + +// debugging functions +/* #include "glf.h" */ +/* uint32_t pysam_glf_depth( glf1_t * g); */ + +/* #include "bam_maqcns.h" */ +/* void pysam_dump_glf( glf1_t * g, bam_maqcns_t * c ); */ + #endif