m6A sites identified on transcript basis

[baibhav-bioinfo]

Hello,
I am using DRS dataset to identify m6A sites in plant dataset using dorado+modkit tools.

If i want to get the m6A sites identified on the transcript locations, can i align the unaligned basecalled output bam files to the transcript.fasta instead of genome.fasta. and then use modkit to get bedmethyl files of m6A sites and do downstream analysis.

and then later find out the genomic locations of the sites for representation and visualisation etc.

[ArtRand]

Hello @baibhav-bioinfo,

Two options I can think of off-the-top are:

1. Align to the transcriptome and the genome, then use their pileups for individual tasks.
2. Align to the genome, then make a BED of the transcripts and use the --include-bed argument to modkit pileup so that you only have m6A positions included. Of course, you'll have to make sure to get the BED file right so that you get the isoforms you want.

You could convert transcript coordinates to genome coordinates with LiftOver or a similar tool, but Modkit doesn't have any such functionality.