smallRNA-seq alignment with bowtie1:what is the argument for handling multi-mapped reads?
0
0
Entering edit mode
3.2 years ago

Hey. I am working on smallRNA-seq data with the final goal of investigating differential expression analysis.

I used bowtie1 for mapping with the following specification: bowtie -p 10 -n 2 -m 1 -k 1 --strata --best -x index -q *fq

The summary statistics I got are as follows:

#reads processed: 2,00,1187
# reads with at least one alignment: 1,559,437 (77.93%)
# reads that failed to align: 441,750 (22.07%)
# reads with alignments suppressed due to -m: 1,368,859 (68.40%)

Given 68% of the reads were suppressed, I feel a huge chunk of the reads were non-unique and I'm not sure if retaining only the best 1 (-m 1 -k 1 --strata --best) is most appropriate. Please any suggestion from experienced smallRNA-seq persons will be appreciated.

Thanks

alignment smallRNA-seq bowtie1 multi-mapped reads • 485 views
ADD COMMENT

Login before adding your answer.

Traffic: 2006 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6