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

Single reads used in error correction but not in assembly #214

Closed
ghost opened this issue Nov 14, 2019 · 1 comment
Closed

Single reads used in error correction but not in assembly #214

ghost opened this issue Nov 14, 2019 · 1 comment

Comments

@ghost
Copy link

ghost commented Nov 14, 2019

I made a test and processed some reads (paired and unpaired reads; (PF = Paired_Forward, PR= Paired_Reverse, UF and UR unpaired forward and reverse) with spades in error_correction mode without assembly and compared the number of reads to the number of corrected reads produced by Unicycler:

Using SPAdes
"spades.py --only-error-correction -1 PF.fastq.gz -2 PR.fastq.gz -s UF.fastq.gz -s UR.fastq.gz --meta -t 6 -m 32 -o error_correction"

I got (counted the number of lines with zgrep):
PF.fastq.00.0_0.cor.fastq.gz:4099780
PR.fastq.00.0_0.cor.fastq.gz:4099780
P_unpaired.00.0_0.cor.fastq.gz:34392
UF.fastq.00.0_1.cor.fastq.gz:211720
UR.fastq.00.0_2.cor.fastq.gz:7672

Unicycler (I joined the two files with unpaired reads before running unicycler).
"unicycler -1 PF.fastq.gz -2 PR.fastq.gz -s SE.fastq.gz -o assembly -t 6"

I got (counted the number of lines with zgrep):
corrected_1.fastq.gz:4099780
corrected_2.fastq.gz:4099780
corrected_u.fastq.gz:34392

Based on this, unicycler ignored all original unpaired reads although I specified them and they were listed in the command while running unicycler.

@rrwick
Copy link
Owner

rrwick commented Jan 22, 2022

I have decided to remove the read-error-correction step from Unicycler's pipeline in the newest version (v0.5.0), as I don't think it's helpful for a QCed read set. Since this issue no longer applies, I'm going to close it now.

@rrwick rrwick closed this as completed Jan 22, 2022
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