Handling Paired-End Reads in ShortBRED

Dear biobakery-Team,

I want to process paired-end reads with ShortBRED. To my knowledge, specifying 2 reads is not mentioned in the user manual. However, ShortBRED runs with 2 parsed reads and it was discussed in this thread as well.

How does ShortBRED subsequently calculate hits? Will it be twice as much hits as compared to single-end reads?


1 Like

I’m hoping for additional clarification on paired-end reads, too. How is shortBRED accounting for the pairs when determining reads aligning to marker sequences?

It seems that you would want to assemble your paired-end reads into contigs and then pass these contigs into shortBRED, but based on the thread plicht references and based on information I’ve found elsewhere, it sounds like biobakery team recommends just passing R1 and R2 in at the same time.

Any suggestions/ideas here?

Hi all,
I wanted to revitalize this topic as this is something I am curious about as well. How should someone interpret the output of two paired-end reads run together? Would it be a different outcome than if the paired-end reads were assembled?
Thank you for your input!