比较两个fasta序列的碱基个数并提取ATG及之后的序列

请求大佬帮忙写个脚本：

比较两个fasta序列的碱基个数，
1、如果fasta1<fasta2，则输出fasta1第一个ATG及之后的碱基序列。同时输出fasta2的第一个ATG及之后的碱基序列。
2、如果fasta1>fasta2，则判断fasta1-fasta2的这段序列中是否存在终止密码子（TAG,TAA,TGA）
如果存在，则输出fasta1文件终止密码子之后的第一个ATG及之后的碱基，同时输出fasta2文件第一个ATG及之后的碱基序列
如果不存在，则输出fasta1第一个ATG及之后的碱基序列，同时输出fasta2文件第一个ATG及之后的碱基序列
3、输出处理后的fasta1、fasta2文件中每条序列的碱基个数，另存为一个文件
4、fasta1、fasta2文件分别有两个，为±方向的序列，＋序列做以上处理，-序列需要先转为+序列然后做以上处理

fasta1+文件示例：

>PGSC0003DMT400004961::chr08:47200220-47203045
GGAATTAGAATTAAAAGAAATTCTTCCTCACTAAAAGTTCTCCATTATGGGCCAAAAAGCTCTTCCTCACTAAAATTATTTATGCTTTTGTAGAGGTAACTCCCACTTGTCCGGTGTGACGAACTCGACCCAGAAGGCAAATCCCTTGCCAACTTGAGACATCCGATTTGAAAATCCCAAGCCCAAACGACTCAGCTTCACAATTTATTAACTTCACGCAGAAACAAACTGATAACAAGACTTGAAAAAAGGTCTGGAGCTTTTTCAATCATGGAGTTAGTACTGCAATACTGTACAAAAAAGCTAGGAAATGAAGCGATTGAGCGCGAATACAGAAGAAGATGGGAAGAAATAAGGAAAGCTTATATGCCCCTATGGAGCAATTCAACTGCAGAGTCTCACTCAAATTAATATTTGAAAATGCAATATAACAAGACAACATCTTACAGGCCTATAGGCCATAATAACATGTAGTATGTCTTATTTAAAGAGAAGAACATTAGCACACTGTAAAGCCAGCTGAAAGTCTTCGATACTATAACTAGTTTACTTCCACTTCTTGAATTTACCTCCGAACATCCCTCCCTTGCCATGCTTGCCATGCTTTCCATGTTTGAATTTCCCTCCCTTGCCATGCTTAAACTTGCCTCCTCCATAATGACCACCTCCATAATGACTGCCGCCTCCAACAAATGGCACACCATGTCCTCCCATGTGAGATCCACCGTGTGACAGCTTATGAGCGCCATACGCAGCAGCTGCAGCTGCAGCACCTCCAGCTAACATTACCCCCATGCCACCATGTCCTCCTGAAAAGTAACAAGCAAACAAACACGAAAAATTAATAATGAGGAAATTGTCAAGATGTACTTGGTCACAAAAGTCGCAACATATTATAAATGTTCATTATTTTCCTGCATCTAGTTGTTCAAGTGATACCATTGCTATCAAAGAAACTTAACTAGACAAATCCAAAAAGGGTGCAACCAATTAGACAGTATCTTATTGGGATGTCCGGGGGTTGTCGACAGCTTATAAATTACTAAATGAAGTCTTTACAAGTTACTTCCAAGAAACTTGGGAAGAAATAATACCTGACTGGTGTGGAGCACCATGATGACCTGCAGGAGGATATCCCTGTGGAGGATATCCTTGTTGTGGAGGATATCCAGCAGGTGGATATCCTCCTTGCGGAGGATAACCCTGTTGTGGATACCCTTGCTGTGGAGGGTACCCTTGTTGTGGAGGGTACCCTTGAGGAGGATACCCTTGCTGCGGAGGGTAACCACCGTGCTGTGGAGGATACTGTCCTGGATGTCCGTGATGACCACCAGCAAGACCACCCATACCTCCCATTATATTGTTAAAGAGTCCCTTATCAGAATCATCCTTACCACCTCCCATATTTCACTAGTACCTAAAACAGATAAATACTTGTCAGACATATTACAACAATGACTAAGTAGGAATGAAATCAAATGAAATTCGTAGATGATAAAAGATTATGAGCCACAAGATATCAATTACTACTCTTAACTTCAAATATGCTTTCTGAAACATCCAATCAATCAAAGAAGAATTATGTGCACAACTAATGATCAAGCAACCACTGATTGAACTTGTAGTACATCATATATAGATAGTGCACATCAAGAAACAAGAAAGGTGGACACCAACATAAAGCTGATTGTTAAGTAGGAGAAGTGTTAAATTTTAACACAACAAGCTGCAGACAGAATTGCAGGAAGCCAGCTAATATGTTCTGAGCAAGTGAAAATGCAAAAAAACAGGAGGGTAGACATAAACATTTTTATAAATTTCTATGCAGTACCCTAGCCTCCATTCCTGAATAAGTTTTCAGATCACCCTAAGAAAAAAAAAAAACCTTCACCTCTCAAAGATTAGAAGCATCAAAAGTTCTTTTTTTGTTTTGGATAAGTTGGAAACAACTTGCATTTTGAGGAAAAAGGAGAGGCCTTTTTTAGAGGGGTGGATAATAAAGAAGGAAAAGTGACTAAATCTAATGAAACAAAAGGAGCTAAAACCAAATATCCATTCTTTTTTCATGATTTTTTTCTTTTTTTGGAATCTAAGTTGCATGACTATTTGTTTATAATACACCTTGAAAAAAGACAAAACTCAATTAAAACAGATCTCTATTCCCTTCTACAGCAAAATCATGCAAGAGGACAAAATCCTCATTTTCTAAAGCATCCAAATTCAGTTCCAGAATAGGACAATTCTCAACCCAAAAACAGAAAAAAGGGAACTAATAGTCTCTCTTTCCAATTCAACAAAACCAAAATTTAAAACCCAATGACTCATCAGCTAAACAGAAAATACAGAGAAAACACAAAATGATTGATCATACAAAAAAAAATCAAAATGGTAAAAACCAAATCTTCTGTTTTTCCATTTAATCAAAACCCAAATCAAAAAAAGTCACAAGGGAAGAAACAATCAAACTTGAATTTCCAAAATCAGTTTTTTTTTTTAAAAAAAAAGCTGATCTAATCTGTAGTAGGTGGTCATAGAAACAAAAATCAAAAGGGTAGACCCAAAATCCTTTAATTTTCCCATTAACAAAAACCCAAATCTAAAAAACATCAAAAGGGTACAAAAAGATCAATTTTTTAATTACAAAAGACCAGATCTTGAATTATCAGACTGATGATCATCATCTGGGCAATTCAGATAATTCATGCATGGCTTAACAAATCATCAAACAAGAGAAAAGTATCATCAATGGGGTAAAAGGAGAAGGGGCTCTTACAGATTTATCTCCTTGAGAAGTTGAACGAACGTTTGTTTAACAAACGATATGG
>PGSC0003DMT400073331::chr06:49955574-49958868
AAACATGCAGTGCCATTTAAAACCTTCTCTTGTATCAGCTGATAGAACTTTCTCCATTGGTTTATATTGTGGTAAAGTATCCCCTGAATTTGAGGTTTTAACAAGATGGTTCACGTGGATGCTGGCAAGTTCTTTTAGTGTAACAGGAGCTCCAAAGCTGTAAGTTGTAACTCATGGATGTTTGATTCCCTTCAGTGTGGTTTATTCCGAGATAACTGAAGCGATGCAATTTTCTATGATTGATGGAAGTGGAAGTTCAACAAATGTGGGAGGCTTAGCTACCTGGTGTAGGGGTGTGATACCATGTTTAAGTCAGCTAAGCTACTATAATTCCACATTTAATGCATTCATCCTGAAGAGGTCAGTGTACGTTATTTTCCGTATGCTCTCTTTTACATCTTAACTTTGTGTGGTAGTGCCATAGCTCCCTATTCATCAAGCTTTGCACATTCTACCCTTATATATGTAATGAATTTTTTGAGTAATTCTAGTCATTGCTGAATAACTTATCGCCTTATCTAGGACTACCTCCTTTTTTACCCAAAAAAATTAAATTGAGTGCTATTCATTTGTTTTTGTTGCTTTTTCCGGGTTATGTATCCTATCCTTGTAAATTATTCAAAACTTGTGATGTGAAGGAGAAGGGTAGTCATTTATCTTGGATTCCAGTCATTGTTTTTGGAAAAAATTGAAAATATAGGATAAATGGTATGAATGGAAATATCAACAGCACGAAAAGTGTGGTTGGAACTTTTGAATACATTGATTGATGATATCTTTTGGTCAATTATAGTCACCTCTCAGGTATAGTATTGGCCTTTGGGACTACGAATTGGTTCACTAGAATAAAGATTGGGGACAAATGTGACGTGTCTGTATATATCAAAAGAAAAAAAGTGATCTGATTTGTCCAAAGCTGTTGAAAGGCAAATAATAGCTCTTTGAACCCAAAAAAAAAACACTGTAAAAGGAAAAGCACATAACCTTCACAGGCTTTCAATATCATACTATAAAGCTTTGTTTTTCTCTTCACCCTTTTTCAGAAGGATAAACAGTGACTTCCTCATCAGTGTGACCACATTGTGAAAGCTCCCCTTCCTTTCAAGTGTTGTCCTGGTTGGAACAACTCTGTATTCTCCCCAACATAACCCATAATCAAGAGAATTGTTGAACTAAGTTGTGGGAACTCTTCACTTTCGATGCCGCACCCATGTCGGATTCTCCAAAAATAGTGTATTTTTGGGGAATATGACATGCACCCAATTGACTTTTTTGAAGAATCCGAGCAACATAGTTGTTGAGTAGTTAATCAAACAATAATGTCCTAATTCTATATTAGTGTGATTTGCAGTATATTTCTCTGTATCCATTCCGCTCAGTTCAGGTTCAATATCACGGTTCAACTAATTTTACGATGATTGTAATAAGCCTATTAATAGATAAAAGATTAACTTCCATTTGCTTGTTGCTGCATATTTTAGATTTTCATTTCCACTGAAAGCATTTGCATATACAATATTTGTTAGCATTATGAAGAAACGAGCAAAAGAATAAAGGTGAATCAAATAAAGCATTGGCCTTTCTGCATAAGGAAATTCTAGGATGCATTCGAATCTGACGTCAACGGCGTAGATGGAGGTAGTGGGCTTCTGGTTTCGTTACTACTTTTGGATTTGGCAATTGGAGGAGACTGCAAAGAAACTGACGATGGATTTGAAGGAGTCGTCTTATTCGAGCTACTCTGAGGTTTCTTGTTGTTTGATGGATTTGAAGATGTAGTCGGCTTATTCGAGCTGCTCTGCGGTTTCTTCTTGTCCCCTTCACCCGATGGATTCGAAGACGGAGTCTGGTTATTCAAGCTGCTCTGAGGTTTCTTCTTGTCCCCTTCACCTGATGGATTCGAAGACGGAGTCTGGTTATTCAAGCTGCTCTGAGGTTTCTTGCTGTTCCCTTCACTTGCCTTGTCCTCCATCCGATGGTGGCGAACTTTGGTAGACTTTAAGCTCATTGGAACAATACAATGGCAAACGAAACCTACACACACCAACAATTACCAAGTTAATAAATTTCTCATATCAACAAGTTGCAGAACACTAACTATTCAAACTTGAGAACTTTATAACCCCTTCAATGACTGCAACAAGCTCATGTTGAAGATTCCTATCGTGTTAAAGGCTACGAAAAACTTTCAACGCAGCCACTGAACTCGGTAGGAATCTTCTGCAATGACTGCAAAAAATCATTGCCAAGAACAAACATTGAGACTTCTATGAAACTACCCAGAAGCTTCCTGATTCAACATGTCTAGAAAACAAGGTACTCTCTCCGTCCCAATTTATGTAACACTTAACTATTCAAAGTGTTTGAATTTTGACCAAACATTTTAAAATGTATTTTCCGAATATTCAAGCATGTGAAACAATAAATCTGACTGAAAATCTTACCAATTCTTGCTAGTCGATTAACCCAACTGGGAATTGGATTGTTAGTGAGTTTAACACAAGCATCACTACAGAAATGGTTACAATTCTTAGTAATAAGATTGTAAGCATTTCCTTTATACTTATCAGCAAGTTCTTCCATAACTCTTCTTACTTCTCTAGGTGTCTTCTCCGTCCATCCGATCAAAATTGTCTTCCTGAAAATAAATCCTTCACACTTTTTCGGCTCCCCTTCGAAAATTCCAGTACTTGGATACTCATGAGCTCCAAATGCATACTCAACTCCATGAACTAGCTACCATAACCCACCACCATCAACAAATAATTAGTATAAACTTTGCTACAGTTAATTAAGAAAAAACCAGATCTGAATAACAATGAAAAGTTTGAAGATTGAGTAACTGACCTTGAACGCCGGAATGGTAAACTCCGAGACCCAGCCAATAAGCATAACCATTCATCGGAGTTAGATCGTATACGTTTAAGTACACCGGCACCGATCCACGGCCATTATTAACTGATTTTTTGCAGCAGAGCATCTTGAAATAAAAACAGAGGAACTTTTCAGGTTCCTCGCATTGTGAAGATGTAAAGAAAGAGAATTCTGAATTTTCCGATTGATTGAGTTGCATTGAATTTCTGGTATATGTGGGCAGTACTTGGTAGGTAAATGAAAGAGGAAATGTGGGTCCTCTTATTTTGGAGTTGGCGGTGGAGTGGGGACTGAGGGTACCACCGGCGGCGATACAGTGGCCGTACCGGTAGGTATAAAAGGAGGAGAAGAGATTAGTTGGAAAATTGGAATCTTTCAGCCGAAGAGTGTATTATAATGGTGGCTGAGAGTGTATTGTGG

