Reads number in bam/sam resluts

[guangzhaocs]

In my original pod5 files, I have 63886 reads. I run dorado basecaller and get the bam file. Then I use samtools to convert it to sam file. I counted the reads number in the sam file and found that contains 67690 reads. I do not know why.

I compared the read id in the original pod5 files and the sam file, and found that the sam file has some reads that do not exist in the pod5 file.

Run environment:

• Dorado version: 0.5.3
• Dorado command: dorado basecaller models/rna002_70bps_hac@v3/ multi_pod5 --estimate-poly-a > res.bam
• samtools view -h -o res.sam res.bam

[HalfPhoton]

@guangzhaocs - These new read ids from the read splitting of concatemers.