python實現根據序列ID從提取fasta文件序列

當序列少的時候，我習慣用 grep -A 1 -f seq.lst seq.fas | sed ‘/^–$/d’ > out.fas提取，可是此次遇到了一個大文件，用grep就太費時了，而後又試了一下TBtools的提取序列功能，發現時間也很長，因此就寫了個python。提取將近100萬條reads耗時也就須要10s左右python #!/usr/bin/python # -*- coding:

>>阅读原文<<