mothur

unique.seqs error


#1

When running unique.seqs command, I am getting the error below:

mothur > unique.seqs(fasta=uic.trim.contigs.good.fasta )


[ERROR]: You already have a sequence named M02127_258_000000000-B4RL2_1_1101_16099_2069 in your fasta file, sequence names must be unique, please correct. This error is repeated for all of the sequences that I will not add to this post. This seems odd that this command would require unique sequences.

The details for the mothur version used are:
Linux version
Using ReadLine
Running 64Bit Version
mothur v.1.39.5

Thanks for your assistance,
Stephanie


#2

What are the commands that you are running upstream of this command?


#3

The previous commands were:

make.contigs(file=uic.files, processors=8)

summary.seqs(fasta=current)

screen.seqs(fasta=current, group=current, summary=current, maxambig=0, maxlength=275)

summary.seqs(fasta=current)

get.current()

unique.seqs(fasta=current)
Using uic.trim.contigs.good.fasta as input file for the fasta parameter.

Best,
Stephanie


#4

I also have run these same set of commands on the newest version for Mac using 1 processor, in which I still get the same error, except for only a large subset of sequences, and not every sequence.


#5

Could you have duplicate files in your uic.files file? If not, can you send your input files and log file to mothur.bugs@gmail.com?


#6

I have been getting this same error recently. When I ran my data using v 1.39.5 the error occurred after make.contigs and identified 68 repeated sequence names. However, when I updated to v.1.40.5 the error didn’t occur until the unique.seqs command and identified thousands of repeated sequence names.


Thank you for the help,

Kelsy


#7

Could you send your input files and logfile to mothur.bugs@gmail.com?


#8

Hi there,

Has this issue beeen fixed? I’m getting the same error with v. 1.41.1 after running unique.seqs. These are the code and errors:

make.contigs(file=dio.files, checkorient=T)

summary.seqs(fasta=dio.trim.contigs.fasta)

screen.seqs(fasta=dio.trim.contigs.fasta, group=dio.contigs.groups, summary=dio.trim.contigs.summary, maxambig=0, maxlength=301, maxhomop=8)

summary.seqs(fasta=dio.trim.contigs.good.fasta)

unique.seqs(fasta=dio.trim.contigs.good.fasta)
[ERROR]: You already have a sequence named M03613_79_000000000-BN6L2_1_2109_20418_2415 in your fasta file, sequence names must be unique, please correct.
[ERROR]: You already have a sequence named M03613_79_000000000-BN6L2_1_1101_5910_25285 in your fasta file, sequence names must be unique, please correct.
[ERROR]: You already have a sequence named M03613_79_000000000-BN6L2_1_2103_11121_6885 in your fasta file, sequence names must be unique, please correct.
[ERROR]: You already have a sequence named M03613_79_000000000-BN6L2_1_2110_3694_9213 in your fasta file, sequence names must be unique, please correct.
[ERROR]: You already have a sequence named M03613_79_000000000-BN6L2_1_1110_12330_22908 in your fasta file, sequence names must be unique, please correct.
[ERROR]: You already have a sequence named M03613_79_000000000-BN6L2_1_1111_24820_9625 in your fasta file, sequence names must be unique, please correct.
[ERROR]: You already have a sequence named M03613_79_000000000-BN6L2_1_2111_22903_5708 in your fasta file, sequence names must be unique, please correct.
[ERROR]: You already have a sequence named M03613_79_000000000-BN6L2_1_2114_7703_9852 in your fasta file, sequence names must be unique, please correct.
[ERROR]: You already have a sequence named M03613_79_000000000-BN6L2_1_1108_15091_17664 in your fasta file, sequence names must be unique, please correct.
[ERROR]: You already have a sequence named M03613_79_000000000-BN6L2_1_2108_5225_20067 in your fasta file, sequence names must be unique, please correct.
[ERROR]: You already have a sequence named M03613_79_000000000-BN6L2_1_2102_12063_5957 in your fasta file, sequence names must be unique, please correct.

I indeed checked the output after make.contigs (dio.trim.contigs.fasta) and there are indeed some duplicates. However, there are no duplicate samples in my dio.files file.

Thanks in advance for your help.


#9

Following this issue further, it seems that there are output differences between Mothur versions.

I ran the following command using raw fastq files (not fastq.gz files):

make.contigs(file=dio.files, checkorient=T)

When I run this command in v. 1.39.0, the number of sequences in the dio.trim.contigs.fasta and dio.contigs.groups files is the same (in my case: 8,131,113). Now, if I run the same command in v. 1.41.1, the number of sequences in the dio.trim.contigs.fasta file is 8,131,113, whereas it outputs 8,129,290 sequences in the dio.contigs.groups file. Clearly, there is something wrong with this command in v. 1.41.1.

I went into the SOP and next ran screen.seqs with the two Mothur versions, using the outputs of v. 1.39.0 above (because those of v. 1.41.1 were wrong), as follows:

screen.seqs(fasta=dio.trim.contigs.fasta, group=dio.contigs.groups, summary=dio.trim.contigs.summary, maxambig=0, maxlength=301, maxhomop=8)

However, both v. 1.39.0 and v. 1.41.1 gave me wrong dio.contigs.good.groups files. Let’s check the number of sequences being trimmed and clean:

dio.trim.contigs.bad.accnos: 4,666,759 (I know, too many crapy sequences, but let’s forget that for the sake of this discussion).

dio.trim.contigs.good.fasta: 3,464,354

The sum of these two is 8,131,113, the correct number of total sequences after make.contigs. I expected then 3,464,354 sequences in the dio.contigs.good.groups file, but instead I got 3,465,917 with v. 1.39.0 or 3,464,094 in v. 1.41.1. In any case, I can’t go further with the SOP since in the unique.seqs command it always tells me that the fasta and groups files differ.

Any idea of what’s going on here?


#10

The group file contains only unique names, but it looks like your fasta file contains duplicates. Can you try this?

mothur > list.seqs(fasta=dio.trim.contigs.good.fasta) - list all the unique names in the fasta file
mothur > get.seqs(fasta=current, accnos=current) - remove any duplicate names, you will see a bunch of warnings
mothur > unique.seqs(fasta=current) - should process without duplicate sequences error

If the resulting fasta file does not match the group file could you send the dio.trim.contigs.fasta and did.contigs.groups files to mothur.bugs@gmail.com?