8 days ago by
You probably do not want to use HEADCROP:30 - that will remove the first 30 nucleotides (!) whereas LEADNING will only remove low quality bases in the 5' end.
In regards to using TRAILING the other option is to filter on the average quality in a window via the SLIDINGWINDOW command. According to this article is better than TRAILING (due to there sometimes being a single semi-low quality base followed by additional high quality bases - this can easily be mapped (as mappers do not require perfect matches) but will be thrown away if using TRAILING).
Other things to pay attention to are:
- The order of arguments matters as that is the order in which the trimming is performed so you always want to end with MINLEN
- You can also remove (any) adapter contamination via ILLUMINACLIP.
- You can find more info on how to use Trimmomatic here.