[Biopython] Help modify this code so it can do what I want it to do

Edson Ishengoma ishengomae at nm-aist.ac.tz
Wed Feb 5 17:52:17 UTC 2014


Hi Peter,

Woow, that made my day. Thank you very much and keep up the good work.

Regards,

Edson

On Wed, Feb 5, 2014 at 7:07 PM, Peter Cock <p.j.a.cock at googlemail.com>wrote:

> Hi Edson,
>
> I can see where the problem stems from now - it did puzzle me for a while.
> For this part to make sense, query and sbjct need to be the FULL sequence
> of the query and the subject (as given to BLAST as input):
>
>     complete_query_seq += str(query[q_start-1:q_end])
>     complete_sbjct_seq += str(sbjct[sb_start-1:sb_end])
>
> (I had assumed these variables were setup at the beginning of the file,
> which I partly why I asked for the full script.)
>
> However, via the for loop, you are using hsp.query, hsp.sbjct as query
> and sbjct, This are the PARTIAL sequences aligned with gap characters.
> This might do what you seemed to want:
>
>     complete_query_seq += query.replace("-", "")
>     complete_sbjct_seq += sbjct.replace("-", "")
>
> However, this will concatenate the fragments with an HSP - any bit of
> the query or subject which did not align will not be included. Any bit
> which appears in more than one HSP will be there twice. And also
> if you're using masking you'll have XXXXX X regions in the sequence
> where the filter said it was low complexity.
>
> I would instead get the original unmodified query/subject sequences
> from the original FASTA files given to BLAST.
>
> Peter
>
>
> On Tue, Feb 4, 2014 at 9:12 AM, Edson Ishengoma
> <ishengomae at nm-aist.ac.tz> wrote:
> > Hi Peter,
> >
> > My apology, I have updated the code at
> > https://gist.github.com/EBIshengoma/efc4ad3e32427891931d to appear
> exactly
> > how I run it from my computer.
> >
> > Thanks.
> >
>



More information about the Biopython mailing list