FelixKrueger / Bismark

A tool to map bisulfite converted sequence reads and determine cytosine methylation states
http://felixkrueger.github.io/Bismark/
GNU General Public License v3.0
366 stars 101 forks source link

Can deduplicate_bismark (the component of bismark) script distinguish single-end/paired-end BAM/SAM Bismark files directly ? #659

Closed changchuanjun closed 4 months ago

changchuanjun commented 4 months ago

@FelixKrueger, In deduplicate_bismark script instruction part, It persents can detect single-end/paired-end BAM/SAM Bismark files automatically. image

But when I ran this script using paired-end BAM/SAM Bismark files,It detect it a single-end BAM/SAM Bismark file. It is very weird. image

So I am confused that whether deduplicate_bismark (the component of bismark) script distinguish single-end/paired-end BAM/SAM Bismark files directly.Or when I ran deduplicate_bismark script with the -p parameter?

FelixKrueger commented 4 months ago

Which version of deduplicate_bismark are you using?

changchuanjun commented 4 months ago

@FelixKrueger Deduplicator Version: v0.24.2 image

changchuanjun commented 4 months ago

@FelixKrueger Deduplicator Version: v0.24.2 image

FelixKrueger commented 4 months ago

Ok thanks. It appears that the top message should probably just say Bismark output files, and lose the '--single-end'. I'll have that changed.

The last line of the screen-shot above says: "Treating file as paired-end file". So that is good, this is the important bit.

changchuanjun commented 4 months ago

OK. So Can deduplicate_bismark script distinguish single-end/paired-end BAM/SAM Bismark files directly? Or when I ran deduplicate_bismark script with the parameter --paired ?

FelixKrueger commented 4 months ago

either is fine. The distinction is that the auto-detection work on a per file basis, whereas -p works for all files supplied.