Hi, I am trying to use "plotAvgProf". Is it possible to generate the same plot that shows peak bindings with respect to transcription START SITE and STOP SITE both?

This is the current command I am using:

plotAvgProf2(peak1[], TxDb=txdb, upstream=3000, downstream=3000,
               xlab="Genomic Region (5'->3')", ylab = "Read Count Frequency", resample = 1000)

Thanks -S

# include your problematic code here with any corresponding output 
# please also include the results of running the following in an R session 

sessionInfo( )