I am a little bit confused about the first intron, just using the first one in the introns.bed of each gene or we should also use the *.genepred file and consider more things about it.
Thank you so much, your code is pretty valuable.
I am looking forward to your reply.