pysam: Use read ID to access location
1
0
Entering edit mode
4.1 years ago
bd2810 • 0

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?

pysam python bam file • 1.6k views
ADD COMMENT
0
Entering edit mode

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

ADD REPLY
1
Entering edit mode
4.1 years ago

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.

ADD COMMENT

Login before adding your answer.

Traffic: 1758 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6