Mapping Very Short (<20 Bp) Reads
2
2
Entering edit mode
13.3 years ago
Hcl ▴ 30

Hi,

Can anyone suggest an aligner which can handle mapping very short reads (<20 bp) to a reference genome?

Thank you.

short aligner • 5.5k views
ADD COMMENT
2
Entering edit mode

What kind of experiment? RNA-seq? What genome? bacteria ?

ADD REPLY
0
Entering edit mode

What kind of experiment? RNA-seq? Genome ?

ADD REPLY
0
Entering edit mode

And how many reads?

ADD REPLY
5
Entering edit mode
13.3 years ago
Leszek 4.2k

I think you can use any of the short reads aligners. Just be careful to set number of mismatches to 0 or 1 and maybe disable indels, as @ 20bp sequence length the probability of random match is quite high. That depends of course on the size of your genome.
What is your reference genome size?

ADD COMMENT
0
Entering edit mode

yes try bowtie with -v 0 parameter to set the number of mismatches to 0

ADD REPLY
1
Entering edit mode
13.3 years ago
Hcl ▴ 30

I tried SOAP2 and bowtie. SOAP2 could only align reads >30 bp but bowtie did the job beautifully.

ADD COMMENT
0
Entering edit mode

I'm using GEM-mapper for most aplications - give it a try:)

ADD REPLY

Login before adding your answer.

Traffic: 986 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