0
$\begingroup$

I am using sickle package from bioconda to trim my reads from the Illumina Miseq.

When I tried following code, it failed to run

sickle pe -t illumina -f file1   -r file2  -o output1  -p output2  -s combined_output -q 27 -l 25  

When I replaced illumina with sanger, it works. But my reads are from the illumina. Could anyone clarify if using this option is acceptable ? Also, what is the optimum value of -q and -l

$\endgroup$
2
$\begingroup$

It'd be rare to run into Illumina encoded quality scores any more, I doubt there are machines running that still produce those. Everything produced at the moment uses Sanger encoded (phred + 33) scores.

There is no optimum value for -q and -l, though reasonable values would typically be <=5 for -q and ~5 for -l.

$\endgroup$
0

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.