FastUniq on linux terminal
1
0
Entering edit mode
9.1 years ago

Dear all,

I am trying to use fastUniq to remove duplicated sequences. I am using the following command on linux platform:

~/FastUniq/source/fastuniq -i <seq1.fq> <seq2.fq> -tq -o fu_1.fq -p fu_2.fq  -c 0

But I got an error. I tried also different versions of –t but did not work…

I found the manual a bit spartan, could you please give some tips on how to write the command?

Thank you

next-gen • 6.5k views
ADD COMMENT
3
Entering edit mode

Questions only stating "I got an error" but give no further details about the nature of the error or specific messages, should be closed as "Not specific enough to answer" in the future, to discourage this.

ADD REPLY
0
Entering edit mode

I am getting: $ Error in open left fastq file @M03595:11:000000000-AG58B:1:1101:16029:1738 1:N:0:26 for read!

ADD REPLY
0
Entering edit mode

If you carefully read my answer bellow, you will find out what went wrong.

ADD REPLY
0
Entering edit mode

Thanks, it worked by providing the list of <seq1.fq> <seq2.fq>

ADD REPLY
1
Entering edit mode

What is the error you are getting?

ADD REPLY
5
Entering edit mode
9.1 years ago
h.mon 35k

FastUniq accepts as input a file listing the fastq file names.

$ ~/FastUniq/source/fastuniq -i list.txt ...

list.txt:

$ cat list

seq1.fq

seq2.fq

See here for a trick to using compressed fastq files.

edit: as venu commented, it is always helpful to copy here the exact error message, instead of just saying "it did not work".

ADD COMMENT

Login before adding your answer.

Traffic: 1913 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6