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.

0 answers

No answers yet.

Log in to answer this question.