Upgrade to Pro — share decks privately, control downloads, hide ads and more …

NextBio @Taipei.py

Liang Bo Wang
August 29, 2013
65

NextBio @Taipei.py

NextBio Python Library and Sprint! introduction

Liang Bo Wang

August 29, 2013
Tweet

Transcript

  1. Biopython Tutorial and Cookbook Jeff Chang, Brad Chapman, Iddo Friedberg,

    Thomas Hamelryck, Michiel de Hoon, Peter Cock, Tiago Antao, Eric Talevich, Bartek Wilczyński Last Update – 28 August 2013 (Biopython 1.62) Contents Chapter 1 Introduction 1.1 What is Biopython? 1.2 What can I find in the Biopython package 1.3 Installing Biopython 1.4 Frequently Asked Questions (FAQ) Chapter 2 Quick Start – What can you do with Biopython? 2.1 General overview of what Biopython provides 2.2 Working with sequences 2.3 A usage example 2.4 Parsing sequence file formats 2.4.1 Simple FASTA parsing example 2.4.2 Simple GenBank parsing example 2.4.3 I love parsing – please don’t stop talking about it! 2.5 Connecting with biological databases 2.6 What to do next Chapter 3 Sequence objects 3.1 Sequences and Alphabets 3.2 Sequences act like strings 3.3 Slicing a sequence 3.4 Turning Seq objects into strings 3.5 Concatenating or adding sequences 3.6 Changing case 3.7 Nucleotide sequences and (reverse) complements 3.8 Transcription 3.9 Translation 3.10 Translation Tables 3.11 Comparing Seq objects 3.12 MutableSeq objects 3.13 UnknownSeq objects 3.14 Working with directly strings Chapter 4 Sequence annotation objects 4.1 The SeqRecord object 4.2 Creating a SeqRecord 4.2.1 SeqRecord objects from scratch 4.2.2 SeqRecord objects from FASTA files 4.2.3 SeqRecord objects from GenBank files 4.3 Feature, location and position objects 4.3.1 SeqFeature objects 4.3.2 Positions and locations 4.3.3 Sequence described by a feature or location 4.4 References 4.5 The format method 4.6 Slicing a SeqRecord 4.7 Adding SeqRecord objects 4.8 Reverse-complementing SeqRecord objects Chapter 5 Sequence Input/Output 5.1 Parsing or Reading Sequences 5.1.1 Reading Sequence Files 5.1.2 Iterating over the records in a sequence file 5.1.3 Getting a list of the records in a sequence file 5.1.4 Extracting data 5.2 Parsing sequences from compressed files 5.3 Parsing sequences from the net 5.3.1 Parsing GenBank records from the net 5.3.2 Parsing SwissProt sequences from the net 5.4 Sequence files as Dictionaries 5.4.1 Sequence files as Dictionaries – In memory 5.4.2 Sequence files as Dictionaries – Indexed files 5.4.3 Sequence files as Dictionaries – Database indexed files 5.4.4 Indexing compressed files 5.4.5 Discussion 5.5 Writing Sequence Files 5.5.1 Round trips 5.5.2 Converting between sequence file formats 5.5.3 Converting a file of sequences to their reverse complements 5.5.4 Getting your SeqRecord objects as formatted strings Chapter 6 Multiple Sequence Alignment objects 6.1 Parsing or Reading Sequence Alignments 6.1.1 Single Alignments 6.1.2 Multiple Alignments 6.1.3 Ambiguous Alignments 6.2 Writing Alignments 6.2.1 Converting between sequence alignment file formats 6.2.2 Getting your alignment objects as formatted strings 6.3 Manipulating Alignments 6.3.1 Slicing alignments 6.3.2 Alignments as arrays 6.4 Alignment Tools 6.4.1 ClustalW 6.4.2 MUSCLE 6.4.3 MUSCLE using stdout 6.4.4 MUSCLE using stdin and stdout 6.4.5 EMBOSS needle and water Chapter 7 BLAST 7.1 Running BLAST over the Internet 7.2 Running BLAST locally 7.2.1 Introduction 7.2.2 Standalone NCBI “legacy” BLAST 7.2.3 Standalone NCBI BLAST+ 7.2.4 WU-BLAST and AB-BLAST 7.3 Parsing BLAST output 7.4 The BLAST record class 7.5 Deprecated BLAST parsers 7.5.1 Parsing plain-text BLAST output 7.5.2 Parsing a plain-text BLAST file full of BLAST runs 7.5.3 Finding a bad record somewhere in a huge plain-text BLAST file 7.6 Dealing with PSI-BLAST 7.7 Dealing with RPS-BLAST Chapter 8 BLAST and other sequence search tools (experimental code) 8.1 The SearchIO object model 8.1.1 QueryResult 8.1.2 Hit 8.1.3 HSP 8.1.4 HSPFragment 8.2 A note about standards and conventions 8.3 Reading search output files 8.4 Dealing with large search output files with indexing 8.5 Writing and converting search output files Chapter 9 Accessing NCBI’s Entrez databases 9.1 Entrez Guidelines 9.2 EInfo: Obtaining information about the Entrez databases 9.3 ESearch: Searching the Entrez databases 9.4 EPost: Uploading a list of identifiers 9.5 ESummary: Retrieving summaries from primary IDs 9.6 EFetch: Downloading full records from Entrez 9.7 ELink: Searching for related items in NCBI Entrez 9.8 EGQuery: Global Query - counts for search terms 9.9 ESpell: Obtaining spelling suggestions 9.10 Parsing huge Entrez XML files 9.11 Handling errors 9.12 Specialized parsers 9.12.1 Parsing Medline records 9.12.2 Parsing GEO records 9.12.3 Parsing UniGene records 9.13 Using a proxy 9.14 Examples 9.14.1 PubMed and Medline 9.14.2 Searching, downloading, and parsing Entrez Nucleotide records 9.14.3 Searching, downloading, and parsing GenBank records 9.14.4 Finding the lineage of an organism 9.15 Using the history and WebEnv 9.15.1 Searching for and downloading sequences using the history 9.15.2 Searching for and downloading abstracts using the history 9.15.3 Searching for citations Chapter 10 Swiss-Prot and ExPASy 10.1 Parsing Swiss-Prot files 10.1.1 Parsing Swiss-Prot records 10.1.2 Parsing the Swiss-Prot keyword and category list 10.2 Parsing Prosite records 10.3 Parsing Prosite documentation records 10.4 Parsing Enzyme records 10.5 Accessing the ExPASy server 10.5.1 Retrieving a Swiss-Prot record 10.5.2 Searching Swiss-Prot 10.5.3 Retrieving Prosite and Prosite documentation records 10.6 Scanning the Prosite database Chapter 11 Going 3D: The PDB module 11.1 Reading and writing crystal structure files 11.1.1 Reading a PDB file 11.1.2 Reading an mmCIF file 11.1.3 Reading files in the PDB XML format 11.1.4 Writing PDB files 11.2 Structure representation 11.2.1 Structure 11.2.2 Model 11.2.3 Chain 11.2.4 Residue 11.2.5 Atom 11.2.6 Extracting a specific Atom/Residue/Chain/Model from a Structure 11.3 Disorder 11.3.1 General approach 11.3.2 Disordered atoms
  2. The ability to parse bioinformatics files into Python utilizable data

    structures, including support for the following formats: Blast output – both from standalone and WWW Blast Clustalw FASTA GenBank PubMed and Medline ExPASy files, like Enzyme and Prosite SCOP, including ‘`dom’ and ‘`lin’ files UniGene SwissProt Files in the supported formats can be iterated over record by record or indexed and accessed via a Dictionary interface. Code to deal with popular on-line bioinformatics destinations such as: NCBI – Blast, Entrez and PubMed services ExPASy – Swiss-Prot and Prosite entries, as well as Prosite searches Interfaces to common bioinformatics programs such as: Standalone Blast from NCBI Clustalw alignment program EMBOSS command line tools A standard sequence class that deals with sequences, ids on sequences, and sequence features. Tools for performing common operations on sequences, such as translation, transcription and weight calculations. Code to perform classification of data using k Nearest Neighbors, Naive Bayes or Support Vector Machines. Code for dealing with alignments, including a standard way to create and deal with substitution matrices. Code making it easy to split up parallelizable tasks into separate processes. GUI-based programs to do basic sequence manipulations, translations, BLASTing, etc. Extensive documentation and help with using the modules, including this file, on-line wiki documentation, the web site, and the mailing list. Integration with BioSQL, a sequence database schema also supported by the BioPerl and BioJava projects. We hope this gives you plenty of reasons to download and start using Biopython! 1.3 Installing Biopython All of the installation information for Biopython was separated from this document to make it easier to keep updated. The short version is go to our downloads page (http://biopython.org/wiki/Download), download and install the listed dependencies, then download and install Biopython. Biopython runs on many platforms (Windows, Mac, and on the various flavors of Linux and Unix). For Windows we provide pre-compiled click-and-run installers, while for Unix and other operating systems you must install from source as described in the included README file. This is usually as simple as the standard commands: python setup.py build python setup.py test sudo python setup.py install (You can in fact skip the build and test, and go straight to the install – but its better to make sure everything seems to be working.) The longer version of our installation instructions covers installation of Python, Biopython dependencies and Biopython itself. It is available in PDF (http://biopython.org/DIST/docs/install/Installation.pdf) and HTML formats (http://biopython.org/DIST/docs/install/Installation.html). 1.4 Frequently Asked Questions (FAQ) 1. How do I cite Biopython in a scientific publication? Please cite our application note [1, Cock et al., 2009] as the main Biopython reference. In addition, please cite any publications from the following list if appropriate, in particular as a reference for specific modules within Biopython (more information can be found on our website): For the official project announcement: [13, Chapman and Chang, 2000]; For Bio.PDB: [18, Hamelryck and Manderick, 2003]; For Bio.Cluster: [14, De Hoon et al., 2004]; For Bio.Graphics.GenomeDiagram: [2, Pritchard et al., 2006]; For Bio.Phylo and Bio.Phylo.PAML: [9, Talevich et al., 2012]; For the FASTQ file format as supported in Biopython, BioPerl, BioRuby, BioJava, and EMBOSS: [7, Cock et al., 2010]. 2. How should I capitalize “Biopython”? Is “BioPython” OK? The correct capitalization is “Biopython”, not “BioPython” (even though that would have matched BioPerl, BioJava and BioRuby). 3. How do I find out what version of Biopython I have installed? Use this: >>> import Bio >>> print Bio.__version__ ... If the “import Bio” line fails, Biopython is not installed. If the second line fails, your version is very out of date. If the version string ends with a plus, you don’t have an official release, but a snapshot of the in development code. 4. Where is the latest version of this document? If you download a Biopython source code archive, it will include the relevant version in both HTML and PDF formats. The latest published version of this document (updated at each release) is online: http://biopython.org/DIST/docs/tutorial/Tutorial.html http://biopython.org/DIST/docs/tutorial/Tutorial.pdf If you are using the very latest unreleased code from our repository you can find copies of the in-progress tutorial here: http://biopython.org/DIST/docs/tutorial/Tutorial-dev.html http://biopython.org/DIST/docs/tutorial/Tutorial-dev.pdf 5. Which “Numerical Python” do I need? For Biopython 1.48 or earlier, you needed the old Numeric module. For Biopython 1.49 onwards, you need the newer NumPy instead. Both Numeric and NumPy can be installed on the same machine fine. See also: http://numpy.scipy.org/ 6. Why is the Seq object missing the (back) transcription & translation methods described in this Tutorial? You need Biopython 1.49 or later. Alternatively, use the Bio.Seq module functions described in Section 3.14. 7. Why is the Seq object missing the upper & lower methods described in this Tutorial? You need Biopython 1.53 or later. Alternatively, use str(my_seq).upper() to get an upper case string. If you need a Seq object, try Seq(str(my_seq).upper()) but be careful about blindly re-using the same alphabet. 8. Why doesn’t the Seq object translation method support the cds option described in this Tutorial? You need Biopython 1.51 or later. 9. Why doesn’t Bio.SeqIO work? It imports fine but there is no parse function etc. You need Biopython 1.43 or later. Older versions did contain some related code under the Bio.SeqIO name which has since been removed - and this is why the import “works”. 10. Why doesn’t Bio.SeqIO.read() work? The module imports fine but there is no read function! You need Biopython 1.45 or later. Or, use Bio.SeqIO.parse(...).next() instead. 11. Why isn’t Bio.AlignIO present? The module import fails! You need Biopython 1.46 or later. 12. What file formats do Bio.SeqIO and Bio.AlignIO read and write? Check the built in docstrings (from Bio import SeqIO , then help(SeqIO) ), or see http://biopython.org/wiki/SeqIO and http://biopython.org/wiki/AlignIO on the wiki for the latest listing. 13. Why don’t the Bio.SeqIO and Bio.AlignIO input functions let me provide a sequence alphabet? You need Biopython 1.49 or later. 14. Why won’t the Bio.SeqIO and Bio.AlignIO functions parse , read and write take filenames? They insist on handles! You need Biopython 1.54 or later, or just use handles explicitly (see Section 22.1). It is especially important to remember to close output handles explicitly after writing your data. 15. Why won’t the Bio.SeqIO.write() and Bio.AlignIO.write() functions accept a single record or alignment? They insist on a list or iterator! You need Biopython 1.54 or later, or just wrap the item with [...] to create a list of one element. 16. Why doesn’t str(...) give me the full sequence of a Seq object? You need Biopython 1.45 or later. Alternatively, rather than str(my_seq) , use my_seq.tostring() (which will also work on recent versions of Biopython). 17. Why doesn’t Bio.Blast work with the latest plain text NCBI blast output? The NCBI keep tweaking the plain text output from the BLAST tools, and keeping our parser up to date is/was an ongoing struggle. If you aren’t using the latest version of Biopython, you could try upgrading. However, we (and the NCBI) recommend you use the XML output instead, which is designed to be read by a computer program. 18. Why doesn’t Bio.Entrez.read() work? The module imports fine but there is no read function! You need Biopython 1.46 or later. 19. Why doesn’t Bio.Entrez.parse() work? The module imports fine but there is no parse function! You need Biopython 1.52 or later. 20. Why has my script using Bio.Entrez.efetch() stopped working? This could be due to NCBI changes in February 2012 introducing EFetch 2.0. First, they changed the default return modes - you probably want to add retmode="text" to your call. Second, they are now stricter about how to provide a list of IDs – Biopython 1.59 onwards turns a list into a comma separated string automatically. 21. Why doesn’t Bio.Blast.NCBIWWW.qblast() give the same results as the NCBI BLAST website? You need to specify the same options – the NCBI often adjust the default settings on the website, and they do not match the QBLAST defaults anymore. Check things like the gap penalties and expectation threshold. 22. Why doesn’t Bio.Blast.NCBIXML.read() work? The module imports but there is no read function! You need Biopython 1.50 or later. Or, use Bio.Blast.NCBIXML.parse(...).next() instead. 23. Why doesn’t my SeqRecord object have a letter_annotations attribute? Per-letter-annotation support was added in Biopython 1.50. 24. Why can’t I slice my SeqRecord to get a sub-record? You need Biopython 1.50 or later. 25. Why can’t I add SeqRecord objects together? You need Biopython 1.53 or later. 26. Why doesn’t Bio.SeqIO.convert() or Bio.AlignIO.convert() work? The modules import fine but there is no convert function! You need Biopython 1.52 or later. Alternatively, combine the parse and write functions as described in this tutorial (see Sections 5.5.2 and 6.2.1). 27. Why doesn’t Bio.SeqIO.index() work? The module imports fine but there is no index function! You need Biopython 1.52 or later. 28. Why doesn’t Bio.SeqIO.index_db() work? The module imports fine but there is no index_db function! You need Biopython 1.57 or later (and a Python with SQLite3 support). 29. Where is the MultipleSeqAlignment object? The Bio.Align module imports fine but this class isn’t there! You need Biopython 1.54 or later. Alternatively, the older Bio.Align.Generic.Alignment class supports some of its functionality, but using this is now discouraged. 30. Why can’t I run command line tools directly from the application wrappers? You need Biopython 1.55 or later. Alternatively, use the Python subprocess module directly. 31. I looked in a directory for code, but I couldn’t find the code that does something. Where’s it hidden? One thing to know is that we put code in __init__.py files. If you are not used to looking for code in this file this can be confusing. The reason we do this is to make the imports easier for users. For instance, instead of having to do a “repetitive” import like from Bio.GenBank import GenBank , you can just use from Bio import GenBank . 32. Why does the code from CVS seem out of date? In late September 2009, just after the release of Biopython 1.52, we switched from using CVS to git, a distributed version control system. The old CVS server will remain available as a static and read only backup, but if you want to grab the latest code, you’ll need to use git instead. See our website for more details.
  3. What is FASTA / Q ? >HWI-H248:87:C1NPGACXX:1:1101:1212:2075 1:N:0:GTGAAA CCTCGAAATACTGGACGATCAACTCCAACTCCCATTGCATTAAGCCCATTGTCAACATA >HWI-H248:87:C1NPGACXX:1:1101:1217:2170

    1:N:0:GTGAAA TCTTCATCAGCAGGAGCAGGAATTGCAGTATAAAGAGGCCAATAGTAGGCACGATCATA @HWI-H248:87:C1NPGACXX:1:1101:1212:2075 1:N:0:GTGAAA CCTCGAAATACTGGACGATCAACTCCAACTCCCATTGCATTAAGCCCATTGTCAACATA + ???BDDD73DDDDE93)@FADCFE?EEE@DEDDDDCCDDCEECDDD=@DDCBDCBBACI @HWI-H248:87:C1NPGACXX:1:1101:1217:2170 1:N:0:GTGAAA TCTTCATCAGCAGGAGCAGGAATTGCAGTATAAAGAGGCCAATAGTAGGCACGATCATA + @@@BDDFBFDFHFGGBHDGFBFEGFGID::CFFB<?CDHGDGGCGFGHGEBFHIGHGIH