pysam.pileup give empty results, any clue?

I'm trying to use pysam's pileup() to obtain reads but gives empty results. While for same region pysam's fetch() gives correct results, any clue for this? Thanks in advance

screen shot for pileup not working

python pysam

• 879 views • link

updated 3.6 years ago by barslmn •

• written 3.6 years ago by PeterWu •

It looks like you need to press enter to run that block.

• 0 views • link

written 3.6 years ago by barslmn •

0 answers

No answers yet.

More posts like this

Interpret IGV
written 1.6 years ago by najim3 •

Hello, Biostars community. Does anybody know the genetic event that results in the reads' alignment in the example below? Note the green reads, which usually …
ChIP seq BamCoverage
written 2.8 years ago by lonekaisar72 •

Dear All, We are trying to analyse some ChIP seq data sets on one of the repeated regions of Genome (Ribosomal DNA). Recently one paper …
Split reads along the genome in my samples
written 2.9 years ago by paulaotero.sanchez •

Hello, I am working with DNA sequencing libraries and I have found out that there are a lot of split reads randomly distributed (not clustered) …
Top 3 gene clusters using MCODE in Cytoscape - Protein-protein interaction analysis
written 3.7 years ago by bioinformatics •

Hi everyone, I am trying to find the top 3 gene clusters using MCODE in Cytoscape by performing cluster analysis. I have imported STRING data …
Exomiser not returning any output
written 4.1 years ago by Taku •

Hi everyone, I am trying to run Exomiser on my VCF files but I keep getting empty output files. If I try the test Pfeiffer.vcf …
Covarage X not per base coverage. As a whole
written 4.8 years ago by Edmond •

Hello, Before anyone marks this post as "duplicated" from [this post][1]. This post and plotCoverage from deeptools, gives us the mean coverage ( or coverage) …
PPI-Network igraph for tissue-specific data
written 4.9 years ago by bi.nb21 •

I am new to bioinformatics in general and would really appreciate some help with the project Im working on. My data of protein-protein interactions is …
Stringtie error : discarding overlapping duplicate miRNA feature
written 5.1 years ago by Hojn •

Hi everyone! I am following a pipeline to extract new variants for A. thaliana. I have already conducted this part of the process but I …
Separate colors in immunohistochemical staining
written 5.4 years ago by salvatore.raieli2

Hi everyone, I was trying to replicate this tutorial: https://scikit-image.org/docs/dev/auto_examples/color_exposure/plot_ihc_color_separation.html the results should be this: ![correct results ][1] I used jupyter, spider and google colab, …
How would you trim this RNA-seq data
written 7.7 years ago by williamsbrian5064

Hi, I am working on an RNA-Seq project. We just got back some sequencing and I had a couple questions about it. I figure someone …