※ 引述《Arton0306 (Ar藤)》之铭言： : ※ 引述《consley (je ne sais quoi)》之铭言： : : 各位版大好 小弟最近刚开始学python 是透过putty连上unix : : 是有关基因序列的 : : 我用raw_input读入一字串 : : 然後希望在字串中找到'>'并列印出该行文字(head) : : 然後再抓取下一行的基因序列进行後续处理 : : 我的问题是要用何种语法找到'>'并列印出该整行？ : : 然後如何再把下一行抓到另一个字串阵列里进行後续处理？ : : 或是各位高手有更好的处理方法 也请不吝指教 : : Example input: : : AGGTGGCCAAAACTCGTCTT : : Example output: : : length: : : G+C: : : Entropy: : 基因库抓出来的 ">"应该是固定当做一笔资料一个 : 不知你为什麽要用raw_input : 要用这个就用alist.index找 或是if ">" in astring 我也是初学者 我之前写过类似的东西 我觉得一定有更好的方法 我把input 存在 list 里面用 alist = file("input.txt", "r").readlines() start_number = alist.index('>') 然後下一行就是 seq = alist[start_number+1] 抛砖引玉一下 请高手提供更聪明的方法 : 只是>应该都在第一位吧 而且一笔资料就只有一个> : 我以前遇到的是这样 : 我是直接读档 f.read().split(">") : 这样就把每一笔基因资料拆开了 --

※ 发信站: 批踢踢实业坊(ptt.cc) ◆ From: 205.175.114.61