Python实现 提取fasta中Homosapiens所在行标题和序列

最新推荐文章于 2021-07-22 20:41:34 发布
最新推荐文章于 2021-07-22 20:41:34 发布
阅读量575 收藏 1
点赞数
分类专栏： Python
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_39212193/article/details/116083466
版权 
file_r = open('SwissProt.fasta', 'r')
file_w = open('SwissProtHomo.fasta', 'w')
seq = ''
for line in file_r:
    if line[0] == '>' and seq == '':   #若seq为空，该行只有标题没有序列
        header = line      #赋予行的内容给header
    elif line[0] != '>':   #若该行没有标题只有序列
        seq = seq+line     #添加行的内容到seq序列
    elif line[0] == '>' and seq != '':   #有标题有序列
        if "Homo sapiens" in header:
            file_w.write(header+seq)
        seq =''   #seq初始化为空字符串
        header = line  #重新初始化到当前行 ，初始化目的：每次循环的header和seq的值都不一样
if "Homo sapiens" in header:
        file_w.write(header+seq)
file_w.close()

结果

>sp|Q66K14|TBC9B_HUMAN TBC1 domain family member 9B OS=Homo sapiens OX=9606 GN=TBC1D9B PE=1 SV=3
MWLSPEEVLVANALWVTERANPFFVLQRRRGHGRGGGLTGLLVGTLDVVLDSSARVAPYR
ILHQTQDSQVYWTVACGSSRKEITKHWEWLENNLLQTLSIFDSEEDITTFVKGKIHGIIA
EENKNLQPQGDEDPGKFKEAELKMRKQFGMPEGEKLVNYYSCSYWKGRVPRQGWLYLTVN
HLCFYSFLLGKEVSLVVQWVDITRLEKNATLLFPESIRVDTRDQELFFSMFLNIGETFKL
MEQLANLAMRQLLDSEGFLEDKALPRPIRPHRNISALKRDLDARAKNECYRATFRLPRDE
RLDGHTSCTLWTPFNKLHIPGQMFISNNYICFASKEEDACHLIIPLREVTIVEKADSSSV
LPSPLSISTKSKMTFLFANLKDRDFLVQRISDFLQKTPSKQPGSIGSRKASVVDPSTESS
PAPQEGSEQPASPASPLSSRQSFCAQEAPTASQGLLKLFQKNSPMEDLGAKGAKEKMKEE
SWHIHFFEYGRGVCMYRTAKTRALVLKGIPESLRGELWLLFSGAWNEMVTHPGYYAELVE
KSTGKYSLATEEIERDLHRSMPEHPAFQNELGIAALRRVLTAYAFRNPTIGYCQAMNIVT
SVLLLYGSEEEAFWLLVALCERMLPDYYNTRVVGALVDQGIFEELTRDFLPQLSEKMQDL
GVISSISLSWFLTLFLSVMPFESAVVIVDCFFYEGIKVILQVALAVLDANMEQLLGCSDE
GEAMTMLGRYLDNVVNKQSVSPPIPHLRALLSSSDDPPAEVDIFELLKVSYEKFSSLRAE
DIEQMRFKQRLKVIQSLEDTAKRSVVRAIPVDIGFSIEELEDLYMVFKAKHLASQYWGCS
RTMAGRRDPSLPYLEQYRIDASQFRELFASLTPWACGSHTPLLAGRMFRLLDENKDSLIN
FKEFVTGMSGMYHGDLTEKLKVLYKLHLPPALSPEEAESALEAAHYFTEDSSSEASPLAS
DLDLFLPWEAQEALPQEEQEGSGSEERGEEKGTSSPDYRHYLRMWAKEKEAQKETIKDLP
KMNQEQFIELCKTLYNMFSEDPMEQDLYHAIATVASLLLRIGEVGKKFSARTGRKPRDCA
TEEDEPPAPELHQDAARELQPPAAGDPQAKAGGDTHLGKAPQESQVVVEGGSGEGQGSPS
QLLSDDETKDDMSMSSYSVVSTGSLQCEDLADDTVLVGGEACSPTARIGGTVDTDWCISF
EQILASILTESVLVNFFEKRVDIGLKIKDQKKVERQFSTASDHEQPGVSG
>sp|Q66K14|TBC9B_HUMAN TBC1 domain family member 9B OS=Homo sapiens OX=9606 GN=TBC1D9B PE=1 SV=3
MWLSPEEVLVANALWVTERANPFFVLQRRRGHGRGGGLTGLLVGTLDVVLDSSARVAPYR
ILHQTQDSQVYWTVACGSSRKEITKHWEWLENNLLQTLSIFDSEEDITTFVKGKIHGIIA
EENKNLQPQGDEDPGKFKEAELKMRKQFGMPEGEKLVNYYSCSYWKGRVPRQGWLYLTVN
HLCFYSFLLGKEVSLVVQWVDITRLEKNATLLFPESIRVDTRDQELFFSMFLNIGETFKL
MEQLANLAMRQLLDSEGFLEDKALPRPIRPHRNISALKRDLDARAKNECYRATFRLPRDE
RLDGHTSCTLWTPFNKLHIPGQMFISNNYICFASKEEDACHLIIPLREVTIVEKADSSSV
LPSPLSISTKSKMTFLFANLKDRDFLVQRISDFLQKTPSKQPGSIGSRKASVVDPSTESS
PAPQEGSEQPASPASPLSSRQSFCAQEAPTASQGLLKLFQKNSPMEDLGAKGAKEKMKEE
SWHIHFFEYGRGVCMYRTAKTRALVLKGIPESLRGELWLLFSGAWNEMVTHPGYYAELVE
KSTGKYSLATEEIERDLHRSMPEHPAFQNELGIAALRRVLTAYAFRNPTIGYCQAMNIVT
SVLLLYGSEEEAFWLLVALCERMLPDYYNTRVVGALVDQGIFEELTRDFLPQLSEKMQDL
GVISSISLSWFLTLFLSVMPFESAVVIVDCFFYEGIKVILQVALAVLDANMEQLLGCSDE
GEAMTMLGRYLDNVVNKQSVSPPIPHLRALLSSSDDPPAEVDIFELLKVSYEKFSSLRAE
DIEQMRFKQRLKVIQSLEDTAKRSVVRAIPVDIGFSIEELEDLYMVFKAKHLASQYWGCS
RTMAGRRDPSLPYLEQYRIDASQFRELFASLTPWACGSHTPLLAGRMFRLLDENKDSLIN
FKEFVTGMSGMYHGDLTEKLKVLYKLHLPPALSPEEAESALEAAHYFTEDSSSEASPLAS
DLDLFLPWEAQEALPQEEQEGSGSEERGEEKGTSSPDYRHYLRMWAKEKEAQKETIKDLP
KMNQEQFIELCKTLYNMFSEDPMEQDLYHAIATVASLLLRIGEVGKKFSARTGRKPRDCA
TEEDEPPAPELHQDAARELQPPAAGDPQAKAGGDTHLGKAPQESQVVVEGGSGEGQGSPS
QLLSDDETKDDMSMSSYSVVSTGSLQCEDLADDTVLVGGEACSPTARIGGTVDTDWCISF
EQILASILTESVLVNFFEKRVDIGLKIKDQKKVERQFSTASDHEQPGVSG

