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

seqkit seq -s option is breaking down the single line sequence #132

Closed
ghost opened this issue May 16, 2020 · 1 comment
Closed

seqkit seq -s option is breaking down the single line sequence #132

ghost opened this issue May 16, 2020 · 1 comment

Comments

@ghost
Copy link

ghost commented May 16, 2020

seqkit version: seqkit v0.12.0

Input fasta:

>419759601_1
GTCTTGCATCTATAGTATTTTTACCCTGGTGGATATCTCTATCATTTCAAAAAAGTCTGGAATCTTGGGTTACTGATTCGTGGAATATTAAACAATCTGCAACTTTTTTGAATGATATTCTAGAAACTAGTCTTCTAGAAAAATTCAACGA

issue:

  1. Input fasta copy/pasted above is a single line fasta.
  2. Seqkit -s is supposed to print sequence as a sequence without changing the input format. Otherwise, user has to type in seqkit -sw0 for extracting sequences in a single line from a flattened fasta file.
  3. Seqkit -s is printing sequence in 3 lines (probably at 60 characters per line- not sure of this) as below:
$ seqkit seq -s test.fa 
GTCTTGCATCTATAGTATTTTTACCCTGGTGGATATCTCTATCATTTCAAAAAAGTCTGG
AATCTTGGGTTACTGATTCGTGGAATATTAAACAATCTGCAACTTTTTTGAATGATATTC
TAGAAACTAGTCTTCTAGAAAAATTCAACGA

To check if sequence is indeed is broken into 3 lines, into three lines, following code was run:

seqkit seq -w 0 test.fa | seqkit seq -s | sed -n '3p'
TAGAAACTAGTCTTCTAGAAAAATTCAACGA

Third line is TAGAAACTAGTCTTCTAGAAAAATTCAACGA

Expected behaviour:

$seqkit seq -s test.fa
GTCTTGCATCTATAGTATTTTTACCCTGGTGGATATCTCTATCATTTCAAAAAAGTCTGGAATCTTGGGTTACTGATTCGTGGAATATTAAACAATCTGCAACTTTTTTGAATGATATTCTAGAAACTAGTCTTCTAGAAAAATTCAACGA

sequence should be printed without gaps, in a single line.

@shenwei356
Copy link
Owner

It makes sense, I'll fix this.

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