Question: Visualizing & plotting variant reads within BAM files programmatically?
gravatar for steve
3.2 years ago by
United States
steve2.3k wrote:

I am looking for essentially the functionality mentioned here; I want to automatically create snapshots of images of reads within a BAM file at specified loci. I know that IGV can do this, but in my experience the results from scripting IGV for this are lackluster since there can be a lot of variability as far as the number of reads on screen and resolution needed for each variant. I also saw that samtools tview has the same functionality. The answers in that previous question are a bit old now, I was hoping that something new might be available? Maybe a Bioconductor package? My searches so far have not found anything quite like this. My immediate needs are for snapshots / static plots of the reads, but ultimately I think this would be a nice application for something like Plotly to enable interactive browsing as well.

alignment • 1.4k views
ADD COMMENTlink modified 3.2 years ago by Pierre Lindenbaum124k • written 3.2 years ago by steve2.3k

there is script for this here:

ADD REPLYlink written 5 weeks ago by steve2.3k
gravatar for Pierre Lindenbaum
3.2 years ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum124k wrote:

I've written a BAM2SVG:

enter image description here

and a BAM2raster:

enter image description here

ADD COMMENTlink written 3.2 years ago by Pierre Lindenbaum124k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1872 users visited in the last hour