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

problem for merging 2 pairs of reads by "SeqPrep" tool #45

Open
modarzi opened this issue Oct 30, 2019 · 2 comments
Open

problem for merging 2 pairs of reads by "SeqPrep" tool #45

modarzi opened this issue Oct 30, 2019 · 2 comments

Comments

@modarzi
Copy link

modarzi commented Oct 30, 2019

For Data processing In GSM1827405 author wrote that "Overlapping reads were merged using SeqPrep and split in half with a Awk routine. Paired-ends with one read shorter than 30 bp were discarded."
I run below codein my linux:
seqprep -f SRR2017741_1.fastq.gz -r SRR2017741_2.fastq.gz -1 SRR2017741_1_seqprep.fastq.gz -2 SRR2017741_2_seqprep.fastq.gz
and I get below result:

Pairs Processed: 37404976
Pairs Merged: 0
Pairs With Adapters: 75653
Pairs Discarded: 28621
CPU Time Used (Minutes): 55.881797

I don't know why Pairs Merged is 0? I appreciate if you share your comment with me for finding my problem and completing my seqprep code for this sample.

@aersoares81
Copy link

Your command line is missing "-s".

Something like "seqprep -f SRR2017741_1.fastq.gz -r SRR2017741_2.fastq.gz -1 SRR2017741_1_seqprep.fastq.gz -2 SRR2017741_2_seqprep.fastq.gz -s SRR2017741_merged.fastq.gz" should probably work.

@Ellmen
Copy link

Ellmen commented Feb 23, 2021

As a newcomer, I think this is kind of a funny UI decision given that the first line in the description is "SeqPrep is a program to merge paired end Illumina reads that are overlapping into a single longer read"

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

3 participants