Question:
I’m trying to make a single gene coverage plot like the one shown here.
I have already plotted the coverage but I still need to insert the exons and introns in the corresponding coordinates. Thank you! My code:

Answer:
Without a minimal reproducible example (MRE) containing some input data and expected output it’s difficult to know how to help you. For instance, here is an MRE with some ‘fake’ data:
Created on 2022-06-16 by the reprex package (v2.0.1)
Would this approach be suitable for your single gene coverage plot? If not, what would you want to change?
If you have better answer, please add a comment about this, thank you!