Query about adapting ployAlength script for gene-level polyA length analysis

[aaaaaaaaaayong]

Hi there,

First, thank you for sharing the ployAlength script—it’s been incredibly useful for analyzing polyA length differences at the transcript level!

I’m wondering if the current script could also be applied to analyze polyA length variation at the gene level (e.g., aggregating results across transcripts of the same gene). If so, are there recommended parameters or modifications to accommodate this? If not, would you have any suggestions for adapting the workflow?

I’d greatly appreciate any guidance! Happy to provide more context if needed.

Thanks again for your excellent tool and for making it openly available!