Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- from Bio import SeqIO
- filename = "rosalind_gc.txt"
- count = 0
- for record in SeqIO.parse(filename, "fasta"):
- seqs = record
- a = seqs[count].count('A')
- t = seqs[count].count('T')
- g = seqs[count].count('G')
- c = seqs[count].count('C')
- Q = (g+c)/(a+g+t+c)
- print("Record " + record.id + ", GC content " + str(Q))
- count = count + 1
Add Comment
Please, Sign In to add comment