<?xml version="1.0" encoding="ascii"?> <!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "DTD/xhtml1-transitional.dtd"> <html xmlns="http://www.w3.org/1999/xhtml" xml:lang="en" lang="en"> <head> <title>Bio.Sequencing.Ace</title> <link rel="stylesheet" href="epydoc.css" type="text/css" /> <script type="text/javascript" src="epydoc.js"></script> </head> <body bgcolor="white" text="black" link="blue" vlink="#204080" alink="#204080"> <!-- ==================== NAVIGATION BAR ==================== --> <table class="navbar" border="0" width="100%" cellpadding="0" bgcolor="#a0c0ff" cellspacing="0"> <tr valign="middle"> <!-- Tree link --> <th> <a href="module-tree.html">Trees</a> </th> <!-- Index link --> <th> <a href="identifier-index.html">Indices</a> </th> <!-- Help link --> <th> <a href="help.html">Help</a> </th> <th class="navbar" width="100%"></th> </tr> </table> <table width="100%" cellpadding="0" cellspacing="0"> <tr valign="top"> <td width="100%"> <span class="breadcrumbs"> <a href="Bio-module.html">Package Bio</a> :: <a href="Bio.Sequencing-module.html">Package Sequencing</a> :: Module Ace </span> </td> <td> <table cellpadding="0" cellspacing="0"> <!-- hide/show private --> <tr><td align="right"><span class="options">[<a href="javascript:void(0);" class="privatelink" onclick="toggle_private();">hide private</a>]</span></td></tr> <tr><td align="right"><span class="options" >[<a href="frames.html" target="_top">frames</a >] | <a href="Bio.Sequencing.Ace-module.html" target="_top">no frames</a>]</span></td></tr> </table> </td> </tr> </table> <!-- ==================== MODULE DESCRIPTION ==================== --> <h1 class="epydoc">Module Ace</h1><p class="nomargin-top"><span class="codelink"><a href="Bio.Sequencing.Ace-pysrc.html">source code</a></span></p> <pre class="literalblock"> Parser for (new) ACE files output by PHRAP. version 1.3, 05/06/2004 Written by Frank Kauff (fkauff@duke.edu) and Cymon J. Cox (cymon@duke.edu) Uses the Biopython Parser interface: ParserSupport.py Usage: There are two ways of reading an ace file: The ACEParser() reads the whole file at once and the RecordParser() reads contig after contig. 1) Parse whole ace file at once: from Bio.Sequencing import Ace aceparser=Ace.ACEParser() acefilerecord=aceparser.parse(open('my_ace_file.ace','r')) This gives you: acefilerecord.ncontigs (the number of contigs in the ace file) acefilerecord.nreads (the number of reads in the ace file) acefilerecord.contigs[] (one instance of the Contig class for each contig) The Contig class holds the info of the CO tag, CT and WA tags, and all the reads used for this contig in a list of instances of the Read class, e.g.: contig3=acefilerecord.contigs[2] read4=contig3.reads[3] RD_of_read4=read4.rd DS_of_read4=read4.ds CT, WA, RT tags from the end of the file can appear anywhere are automatically sorted into the right place. see _RecordConsumer for details. 2) Or you can iterate over the contigs of an ace file one by one in the ususal way: from Bio.Sequencing import Ace recordparser=Ace.RecordParser() iterator=Ace.Iterator(open('my_ace_file.ace','r'),recordparser) for contig in iterator : print contig.name ... Please note that for memory efficiency, when using the iterator approach, only one contig is kept in memory at once. However, there can be a footer to the ACE file containing WA, CT, RT or WR tags which contain additional meta-data on the contigs. Because the parser doesn't see this data until the final record, it cannot be added to the appropriate records. Instead these tags will be returned with the last contig record. Thus an ace file does not entirerly suit the concept of iterating. If WA, CT, RT, WR tags are needed, the ACEParser instead of the RecordParser might be appropriate. </pre> <!-- ==================== CLASSES ==================== --> <a name="section-Classes"></a> <table class="summary" border="1" cellpadding="3" cellspacing="0" width="100%" bgcolor="white"> <tr bgcolor="#70b0f0" class="table-header"> <td colspan="2" class="table-header"> <table border="0" cellpadding="0" cellspacing="0" width="100%"> <tr valign="top"> <td align="left"><span class="table-header">Classes</span></td> <td align="right" valign="top" ><span class="options">[<a href="#section-Classes" class="privatelink" onclick="toggle_private();" >hide private</a>]</span></td> </tr> </table> </td> </tr> <tr> <td width="15%" align="right" valign="top" class="summary"> <span class="summary-type"> </span> </td><td class="summary"> <a href="Bio.Sequencing.Ace.rd-class.html" class="summary-name">rd</a><br /> RD (reads), store a read with its name, sequence etc. </td> </tr> <tr> <td width="15%" align="right" valign="top" class="summary"> <span class="summary-type"> </span> </td><td class="summary"> <a href="Bio.Sequencing.Ace.qa-class.html" class="summary-name">qa</a><br /> QA (read quality), including which part if any was used as the consensus. </td> </tr> <tr> <td width="15%" align="right" valign="top" class="summary"> <span class="summary-type"> </span> </td><td class="summary"> <a href="Bio.Sequencing.Ace.ds-class.html" class="summary-name">ds</a><br /> DS lines, include file name of a read's chromatogram file. </td> </tr> <tr> <td width="15%" align="right" valign="top" class="summary"> <span class="summary-type"> </span> </td><td class="summary"> <a href="Bio.Sequencing.Ace.af-class.html" class="summary-name">af</a><br /> AF lines, define the location of the read within the contig. </td> </tr> <tr> <td width="15%" align="right" valign="top" class="summary"> <span class="summary-type"> </span> </td><td class="summary"> <a href="Bio.Sequencing.Ace.bs-class.html" class="summary-name">bs</a><br /> "BS (base segment), which read was chosen as the consensus at each position. </td> </tr> <tr> <td width="15%" align="right" valign="top" class="summary"> <span class="summary-type"> </span> </td><td class="summary"> <a href="Bio.Sequencing.Ace.rt-class.html" class="summary-name">rt</a><br /> RT (transient read tags), generated by crossmatch and phrap. </td> </tr> <tr> <td width="15%" align="right" valign="top" class="summary"> <span class="summary-type"> </span> </td><td class="summary"> <a href="Bio.Sequencing.Ace.ct-class.html" class="summary-name">ct</a><br /> CT (consensus tags). </td> </tr> <tr> <td width="15%" align="right" valign="top" class="summary"> <span class="summary-type"> </span> </td><td class="summary"> <a href="Bio.Sequencing.Ace.wa-class.html" class="summary-name">wa</a><br /> WA (whole assembly tag), holds the assembly program name, version, etc. </td> </tr> <tr> <td width="15%" align="right" valign="top" class="summary"> <span class="summary-type"> </span> </td><td class="summary"> <a href="Bio.Sequencing.Ace.wr-class.html" class="summary-name">wr</a><br /> WR lines. </td> </tr> <tr> <td width="15%" align="right" valign="top" class="summary"> <span class="summary-type"> </span> </td><td class="summary"> <a href="Bio.Sequencing.Ace.Reads-class.html" class="summary-name">Reads</a><br /> Holds information about a read supporting an ACE contig. </td> </tr> <tr> <td width="15%" align="right" valign="top" class="summary"> <span class="summary-type"> </span> </td><td class="summary"> <a href="Bio.Sequencing.Ace.Contig-class.html" class="summary-name">Contig</a><br /> Holds information about a contig from an ACE record. </td> </tr> <tr> <td width="15%" align="right" valign="top" class="summary"> <span class="summary-type"> </span> </td><td class="summary"> <a href="Bio.Sequencing.Ace.Iterator-class.html" class="summary-name">Iterator</a><br /> Iterates over an ACE-file with multiple contigs. </td> </tr> <tr> <td width="15%" align="right" valign="top" class="summary"> <span class="summary-type"> </span> </td><td class="summary"> <a href="Bio.Sequencing.Ace.RecordParser-class.html" class="summary-name">RecordParser</a><br /> Parses ACE file data into a Record object. </td> </tr> <tr> <td width="15%" align="right" valign="top" class="summary"> <span class="summary-type"> </span> </td><td class="summary"> <a href="Bio.Sequencing.Ace.ACEFileRecord-class.html" class="summary-name">ACEFileRecord</a><br /> Holds data of an ACE file. </td> </tr> <tr> <td width="15%" align="right" valign="top" class="summary"> <span class="summary-type"> </span> </td><td class="summary"> <a href="Bio.Sequencing.Ace.ACEParser-class.html" class="summary-name">ACEParser</a><br /> Parses full ACE file in list of contigs. </td> </tr> <tr class="private"> <td width="15%" align="right" valign="top" class="summary"> <span class="summary-type"> </span> </td><td class="summary"> <a href="Bio.Sequencing.Ace._Scanner-class.html" class="summary-name" onclick="show_private();">_Scanner</a><br /> Scans an ACE-formatted file. </td> </tr> <tr class="private"> <td width="15%" align="right" valign="top" class="summary"> <span class="summary-type"> </span> </td><td class="summary"> <a href="Bio.Sequencing.Ace._RecordConsumer-class.html" class="summary-name" onclick="show_private();">_RecordConsumer</a><br /> Reads the ace tags into data records. </td> </tr> </table> <!-- ==================== VARIABLES ==================== --> <a name="section-Variables"></a> <table class="summary" border="1" cellpadding="3" cellspacing="0" width="100%" bgcolor="white"> <tr bgcolor="#70b0f0" class="table-header"> <td colspan="2" class="table-header"> <table border="0" cellpadding="0" cellspacing="0" width="100%"> <tr valign="top"> <td align="left"><span class="table-header">Variables</span></td> <td align="right" valign="top" ><span class="options">[<a href="#section-Variables" class="privatelink" onclick="toggle_private();" >hide private</a>]</span></td> </tr> </table> </td> </tr> <tr> <td width="15%" align="right" valign="top" class="summary"> <span class="summary-type"> </span> </td><td class="summary"> <a name="StringTypes"></a><span class="summary-name">StringTypes</span> = <code title="(<type 'str'>, <type 'unicode'>)"><code class="variable-group">(</code><type 'str'><code class="variable-op">, </code><type 'unicode'><code class="variable-group">)</code></code> </td> </tr> <tr> <td width="15%" align="right" valign="top" class="summary"> <span class="summary-type"> </span> </td><td class="summary"> <a name="xml_support"></a><span class="summary-name">xml_support</span> = <code title="1">1</code> </td> </tr> </table> <!-- ==================== NAVIGATION BAR ==================== --> <table class="navbar" border="0" width="100%" cellpadding="0" bgcolor="#a0c0ff" cellspacing="0"> <tr valign="middle"> <!-- Tree link --> <th> <a href="module-tree.html">Trees</a> </th> <!-- Index link --> <th> <a href="identifier-index.html">Indices</a> </th> <!-- Help link --> <th> <a href="help.html">Help</a> </th> <th class="navbar" width="100%"></th> </tr> </table> <table border="0" cellpadding="0" cellspacing="0" width="100%%"> <tr> <td align="left" class="footer"> Generated by Epydoc 3.0.1 on Mon Sep 15 09:26:25 2008 </td> <td align="right" class="footer"> <a target="mainFrame" href="http://epydoc.sourceforge.net" >http://epydoc.sourceforge.net</a> </td> </tr> </table> <script type="text/javascript"> <!-- // Private objects are initially displayed (because if // javascript is turned off then we want them to be // visible); but by default, we want to hide them. So hide // them unless we have a cookie that says to show them. checkCookie(); // --> </script> </body> </html>