RNA-seqデータのTrimmomaticによるトリミングがうまくいかずに
困っています。お知恵をお貸しください。
使用したコマンドは下記になります。
java -jar /path/to/trimmomatic-0.39.jar SE -threads 12 \
-phred33 -trimlog log.txt /path/to/input.fastq \
/path/to/output.fastq MINLEN:30 TRAILING:20 \
SLIDINGWINDOW:4:15
トリミング後のFasrQC結果を見ると、TRAILINGオプションはworkしているようなのですが、
リード長30塩基以下のものが除去されていませんでした。何故でしょうか?
原因や問題解決の方法について何かありましたご教示願います。
よろしくお願いいたします。
(ちなみに同じデータをprinseqを用いて同様のトリミングを行なった場合は、特に問題ありませんでした。) |
|