This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Plot reads from .BAM file above reference sequence in R?

I want to plot something like this:

enter image description here

This was produced with Gviz. My problem with Gviz is the empty space added at the bottom of the plot (below the reads). I posted a question about that here.

Do you have any suggestions on how I can plot reads aligned to a reference sequence in R, and not use Gviz?

Edit: The observant might have noticed that the "empty space problem" with the Gviz plot was solved in the post I linked to. I'm still curious if there are other solutions, though.

r plot
• 6,521 views • link
written by stianlagstad stianlagstad

I surmise it might be easier to do it manually in inkscape or gimp or semi-manually with imagemagick convert. Surely it is possible in pure R but it might be a pain... (I'd be happy to be proved wrong though...)

• 0 views • link
written by dariober dariober

Thank you for the answer! For a single case it would be easier to do it manually, but I have to do this automatically in R :) I have a solution that works with Gviz, but that empty space is giving me problems when I export the plot in scalable formats.

• 0 views • link
written by stianlagstad stianlagstad

Thank you for the answer! I have looked at the Sushi manual before but I couldn't see anything that would help me. Could you point me in the right direction?

• 0 views • link
written by stianlagstad stianlagstad

Thank you, but rest assured that I've done more than a few searches and I've looked at more than a couple of packages before coming here to ask the question. If you have a specific suggestion that would be helpful :)

• 0 views • link
written by stianlagstad stianlagstad

Just for the time being, if you want you can use the illustrator trick, if you export the pdf plot from Gviz. Good Luck!

• 0 views • link
written by Sukhi Singh Sukhi Singh

0 answers

No answers yet.

Log in to answer this question.

More posts like this