I have Illumina smallRNA-RNAseq reads that have already had their adapters removed and have been filtered for tRNA/rRNA.
I'm trying to use SeqTrimMap to align these reads to my reference genome with the following command:
The program runs well and gives me output in SAM format as requested, but the program will only map reads at lengths 20-24 nt. Does anyone have a suggestion that might allow me to map reads at additional lengths?
I'm trying to use SeqTrimMap to align these reads to my reference genome with the following command:
Code:
$SeqTrimMap -b -S -i -v 1 -o ~/seqtrimmap_B_output_496_1 -s -u 18 -t 28 ~/496_1_filtered_18_28nt.fa ~/zma_b.fa
Comment