I've been testing SeqPrep (http://seqanswers.com/wiki/SeqPrep) for trimming adapters and merging reads. However I'm unsure why some reads get discarded. The command line is as follows:
SeqPrep -f s_5_1_sequence.txt.gz -r s_5_2_sequence.txt.gz -1 s_5_1-trimmed.fastq.gz -2 s_5_2-trimmed.fastq.gz -A AGATCGGAAGAGCACACGTCTGAACTCCAGTCAC -B AGATCGGAAGAGCGTCGTGTAGGGAAAGAGT -s merged.fastq.gz -L 1 -6 -E alignment.txt.gz
I'm assuming -L 1 would prevent removal of any reads (the help on command line differs abit from the wiki page). However reads still get removed. What causes this and what should I do to affect this?
SeqPrep -f s_5_1_sequence.txt.gz -r s_5_2_sequence.txt.gz -1 s_5_1-trimmed.fastq.gz -2 s_5_2-trimmed.fastq.gz -A AGATCGGAAGAGCACACGTCTGAACTCCAGTCAC -B AGATCGGAAGAGCGTCGTGTAGGGAAAGAGT -s merged.fastq.gz -L 1 -6 -E alignment.txt.gz
I'm assuming -L 1 would prevent removal of any reads (the help on command line differs abit from the wiki page). However reads still get removed. What causes this and what should I do to affect this?
Comment