compmetagen / micca

micca - MICrobial Community Analysis
http://compmetagen.github.io/micca
GNU General Public License v3.0
20 stars 9 forks source link

various vsearch errors in mergepairs and otu #15

Closed Phylloxera closed 4 years ago

Phylloxera commented 5 years ago

micca v1.7.2 run via docker on a mac. We had been working through the default paired-end workflow but realized or 'R2' files had forward primer and vice-versa so we changed the i,p&e flags in micca mergepairs as follows: -i <path>/*_R2*.fastq -p _R2 -e _R1 and got the following error: VSEARCH error: vsearch v 2.7.1 Fatal error: Invalid line 172 in FastQ Illegal character ' ' However, the analysis continued to run and produced a merged fastq around the expected size and this was used in successful trim, filterstats, and filter commands subsequently. However, at the micca otu stage, we again encountered a VSEARCH error:

root@7edcd6962efa:/micca# micca otu -i filtered_R2ASforward.fasta -o Masstreat_otus -d 0.97 -c -t 4 Error: vsearch v2.7.1_linux_x86_64, 2.0GB RAM, 4 cores https://github.com/torognes/vsearch Reading file filtered_R2ASforward.fasta

This command produced a folder, but no files. I note that the vsearch errors started appearing when we changed up the forward pattern and my colleague wonders if she might be running into storage/processing limitations of her machine. What do you think? Have you heard of this before? The input filtered.fasta file looks similar to others we have used successfully. Thanks