Python | 統(tǒng)計(jì)FAST文件里多個(gè)基因的長(zhǎng)度

假如以下這個(gè)腳本為Gene_length.py, 測(cè)試的fasta文件為test.fasta,輸出文件為outfile

代碼的運(yùn)行命令（在linux中）：python Gene_length.py test.fasta outfile

* 其中 ID = re.split(r'\.',line)[2] 這一步是因?yàn)樵糵asta文件里ID名稱是這樣的>evm.model.Contig62652.1.我們需要的是Contig62652這一部分，這個(gè)ID是用.來(lái)分割，而我們需要的是第三部分，在列表中是屬于【2】

代碼實(shí)例
最后編輯于 ：
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
【社區(qū)內(nèi)容提示】社區(qū)部分內(nèi)容疑似由AI輔助生成，瀏覽時(shí)請(qǐng)結(jié)合常識(shí)與多方信息審慎甄別。
平臺(tái)聲明：文章內(nèi)容（如有圖片或視頻亦包括在內(nèi)）由作者上傳并發(fā)布，文章內(nèi)容僅代表作者本人觀點(diǎn)，簡(jiǎn)書系信息發(fā)布平臺(tái)，僅提供信息存儲(chǔ)服務(wù)。

友情鏈接更多精彩內(nèi)容