fasta

>sp|Q664P8|TAUB_YERPS Taurine import ATP-binding protein TauB OS=Yersinia pseudotuberculosis serotype I (strain IP32953) OX=273123 GN=tauB PE=3 SV=1
MLNVSGLWAEYQGKPALQDVSLQIASGQLVVVLGPSGCGKTTLLNLIAGFMTPSAGVITL
DNIPVSGPSAERGVVFQNEGLLPWRDVVSNVEFGLQLAGMSKEQRRVTALKMLNRVGLAG
FEHHFIWQLSGGMRQRVGIARALAVDPRLLLLDEPFGALDAFTREQMQELLLTIWRDTGK
QILLITHDIEEAVFLASELLLLSPGPGQVVERLSLNFGQRYAEGEPCRAIKSDPEFIARR
EYVLGKVFQQREVLI
>sp|Q66K14|TBC9B_HUMAN TBC1 domain family member 9B OS=Homo sapiens OX=9606 GN=TBC1D9B PE=1 SV=3
MWLSPEEVLVANALWVTERANPFFVLQRRRGHGRGGGLTGLLVGTLDVVLDSSARVAPYR
ILHQTQDSQVYWTVACGSSRKEITKHWEWLENNLLQTLSIFDSEEDITTFVKGKIHGIIA
EENKNLQPQGDEDPGKFKEAELKMRKQFGMPEGEKLVNYYSCSYWKGRVPRQGWLYLTVN
HLCFYSFLLGKEVSLVVQWVDITRLEKNATLLFPESIRVDTRDQELFFSMFLNIGETFKL
MEQLANLAMRQLLDSEGFLEDKALPRPIRPHRNISALKRDLDARAKNECYRATFRLPRDE
RLDGHTSCTLWTPFNKLHIPGQMFISNNYICFASKEEDACHLIIPLREVTIVEKADSSSV
LPSPLSISTKSKMTFLFANLKDRDFLVQRISDFLQKTPSKQPGSIGSRKASVVDPSTESS
PAPQEGSEQPASPASPLSSRQSFCAQEAPTASQGLLKLFQKNSPMEDLGAKGAKEKMKEE
SWHIHFFEYGRGVCMYRTAKTRALVLKGIPESLRGELWLLFSGAWNEMVTHPGYYAELVE
KSTGKYSLATEEIERDLHRSMPEHPAFQNELGIAALRRVLTAYAFRNPTIGYCQAMNIVT
SVLLLYGSEEEAFWLLVALCERMLPDYYNTRVVGALVDQGIFEELTRDFLPQLSEKMQDL
GVISSISLSWFLTLFLSVMPFESAVVIVDCFFYEGIKVILQVALAVLDANMEQLLGCSDE
GEAMTMLGRYLDNVVNKQSVSPPIPHLRALLSSSDDPPAEVDIFELLKVSYEKFSSLRAE
DIEQMRFKQRLKVIQSLEDTAKRSVVRAIPVDIGFSIEELEDLYMVFKAKHLASQYWGCS
RTMAGRRDPSLPYLEQYRIDASQFRELFASLTPWACGSHTPLLAGRMFRLLDENKDSLIN
FKEFVTGMSGMYHGDLTEKLKVLYKLHLPPALSPEEAESALEAAHYFTEDSSSEASPLAS
DLDLFLPWEAQEALPQEEQEGSGSEERGEEKGTSSPDYRHYLRMWAKEKEAQKETIKDLP
KMNQEQFIELCKTLYNMFSEDPMEQDLYHAIATVASLLLRIGEVGKKFSARTGRKPRDCA
TEEDEPPAPELHQDAARELQPPAAGDPQAKAGGDTHLGKAPQESQVVVEGGSGEGQGSPS
QLLSDDETKDDMSMSSYSVVSTGSLQCEDLADDTVLVGGEACSPTARIGGTVDTDWCISF
EQILASILTESVLVNFFEKRVDIGLKIKDQKKVERQFSTASDHEQPGVSG
>sp|Q8K9I1|SYV_BUCAP Valine--tRNA ligase OS=Buchnera aphidicola subsp. Schizaphis graminum (strain Sg) OX=198804 GN=valS PE=3 SV=1
MKKNYNPKDIEEHLYNFWEKNGFFKPNNNLNKPAFCIMMPPPNITGNLHMGHAFQQTIMD
ILIRYNRMQGKNTLWQVGTDHAGIATQILIERQIFSEERKTKKDYSRNDFIKKIWKWKKK
SNFSVKKQMKRLGNSVDWDREKFTLDPDISNSVKEAFIILYKNNLIYQKKRLVHWDSKLE
TVISDLEVEHRLIKSKKWFIRYPIIKNIKNINIEYLLVATTRPETLLGDTALAINPKDDK
YNHLIGQSVICPIVNRIIPIIADHYADMNKDTGCVKITPGHDFNDYEVGQRHKLPMINIF
TFNGKIKSNFSIYDYQGSKSNFYDSSIPTEFQNLDILSARKKIIYEIEKLGLLEKIEECN
FFTPYSERSGVIIQPMLTNQWYLKTSHLSQSAIDVVREKKIKFIPNQYKSMYLSWMNNIE
DWCISRQLWWGHQIPVWYDDKKNIYVGHSEKKIREEYNISDDMILNQDNDVLDTWFSSGL
WTFSTLGWPEKTEFLKIFHSTDVLVSGFDIIFFWIARMIMLTMYLVKDSYGNPQIPFKDV
YITGLIRDEEGKKMSKSKGNVIDPIDMIDGISLNELIEKRTSNLLQPHLSQKIRYHTIKQ
FPNGISATGTDALRFTFSALASNTRDIQWDMNRLKGYRNFCNKLWNASRFVLKNTKDHDY
FNFSVNDNMLLINKWILIKFNNTVKSYRNSLDSYRFDIAANILYDFIWNVFCDWYLEFVK
SVIKSGSYQDIYFTKNVLIHVLELLLRLSHPIMPFITEAIWQRVKIIKHIKDRTIMLQSF
PEYNDQLFDKSTLSNINWIKKIIIFIRNTRSKMNISSTKLLSLFLKNINSEKKKVIQENK
FILKNIASLEKISILSKQDDEPCLSLKEIIDGVDILVPVLKAIDKEIELKRLNKEIEKIK
SKMLISEKKMSNQDFLSYAPKNIIDKEIKKLKSLNEIYLTLSQQLESLHDAFCKKNKIFN
>sp|Q664P8|TAUB_YERPS Taurine import ATP-binding protein TauB OS=Yersinia pseudotuberculosis serotype I (strain IP32953) OX=273123 GN=tauB PE=3 SV=1
MLNVSGLWAEYQGKPALQDVSLQIASGQLVVVLGPSGCGKTTLLNLIAGFMTPSAGVITL
DNIPVSGPSAERGVVFQNEGLLPWRDVVSNVEFGLQLAGMSKEQRRVTALKMLNRVGLAG
FEHHFIWQLSGGMRQRVGIARALAVDPRLLLLDEPFGALDAFTREQMQELLLTIWRDTGK
QILLITHDIEEAVFLASELLLLSPGPGQVVERLSLNFGQRYAEGEPCRAIKSDPEFIARR
EYVLGKVFQQREVLI
>sp|Q66K14|TBC9B_HUMAN TBC1 domain family member 9B OS=Homo sapiens OX=9606 GN=TBC1D9B PE=1 SV=3
MWLSPEEVLVANALWVTERANPFFVLQRRRGHGRGGGLTGLLVGTLDVVLDSSARVAPYR
ILHQTQDSQVYWTVACGSSRKEITKHWEWLENNLLQTLSIFDSEEDITTFVKGKIHGIIA
EENKNLQPQGDEDPGKFKEAELKMRKQFGMPEGEKLVNYYSCSYWKGRVPRQGWLYLTVN
HLCFYSFLLGKEVSLVVQWVDITRLEKNATLLFPESIRVDTRDQELFFSMFLNIGETFKL
MEQLANLAMRQLLDSEGFLEDKALPRPIRPHRNISALKRDLDARAKNECYRATFRLPRDE
RLDGHTSCTLWTPFNKLHIPGQMFISNNYICFASKEEDACHLIIPLREVTIVEKADSSSV
LPSPLSISTKSKMTFLFANLKDRDFLVQRISDFLQKTPSKQPGSIGSRKASVVDPSTESS
PAPQEGSEQPASPASPLSSRQSFCAQEAPTASQGLLKLFQKNSPMEDLGAKGAKEKMKEE
SWHIHFFEYGRGVCMYRTAKTRALVLKGIPESLRGELWLLFSGAWNEMVTHPGYYAELVE
KSTGKYSLATEEIERDLHRSMPEHPAFQNELGIAALRRVLTAYAFRNPTIGYCQAMNIVT
SVLLLYGSEEEAFWLLVALCERMLPDYYNTRVVGALVDQGIFEELTRDFLPQLSEKMQDL
GVISSISLSWFLTLFLSVMPFESAVVIVDCFFYEGIKVILQVALAVLDANMEQLLGCSDE
GEAMTMLGRYLDNVVNKQSVSPPIPHLRALLSSSDDPPAEVDIFELLKVSYEKFSSLRAE
DIEQMRFKQRLKVIQSLEDTAKRSVVRAIPVDIGFSIEELEDLYMVFKAKHLASQYWGCS
RTMAGRRDPSLPYLEQYRIDASQFRELFASLTPWACGSHTPLLAGRMFRLLDENKDSLIN
FKEFVTGMSGMYHGDLTEKLKVLYKLHLPPALSPEEAESALEAAHYFTEDSSSEASPLAS
DLDLFLPWEAQEALPQEEQEGSGSEERGEEKGTSSPDYRHYLRMWAKEKEAQKETIKDLP
KMNQEQFIELCKTLYNMFSEDPMEQDLYHAIATVASLLLRIGEVGKKFSARTGRKPRDCA
TEEDEPPAPELHQDAARELQPPAAGDPQAKAGGDTHLGKAPQESQVVVEGGSGEGQGSPS
QLLSDDETKDDMSMSSYSVVSTGSLQCEDLADDTVLVGGEACSPTARIGGTVDTDWCISF
EQILASILTESVLVNFFEKRVDIGLKIKDQKKVERQFSTASDHEQPGVSG
>sp|Q8E4B4|TARI_STRA3 Ribitol-5-phosphate cytidylyltransferase OS=Streptococcus agalactiae serotype III (strain NEM316) OX=211110 GN=tarI PE=3 SV=1
MNIGVIFAGGVGRRMNTKGKPKQFLEVHGKPIIVHTIDIFQNTEAIDAVVVVCVSDWLDY
MNNLVERFNLTKVKAVVAGGETGQMSIFKGLEAAEQLATDDAVVLIHDGVRPLINEEVIN
ANIKSVKETGSAVTSVRAKETVVLVNDSSKISEVVDRTRSFIAKAPQSFYLSDILSVERD
AISKGITDAIDSSTLMGMYNRELTIVEGPYENIKITTPDDFYMFKALYDARENEQIYGM
>sp|B3CQ06|SYS_WOLPP Serine--tRNA ligase OS=Wolbachia pipientis subsp. Culex pipiens (strain wPip) OX=570417 GN=serS PE=3 SV=1
MHDIEHIRKNPKGFEKAIKSRGVKEFTAKEILEIDHKKRSLTTKLQALNKQRNEVTEEIK
RLKMNKSPCEEQVKLSKSITSEIETISLKEQTEKNKLVDILSNLPNISAQNVPIGEDESS
NVEIRKYGKKRKFDFTPKFHYELGERLGLMDFEQAAKISGSRFTILKGQLAKLGRALINF
MLETHVNEFAYTEVYHPALVKNEAMYNVGQLPKFSDDSYLTTDKLRLIPTSEVVLTNLVA
DKIIEEKELPIRFTAYSECFRKEAGSAGRDTRGMIRQHQFGKVELVSITTEDQSKDELER
MTNAAEEILKKLELPYRIMLLCSGDMGFAAQKTYDIEVWLPEQNKYREISSCSNCGTFQA
RRMNTKYFLETDRKTKKYVHTLNGSALAIGRTIVAIMENYQNSDGSVTIPNVLQRYMSND
TVISKQ
>sp|Q9ATB4|TAD2B_ARATH Transcriptional adapter ADA2b OS=Arabidopsis thaliana OX=3702 GN=ADA2B PE=1 SV=1
MGRSRGNFQNFEDPTQRTRKKKNAANVENFESTSLVPGAEGGGKYNCDYCQKDITGKIRI
KCAVCPDFDLCIECMSVGAEITPHKCDHPYRVMGNLTFPLICPDWSADDEMLLLEGLEIY
GLGNWAEVAEHVGTKSKEQCLEHYRNIYLNSPFFPLPDMSHVAGKNRKELQAMAKGRIDD
KKAEQNMKEEYPFSPPKVKVEDTQKESFVDRSFGGKKPVSTSVNNSLVELSNYNQKREEF
DPEYDNDAEQLLAEMEFKENDTPEEHELKLRVLRIYSKRLDERKRRKEFIIERNLLYPNP
FEKDLSQEEKVQCRRLDVFMRFHSKEEHDELLRNVVSEYRMVKRLKDLKEAQVAGCRSTA
EAERYLGRKRKRENEEGMNRGKESGQFGQIAGEMGSRPPVQASSSYVNDLDLIGFTESQL
LSESEKRLCSEVKLVPPVYLQMQQVMSHEIFKGNVTKKSDAYSLFKIDPTKVDRVYDMLV
KKGIAQL
>sp|Q83JA5|SYW_SHIFL Tryptophan--tRNA ligase OS=Shigella flexneri OX=623 GN=trpS PE=3 SV=1
MTKPIVFSGAQPSGELTIGNYMGALRQWVNMQDDYHCIYCIVDQHAITVRQDAQKLRKAT
LDTLALYLACGIDPEKSTIFVQSHVPEHAQLGWALNCYTYFGELSRMTQFKDKSARYAEN
INAGLFGYPVLMAADILLYQTNLVPVGEDQKQHLELSRDIAQRFNALYGEIFKVPEPFIP
KSGARVMSLLEPTKKMSKSDDNRNNVIGLLEDPKSVVKKIKRAVTDSDEPPVVRYDVQNK
AGVSNLLDILSAVTGQSIPELEKQFEGKMYGHLKGEVADAVSGMLTELQERYHRFRNDEA
FLQQVMKDGAEKASVHASRTLKAVYEAIGFVAKP
>sp|P14213|TAC1_TACTR Tachyplesin-1 OS=Tachypleus tridentatus OX=6853 PE=1 SV=2
MKKLVIALCLMMVLAVMVEEAEAKWCFRVCYRGICYRRCRGKRNEVRQYRDRGYDVRAIP
EETFFTRQDEDEDDDEE
>sp|Q7SZM9|TB1RA_XENLA F-box-like/WD repeat-containing protein TBL1XR1-A OS=Xenopus laevis OX=8355 GN=tbl1xr1-a PE=1 SV=1
MSISSDEVNFLVYRYLQESGFSHSAFTFGIESHISQSNINGALAPPAALISIIQKGLQYV
EAEVSINEDGTLFDGRPIESLSLIDAVMPDVVQTRQQAYRDKLAQQQTAAAAAAAAAAAA
TPNNQQPPAKNGENTANGEENGGHALANNHTDMMEVDGDVEIPSSKAVVLRGHESEVFIC
AWNPVSDLLASGSGDSTARIWNLSENSTSGSTQLVLRHCIREGGQDVPSNKDVTSLDWNS
EGTLLATGSYDGFARIWTKDGNLASTLGQHKGPIFALKWNKKGNFILSAGVDKTTIIWDA
HTGEAKQQFPFHSAPALDVDWQSNNTFASCSTDMCIHVCKLGQDRPIKTFQGHTNEVNAI
KWDPTGNLLASCSDDMTLKIWSMKHDTCVHDLQAHNKEIYTIKWSPTGPGTNNPNANLML
ASASFDSTVRLWDVDRGICIHTLTKHQEPVYSVAFSPDGRYLASGSFDKCVHIWNTQTGA
LVHSYRGTGGIFEVCWNAAGDKVGASASDGSVCVLDLRK
>sp|Q9FGE9|TBL12_ARATH Protein trichome birefringence-like 12 OS=Arabidopsis thaliana OX=3702 GN=TBL12 PE=2 SV=1
MELGSRRIYTTMPSKLRSSSSLLPRILLLSLLLLLFYSLILRRPITSNIASPPPCDLFSG
RWVFNPETPKPLYDETCPFHRNAWNCLRNKRDNMDVINSWRWEPNGCGLSRIDPTRFLGM
MRNKNVGFVGDSLNENFLVSFLCILRVADPSAIKWKKKKAWRGAYFPKFNVTVAYHRAVL
LAKYQWQARSSAEANQDGVKGTYRVDVDVPANEWINVTSFYDVLIFNSGHWWGYDKFPKE
TPLVFYRKGKPINPPLDILPGFELVLQNMVSYIQREVPAKTLKFWRLQSPRHFYGGDWNQ
NGSCLLDKPLEENQLDLWFDPRNNGVNKEARKINQIIKNELQTTKIKLLDLTHLSEFRAD
AHPAIWLGKQDAVAIWGQDCMHWCLPGVPDTWVDILAELILTNLKTE
lengedd
关注
  • 0
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
extract_fasta_序列提取_
10-03
给定物种名，批量分子矩阵中提取给定的物种名称对应的序列，并生成新矩阵用于系统发育分析。
python 行和列_Python：区分行和列向量
weixin_35950531的博客
01-15 1644
您可以通过向数组添加另一个维度来明确区分。>>> a = np.array([1, 2, 3])>>> aarray([1, 2, 3])>>> a.transpose()array([1, 2, 3])>>> a.dot(a.transpose())14现在强制它成为一个列向量：>>> a.shape = ...
一个处理fasta标题的小程序
03-04 482
问题：有一个fasta文件，其中>后面只有名字，现在要依据另外一个txt文件将附属说明添加到对应名字中。 解决：用perl很方便的可以解决。 思路： 1、先打开txt文件，读取关键字和附加说明。 2、再打开fasta文件，将关键字和>后面的名字对比。如果相同就将附加说明添加上去。并将其后面的序列重新写入新的文件中。 这一块用了flag标识。0代表未匹配状态，1代表已经匹配，2代表链接完后
用python提取fasta序列,如何使用Python随机提取FASTA序列？
weixin_39716088的博客
12-16 404
I have the following sequences which is in a fasta format with sequence header and its nucleotides. How can I randomly extract the sequences. For example I would like to randomly select 2 sequences ou...
python实现根据序列ID从提取fasta文件序列
little^raccoon
11-21 8436
当序列少的时候，我习惯用 grep -A 1 -f seq.lst seq.fas | sed ‘/^–$/d’ > out.fas提取，但是这次遇到了一个大文件，用grep就太费时了，然后又试了一下TBtools的提取序列功能，发现时间也很长，所以就写了个python。提取将近100万条reads耗时也就需要10s左右 #!/usr/bin/python # -*- coding: u...
DNA-FASTA-Python:使用Python分析Multi-Fasta格式的DNA序列
05-04
在第一列中，标题行与序列数据之间用大于号（“>”）隔开。 “>”符号后的单词是序列的标识符，该行的其余部分是该条目的可选描述。 “>”和标识符的首字母之间不应有空格。 文件中序列的长度是多少？ 最长的序列...
Remove-duplicate-fasta:Python脚本删除重复的Fasta序列
05-13
-Python script to remove whole duplicate fasta sequences i.e identical sequence and header -input file must be in fasta format usage: python remove_duplicate_fasta.py inputfile outputfile 例子： ...
氨基酸符号序列转换为FASTA格式的蛋白质序列
05-18
氨基酸符号序列转换为FASTA格式的蛋白质序列，直接编译，将英文氨基酸序列粘贴到窗口，按回车即可输出转换结果。
分割fasta文件的python脚本
08-21
文件脚本可将大的fasta文件中的序列，按照个数均分，分割成多个fasta文件，便于对各个小文件中的序列进行后续操作
python绘制直方图
热门推荐
qq_39212193的博客
07-01 1万+
用两列数据绘制直方图 #coding=gbk import xlwings as xw import pandas as pd import matplotlib.pyplot as plt #pd.set_option('display.max_columns', None) #解决表格多列时中间省略显示问题 #pd.set_option('display.max_rows', None) #解决表格多行时中间省略显示问题 #读取excel文件中的数据 app = xw.App(visi
python绘制热图
qq_39212193的博客
07-01 1万+
# coding=gbk import matplotlib.pyplot as plt import pandas as pd import seaborn as sns import numpy as np #结果图编码从0开始，编号依次加一对应蛋白的编号,其中50为m7 file_name = pd.read_excel('lactose-lactuse.xlsx',engine='openpyxl') data = file_name.pivot("lactose","lactulose","l
python爬虫报错requests.exceptions.MissingSchema: Invalid URL
qq_39212193的博客
07-09 6806
把开头的 #encoding=gbk改为utf-8
用seaborn作热图报错：urllib.error.URLError: ＜urlopen error [Errno 11004] getaddrinfo failed＞
qq_39212193的博客
02-04 3211
seaborn-data数据下载网址 https://codechina.csdn.net/mirrors/mwaskom/seaborn-data?utm_source=csdn_github_accelerator import matplotlib.pyplot as plt import seaborn as sns sns.set() #用行和列标签绘制 flights_long = sns.load_dataset("flights") flights = flights_long.pivot(
计算标准差和中位数
qq_39212193的博客
04-16 2724
计算标准差 计算标准差：要先算出平均值，—>（value-aver）**2,相加后除以长度，再开平方跟 import math data = [3.2, 4.5, 2.4, 9.2, 8.3] aver = sum(data) / len(data) total = 0.0 for value in data: total +=(value - aver) ** 2 pfc = math.sqrt(total/len(data)) print(pfc) 计算中位数 data =
用seaborn作聚类热图
qq_39212193的博客
02-04 2167
import matplotlib.pyplot as plt import seaborn as sns import pandas as pd sns.set() data = pd.read_csv("fs.csv") #print(data.head()) flights = data.pivot("month","year","passengers") ## 绘制x-y-z的热力图，比如 年-月-销量 的聚类热图 g = sns.clustermap(flights, fmt='d', c
python通过txt文件将指定文件复制到另一个文件夹
qq_39212193的博客
07-22 1999
import shutil import numpy as np data = [] for line in open("test5.txt", "r"): # 设置文件对象并读取每一行文件 data.append(line) # print(data) for a in data: srcfile_path = 'pdbname_list/{}'.format(a[:-1]) print(srcfile_path) tarfile_path = 'targt1'
python绘制散点聚类图及标记序号
qq_39212193的博客
07-01 1945
# coding=gbk import pandas as pd import matplotlib.pyplot as plt from sklearn.cluster import KMeans # data = pd.read_excel('lactose-lactuse-scatter.xlsx', engine='openpyxl') # #遍历所有数据 # x = data.iloc[:] # # #5个簇，1个自由点 # mod = KMeans(n_clusters=5, random_s
python实现pdb转fasta
qq_39212193的博客
07-16 1683
import sys if len(sys.argv) <= 1: print('python pdb2fasta.py desgin_3.txt > desgin_3.fasta') exit() input_file = open(sys.argv[1]) letters = {'ALA': 'A', 'ARG': 'R', 'ASN': 'N', 'ASP': 'D', 'CYS': 'C', 'GLU': 'E', 'GLN': 'Q', 'GLY': 'G', '
用awk提取fasta文件的序列ID和该序列中的N
最新发布
05-10
不好意思，我作为语言模型，无法与您进行实时的计算和代码操作。但是，我可以为您提供相关的代码示例来...这个awk命令会从fasta文件中过滤掉序列ID行，并且计算所有剩余行中"N"的数量。 希望这些示例对您有所帮助！

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值