Entering edit mode
1 day ago
baibhu1234
▴
20
Hello everyone,
I have identified m6A sites in DRACH context using plants DRS dataset using dorado + Modkit pipeline.
Now i am trying to get the Differentially modified sites between conditions using a number of modified reads/total number of reads at that position. I am making a master list of sites in two conditions by taking union of all sites in them.
This way some sites are there which are having no modification in one condition but some in other. so there is #N/A showing for that condition (See attached screenshot for reference). To have comparison i replaced the #N/A with 0 and then added 1 in all the count or reads.
is this a correct way to do the analysis?
or it will lead to any biases towards to some sites?
