Guest User

Untitled

a guest
Jun 24th, 2018
118
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.35 KB | None | 0 0
  1. from Bio import SeqIO
  2. filename = "rosalind_gc.txt"
  3. count = 0
  4.  
  5. for record in SeqIO.parse(filename, "fasta"):
  6. seqs = record
  7. a = seqs[count].count('A')
  8. t = seqs[count].count('T')
  9. g = seqs[count].count('G')
  10. c = seqs[count].count('C')
  11.  
  12. Q = (g+c)/(a+g+t+c)
  13.  
  14. print("Record " + record.id + ", GC content " + str(Q))
  15. count = count + 1
Add Comment
Please, Sign In to add comment