Pysam usage problem · Biostar

This is a test version of Biostars. For the public version, visit https://www.biostars.org.

Pysam usage problem

Hi BioStars,

I am using pysam to count read from bam file. I use pysam.AlignmentFile to read bam file and pysam.AlignmentFile.count to do read count, but it requires bai file which is an index should in the same path where bam file is. However I have no writing permission in bam file path. Is it possible to input index bai file from other path, I mean save bam file and its bai file in different path? Thanks very much if some guys can help this!

software error

• 2,331 views • link

updated 7.3 years ago by Biostar

• written 9.6 years ago by Rongbin •

1 answer

Easiest is probably to create a symbolic link to the bam (ln - s) in a directory where you can write, or /tmp/, and write the index there. You can probably fix this using the os module

• 0 views • link

updated 7.3 years ago by Ram

• written 9.6 years ago by WouterDeCoster

Log in to answer this question.

More posts like this

Nextflow: how to pass Bam and Bam index as Input Channel?
written 3.1 years ago by Eliveri •

I would like to pass in bam files `pair_id.sorted.bam` and their corresponding index files `pair_id.sorted.bam.csi` into a nextflow workflow. However I am having trouble passing …
[E::idx_find_and_load] Could not retrieve index file for Singularity /NextFlow
written 3.1 years ago by Kiran

Hi I am writing my first code is dsl2 Nextflow using Singularity container as below main.nf process pbc_varicall { publishDir "/data/shared/clinical/LongRead/Data/resources/" container 'docker://google/deepvariant:1.5.0' input: path …
pysam error when reading .bam file ValueError: file has no sequences defined (mode='rb') - is it SA…
written 5.5 years ago by daewowo •

Error: f = pysam.AlignmentFile("SRA_sorted.bam","rb") File "pysam/libcalignmentfile.pyx", line 991, in pysam.libcalignmentfile.AlignmentFile._open ValueError: file has no sequences defined (mode='rb') - is it SAM/BAM format? Consider opening with …
suppress pysam index warnings.
written 6.8 years ago by b10hazard •

When I read a bam file with pysam I get a lot of these warnings... [W::hts_idx_load2] The index file is older than the data file: …
pysam alignment error
written 6.8 years ago by rthapa •

Hi, I am trying to read bam file using a command: samfile = pysam.AlignmentFile(os.path.join(path,bam),'rb') but I keep on getting an error: Traceback (most recent call …
pysam + multiple_iterators = Type Error during multiprocessing
written 8.8 years ago by QVINTVS_FABIVS_MAXIMVS

I would like to parse a BAM file in parallel using pysam and multiple_iterators Here is my code import pysam import sys from multiprocessing import …
secondaryFiles for bam: how to handle .bai or/and ^.bai ?
written 9.0 years ago by liuxf09 •

Bam file could have two different file extensions depending on the tools. Given bam file test.bam, - test.bam.bai secondaryFiles = [.bai] - test.bai secondaryFiles = …
Reading CRAM files in pysam
written 10.6 years ago by QVINTVS_FABIVS_MAXIMVS

I'm new to CRAM format and I am wondering if I can read a CRAM file in pysam or samtools. Or generally speaking how do …
Pysam cannot find index
written 11.0 years ago by eric.kern13

Hi Biostars, This is my first question on the site, so thanks for all the help I've gotten from lurking, and thanks for the help …
Java Library For Reading Bai Files And Understanding Index With Respect To Bam
written 14.6 years ago by User 9906 •

<p>Hi,</p> <ol> <li><p>Is there a library in java which allows me to read BAM index files(.bai) ? </p></li> <li><p>I want to understand how the index …