Mapping variable read length

Hi There,

Does anyone have a recommendation on tools that allow to map reads of variable length (between 18-300 nt) as can be seen here

https://www.dropbox.com/s/10dn098492zst9i/Length.png?dl=0

I had previously (6 years back) used blat to map such cases, and use percent identity to get mapped reads. Are there some new tools that can efficiently map sequence reads of variable length. thanks !!

alignment

• 1,774 views • link

updated 5.9 years ago by Biostar

• written 8.5 years ago by Chirag Nepal

2 answers

Any of the popular mapping tools, e.g. bowtie2 or bwa should be able to handle this. What kind of data is this and how certain are you that the small fragments are indeed meaningful and not any noisy byproducts?

• 0 views • link

written 8.5 years ago by ATpoint

bbmap.sh from BBMap suite should be able to do this.

• 0 views • link

written 8.5 years ago by GenoMax

More posts like this

Very low average coverage obtained
written 3.1 years ago by mail2steff

Dear All, I got the viral genome assembly using metaviral spade. Since this is a relatively new genome, I don't have any closest proper reference. …
What could be the rage for Average Coverage after Assembly
written 3.1 years ago by mail2steff

Dear all, I performed viral genome assembly using metaviral spade and I wanted to map my reads back to contigs to know the coverage. I …
Tutorial: Yes .. BBMap can do that! - Part II randomreads (read simulation), demuxbyname/filterbyna…
written 5.6 years ago by GenoMax

NOTE: This collection was originally posted at SeqAnswers.com. Creating a copy here to preserve the information. Part I is available here: https://www.biostars.org/p/483592/ Part III is …
Tutorial: Yes .. BBMap can do that! - Part I : bbmap (aligner), bbduk (scan/trim), repair (fix PE r…
written 5.6 years ago by GenoMax

NOTE: This collection was originally posted at SeqAnswers.com. Creating a copy here to preserve the information. Part II is available here: https://www.biostars.org/p/483593/ Part III is …
Best/right way to quantify small RNA transcripts
written 8.5 years ago by Chirag Nepal

Hi there, I want to quantify read counts of small RNAs into reads per million (RPM). Which is the correct/best way to do ? small …
Pair-end reads are merged into single file: need to separate
written 11.9 years ago by Chirag Nepal

Hi there, I downloaded publicly available datasets. which is 100bp long pair-end reads (with 200 nt of insert size). I downloaded fastq to map to …
Tool: Bedops V2.2 Released
written 13.3 years ago by Alex Reynolds

<p><strong><a href='http://code.google.com/p/bedops/'>BEDOPS</a></strong> is a suite of tools to address common questions raised in genomic studies — mostly with regard to overlap and proximity relationships between …
Tutorial: Mapping Reads With Bwa And Bowtie
written 14.3 years ago by MSU course 2011 •

**Edit in 2019** by https://www.biostars.org/u/25721/ : For users who stumble over this in 2019 (or later), please note that many of the commands and options …