Re: Adding overall coverage to blast-xml
John T <[email protected]> Mon, 11 Sep 2017 22:02:47 +0100
| Newsgroups | gmane.comp.python.bio.devel |
|---|---|
| Message-ID | <CAL0Rc+zgimuf4niGymbACBLOest25c3sSaW6ttGk=gHFnFdXjA@mail.gmail.com> |
--===============2704506561088095986== Content-Type: multipart/alternative; boundary="f403045c8feab7446f0558f03fec" --f403045c8feab7446f0558f03fec Content-Type: text/plain; charset="UTF-8" It is indeed qcovs, I was stuck using XML, didn't investigate that before. I haven't looked into the SearchIO version yet, but I assume it won't be more difficult to implement there so I'm happy to commit to that. John On 11 September 2017 at 21:14, Peter Cock <[email protected]> wrote: > Hi John, CC Bow, > > Is this information in any of the BLAST output formats? > > I suspect you are talking about recreating one or more > of the coverage fields which can be requested in the > BLAST csv or tabular output files formats? i.e. qcovs > or qcovhsp perhaps? > > If so, I can see a rational for computing this for the > BLAST XML parsing... although given we'd like to > push people towards SearchIO (which supports the > BLAST tabular format as well), doing it there might > make more sense than in NCBIXML? > > Peter > > On Mon, Sep 11, 2017 at 5:40 PM, John T <[email protected]> wrote: > > Hi all, > > > > I've written a script that calculates the overall coverage for all HSP > in an > > alignment, like that given on the BLAST homepage. Would it make sense to > add > > this to Bio.Blast.Record.Alignment (and/or the BlastIO equivalent)? It's > not > > something that's explicitly defined in the XML and so may be out of > scope, > > but it is something that's given in the website version of BLAST results > so > > people might reasonably expect it to be there. I'd add it as a @property > so > > it as only calculated when required. > > > > I've never tried contributing to an open source project before so I > didn't > > want to dive in and do something stupid. > > > > Cheers, > > John > > > > > > _______________________________________________ > > Biopython-dev mailing list > > [email protected] > > http://mailman.open-bio.org/mailman/listinfo/biopython-dev > --f403045c8feab7446f0558f03fec Content-Type: text/html; charset="UTF-8" Content-Transfer-Encoding: quoted-printable <div dir=3D"ltr"><div>It is indeed qcovs, I was stuck using XML, didn't= investigate that before. I haven't looked into the SearchIO version ye= t, but I assume it won't be more difficult to implement there so I'= m happy to commit to that.<br></div><div><br></div><div>John</div></div><di= v class=3D"gmail_extra"><br><div class=3D"gmail_quote">On 11 September 2017= at 21:14, Peter Cock <span dir=3D"ltr"><<a href=3D"mailto:p.j.a.cock@go= oglemail.com" target=3D"_blank">[email protected]</a>></span> wr= ote:<br><blockquote class=3D"gmail_quote" style=3D"margin:0 0 0 .8ex;border= -left:1px #ccc solid;padding-left:1ex">Hi John, CC Bow,<br> <br> Is this information in any of the BLAST output formats?<br> <br> I suspect you are talking about recreating one or more<br> of the coverage fields which can be requested in the<br> BLAST csv or tabular output files formats? i.e. qcovs<br> or qcovhsp perhaps?<br> <br> If so, I can see a rational for computing this for the<br> BLAST XML parsing... although given we'd like to<br> push people towards SearchIO (which supports the<br> BLAST tabular format as well), doing it there might<br> make more sense than in NCBIXML?<br> <br> Peter<br> <div><div class=3D"h5"><br> On Mon, Sep 11, 2017 at 5:40 PM, John T <<a href=3D"mailto:jaytee00@gmai= l.com">[email protected]</a>> wrote:<br> > Hi all,<br> ><br> > I've written a script that calculates the overall coverage for all= HSP in an<br> > alignment, like that given on the BLAST homepage. Would it make sense = to add<br> > this to Bio.Blast.Record.Alignment (and/or the BlastIO equivalent)? It= 's not<br> > something that's explicitly defined in the XML and so may be out o= f scope,<br> > but it is something that's given in the website version of BLAST r= esults so<br> > people might reasonably expect it to be there. I'd add it as a @pr= operty so<br> > it as only calculated when required.<br> ><br> > I've never tried contributing to an open source project before so = I didn't<br> > want to dive in and do something stupid.<br> ><br> > Cheers,<br> > John<br> ><br> ><br> </div></div>> ______________________________<wbr>_________________<br> > Biopython-dev mailing list<br> > <a href=3D"mailto:[email protected]">Biopython-dev@ma= ilman.open-<wbr>bio.org</a><br> > <a href=3D"http://mailman.open-bio.org/mailman/listinfo/biopython-dev"= rel=3D"noreferrer" target=3D"_blank">http://mailman.open-bio.org/<wbr>mail= man/listinfo/biopython-dev</a><br> </blockquote></div><br></div> --f403045c8feab7446f0558f03fec-- --===============2704506561088095986== Content-Type: text/plain; charset="us-ascii" MIME-Version: 1.0 Content-Transfer-Encoding: 7bit Content-Disposition: inline _______________________________________________ Biopython-dev mailing list [email protected] http://mailman.open-bio.org/mailman/listinfo/biopython-dev --===============2704506561088095986==--