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?
• 894 views
•
link
0 answers
No answers yet.
Log in to answer this question.
As I recall
sabretransfers 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