Pysam--how to check for empty results in fetch(region) of bam file
I'm reading a bam file using Pysam, but I couldn't find one file for the whole genome, so I'm iterating over the chromosomes. When I fetch a region by chromosome, only the file-specific chromosome will get results (ie if I fetch "21" in the chrom 1 file, I shouldn't get any results). Is there a way to determine at the fetch that there's no data so I can move on? Thanks.
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Invalid Contig error when using bamfile.fetch
written by bxxxmxxxx •
-
Discrepancy in pysam fetching region
written by mhasa006
-
Repeats that affects NIPT results
written by khanhlpbao •
-
Pysam pileup method gives more columns than expected
written by Didrik.Olofsson •
-
Parsing tabix with pysam
written by Alex Reynolds
-
pulling out a single pseudochromosome from my vcf file
written by truebeliever24 •
-
Changing Chromosome Notation in Bam Files to Include Sample ID
written by dthorbur
-
Is there a way to fetch genomic sequences at given coordinates without downloading fasta files?
written by ericbrenner •
-
Pysam - Different amount of mapped and aligned reads
written by AlchemistMoz •
-
Reading CRAM files in pysam
written by QVINTVS_FABIVS_MAXIMVS
I found count(region), which does more or less what I want (ie it tells me if the region returns nothing), but running it on a non-empty region takes a while and I'm hoping there's a better way.
The most simply way would be to name the files that way, that you know which chromosome is in.
Another way would be to use pysam's get_index_statistics() which will give you the number of mapped, unmapped and total reads per contig.