This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Aligner For Allowing Multiple Matches (#N) For A Read

Hi All

Based on your experience could you please recommend an aligner which can accommodate multiple read mapping positions. I have already tried BWA and Mosaik with the option to allow for multiple mappings.

With BWA I dont get the Mapping quality for all mappings and Mosaik gives me low % of reads aligned. I was wondering may be I could try something else.

Thanks! -Abhi

short aligner
• 3,710 views • link
updated by Ram Ram • written by Abhi Abhi

Mapping quality can only be calculated for the top hit.

• 0 views • link
written by lh3 lh3

3 answers

try gem-mapper

• 0 views • link
written by Leszek Leszek

Hi, recently, I have a problem about GEM. When I want build index for human using GEM, it reported

Creating sequence and location files...
Fatal error: exception
Failure("Command 'echo "-------gem-indexer_fasta2meta+cont" > human_g1k_v37.fasta.log && gem-indexer_fasta2meta+cont -i human_g1k_v37.fasta -c dna --filter-function iupac-dna --strip-unknown-bases-threshold 50 --complement-size-threshold 2000000000 --comp.

The command I used is

./gem-indexer -i /devdata/chhy/masai/bin/human_g1k_v37.fasta -o hum.index -T 4

Could you help me to solve this problem?

Thanks!

• 0 views • link
updated by Ram Ram • written by chhylp123 • chhylp123

You can try the MUMmer suite with the maxmatch option. So allowing multiple matches. Another alternative is to simply use the fast BLAT (not BLAST) aligner.

• 0 views • link
written by ALchEmiXt ALchEmiXt

bowtie does this beautifully.

"With the -k, -a, -m, -M, --best and --strata options, the user can flexibily select which alignments are reported"

http://bowtie-bio.sourceforge.net/manual.shtml#reporting-modes

• 0 views • link
updated by Ram Ram • written by Lee Brooks • Lee Brooks

Log in to answer this question.

More posts like this