m6Anet result differences in same bam file.

[DayeaPark]

Hello.

I am trying to analyze m6A sequencing data using m6Anet. I currently generated data from the Direct RNA sequencing kit. I used guppy and minimap2 in MinKNOW software for base-calling and alignment.

I separated one bam file into two files and run m6Anet with these three bam files. However, I realized the results (n_read count, modi_ratio, probability etc...) from these three bam files are different. I expect at least matched n_read numbers. (If file_A.bam and file_B.bam from original_file.bam, the sum of n_read of certain genes from file_A and file_B has to be the same as the n_read from the original_file.

Do you have any idea why it makes a discrepancy?

Thank you.