This is a test version of Biostars. For the public version, visit https://www.biostars.org.
bioconda sabre reads shorter

I have a bit of a problem. im currently building a snakemake pipeline where i use sabre to demultiplex data. when i look at the length of the reads before demultiplexing, the length is 60, while the reads in the files of reads that match with a barcode are all 52 long. the reads in the -u file are still sixty. i use the pe (paired end) option and the only other parameter i give is -m 1. can anyone tell me why this might happen?

snakemake sabre

As I recall sabre transfers the barcode (in your case it seems to be 8 bp) from the main read to the header. So reads that did not contain the barcode were not trimmed and thus remained full length.

Thank you for your awnser, it solved my problem

0 answers

No answers yet.

Log in to answer this question.