This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Bwa : Allow Low Stringency In Alignments

Hi,

I wonder if it's possible to allow small local alignment using bwa. Per example: with 2x100bp reads, to output all alignment with minimum 10 bases mapped ? So a very loose alignment.

Thanks

bwa
• 3,850 views • link
updated by Chris Fields Chris Fields • written by Nicolas Rosewick Nicolas Rosewick

I think you might be better of with a Smith-Waterman implementation like 'bwa bwasw', but am not sure what the smallest possible match is.

Andreas

• 0 views • link
written by Andreas Andreas

1 answer

Maybe have a look at bwa mem (in newer releases of bwa, post 0.7 I believe). It has a mismatch penalty, so you can try playing with that as well as seed length.

lh3's paper on bwa-mem is in arxiv and is also described elsewhere in biostar.

• 0 views • link
updated by Ram Ram • written by Chris Fields Chris Fields

Log in to answer this question.

More posts like this
  • How Does Bwa Report Multi-Mapped Reads
    written by bbio bbio
    <p>I am mapping some RNA-seq data with <a href='http://bio-bwa.sourceforge.net/'>bwa</a> and would like to do some analysis on where multi-mapped reads fall.</p> <p>I know that I …
  • Bwa On Mac : How To Run It On Terminal
    written by freddy • freddy
    <p>Can any one help me please I want use <a href='http://bio-bwa.sourceforge.net/'>BWA</a> aligner, when I downloaded I cd into the file and made it and it …
  • Remove Soft Clipped Bases
    written by Clare Clare
    <p>I want to conduct some computations using a python script directly on some <a href='http://bio-bwa.sourceforge.net/'>BWA</a> aligned <a href='http://samtools.sourceforge.net/SAM1.pdf'>bam</a> files, and to do this I need …
  • Bwa Index Output And Bwa Aln Problem
    written by Tonyzeng • Tonyzeng
    <p>My output of <a href='http://bio-bwa.sourceforge.net/'>BWA</a> index has only 4 files as </p> <p>mm10.fa.amb mm10.fa.ann mm10.fa.bwt mm10.fq.pac</p> <p>I read a thread that a guy's <a href='http://bio-bwa.sourceforge.net/'>BWA</a> …
  • What Causes The Bwa Error: "Assertion `Re - Rb == Rlen' Failed."
    written by Cacau Cacau
    <p>Neither <a href='http://bio-bwa.sourceforge.net/'>bwa</a> 6.0.2 nor <a href='http://bio-bwa.sourceforge.net/'>bwa</a> 7.0.5 worked for my alignment. The size of my genome is about 1.5M and the number reads are …
  • Bwa Alignments With Gaps
    written by lob • lob
    <p>I am looking for a description of <a href='http://bio-bwa.sourceforge.net/'>bwa</a> algorithm for mapping with gaps. More precisely I wonder how <a href='http://bio-bwa.sourceforge.net/'>bwa</a> handles mismatches in reads …
  • Whole Genome Sequencing: Bwa
    written by mathiasf52 • mathiasf52
    <p>Hi all, </p> <p>I want to use <a href='http://bio-bwa.sourceforge.net/'>BWA</a> for aligning human WGS paired-end data generated by Illumina. I read the manual of <a href='http://bio-bwa.sourceforge.net/'>BWA</a> …
  • Mapping Quality Issues With Bwa/Stampy Hybrid Alignment
    written by oadaily • oadaily
    <p>Hi all,</p> <p>I am new to biostar so I want start by saying thank you all for sharing your expertise. </p> <p>I recently started testing …
  • Bwa Alignment: Gap Vs Soft-Clipped Sequences
    written by michealsmith michealsmith
    <p>Suppose we have a 30bp-50bp deletion, and read length is 100bp. If the read happens to cover the deletion, I'm just wondering, will <a href='http://bio-bwa.sourceforge.net/'>BWA</a> …
  • Clarification For Running Bwa In Parallel‏
    written by Mohamed • Mohamed
    <p>Hi,</p> <p>I tried running <a href='http://bio-bwa.sourceforge.net/'>BWA</a> in parallel by splitting the input FASTQ file into 'n' splits and ran single ended alignment (aln and samse …