Latest
Open
RNA-Seq
ChIP-Seq
SNP
Assembly
Forum
Planet
All »
Home
Log In
Question: Best mapping algorithm
0
gravatar for heather.watkins
6 months ago by
heather.watkins • 0
heather.watkins • 0 wrote:

What is the best mapping algorithm to use? My reference genome is about 70,000 sequences, each 20 bp long. My reads are mostly longer, > 100 bp. Thank you.
assembly mapping • 256 views
ADD COMMENT • link •
modified 6 months ago by Jennifer Hillman Jackson ♦ 25k • written 6 months ago by heather.watkins • 0
1
gravatar for Jennifer Hillman Jackson
6 months ago by
Jennifer Hillman Jackson ♦ 25k
United States
Jennifer Hillman Jackson ♦ 25k wrote:

Hello,

Mapping reads against genome/transcriptome/exon fragments, especially when the reads are longer than the target, is not going to work with NGS mapping tools like HISAT2, BWA, etc.

You could try Blast+ blastn instead. I am not sure if makeblastdb will work with this many target sequences as implemented at Galaxy Main https://usegalaxy.org, and using a fasta from the history will almost certainly be problematic.

You'll need to test it out, tune parameters, retest. Tools can run into memory problems with this type of target database with the resources available at Galaxy Main. If you do run into memory errors, those do not need to be reported as a bug - it would be an expected failure.

Galaxy tutorials: https://galaxyproject.org/learn/

Thanks! Jen, Galaxy team
ADD COMMENT • link written 6 months ago by Jennifer Hillman Jackson ♦ 25k

Thank you Jen. This is just what I needed.

Heather

ADD REPLY • link written 6 months ago by heather.watkins • 0
Please log in to add an answer.

Similar posts • Search »

Content
Help
Access
Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 183 users visited in the last hour