<div dir="ltr">Hi Lior,<div><br></div><div>I don't know if you've resolved this yet - it's been several months! But I am not aware of any way to do this in Biopython.</div><div><br></div><div>I was doing a similar task so I made a gist using your example:</div><div><br></div><div><a href="https://gist.github.com/lennax/9925c08804a57bb5ef72">https://gist.github.com/lennax/9925c08804a57bb5ef72</a><br></div><div><br></div><div>Cheers,</div><div><br></div><div>Lenna</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Jun 16, 2015 at 3:35 AM, Lior Glick <span dir="ltr"><<a href="mailto:liorglic@mail.tau.ac.il" target="_blank">liorglic@mail.tau.ac.il</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="rtl"><div dir="ltr"><font face="arial, helvetica, sans-serif">Hi all,</font></div><div dir="ltr"><font face="arial, helvetica, sans-serif"><br></font></div><div dir="ltr"><font face="arial, helvetica, sans-serif">I have a pairwise alignment of protein sequences, created with <span style="font-size:13px;line-height:10.8333330154419px;white-space:pre-wrap;color:rgb(102,0,102)">Bio</span><span style="font-size:13px;line-height:10.8333330154419px;white-space:pre-wrap;color:rgb(102,102,0)">.</span><span style="font-size:13px;line-height:10.8333330154419px;white-space:pre-wrap;color:rgb(0,0,0)">pairwise2, something like this:</span></font></div><div dir="ltr"><span style="font-family:Monaco,Menlo,Consolas,'Courier New',monospace;font-size:13px;line-height:10.8333330154419px;white-space:pre-wrap;color:rgb(0,0,0)"><br></span></div><div dir="ltr"><span style="font-family:Monaco,Menlo,Consolas,'Courier New',monospace;font-size:13px;line-height:10.8333330154419px;white-space:pre-wrap;color:rgb(0,0,0)"><pre style="padding:10px;font-family:Monaco,Menlo,Consolas,'Courier New',monospace;color:rgb(51,51,51);border-radius:5px;margin:8px;line-height:10.8333330154419px;word-break:break-all;word-wrap:break-word;white-space:pre-wrap;border:1px solid rgb(198,201,204);background-image:-webkit-linear-gradient(top,rgb(248,248,249) 0%,rgb(236,240,243) 100%);background-repeat:initial"><span style="color:rgb(0,0,0)">C</span><span style="color:rgb(102,102,0)">-</span><span style="color:rgb(0,0,0)">GGLK</span><span style="color:rgb(102,102,0)">-</span><span style="color:rgb(0,0,0)">GAPFLAC
</span><span style="color:rgb(102,102,0)">||||||||||||||</span><span style="color:rgb(0,0,0)">
CHG</span><span style="color:rgb(102,102,0)">--</span><span style="color:rgb(0,0,0)">KTGASFLQC</span></pre></span></div><div dir="ltr"><font face="arial, helvetica, sans-serif">I would like to <span style="color:rgb(51,51,51);font-size:13px;line-height:10.8333330154419px">know which position in the subject sequence is aligned to which position in the query sequence. </span><span style="color:rgb(51,51,51);font-size:13px;line-height:10.8333330154419px">In this example if the top sequence is the subject than I would like to know that position 1 in subject is aligned to position 1 in query, 2 is aligned to 3, 5 to 4, 6 to 6 etc...</span></font></div><div dir="ltr"><span style="font-size:13px;line-height:10.8333330154419px;white-space:pre-wrap;color:rgb(0,0,0)"><font face="arial, helvetica, sans-serif"><br></font></span></div><div dir="ltr"><span style="font-size:13px;line-height:10.8333330154419px;white-space:pre-wrap;color:rgb(0,0,0)"><font face="arial, helvetica, sans-serif">Does Biopython have this functionality, or do I have to write my own function? Any suggestions?</font></span></div><div dir="ltr"><span style="font-size:13px;line-height:10.8333330154419px;white-space:pre-wrap;color:rgb(0,0,0)"><font face="arial, helvetica, sans-serif"><br></font></span></div><div dir="ltr"><span style="font-size:13px;line-height:10.8333330154419px;white-space:pre-wrap;color:rgb(0,0,0)"><font face="arial, helvetica, sans-serif">Thanks!</font></span></div><div dir="ltr"><span style="font-family:Monaco,Menlo,Consolas,'Courier New',monospace;font-size:13px;line-height:10.8333330154419px;white-space:pre-wrap;color:rgb(0,0,0)"><br></span></div><div dir="ltr"><span style="font-family:Monaco,Menlo,Consolas,'Courier New',monospace;font-size:13px;line-height:10.8333330154419px;white-space:pre-wrap;color:rgb(0,0,0)"><br></span></div><div dir="ltr"><span style="font-family:Monaco,Menlo,Consolas,'Courier New',monospace;font-size:13px;line-height:10.8333330154419px;white-space:pre-wrap;color:rgb(0,0,0)"><br></span></div><div dir="ltr"><span style="font-family:Monaco,Menlo,Consolas,'Courier New',monospace;font-size:13px;line-height:10.8333330154419px;white-space:pre-wrap;color:rgb(0,0,0)"><br></span></div><div dir="ltr"><span style="font-family:Monaco,Menlo,Consolas,'Courier New',monospace;font-size:13px;line-height:10.8333330154419px;white-space:pre-wrap;color:rgb(0,0,0)"><br></span></div></div>
<br>_______________________________________________<br>
Biopython mailing list  -  <a href="mailto:Biopython@mailman.open-bio.org">Biopython@mailman.open-bio.org</a><br>
<a href="http://mailman.open-bio.org/mailman/listinfo/biopython" rel="noreferrer" target="_blank">http://mailman.open-bio.org/mailman/listinfo/biopython</a><br></blockquote></div><br></div>