比较两个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条)

报告相同问题？

关注问题

比较两个fasta序列的碱基个数并提取ATG及之后的序列 python 有问必答
2021-05-19 15:51

回答 2 已采纳参考一下：https://blog.csdn.net/david10251025/article/details/102093869
怎么用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 文字教程真香！！！
2021-01-11 21:39

生信宝典的博客欢迎来到Python的世界，本教程将带你遨游Python，领悟Python的魅力。本教程专注于帮助初学者，尤其是生物信息分析人员快速学会Python的常用功能和使用方式，因此只精选了部分...
怎么用Python将Excel转化为fasta？ python 有问必答
2022-04-11 12:04

回答 3 已采纳在github上有处理代码，可参考一下：https://github.com/kuod/cctf/blob/master/parser.py
如何用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了
序列和序列对象
2019-12-09 22:44

wangchuang2017的博客序列和序列对象 Seq 类 Seq类是Biopython最基础的一类, 储存序列信息. from Bio.Seq import Seq. 该类基本格式是Seq(self, data, alphabet=Alphabet()). 类似于字符串, 能够储存蛋白, DNA, RNA序列信息. 该类是不...
初学python，批量选取文件中某一部分内容，然后将内容逐条存于一个文件中 python
2019-06-27 20:46

回答 2 已采纳 ``` # -*- coding:utf-8 -*- import os from Bio import SeqIO # root_dir为要读取文件的根目录 root_dir =
为什么老师的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
python for bioinformatics相关题目
2019-12-09 11:21

wangchuang2017的博客题目完整版来自：... ... 每个题可能有多种解法，不同解法用分别用## 1/2/3表示，通常## 1 是我自己脚本，而其他解法是参考其他网友的脚本。如果你也同我一样刚用python处理生信数据的话，请务必先自...
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
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的常用功能和使用方式，因此只精选了部分...
第八章遗传密码
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 ...
biopython有什么用_BioPython学习笔记
2020-12-04 09:24

weixin_39551188的博客序列和序列对象Seq 类Seq类是Biopython最基础的一类, 储存序列信息. from Bio.Seq import Seq. 该类基本格式是Seq(self, data, alphabet=Alphabet()). 类似于字符串, 能够储存蛋白, DNA, RNA序列信息. 该类是不可变...
含有多段rna的.fasta文件读取，并打印最长的阅读框
2022-11-21 09:07

摸鱼小李的博客阅读框指基因序列中的一段无终止序列打断的碱基序列，可编码相应的蛋白，包含起始密码子（ATG）和终止密码子（TAA，TAG或TGA），如ATGAAATAG 是长度为9的一个阅读框。分别为阅读框1，2，和3。
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 ] ...
没有解决我的问题, 去提问

悬赏问题

¥15 AT89C51控制8位八段数码管显示时钟。
¥15 真我手机蓝牙传输进度消息被关闭了，怎么打开？(关键词-消息通知)
¥15 下图接收小电路，谁知道原理
¥15 装 pytorch 的时候出了好多问题，遇到这种情况怎么处理？
¥20 IOS游览器某宝手机网页版自动立即购买JavaScript脚本
¥15 手机接入宽带网线，如何释放宽带全部速度
¥30 关于#r语言#的问题：如何对R语言中mfgarch包中构建的garch-midas模型进行样本内长期波动率预测和样本外长期波动率预测
¥15 ETLCloud 处理json多层级问题
¥15 matlab中使用gurobi时报错
¥15 这个主板怎么能扩出一两个sata口

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

2条回答 默认 最新

悬赏问题

2条回答默认最新