Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • Proteos
    Junior Member
    • Dec 2011
    • 6

    #1

    HTSeq - how do I cimport it for Cython?

    TL;DR
    How do I `cimport HTSeq` in Cython module?

    --- Details ---
    I am trying to use HTSeq to read BAM files that are quite large but it takes days to process them in pure Python. Therefore I decided to use Cython and cythonize reading BAM file part.
    Here is example code:
    Code:
    import HTSeq
    
    bam_file = '.../test.bam'
    bam = HTSeq.BAM_Reader(bam_file)
    for aln in bam:
        # process alignment
    Now I am trying to replace `import HTSeq` with `cimport HTSeq` but the Cython can not find `HTSeq.pxd` file. I found that the header file is actually `src/HTSeq/_HTSeq.pxd`(full link: https://github.com/simon-anders/htse...Seq/_HTSeq.pxd )
    as also described here: http://htseq.readthedocs.io/en/master/contrib.html
    So, the file starts with underscore and is in the htseq repo but `pip install` does not copy it to any of Includes directories. Therefore, I copied it manually to the root directory of my package so that `setup.py` can see it. Then I added `cimport _HTSeq as HTSeq` to my `*.pyx` file and it got compiled to an `*.so` file but when I ran the app it throws error:
    Code:
    ImportError: No module named _HTSeq
    My system: macOS 10.13.3 High Sierra; Python 2.7.14 (will port the code to 3.6 later); HTSeq 0.6.0;

    Any suggestion about solving this problem OR speeding up reading BAM files is appreciated. (One more trick I want to try next is to extract chromosome info from BAM header file and run them in parallel using multiprocessing/cython_nogil/openMPI or something; still don't know what will work best with Cython.)

    Thanks!
    Last edited by Proteos; 03-26-2018, 01:28 PM.

Latest Articles

Collapse

  • SEQadmin2
    Beyond CRISPR/Cas9: Understand, Choose, and Use the Right Genome Editing Tool
    by SEQadmin2



    CRISPR/Cas9 sparked the gene editing revolution for both research and therapeutics.1 But this system still showed severe issues that limited its applications. The most prominent were the heavy reliance on PAM sequences, delivery limitations, double-stranded breaks that prompt unintended edits and cell death, and editing inefficiency (both in targeting and in knock-in reliability).

    Despite this, “CRISPR helped turn genome editing from a specialized technique into
    ...
    07-31-2026, 11:01 AM
  • SEQadmin2
    Proteomic Platforms: How to Choose the Right Analytical Strategy to Improve Detection and Clinical Applications
    by SEQadmin2


    Proteomics platforms are evolving rapidly, with advances in mass spectrometry and affinity-based approaches expanding what researchers can detect and at what scale. As the field moves toward deeper proteome coverage and clinical applications, scientists face an increasingly complex landscape of tools. This article will explore how researchers are navigating these choices to find the right platform for their work.

    The systematic characterization of the human proteome has
    ...
    07-20-2026, 11:48 AM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by SEQadmin2, 08-06-2026, 07:41 AM
0 responses
14 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 08-03-2026, 10:13 AM
0 responses
31 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 07-31-2026, 02:55 AM
0 responses
41 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 07-24-2026, 12:17 PM
0 responses
26 views
0 reactions
Last Post SEQadmin2  
Working...