Question: Bowtie1 Options For Mirna
1
gravatar for Varun Gupta
5.0 years ago by
Varun Gupta1.1k
United States
Varun Gupta1.1k wrote:

Hi I am using a data set with 36 bp read sequence. I am using bowtie1 to align it to human genome. Since I am interested in miRNA and those are generally 21-23 bp long, is there a way in bowtie1 that if I specify seed length -l 25 with n option as 2 then if the first 25 bp matches the genomic seq with 2 mis matches and remaining bases does not match, just drop the remaining bases and report the matched 25bp.

If not is there a way I can perform such alignment.

Hope to hear from you soon. Regards Varun

mirna • 1.1k views
ADD COMMENTlink modified 24 months ago by Biostar ♦♦ 20 • written 5.0 years ago by Varun Gupta1.1k

my advice would be to use one of the existing tools for miRNA (http://omictools.com/srna-seq/), like miRDeep2 (https://www.mdc-berlin.de/rajewsky/miRDeep)

ADD REPLYlink written 5.0 years ago by Martombo2.4k
0
gravatar for Ido Tamir
5.0 years ago by
Ido Tamir4.9k
Austria
Ido Tamir4.9k wrote:

If you insist on doing your own pipeline you could use bowtie2 with local alignment or remove the adaptor with one of the many tools that do this.

ADD COMMENTlink written 5.0 years ago by Ido Tamir4.9k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2543 users visited in the last hour