fasta2+文件示例：

>PGSC0003DMT400004961::chr08:47200220-47203144
GGAATTAGAATTAAAAGAAATTCTTCCTCACTAAAAGTTCTCCATTATGGGCCAAAAAGCTCTTCCTCACTAAAATTATTTATGCTTTTGTAGAGGTAACTCCCACTTGTCCGGTGTGACGAACTCGACCCAGAAGGCAAATCCCTTGCCAACTTGAGACATCCGATTTGAAAATCCCAAGCCCAAACGACTCAGCTTCACAATTTATTAACTTCACGCAGAAACAAACTGATAACAAGACTTGAAAAAAGGTCTGGAGCTTTTTCAATCATGGAGTTAGTACTGCAATACTGTACAAAAAAGCTAGGAAATGAAGCGATTGAGCGCGAATACAGAAGAAGATGGGAAGAAATAAGGAAAGCTTATATGCCCCTATGGAGCAATTCAACTGCAGAGTCTCACTCAAATTAATATTTGAAAATGCAATATAACAAGACAACATCTTACAGGCCTATAGGCCATAATAACATGTAGTATGTCTTATTTAAAGAGAAGAACATTAGCACACTGTAAAGCCAGCTGAAAGTCTTCGATACTATAACTAGTTTACTTCCACTTCTTGAATTTACCTCCGAACATCCCTCCCTTGCCATGCTTGCCATGCTTTCCATGTTTGAATTTCCCTCCCTTGCCATGCTTAAACTTGCCTCCTCCATAATGACCACCTCCATAATGACTGCCGCCTCCAACAAATGGCACACCATGTCCTCCCATGTGAGATCCACCGTGTGACAGCTTATGAGCGCCATACGCAGCAGCTGCAGCTGCAGCACCTCCAGCTAACATTACCCCCATGCCACCATGTCCTCCTGAAAAGTAACAAGCAAACAAACACGAAAAATTAATAATGAGGAAATTGTCAAGATGTACTTGGTCACAAAAGTCGCAACATATTATAAATGTTCATTATTTTCCTGCATCTAGTTGTTCAAGTGATACCATTGCTATCAAAGAAACTTAACTAGACAAATCCAAAAAGGGTGCAACCAATTAGACAGTATCTTATTGGGATGTCCGGGGGTTGTCGACAGCTTATAAATTACTAAATGAAGTCTTTACAAGTTACTTCCAAGAAACTTGGGAAGAAATAATACCTGACTGGTGTGGAGCACCATGATGACCTGCAGGAGGATATCCCTGTGGAGGATATCCTTGTTGTGGAGGATATCCAGCAGGTGGATATCCTCCTTGCGGAGGATAACCCTGTTGTGGATACCCTTGCTGTGGAGGGTACCCTTGTTGTGGAGGGTACCCTTGAGGAGGATACCCTTGCTGCGGAGGGTAACCACCGTGCTGTGGAGGATACTGTCCTGGATGTCCGTGATGACCACCAGCAAGACCACCCATACCTCCCATTATATTGTTAAAGAGTCCCTTATCAGAATCATCCTTACCACCTCCCATATTTCACTAGTACCTAAAACAGATAAATACTTGTCAGACATATTACAACAATGACTAAGTAGGAATGAAATCAAATGAAATTCGTAGATGATAAAAGATTATGAGCCACAAGATATCAATTACTACTCTTAACTTCAAATATGCTTTCTGAAACATCCAATCAATCAAAGAAGAATTATGTGCACAACTAATGATCAAGCAACCACTGATTGAACTTGTAGTACATCATATATAGATAGTGCACATCAAGAAACAAGAAAGGTGGACACCAACATAAAGCTGATTGTTAAGTAGGAGAAGTGTTAAATTTTAACACAACAAGCTGCAGACAGAATTGCAGGAAGCCAGCTAATATGTTCTGAGCAAGTGAAAATGCAAAAAAACAGGAGGGTAGACATAAACATTTTTATAAATTTCTATGCAGTACCCTAGCCTCCATTCCTGAATAAGTTTTCAGATCACCCTAAGAAAAAAAAAAAACCTTCACCTCTCAAAGATTAGAAGCATCAAAAGTTCTTTTTTTGTTTTGGATAAGTTGGAAACAACTTGCATTTTGAGGAAAAAGGAGAGGCCTTTTTTAGAGGGGTGGATAATAAAGAAGGAAAAGTGACTAAATCTAATGAAACAAAAGGAGCTAAAACCAAATATCCATTCTTTTTTCATGATTTTTTTCTTTTTTTGGAATCTAAGTTGCATGACTATTTGTTTATAATACACCTTGAAAAAAGACAAAACTCAATTAAAACAGATCTCTATTCCCTTCTACAGCAAAATCATGCAAGAGGACAAAATCCTCATTTTCTAAAGCATCCAAATTCAGTTCCAGAATAGGACAATTCTCAACCCAAAAACAGAAAAAAGGGAACTAATAGTCTCTCTTTCCAATTCAACAAAACCAAAATTTAAAACCCAATGACTCATCAGCTAAACAGAAAATACAGAGAAAACACAAAATGATTGATCATACAAAAAAAAATCAAAATGGTAAAAACCAAATCTTCTGTTTTTCCATTTAATCAAAACCCAAATCAAAAAAAGTCACAAGGGAAGAAACAATCAAACTTGAATTTCCAAAATCAGTTTTTTTTTTTAAAAAAAAAGCTGATCTAATCTGTAGTAGGTGGTCATAGAAACAAAAATCAAAAGGGTAGACCCAAAATCCTTTAATTTTCCCATTAACAAAAACCCAAATCTAAAAAACATCAAAAGGGTACAAAAAGATCAATTTTTTAATTACAAAAGACCAGATCTTGAATTATCAGACTGATGATCATCATCTGGGCAATTCAGATAATTCATGCATGGCTTAACAAATCATCAAACAAGAGAAAAGTATCATCAATGGGGTAAAAGGAGAAGGGGCTCTTACAGATTTATCTCCTTGAGAAGTTGAACGAACGTTTGTTTAACAAACGATATGGTCTCTGAATGGTATTTGAAGGGCTATTTGTGAATGAGATATTTCAATGATTCCACGCGTTACTTAAACGCGTACATGACGTGTGCTACTTTCAATAAGG
>PGSC0003DMT400073331::chr06:49955574-49958966
AAACATGCAGTGCCATTTAAAACCTTCTCTTGTATCAGCTGATAGAACTTTCTCCATTGGTTTATATTGTGGTAAAGTATCCCCTGAATTTGAGGTTTTAACAAGATGGTTCACGTGGATGCTGGCAAGTTCTTTTAGTGTAACAGGAGCTCCAAAGCTGTAAGTTGTAACTCATGGATGTTTGATTCCCTTCAGTGTGGTTTATTCCGAGATAACTGAAGCGATGCAATTTTCTATGATTGATGGAAGTGGAAGTTCAACAAATGTGGGAGGCTTAGCTACCTGGTGTAGGGGTGTGATACCATGTTTAAGTCAGCTAAGCTACTATAATTCCACATTTAATGCATTCATCCTGAAGAGGTCAGTGTACGTTATTTTCCGTATGCTCTCTTTTACATCTTAACTTTGTGTGGTAGTGCCATAGCTCCCTATTCATCAAGCTTTGCACATTCTACCCTTATATATGTAATGAATTTTTTGAGTAATTCTAGTCATTGCTGAATAACTTATCGCCTTATCTAGGACTACCTCCTTTTTTACCCAAAAAAATTAAATTGAGTGCTATTCATTTGTTTTTGTTGCTTTTTCCGGGTTATGTATCCTATCCTTGTAAATTATTCAAAACTTGTGATGTGAAGGAGAAGGGTAGTCATTTATCTTGGATTCCAGTCATTGTTTTTGGAAAAAATTGAAAATATAGGATAAATGGTATGAATGGAAATATCAACAGCACGAAAAGTGTGGTTGGAACTTTTGAATACATTGATTGATGATATCTTTTGGTCAATTATAGTCACCTCTCAGGTATAGTATTGGCCTTTGGGACTACGAATTGGTTCACTAGAATAAAGATTGGGGACAAATGTGACGTGTCTGTATATATCAAAAGAAAAAAAGTGATCTGATTTGTCCAAAGCTGTTGAAAGGCAAATAATAGCTCTTTGAACCCAAAAAAAAAACACTGTAAAAGGAAAAGCACATAACCTTCACAGGCTTTCAATATCATACTATAAAGCTTTGTTTTTCTCTTCACCCTTTTTCAGAAGGATAAACAGTGACTTCCTCATCAGTGTGACCACATTGTGAAAGCTCCCCTTCCTTTCAAGTGTTGTCCTGGTTGGAACAACTCTGTATTCTCCCCAACATAACCCATAATCAAGAGAATTGTTGAACTAAGTTGTGGGAACTCTTCACTTTCGATGCCGCACCCATGTCGGATTCTCCAAAAATAGTGTATTTTTGGGGAATATGACATGCACCCAATTGACTTTTTTGAAGAATCCGAGCAACATAGTTGTTGAGTAGTTAATCAAACAATAATGTCCTAATTCTATATTAGTGTGATTTGCAGTATATTTCTCTGTATCCATTCCGCTCAGTTCAGGTTCAATATCACGGTTCAACTAATTTTACGATGATTGTAATAAGCCTATTAATAGATAAAAGATTAACTTCCATTTGCTTGTTGCTGCATATTTTAGATTTTCATTTCCACTGAAAGCATTTGCATATACAATATTTGTTAGCATTATGAAGAAACGAGCAAAAGAATAAAGGTGAATCAAATAAAGCATTGGCCTTTCTGCATAAGGAAATTCTAGGATGCATTCGAATCTGACGTCAACGGCGTAGATGGAGGTAGTGGGCTTCTGGTTTCGTTACTACTTTTGGATTTGGCAATTGGAGGAGACTGCAAAGAAACTGACGATGGATTTGAAGGAGTCGTCTTATTCGAGCTACTCTGAGGTTTCTTGTTGTTTGATGGATTTGAAGATGTAGTCGGCTTATTCGAGCTGCTCTGCGGTTTCTTCTTGTCCCCTTCACCCGATGGATTCGAAGACGGAGTCTGGTTATTCAAGCTGCTCTGAGGTTTCTTCTTGTCCCCTTCACCTGATGGATTCGAAGACGGAGTCTGGTTATTCAAGCTGCTCTGAGGTTTCTTGCTGTTCCCTTCACTTGCCTTGTCCTCCATCCGATGGTGGCGAACTTTGGTAGACTTTAAGCTCATTGGAACAATACAATGGCAAACGAAACCTACACACACCAACAATTACCAAGTTAATAAATTTCTCATATCAACAAGTTGCAGAACACTAACTATTCAAACTTGAGAACTTTATAACCCCTTCAATGACTGCAACAAGCTCATGTTGAAGATTCCTATCGTGTTAAAGGCTACGAAAAACTTTCAACGCAGCCACTGAACTCGGTAGGAATCTTCTGCAATGACTGCAAAAAATCATTGCCAAGAACAAACATTGAGACTTCTATGAAACTACCCAGAAGCTTCCTGATTCAACATGTCTAGAAAACAAGGTACTCTCTCCGTCCCAATTTATGTAACACTTAACTATTCAAAGTGTTTGAATTTTGACCAAACATTTTAAAATGTATTTTCCGAATATTCAAGCATGTGAAACAATAAATCTGACTGAAAATCTTACCAATTCTTGCTAGTCGATTAACCCAACTGGGAATTGGATTGTTAGTGAGTTTAACACAAGCATCACTACAGAAATGGTTACAATTCTTAGTAATAAGATTGTAAGCATTTCCTTTATACTTATCAGCAAGTTCTTCCATAACTCTTCTTACTTCTCTAGGTGTCTTCTCCGTCCATCCGATCAAAATTGTCTTCCTGAAAATAAATCCTTCACACTTTTTCGGCTCCCCTTCGAAAATTCCAGTACTTGGATACTCATGAGCTCCAAATGCATACTCAACTCCATGAACTAGCTACCATAACCCACCACCATCAACAAATAATTAGTATAAACTTTGCTACAGTTAATTAAGAAAAAACCAGATCTGAATAACAATGAAAAGTTTGAAGATTGAGTAACTGACCTTGAACGCCGGAATGGTAAACTCCGAGACCCAGCCAATAAGCATAACCATTCATCGGAGTTAGATCGTATACGTTTAAGTACACCGGCACCGATCCACGGCCATTATTAACTGATTTTTTGCAGCAGAGCATCTTGAAATAAAAACAGAGGAACTTTTCAGGTTCCTCGCATTGTGAAGATGTAAAGAAAGAGAATTCTGAATTTTCCGATTGATTGAGTTGCA

