python中 pysam包FastxFile函数

2022-08-20 01:32:39 阅读：126 来源： 互联网

标签：... pysam python ## Rosalind FastxFile file fasta

001、读取fasta文件

root@PC1:/home/test# ls
a.fasta
root@PC1:/home/test# cat a.fasta                         ## 测试数据
>Rosalind_1
ATCCAGCT
>Rosalind_2
GGGCAACT
>Rosalind_3
ATGGATCT
>Rosalind_4
AAGCAACC
>Rosalind_5
TTGGAACT
>Rosalind_6
ATGCCATT
>Rosalind_7
ATGGCACT
root@PC1:/home/test# python                                       ## 启动python
Python 3.10.4 (main, Mar 31 2022, 08:41:55) [GCC 7.5.0] on linux
Type "help", "copyright", "credits" or "license" for more information.
>>> import pysam                                                  ## 导入pysam包
>>> in_file = pysam.FastxFile("a.fasta")                          ## 读取数据，生成可迭代对象
>>> for i in in_file:
...     print(i)
...
>Rosalind_1
ATCCAGCT
>Rosalind_2
GGGCAACT
>Rosalind_3
ATGGATCT
>Rosalind_4
AAGCAACC
>Rosalind_5
TTGGAACT
>Rosalind_6
ATGCCATT
>Rosalind_7
ATGGCACT

>>> in_file = pysam.FastxFile("a.fasta")
>>> for i in in_file:
...     print(i.name)                     ## 输出scaffold
...
Rosalind_1
Rosalind_2
Rosalind_3
Rosalind_4
Rosalind_5
Rosalind_6
Rosalind_7
>>> in_file = pysam.FastxFile("a.fasta")
>>> for i in in_file:
...     print(i.sequence)                ## 输出序列
...
ATCCAGCT
GGGCAACT
ATGGATCT
AAGCAACC
TTGGAACT
ATGCCATT
ATGGCACT

002、

(base) root@PC1:/home/test# ls
a.fasta  b.fasta
(base) root@PC1:/home/test# cat b.fasta                        ## 测试fasta文件
>Rosalind_1
ATCCAGCT
ATCCAGCT
TTTTT
>Rosalind_2
GGGCAACT
GGGCAACT
GGGCAACT
>Rosalind_3
ATGGATCT
ATGGATCT
ATGGATCT
(base) root@PC1:/home/test# python                                   ## 启动python
Python 3.10.4 (main, Mar 31 2022, 08:41:55) [GCC 7.5.0] on linux
Type "help", "copyright", "credits" or "license" for more information.
>>> import pysam                                                     ## 导入pysam包                          
>>> in_file = pysam.FastxFile("b.fasta")
>>> for i in in_file:
...     print(i)                                               ## 将多行序列合并为一行
...
>Rosalind_1
ATCCAGCTATCCAGCTTTTTT
>Rosalind_2
GGGCAACTGGGCAACTGGGCAACT
>Rosalind_3
ATGGATCTATGGATCTATGGATCT

参考：https://mp.weixin.qq.com/s?__biz=MzIxMjQxMDYxNA==&mid=2247484172&idx=1&sn=d8dec9ae5ffea81ef02e8f0d7ea4672b&chksm=9747ca95a030438313f483f6c62c9c32551e23682f98be6868edf423ea88180165e21c5dedc8&scene=178&cur_album_id=1635727573621997580#rd

标签：...,pysam,python,##,Rosalind,FastxFile,file,fasta
来源： https://www.cnblogs.com/liujiaxin2018/p/16607004.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

ICode9

python中 pysam包FastxFile函数