def files(self, labels, repertoire): ''' Provide a list of list of sequences ''' sequences = [] for filenames in repertoire: current_sequences = [] for filename in filenames: current_sequences += list(fasta.parse_as_Fasta(open(filename))) sequences.append(current_sequences) return self(labels, sequences)
''' python revcomp-fasta.py < in.fa > out.fa ''' import fasta import sys for f in fasta.parse_as_Fasta(sys.stdin): f.revcomp() sys.stdout.write(str(f))