BioPython的安装和使用

BioPython 是一个用来处理序列和生物信息的python包,里面包含了很多的工具,可以用来直接读取fasta格式。安装可以通过两种方式,pip方式:

1. pip 方式

pip3 install biopython

打开python终端,

>>> import Bio
>>> from Bio import SeqIO

报错:

---------------------------------------------------------------------------
ImportError                               Traceback (most recent call last)
<ipython-input-5-5d24a7c49c42> in <module>
----> 1 from Bio import SeqIO

ImportError: cannot import name 'SeqIO' from 'Bio' (unknown location)

Bio可以导入,但package显示未知的位置,查看biopython包的位置:

>>> import Bio
>>> help(Bio)

输出:

Help on package Bio:

NAME
    Bio

PACKAGE CONTENTS


FILE
    (built-in)

没有发现包,由于我用的是jupyter notebook,大家可以关闭kernel然后重启,或者关闭重新打开了一个页面,再测试。如果还有问题,尝试使用第二个方式:

2. 使用conda安装

确保计算机安装了anaconda 或者miniconda
直接在终端中输入:

conda install -c conda-forge biopython

出现 Proceed ([y]/n)?时按y,安装成功。
此时进入python 终端,测试

Python 3.7.6 | packaged by conda-forge | (default, Jan  7 2020, 22:33:48) 
[GCC 7.3.0] on linux
Type "help", "copyright", "credits" or "license" for more information.
>>> from Bio import SeqIO
>>> 

没有问题。如果使用的jupernotebook,要重启kernel,然后测试。观察包的位置

>>> from Bio import SeqIO
>>> import Bio
>>> help(Bio)
>

out:

Help on package Bio:

NAME
    Bio - Collection of modules for dealing with biological data in Python.

DESCRIPTION
    The Biopython Project is an international association of developers
    of freely available Python tools for computational molecular biology.
    
    http://biopython.org

PACKAGE CONTENTS
    Affy (package)
    Align (package)
    AlignIO (package)
    Alphabet (package)
    Application (package)
    Blast (package)
    CAPS (package)
    Cluster (package)
    Compass (package)
    Crystal (package)
    Data (package)
    Emboss (package)
    Entrez (package)
    ExPASy (package)
    FSSP (package)
    File
    GenBank (package)
    Geo (package)
    Graphics (package)
    HMM (package)
    Index
    KDTree (package)
    KEGG (package)
    LogisticRegression
    MarkovModel
    MaxEntropy
    Medline (package)
    NMR (package)
    NaiveBayes
    Nexus (package)
    PDB (package)
    Pathway (package)
    Phylo (package)
    PopGen (package)
    Restriction (package)
    SCOP (package)
    .....

安装成功

3. 使用biopython读取fasta格式

from Bio import SeqIO
for seq_record in SeqIO.parse('../input/example.fa', "fasta"):
    print(seq_record.id)
    print(seq_record.seq)

输出:

ENST00000435737.5
ATGTTTCGCATCACCAACATTGAGTTTCTTCCCGAATACCGACAAAAGGAGTCCAGGGAATTTCTTTCAGTGTCACGGACTGTGCAGCAAGTGATAAACCTGGTTTATACAACATCTGCCTTCTCCAAATTTTATGAGCAGTCTGTTGTTGCAGATGTCAGCAACAACAAAGGCGGCCTCCTTGTCCACTTTTGGATTGTTTTTGTCATGCCACGTGCCAAAGGCCACATCTTCTGTGAAGACTGTGTTGCCGCCATCTTGAAGGACTCCATCCAGACAAGCATCATAAACCGGACCTCTGTGGGGAGCTTGCAGGGACTGGCTGTGGACATGGACTCTGTGGTACTAAATGAAGTCCTGGGGCTGACTCTCATTGTCTGGATTGACTGA
ENST00000419127.5
ATGTTTCGCATCACCAACATTGAGTTTCTTCCCGAATACCGACAAAAGGAGTCCAGGGAATTTCTTTCAGTGTCACGGACTGTGCAGCAAGTGATAAACCTGGTTTATACAACATCTGCCTTCTCCAAATTTTATGAGCAGTCTGTTGTTGCAGATGTCAGCAACAACAAAGGCGGCCTCCTTGTCCACTTTTGGATTGTTTTTGTCATGCCACGTGCCAAAGGCCACATCTTCTGTGAAGACTGTGTTGCCGCCATCTTGAAGGACTCCATCCAGACAAGCATCATAAACCGGACCTCTGTGGGGAGCTTGCAGGGACTGGCTGTGGACATGGACTCTGTGGTACTAAATGACAAAGGCTGCTCTCAGTACTTCTATGCAGAGCATCTGTCTCTCCACTACCCGCTGGAGATTTCTGCAGCCTCAGGGAGGCTGATGTGTCACTTCAAGCTGGTGGCCATAGTGGGCTACCTGATTCGTCTCTCAATCAAGTCCATCCAAATCGAAGCCGACAACTGTGTCACTGACTCCCTGACCATTTACGACTCCCTTTTGCCCATCCGGAGCAGCATCT
....
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值