Fastest way to use pysam view in order to get the reads in different regions
Dear all,
I am trying to retrieve all reads in several different regions for several bam files. I am writing my code in python, using the pysam library as follows,
pysam.view("-F0X100", bam_file, region_to_query)
Everything works, but I am concerned about the time it takes to process each bam file. It seems that the I/O takes a long time since I am querying each region at a time.
Would you have any suggestions on how I can speed up the querying of each bam file?
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Pysam and CRAM file
written by MarVi •
-
Retrieve initial query in getBM()
written by helia.ziaei •
-
Quicker way to know if a SNV is involved in a mismatch alignment
written by MarVi •
-
Extracting CIGAR string from a specific region of a read in a bam file using pysam?
written by kiki4133 •
-
Pysam: Is there a way to check if there are inserts when looking at read overalps.
written by ahmed.kidwai •
-
Get the regions mapped by the same reads
written by Corentin
-
Problem with samtools view - retrieving reads mapped to a region
written by MarVi •
-
Gatk Indel Realignment Taking Forever! Help
written by Ds89 •
-
Fetching Read By Its Id From A Bam File With Pysam In Python
written by User 9996
-
How To Get Number Of Reads In Bam File Efficiently In Python?
written by User 9996
Have you tried "bedtools intersect"? I guess you can easily intersect each bam file with the bed file (regions) in parallel by bedtools, then merge the output bam files finally.