This is a test version of Biostars. For the public version, visit https://www.biostars.org.
pysam: Use read ID to access location

I know the read_ids of my reads, but I want to know, i.e. "chr1:100-120". Therefore I'm trying to access the read from a bam/sam file using the read id but I have not had any luck with pysam so far. Does anyone know how I could do this?

bam python pysam
• 2,813 views • link
updated by Ram Ram • written by bd2810 • bd2810

It's always clearer with an example. Do you want to access the sequence of your read using chr1:100-120 or the possible gene at this position ?

Could you detailed this please :

I want to know the location

• 1 views • link
written by Bastien Hervé Bastien Hervé

1 answer

You can't do that, you can only iterate over all of the reads and check if they match your read_ids. Note that this is not a limitation of pysam, it's implicitly not meant to work in the BAM and CRAM formats.

• 1 views • link
written by Devon Ryan Devon Ryan

Log in to answer this question.

More posts like this