Exemplo n.º 1
0
    def _extract_fastas_from_fast5(self):
        """
		Return the sequence in the FAST5 file in FASTA format
		"""
        for id, h5path in fastq_paths[self.version].iteritems():
            try:
                table = self.hdf5file[h5path % self.group]
                fa = formats.Fasta(table['Fastq'][()])
                fa.name += " " + self.filename
                self.fastas[id] = fa
            except Exception, e:
                pass
Exemplo n.º 2
0
	def _extract_fastas_from_fast5(self):
		"""
		Return the sequence in the FAST5 file in FASTA format
		"""
		for id, h5path in fastq_paths.iteritems(): 
			try:
				table = self.hdf5file.getNode(h5path)
				fa = formats.Fasta(table.Fastq[()])
				fa.name += "_" + id + " " + self.filename
				self.fastas[id] = fa
			except Exception, e:
				pass