Go Back   SEQanswers > General

Similar Threads
Thread Thread Starter Forum Replies Last Post
Reformatting BLAST+ pairwise alignment? Genomics101 Bioinformatics 2 12-05-2013 02:53 AM
Merging Pairwise Sequence Alignment loba17 Bioinformatics 10 12-17-2012 04:59 AM
fastest pairwise alignment method Trilliput Bioinformatics 0 10-19-2012 11:59 AM
R mismatch position in pairwise alignment NicoBxl Bioinformatics 4 10-27-2010 08:36 AM

Thread Tools
Old 03-13-2014, 08:59 AM   #1
Junior Member
Location: newyork

Join Date: Mar 2014
Posts: 1
Default pairwise sequence alignment

This is a primitive question but somehow I could not find a solution to it. I have two sequences 'large' and 'small' as given below.



I need to align the 'small' sequence to the 'large' sequence. Clearly there are two places where it can be aligned. I need to get indices of both the locations. I was trying BioPython's "pairwise2.align.globalms" function but it is only able to align to the second position.

pairwise2.align.globalms(largeStr, smallStr, 2, -1, -1, 0, penalize_end_gaps=False)

Which parameter can I change here or which other pachage/lightweight free software can compute this?

marykindall is offline   Reply With Quote
Old 03-13-2014, 12:32 PM   #2
Rick Westerman
Location: Purdue University, Indiana, USA

Join Date: Jun 2008
Posts: 1,104

bowtie or bowtie2, although overkill for a single reference, will report all alignments. For a short query sequence bowtie is probably better.
westerman is offline   Reply With Quote

pairwise, python, sequence alignments

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 03:01 AM.

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