Biopython Tutorial and Cookbook - Biopython · Biopython

1 Biopython Tutorial and CookbookJeff Chang, Brad Chapman, Iddo Friedberg, Thomas Hamelryck,Michiel de Hoon, Peter Cock, Tiago Antao, Eric Talevich, Bartek Wilczy nskiLast Update 25 May 2020 ( Biopython )Contents1 What is Biopython ? .. What can I find in the Biopython package .. Installing Biopython .. Frequently Asked Questions (FAQ) ..102 Quick Start What can you do with Biopython ? General overview of what Biopython provides .. Working with sequences .. A usage example .. Parsing sequence file formats .. FASTA parsing example .. GenBank parsing example .. love parsing please don t stop talking about it!

2 Connecting with biological databases .. What to do next ..183 Sequence Sequences and Alphabets .. Sequences act like strings .. Slicing a sequence .. Turning Seq objects into strings .. Concatenating or adding sequences .. Changing case .. Nucleotide sequences and (reverse) complements .. Transcription .. Translation .. Translation Tables .. Comparing Seq objects .. MutableSeq objects .. UnknownSeq objects .. Working with strings directly ..334 Sequence annotation The SeqRecord object .. Creating a SeqRecord .. objects from scratch.

3 Objects from FASTA files .. objects from GenBank files .. Feature, location and position objects .. objects .. and locations .. described by a feature or location .. Comparison .. References .. The format method .. Slicing a SeqRecord .. Adding SeqRecord objects .. Reverse-complementing SeqRecord objects ..495 Sequence Parsing or Reading Sequences .. Sequence Files .. over the records in a sequence file .. a list of the records in a sequence file .. data .. data .. Parsing sequences from compressed files .. Parsing sequences from the net.

4 GenBank records from the net .. SwissProt sequences from the net .. Sequence files as Dictionaries .. files as Dictionaries In memory .. files as Dictionaries Indexed files .. files as Dictionaries Database indexed files .. compressed files .. Writing Sequence Files .. trips .. between sequence file formats .. a file of sequences to their reverse complements .. your SeqRecord objects as formatted strings .. Low level FASTA and FASTQ parsers ..706 Multiple Sequence Alignment Parsing or Reading Sequence Alignments .. Alignments .. Alignments .. Alignments.

5 Writing Alignments .. between sequence alignment file formats .. your alignment objects as formatted strings .. Manipulating Alignments .. alignments .. as arrays .. Alignment Tools .. using stdout .. using stdin and stdout .. needle and water .. Pairwise sequence alignment .. Substitution matrices .. 1127 Running BLAST over the Internet .. Running BLAST locally .. NCBI BLAST+ .. versions of BLAST .. Parsing BLAST output .. The BLAST record class .. Dealing with PSI-BLAST .. Dealing with RPS-BLAST .. 1268 BLAST and other sequence search The SearchIO object model.

6 A note about standards and conventions .. Reading search output files .. Dealing with large search output files with indexing .. Writing and converting search output files .. 1449 Accessing NCBI s Entrez Entrez Guidelines .. EInfo: Obtaining information about the Entrez databases .. ESearch: Searching the Entrez databases .. EPost: Uploading a list of identifiers .. ESummary: Retrieving summaries from primary IDs .. EFetch: Downloading full records from Entrez .. ELink: Searching for related items in NCBI Entrez .. EGQuery: Global Query - counts for search terms.

7 ESpell: Obtaining spelling suggestions .. Parsing huge Entrez XML files .. HTML escape characters .. Handling errors .. Specialized parsers .. Parsing Medline records .. Parsing GEO records .. Parsing UniGene records .. Using a proxy .. Examples .. PubMed and Medline .. Searching, downloading, and parsing Entrez Nucleotide records .. Searching, downloading, and parsing GenBank records .. Finding the lineage of an organism .. Using the history and WebEnv .. Searching for and downloading sequences using the history .. Searching for and downloading abstracts using the history.

8 Searching for citations .. 17310 Swiss-Prot and Parsing Swiss-Prot files .. Parsing Swiss-Prot records .. Parsing the Swiss-Prot keyword and category list .. Parsing Prosite records .. Parsing Prosite documentation records .. Parsing Enzyme records .. Accessing the ExPASy server .. Retrieving a Swiss-Prot record .. Searching Swiss-Prot .. Retrieving Prosite and Prosite documentation records .. Scanning the Prosite database .. 18211 Going 3D: The PDB Reading and writing crystal structure files .. Reading an mmCIF file .. Reading files in the MMTF format.

9 Reading a PDB file .. Reading a PQR file .. Reading files in the PDB XML format .. Writing mmCIF files .. Writing PDB files .. Writing PQR files .. Writing MMTF files .. Structure representation .. Structure .. Model .. Chain .. Residue .. Atom .. Extracting a specificAtom/Residue/Chain/Modelfrom a Structure .. Disorder .. General approach .. Disordered atoms .. Disordered residues .. Hetero residues .. Associated problems .. Water residues .. Other hetero residues .. Navigating through a Structure object.

10 Analyzing structures .. Measuring distances .. Measuring angles .. Measuring torsion angles .. Internal coordinates for standard residues .. Determining atom-atom contacts .. Superimposing two structures .. Mapping the residues of two related structures onto each other .. Calculating the Half Sphere Exposure .. Determining the secondary structure .. Calculating the residue depth .. Common problems in PDB files .. Examples .. Automatic correction .. Fatal errors .. Accessing the Protein Data Bank .. Downloading structures from the Protein Data Bank.