-
Notifications
You must be signed in to change notification settings - Fork 110
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
Make.contigs - name mismatch BUG #167
Comments
http://www.mothur.org/forum/viewtopic.php?f=3&t=4087 make.contigs(file=stability.files, processors=8) I get an error saying that a sequence is in my R1 file that is not in my R2 file and to remove these sequences using remove.seqs. The error message is from the splitting process. Not sure if it matters but the sequence names contain '-' characters. |
@pschloss, how do you want to handle this bug? Coding around it will slow the command. We currently look for an exact match, then look for a match with the last char dropped to account for mismatches like: /**
If both of these fail we then start the skipping reads make matches and process. This seems like a rare case. Should we ignore it or should I dive in? |
Can you remind me what they should look like? |
Completed with b55bcc8. |
http://www.mothur.org/forum/viewtopic.php?f=4&t=3892&p=11629#p11629
@A1
CCAGC
+
ABBCB
@a2
CCAGC
+
CCCCC
@A3
CCAGC
+
CCCCB
$ mothur '#make.contigs(ffastq=not_working_R1.fastq, rfastq=not_working_R2.fastq, processors=2)'
Using 2 processors.
Making contigs...
3
[WARNING]: name mismatch in forward and reverse fastq file. Ignoring, A4.
1
However, changing the sequence names in the two files to "A1", "B2", "C3", "D4", "E5" allows the file to be properly processed. I can provide working and non-working example files if needed.
The text was updated successfully, but these errors were encountered: