<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
</head>
<body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
<div class="">Is there any particular standard for these comment types? I see these within WGS master records all the time (denotes basic metadata on assembly).</div>
<div class=""><br class="">
</div>
<div class=""><a href="http://www.ncbi.nlm.nih.gov/nuccore/635626163" class="">http://www.ncbi.nlm.nih.gov/nuccore/635626163</a></div>
<div class=""><a href="http://www.ncbi.nlm.nih.gov/nuccore/AFTI00000000.1" class="">http://www.ncbi.nlm.nih.gov/nuccore/AFTI00000000.1</a></div>
<div class=""><br class="">
</div>
<div class="">(you can’t tell I’m working on a mussel assembly, right?)</div>
<div class=""><br class="">
</div>
<div class="">chris</div>
<div class=""><br class="">
<div>
<blockquote type="cite" class="">
<div class="">On Sep 9, 2015, at 12:08 PM, Peter Cock <<a href="mailto:p.j.a.cock@googlemail.com" class="">p.j.a.cock@googlemail.com</a>> wrote:</div>
<br class="Apple-interchange-newline">
<div class="">That's ... nasty. How about record.features["comment"] for the old style<br class="">
plain text (as a Python string) and record.features["structured_comment"]<br class="">
as a (sorted) Python dict?<br class="">
<br class="">
This might make GenBank output easier too...<br class="">
<br class="">
Peter<br class="">
<br class="">
On Wed, Sep 9, 2015 at 6:01 PM, Brian Osborne <<a href="mailto:bosborne11@verizon.net" class="">bosborne11@verizon.net</a>> wrote:<br class="">
<blockquote type="cite" class="">Peter,<br class="">
<br class="">
That is an interesting idea. What would be returned if the COMMENT has both plain and “structured comments” in it? Here’s one:<br class="">
<br class="">
<a href="http://www.ncbi.nlm.nih.gov/nuccore/FJ966082" class="">http://www.ncbi.nlm.nih.gov/nuccore/FJ966082</a><br class="">
<br class="">
Thanks again,<br class="">
<br class="">
Brian O.<br class="">
<br class="">
<br class="">
<br class="">
<br class="">
<blockquote type="cite" class="">On Sep 9, 2015, at 7:27 AM, Peter Cock <p.j.a.cock@googlemail.com> wrote:<br class="">
<br class="">
This sounds good - would you turn these into a Python dict?<br class="">
<br class="">
Peter<br class="">
<br class="">
On Wed, Sep 9, 2015 at 2:56 PM, Brian Osborne <bosborne11@verizon.net> wrote:<br class="">
<blockquote type="cite" class="">All,<br class="">
<br class="">
I noticed that BioPython, like the versions of BioPerl in CPAN, does not<br class="">
handle GenBank structured comments<br class="">
(http://www.ncbi.nlm.nih.gov/genbank/structuredcomment) in the ideal way.<br class="">
Here’s an example structured comment:<br class="">
<br class="">
COMMENT ##FluData-START##<br class="">
EPI_ISOLATE_ID :: EPI_ISL_77637<br class="">
NAME :: A/California/07/2009<br class="">
TYPE :: H1N1<br class="">
Segment_name :: M'<br class="">
HOST_AGE :: 54<br class="">
HOST_GENDER :: F'<br class="">
PASSAGE :: M1/C1 (2009-04-24)<br class="">
LOCATION :: United States / California'<br class="">
COLLECT_DATE :: 09-Apr-2009<br class="">
Lineage :: A(H1N1)pdm09<br class="">
RESIST_TO_ADAMANTANES :: Resistant'<br class="">
RESIST_TO_OSELTAMIVIR :: Sensitive'<br class="">
RESIST_TO_ZANAMVIR :: Sensitive'<br class="">
SPECIMEN_ID :: H13596<br class="">
SENDER_LAB :: Naval Health Research Center'<br class="">
SEQLAB_SAMPLE_ID :: 2009712111<br class="">
EPI_SEQUENCE_ID :: EPI273604<br class="">
##FluData-END##<br class="">
<br class="">
Or here: http://www.ncbi.nlm.nih.gov/nuccore/291609868<br class="">
<br class="">
A table, with tag/value pairs. A fair number of bacterial genomes in GenBank<br class="">
use the structured comment to hold MIGS/MIMS data. The comment() method<br class="">
should return something like this, which is easily parsed:<br class="">
<br class="">
##FluData-START##<br class="">
EPI_ISOLATE_ID :: EPI_ISL_77637<br class="">
NAME :: A/California/07/2009<br class="">
TYPE :: H1N1<br class="">
Segment_name :: M'<br class="">
HOST_AGE :: 54<br class="">
HOST_GENDER :: F'<br class="">
PASSAGE :: M1/C1 (2009-04-24)<br class="">
LOCATION :: United States / California'<br class="">
COLLECT_DATE :: 09-Apr-2009<br class="">
Lineage :: A(H1N1)pdm09<br class="">
RESIST_TO_ADAMANTANES :: Resistant'<br class="">
RESIST_TO_OSELTAMIVIR :: Sensitive'<br class="">
RESIST_TO_ZANAMVIR :: Sensitive'<br class="">
SPECIMEN_ID :: H13596<br class="">
SENDER_LAB :: Naval Health Research Center'<br class="">
SEQLAB_SAMPLE_ID :: 2009712111<br class="">
EPI_SEQUENCE_ID :: EPI273604<br class="">
##FluData-END##<br class="">
<br class="">
Rather than this, which is what it currently returns:<br class="">
<br class="">
##FluData-START## EPI_ISOLATE_ID :: EPI_ISL_77637 NAME<br class="">
:: A/California/07/2009 TYPE :: H1N1 Segment_name<br class="">
:: M' HOST_AGE :: 54 HOST_GENDER :: F' PASSAGE<br class="">
:: M1/C1 (2009-04-24) LOCATION :: United States / California'<br class="">
COLLECT_DATE :: 09-Apr-2009 Lineage :: A(H1N1)pdm09<br class="">
RESIST_TO_ADAMANTANES :: Resistant' RESIST_TO_OSELTAMIVIR :: Sensitive'<br class="">
RESIST_TO_ZANAMVIR :: Sensitive' SPECIMEN_ID :: H13596<br class="">
SENDER_LAB :: Naval Health Research Center' SEQLAB_SAMPLE_ID<br class="">
:: 2009712111 EPI_SEQUENCE_ID :: EPI273604 ##FluData-END##<br class="">
<br class="">
Are there any objections to me putting in a pull request with this change? I<br class="">
made this same fix in BioPerl. Of course, if the comment is a “normal” one,<br class="">
it will be treated the same as it is treated now. Another words, the vast<br class="">
majority of comments stay the same.<br class="">
<br class="">
I’ll also add tests.<br class="">
<br class="">
Thanks again,<br class="">
<br class="">
Brian O.<br class="">
<br class="">
_______________________________________________<br class="">
Biopython-dev mailing list<br class="">
Biopython-dev@mailman.open-bio.org<br class="">
http://mailman.open-bio.org/mailman/listinfo/biopython-dev<br class="">
</blockquote>
</blockquote>
<br class="">
</blockquote>
<br class="">
_______________________________________________<br class="">
Biopython-dev mailing list<br class="">
<a href="mailto:Biopython-dev@mailman.open-bio.org" class="">Biopython-dev@mailman.open-bio.org</a><br class="">
http://mailman.open-bio.org/mailman/listinfo/biopython-dev</div>
</blockquote>
</div>
<br class="">
</div>
</body>
</html>