Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- peaks = []
- with open("His.PSC.50.AllAg.AllCell.bed")as f:
- for line in f:
- peaks.append(line.strip().split())
- f2 = open('H3K36me3Peaks.bed', 'w')
- for item in peaks:
- if item[3].find('H3K36me3') > 0:
- d.write(item[0] + ' ' + item[1] + ' ' + item[2] + '\n')
- f2.close()
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement