hildebra
commented
3 months ago Hey, could you check how many reads passed the initial sdm filter? This should be shown on the console, or in the Log dir (otuS_LSU_out/LotuSLogs/) there should be sdm named files. My first guess would be that no read was dereplicated, because no read passed the quality controls (or too few reads). If this is the case: you need to further lower the qual filter, and the log file would guide which qual filter caused most reads being removed. Further, you can try to lower the dereplication parameters, by setting "-derepMin 0" or similar. However, note that LotuS2 was never programmed to work with ONT reads, ie many assumptions of the read clustering will be broken by the (usually) really low quality of ONT reads. hth, Falk
SebaZambrano
4less
I'm trying to cluster with VSEARCH and assign taxonomy to a set of samples using partial LSU amplicons (in fasta format) PREVIOUSLY extracted from ONT long reads using ITSx. I used the following parameters:
lotus2 -m LSUmap.txt -i LSU_lotus/ -o lotuS_LSU_out/ -CL VSEARCH -amplicon_type LSU -ITSx 0 -s sdm_ONT_LSSU.txt -refDB UNITE -taxAligner blast -tax_group fungiThe pipeline stopped due to a dereplication error (it failed to identify unique sequences). The error was reported as:
The sdm dereplicated output file was either empty or not existing, aborting lotus. lotuS_LSU_out//tmpFiles//derep.fasNote that I modified the sdm file, in particular the max and min length parameters, nothing else. I'm attaching the sdm file. sdm_ONT_LSSU.txt
Any help would be very much appreciated!