diff --git a/source/_static/code/fasta_object.py b/source/_static/code/fasta_object.py index a1a32c8ad3dfa247a16386ce45e9df5503d07f8d..abb474a9f34b3f34295c3057977ccfd213b24757 100644 --- a/source/_static/code/fasta_object.py +++ b/source/_static/code/fasta_object.py @@ -5,6 +5,9 @@ class Sequence(object): self.comment = comment self.sequence = sequence + def gc_percent(self): + seq = self.sequence.upper() + return float(seq.count('G') + seq.count('C')) / float(len(seq)) class FastaParser(object): @@ -73,4 +76,5 @@ if __name__ == '__main__': fasta_parser = FastaParser(fasta_path) for sequence in fasta_parser: print "----------------" - print sequence.id \ No newline at end of file + print "{seqid} = {gc:.3%}".format(gc=sequence.gc_percent(), + seqid = sequence.id) \ No newline at end of file