aligning reads against shorter region
1
0
Entering edit mode
3.5 years ago
Juke34 8.9k

I would like to align 150bp reads against a reference region slightly shorter e.g 120bp from amplicon data. Would the alignment perform well or should we necessarily use a longer reference?

alingment Bowtie2 • 778 views
ADD COMMENT
1
Entering edit mode
3.5 years ago
pbpanigrahi ▴ 430

Manual says

"Bowtie 2 is geared toward aligning relatively short sequencing reads to long genomes. That said, it handles arbitrarily small reference sequences (e.g. amplicons) and very long reads (i.e. upwards of 10s or 100s of kilobases), though it is slower in those settings. It is optimized for the read lengths and error modes yielded by typical Illumina sequencers."

You can explore --local option

ADD COMMENT
0
Entering edit mode

Thanks sounds the way to go with the --local.

ADD REPLY

Login before adding your answer.

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