写回答
好问题 0 提建议
关注问题
分享
邀请回答
编辑收藏删除
收藏举报

2条回答默认最新

关注

码龄粉丝数原力等级 --

被采纳

被点赞

采纳率
CSDN专家-黄老师 2021-05-20 12:49
关注
参考一下：https://blog.csdn.net/david10251025/article/details/102093869

本回答被题主选为最佳回答 , 对您是否有帮助呢?

解决无用
评论打赏
分享
举报

评论

按下Enter换行，Ctrl+Enter发表内容

查看更多回答(1条)

报告相同问题？

关注问题

怎么用python修改fasta序列的ID名字 python
2022-04-19 08:29

回答 1 已采纳 import re def run(): fasta_str = """WP_018731760.1 NAD(P)-binding domain-containing protein [Sal
Python获取基因库序列并排序 python 有问必答
2021-05-29 00:51

回答 2 已采纳 list index out of range 错误出现的原因主要有两个，一个可能是下标超出范围，一个可能是list是空的，没有一个元素比如 a = random.randint(0,
怎么用Python将Excel转化为fasta？ python 有问必答
2022-04-11 12:04

回答 3 已采纳在github上有处理代码，可参考一下：https://github.com/kuod/cctf/blob/master/parser.py
串联fasta文件中的序列，去除ATG、TAA
2024-07-23 19:04

简若573的博客和 Shell 脚本位于同一目录，并正确设置文件路径和权限。目的：根据提供的模式修改序列 ID 并修剪特定碱基。
如何用python写一个循环，或者其他方式，完成批量选取内容并保存。 python
2019-06-27 21:38

回答 3 已采纳 from Bio import SeqIO for xx in range(1, 68): xn = "split__" + str(xx) + ".fasta" yn = "selec
如何实现读取多个文件夹下的fasta文件，大脑缺氧，求救命！ python
2019-07-03 21:26

回答 2 已采纳使用os.listdir('目录路径') 列表形式返回该目录下的所有文件名，然后你再依次读取和处理就ok了
初学python，批量选取文件中某一部分内容，然后将内容逐条存于一个文件中 python
2019-06-27 20:46

回答 2 已采纳 ``` # -*- coding:utf-8 -*- import os from Bio import SeqIO # root_dir为要读取文件的根目录 root_dir =
DNA-蛋白翻译过程的Python实现
2021-04-08 15:40

EmmettPeng的博客最近为了给平台上加上一个将DNA序列翻译为蛋白序列的工具，写了一个任何生信玩家初学时都会写的代码。看了一些别人的翻译工具，我也想尽量把代码写的完整一点，在这个过程中首次接触并使用了BioPython，目前看起来...
为什么老师的Python2能产生一个TXT文件而我的python3却不能？ python
2017-04-25 05:44

回答 3 已采纳 try this ``` if __name__=="__main__": main(args) ```
shell或python脚本批量处理txt文件 bash linux python 有问必答
2022-03-20 04:01

回答 6 已采纳 # -*- utf8 -*- """ pyCode : test.py Created by: Taro Description: 2022/3/20 0020 """ import os impor
re.sub()对A文件中每个以>开头的字符串，在B文件查找并以该字符串在B所在列名替换A中的该字符串 python 大数据有问必答正则表达式
2021-10-26 10:35

回答 1 已采纳是不是类似这样操作，循环替换。 import re s=['gene>ggcctt>tggc','gene>ctct'] s1=[] for x in s: x1=re.s
python for bioinformatics相关题目
2019-12-09 11:21

wangchuang2017的博客题目完整版来自：... ... 每个题可能有多种解法，不同解法用分别用## 1/2/3表示，通常## 1 是我自己脚本，而其他解法是参考其他网友的脚本。如果你也同我一样刚用python处理生信数据的话，请务必先自...
这个为生信学习打造的开源 Python 文字教程真香！！！
2021-01-11 21:39

生信宝典的博客欢迎来到Python的世界，本教程将带你遨游Python，领悟Python的魅力。本教程专注于帮助初学者，尤其是生物信息分析人员快速学会Python的常用功能和使用方式，因此只精选了部分...
序列和序列对象
2019-12-09 22:44

wangchuang2017的博客序列和序列对象 Seq 类 Seq类是Biopython最基础的一类, 储存序列信息. from Bio.Seq import Seq. 该类基本格式是Seq(self, data, alphabet=Alphabet()). 类似于字符串, 能够储存蛋白, DNA, RNA序列信息. 该类是不...
biopython：2：序列组成
2018-03-04 23:20

genome_denovo的博客生物信息当中的序列（DNA、RNA、amino acid）与python当中的字符串有两点不同（其余属性大致相同，比如序列长度，可以进行切片等）第一点不同：序列有 translate() ，即翻译；序列有reverse_complement() ，即反向...
送书 | 令附生信专用简明 Python 文字和视频教程
2020-12-13 22:04

生信宝典的博客欢迎来到Python的世界，本教程将带你遨游Python，领悟Python的魅力。本教程专注于帮助初学者，尤其是生物信息分析人员快速学会Python的常用功能和使用方式，因此只精选了部分...
pyfasta-0.4.1
2011-08-25 19:17

《Pyfasta：高效访问Fasta序列文件的Python库》 Pyfasta是一款专为生物信息学设计的Python库，它的核心功能是提供快速、内存高效且符合Python编程风格的接口，以便于处理和操作Fasta格式的序列文件。Fasta是一种...
第八章遗传密码
2018-09-08 16:36

weixin_30698527的博客到目前为止，我们已经使用Python搜索模体、模拟DNA突变、生成随机序列、将DNA转录为RNA，这些都是重要的活动，它们可以很好地介绍你可以用来研究生物系统的计算技术。在本章中，我们将编写Python程序来模拟遗传...
Perl的一些初级练习题
2020-05-15 23:25

bioittang的博客 1.读入一个文件，将序列变成一行输出 #!/usr/bin/perl -w use strict; use 5.026; open IN,"<$ARGV[0]",or die $!; open OUT,">out.txt",or die $!; while(<IN>){ chomp; print OUT $_; } print OUT ...
0049-【宏基因组】-16S分析qiime1极简教程2
2018-06-22 13:26

leadingsci的博客 : /Users/caporaso/.virtualenvs/qiime/lib/python2. 7 /site-packages/qiime_default_reference/gg_13_8_otus/rep_set/ 97 _otus.fasta] --sortmerna_e_value Maximum E-value when clustering [ default = 1 ] ...
没有解决我的问题, 去提问

比较两个fasta序列的碱基个数并提取ATG及之后的序列

2条回答 默认 最新

2条回答默认最新