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

Empty sequence in output file #36

Open
Isoobel opened this issue Jun 9, 2016 · 0 comments
Open

Empty sequence in output file #36

Isoobel opened this issue Jun 9, 2016 · 0 comments

Comments

@Isoobel
Copy link

Isoobel commented Jun 9, 2016

Hello everyone,

I ran a short test on skewer and I got an empty sequence in the output file.
I know this case is not supposed to happen since I tried to remove an adapter supposed to be located at the read head, but in fact it is at the tail.
But maybe it would be nice remove the header too?

Here is the input file I used to test the adapters filtering:

seq1
AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAATTTTTTTTTTTTTTTTTTTT
seq2
AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAATTTTTTTTTTTTTTTTTTTTAAAA
seq3
AAAAAAAAAAAAAAAAAAAAAAATTTTTTTTTTTTTTTTTTTTAAAAAAAAAAAAAAAAAAAAAAAAA
seq4
TTTTTTTTTTTTTTTTTTTTAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA

I launched skewer with the following command line:
skewer -x TTTTTTTTTTTTTTTTTTTT -1 -l 0 -k 15 -t 16 -m head input.fasta >output.fasta

Here is the output fasta file I obtained:

seq1

seq2
AAAA
seq3
AAAAAAAAAAAAAAAAAAAAAAAAA
seq4
AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA

Thanks for your help!

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