生物信息基因处理工具 Biopython、pysam

    注意：
    pysam暂时不支持windows，在linux直接安装pip安装

    ## 读取
    from Bio import SeqIO
    
    for record in SeqIO.parse("example.fasta", "fasta"):
    print(record.id)
    ## SeqIO.read("NC_006581.gbk", "genbank")

    ##保存
    from Bio import SeqIO
    
    sequences = ...  # add code here
    SeqIO.write(sequences, "example.fasta", "fasta")

    ##序列转换，gb格式转fasta格式
    from Bio import SeqIO
    
    count = SeqIO.convert("cor6_6.gb", "genbank", "cor6_6.fasta", "fasta")
    print("Converted %i records" % count)

    from Bio.Seq import Seq
    
    # 新建一个DNA序列对象
    dna_seq = Seq("GGATGGTTGTCTATTAACTTGTTCAAAAAAGTATCAGGAGTTGTCAAGGCAGAGAAGAGAGTGTTTGCA")
    # # 新建一个RNA序列对象
    rna_seq = Seq("UGCAAACACUCUCUUC")
    # # # 新建一个蛋白质序列对象
    protein_seq = Seq("GWLSINLFKKVSGVVKAEKRVFA")
    
    
    dna_seq.translate(),dna_seq.transcribe(),dna_seq.reverse_complement(),dna_seq.reverse_complement_rna()
    
    
    rna_seq,rna_seq.back_transcribe()

    from Bio import PDB
    parser = PDB.PDBParser()
    
    
    pdb1 =r"C:\***ata\7b85.pdb"
    structure = parser.get_structure("7b85", pdb1) 
    model = structure[0]
    res_no = 0
    non_resi = 0
    
    for model in structure:
    for chain in model:
        for r in chain.get_residues():
            if r.id[0] == ' ':
                res_no +=1
            else:
                non_resi +=1
    
    print ("Residues:  %i" % (res_no))
    print ("Other:     %i" % (non_resi))