Worked perfectly! Thank you.
Hi everyone,
I am attempting to use FastUniq to remove PCR duplicates from my trimmed and paired sequencing reads (150 BP PE done on the Illumina HiSeq X).
I first created an input file:
ls site1_R1_P.fq site1_R2_P.fq > site1_input.txt
Everything looks fine:
$ cat site1_input.txt
site1_R1_P.fq
site1_R2_P.fq
But I am getting the error: "Error in open left fastq file site1_R1_P_qtrim.fq for read!"
I searched through the previous questions that relate to this error but the resolutions to those don't seem to apply - ex. misplaced spaces in the input.txt file (made sure there is no spaces), unusual symbols (tried running without any underscores in case that counted as "unusual").
At a bit of a loss and I bet it's something silly so I appreciate any help!
Thank you!
2 answers
Not a direct answer but I recommend that you use clumpify.sh from BBMap suite for this application (Introducing Clumpify: Create 30% Smaller, Faster Gzipped Fastq Files. And remove duplicates. ). It will be fast, easy to use and will allow you to separate optical dups from PCR dups.
Maybe you can open your "site1_input.txt" with vim,if behind "site1_R1_P.fq" have space except "\n",please delete space and try again!
Log in to answer this question.
If you have
dos2unixcommand available on your machine just pass yourinput.txtfile through it to ensure that line endings are proper.Thank you for the suggestion! I read through the material and will likely go that way.