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

How to interlace pair-end reads via seqtk and make the headers end in /1 and /2 corresponding to its read orientation (R1 and R2, respectively)? #56

Open
yoyohashao opened this issue May 20, 2015 · 1 comment

Comments

@yoyohashao
Copy link

Hello, sorry for the empty previous post cuz I tapped Enter key so fast...
I am preparing datasets for kmernorm which requires the files look like
@seq_1/1
AAAAAAAACCCCCCCTTTTTTTTTGGGGGGGG
+
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
@seq_1/2
AAAAAAAACCCCCCCTTTTTTTTTGGGGGGGG
+
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
@seq_2/1
AAAAAAAACCCCCCCTTTTTTTTTGGGGGGGG
+
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
@seq_2/2
AAAAAAAACCCCCCCTTTTTTTTTGGGGGGGG
+
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&

And a bioinformatician told me seqtk could archive this. But after screening the examples of seqtk, I didn't find the command. Could someone tell me how to do this in seqtk? Many thanks!

@tseemann
Copy link

To interleave PE reads do this:

seqtk mergepe R1.fq,gz R2.fq.gz | gzip > Interleaved.fq.gz

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

2 participants