Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Quality trimming on the 3' end #45

Open
royfrancis opened this issue Oct 30, 2017 · 0 comments
Open

Quality trimming on the 3' end #45

royfrancis opened this issue Oct 30, 2017 · 0 comments

Comments

@royfrancis
Copy link

I start with this (showing only pair 1):
screenshot 2017-10-30 12 22 44

I would like to have to bases in any read below quality 28. So I need all reads with bases below 28 to be trimmed from the 3' end till a base where the quality is >=28.

I used the -q argument like below.

skewer \
-x "dna-adapters_1.fa" \
-y "dna-adapters_2.fa" \
-m pe \
-u \
-n \
-k 6 \
-l 20 \
-t 12 \
-q 28 \
-o "sample" \
$1 $2

I get this (showing only pair 1):

screenshot 2017-10-30 12 22 50

It seems like there is some trimming on the last 5 bases on the 3' end probably from removing 3' adapters. But there are still bases in many reads below quality 28 on the 3' end.

How do I remove these?

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant