Hiya, I'm trying to make contigs in mothur (v.1.41.1) from my mrDNA-sequenced samples. I've used MrDNA:s Fastq Processor program to get the reads indexed, and made file with mothur make.file command. Everything looks ok at this point, however, I don't have the same number of reads in forward and reverse, as FastQProcessor result.file tells me. For example, Barcodes R1 Seqs R2 Seqs ATTCCGAC 37762 34504
Therefore, when I try make.contigs using the file file, I get an error message "xxx is in your forward fastq file and not in your reverse file, please remove it using the remove.seqs command before proceeding." and mothur crashes because of a segmentation fault.
I've been trying to find a solution by reading help forums, and remove.seqs could be the answer, but I don't understand how to use the command to remove unpaired reads from both R1 and R2. The reads are at this point in fastq-format. Cheers, Lotta
I think you dont need to look for a solution to remove reads to make R1 and R2 equal again. You should look at why they dont have the same amount of reads anyways. Maybe something went rong with that Fastq Processor program. Were they equal before you processed them with Fastq Processor?
At least the file size is not the same. Any suggestions on how to check whether the raw files R1 and R2 have equal amount of reads?
How to count fastq reads
Thanks! So in the raw files (where all samples are in one file, before using fastqProcessor by MrDNA), same amount of reads, but after separating the reads with the program, different number. So likely there is something fishy going on with fastqProcessor...
Exactly, or you need to change settings or something. Just as far as I know you should not remove reads afterwards to make them the same. That is asking for problems. Good luck