Re: Writing protein records in Swiss-Prot format
Peter Cock <[email protected]> Mon, 25 Dec 2017 22:23:18 +0000
| Newsgroups | gmane.comp.python.bio.devel |
|---|---|
| Message-ID | <CAKVJ-_7Fjkzm9naUv+w-5iUpcdaob3U8+CdXdv5XuASLux6egg@mail.gmail.com> |
--===============3261069272977607231== Content-Type: multipart/alternative; boundary="94eb2c0599400328730561319d61" --94eb2c0599400328730561319d61 Content-Type: text/plain; charset="UTF-8" Content-Transfer-Encoding: quoted-printable Thanks for the explanation. Note I have sometimes needed to *filter* a SwissProt format file, and you can do this via SeqIO's index get_raw functionality (grab the original raw record as a string, and save it to a file as is). There's an example of this in our Tutorial, see the section "Getting the raw data for a record": http://biopython.org/DIST/docs/tutorial/Tutorial.html Peter On Mon, Dec 25, 2017 at 4:05 PM, Adam Sj=C3=B8gren <[email protected]> wrot= e: > Peter writes: > > > I don't recall anyone asking about this, so there doesn't seem to be > > much need for it. And perhaps because of this, no one has ever > > written the code to do it. I've only ever needed to read these files. > > > > Why do you need Swiss aka SwissProt aka UniProt DAT format > > output? > > I guess we should just use GenBank format. > > Traditionally (which in this case means: using BioPerl) we have been > using SwissProt format when writing protein records, so I was just > trying to do the same in BioPython, without any deeper thought going > into it :-) > > > Thanks for the reply & merry christmas, > > Adam > > -- > "There are forty people in the world and five of them Adam Sj=C3= =B8gren > are hamburgers." [email protected]= k > > _______________________________________________ > Biopython-dev mailing list > [email protected] > http://mailman.open-bio.org/mailman/listinfo/biopython-dev > --94eb2c0599400328730561319d61 Content-Type: text/html; charset="UTF-8" Content-Transfer-Encoding: quoted-printable <div dir=3D"ltr">Thanks for the explanation.<div><br></div><div>Note I have= sometimes needed to *filter* a SwissProt</div><div>format file, and you ca= n do this via SeqIO's index</div><div>get_raw functionality (grab the o= riginal raw record</div><div>as a string, and save it to a file as is). The= re's an</div><div>example of this in our Tutorial, see the section</div= ><div>"Getting the raw data for a record":</div><div><br></div><d= iv><a href=3D"http://biopython.org/DIST/docs/tutorial/Tutorial.html">http:/= /biopython.org/DIST/docs/tutorial/Tutorial.html</a><br></div><div><br></div= ><div>Peter</div></div><div class=3D"gmail_extra"><br><div class=3D"gmail_q= uote">On Mon, Dec 25, 2017 at 4:05 PM, Adam Sj=C3=B8gren <span dir=3D"ltr">= <<a href=3D"mailto:[email protected]" target=3D"_blank">[email protected]= k</a>></span> wrote:<br><blockquote class=3D"gmail_quote" style=3D"margi= n:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span class=3D"">= Peter writes:<br> <br> > I don't recall anyone asking about this, so there doesn't seem= to be<br> > much need for it. And perhaps because of this, no one has ever<br> > written the code to do it. I've only ever needed to read these fil= es.<br> ><br> > Why do you need Swiss aka SwissProt aka UniProt DAT format<br> > output?<br> <br> </span>I guess we should just use GenBank format.<br> <br> Traditionally (which in this case means: using BioPerl) we have been<br> using SwissProt format when writing protein records, so I was just<br> trying to do the same in BioPython, without any deeper thought going<br> into it :-)<br> <br> <br> =C2=A0 Thanks for the reply & merry christmas,<br> <br> =C2=A0 =C2=A0 Adam<br> <span class=3D"HOEnZb"><font color=3D"#888888"><br> --<br> =C2=A0"There are forty people in the world and five of them=C2=A0 =C2= =A0 =C2=A0 =C2=A0 Adam Sj=C3=B8gren<br> =C2=A0 are hamburgers."=C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2= =A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 = =C2=A0 =C2=A0 =C2=A0<a href=3D"mailto:[email protected]">[email protected]<= /a><br> </font></span><div class=3D"HOEnZb"><div class=3D"h5"><br> ______________________________<wbr>_________________<br> Biopython-dev mailing list<br> <a href=3D"mailto:[email protected]">Biopython-dev@mailman= .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>mailman/= listinfo/biopython-dev</a></div></div></blockquote></div><br></div> --94eb2c0599400328730561319d61-- --===============3261069272977607231== 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 --===============3261069272977607231==--