SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Finding Perfect matches and mismatches bhargavpbk88 Bioinformatics 2 01-19-2012 07:00 AM
Perfect 454 FLX for sale (with FREE upgrade to long reads FLX+ available) 454FLXforsale Service Providers 0 01-18-2012 04:38 AM
perfect alignment match NicoBxl Bioinformatics 3 10-07-2010 04:28 AM
output all possible match using bwa mingkunli Bioinformatics 2 09-05-2009 08:06 AM
Perfect match disagreement between bowtie and BLAT on human genome danielsbrewer Bioinformatics 3 03-10-2009 09:10 PM

Reply
 
Thread Tools
Old 04-01-2011, 12:13 PM   #1
guil
Junior Member
 
Location: brazil

Join Date: Jul 2009
Posts: 4
Default bwa perfect match only

Hi everyone,

I'm trying to do an aling of short reads (putative micrornas) to putative precursors using bwa...

but I wanna only perfect matchs between mature and precursor....

in bwa manual I haven't found any parameter to map a read with no mismach.. only gap opening an extesion penalties...


Are there another way to do this???


tnks
guil is offline   Reply With Quote
Old 04-01-2011, 12:24 PM   #2
ttnguyen
Member
 
Location: Ireland

Join Date: Mar 2010
Posts: 41
Default

Have you tried this?
-k 0
ttnguyen is offline   Reply With Quote
Old 04-01-2011, 05:31 PM   #3
guil
Junior Member
 
Location: brazil

Join Date: Jul 2009
Posts: 4
Default

thanks ttnguyen,

but i tried the

-k 0 and it doesn't works....

I also tried to extend the mismatch penalty --> 10

and still results in mapping with mismatches...

any ohter idea???

best
guil is offline   Reply With Quote
Old 04-06-2011, 09:06 AM   #4
jiangyuehit
Junior Member
 
Location: toronto

Join Date: Apr 2011
Posts: 2
Default

what about filter out alignment that NM value is not 0?
jiangyuehit is offline   Reply With Quote
Old 08-30-2013, 09:18 AM   #5
yifangt
Member
 
Location: Canada

Join Date: Feb 2011
Posts: 61
Default

I am following the same question. There may be two options:
1) Redo the mapping with bowtie or bwa, what is the option to accomplish this? -k 0 does not work!

2) Filter out the SAM/BAM file. What is the option to filter no-mismatch!
Seems not enough information about the option "H:i:0"
I search over this forum and biostar, many discussion with the unique mapped reads, ot unmapped reads. No discussion about the perfectly mapped reads, i.e. mapped reads without any mismatch.
For column 6 of CIGAR string, I used
Code:
 grep "\s\d+M\s"
to filter out any matched with deletion and insertion, but the problem after filtering is M can stands for mismatch too, so that the reads still have some that contains mismatch!
Can anybody give some suggestion.
Thanks a lot!

Last edited by yifangt; 08-30-2013 at 10:36 AM.
yifangt is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 11:28 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO