| 查看: 820 | 回复: 0 | ||
emanlee木虫 (小有名气)
|
[求助]
从两个不同fasta文件中寻找不重复的序列
|
|
题目：从两个不同fasta文件中寻找不重复的序列 第一个fasta文件aaa.fa中有40000条碱基序列或者氨基酸序列： >gi|118600994|ref|NM_001079530.1| Homo sapiens cripto, FRL-1, cryptic family 1B (CFC1B), mRNA ATGCCAAATACAGCCATGAAGAAAAAGGTGCTGCTGATGGGGAAGAGCGGGTCGGGGAAGACCAGCATGAGGTCGATAATCTT >gi|57863286|ref|NM_006570.4| Homo sapiens Ras-related GTP binding A (RRAGA), mRNA ACGCTCTACAAAGCCTGGTCCAGCATCGTCTACCAGCTGATTCCCAACGTTCAGCAGCTGGAGATGAACCTCAGGAATTTTG >gi|254587897|ref|NM_178495.5| Homo sapiens inositol 1,4,5-trisphosphate receptor CGCCAATTACATTGCTCGCGACACCCGGCGCCTGGGGGCCACCATTGACGTGGAACACTCCCACGTCCGATTCCTAGGGAACC >gi|191252813|ref|NM_001128635.1| Homo sapiens RIMS binding protein 3B (RIMBP3B), mRNA TGGTGCTGAACCTGTGGGACTGTGGCGGTCAGGACACCTTCATGGAAAATTACTTCACCAGCCAGCGAGACAATATCTTCCGTA >gi|61656209|ref|NM_001013355.1| Homo sapiens olfactory receptor, family 2, subfamily G, member 6 (OR2G6), mRNA ACGTGGAAGTTTTGATTTACGTGTTTGACGTGGAGAGCCGCGAACTGGAAAAGGACATGCATTATTACCAGTCGTGTCTGGAGG 第二个fasta文件bbb.fa中有40000条碱基序列或者氨基酸序列： >gi|83267870|ref|NM_080431.4| Homo sapiens actin-related protein T2 (ACTRT2), mRNA CCATCCTCCAGAACTCTCCTGACGCCAAAATCTTCTGCCTGGTGCACAAAATGGATCTGGTTCAGGAGGATCAGCGTGACCTGA >gi|53828675|ref|NM_001001923.1| Homo sapiens olfactory receptor, family 5, subfamily C, member 1 (OR5C1), mRNA TTTTTAAAGAGCGAGAGGAAGACCTGAGGCGTCTGTCTCGCCCGCTGGAGTGTGCTTGTTTTCGAACGTCCATCTGGGATGAG >gi|52627150|ref|NM_001005276.1| Homo sapiens olfactory receptor, family 2, subfamily AE, member 1 (OR2AE1), mRNA TTTTTAAAGAGCGAGAGGAAGACCTGAGGCGTCTGTCTCGCCCGCTGGAGTGTGCTTGTTTTCGAACGTCCATCTGGGATGAG >gi|61656211|ref|NM_001013357.1| Homo sapiens olfactory receptor, family 8, subfamily U, member 9 (OR8U9), mRNA ACGCTCTACAAAGCCTGGTCCAGCATCGTCTACCAGCTGATTCCCAACGTTCAGCAGCTGGAGATGAACCTCAGGAATTTTG >gi|51871366|ref|NM_001004124.1| Homo sapiens olfactory receptor, family 4, subfamily P, member 4 (OR4P4), mRNA ATGCCAAATACAGCCATGAAGAAAAAGGTGCTGCTGATGGGGAAGAGCGGGTCGGGGAAGACCAGCATGAGGTCGATAATCTT 我们想从第二个文件bbb.fa中找出与aaa.fa中的序列不overlapping的序列（overlap-ratio<0.8），请问如何使用blast比对？ 是否有现成的perl或者python，或者C代码可以直接使用？ |
回复此楼» 猜你喜欢
Wnt-C59：PORCN靶向化学探针，Wnt信号通路阻断机制与科研应用全解析
已经有0人回复
-
Fexaramine：肠道限制性FXR化学探针，结构特征与信号通路研究全解析
已经有0人回复
-
化学工程及工业化学论文润色/翻译怎么收费?
已经有280人回复
-
NVS-STG2：靶向STING跨膜域的分子胶类激动剂，结构特征与寡聚化机制研究全解析
已经有0人回复
-
DOTA-GPC3肽（CAS 3066439-94-6）：双功能分子探针前体的结构特征与标记应用解析
已经有0人回复
-
科研用 TLR7/8 双激动剂怎么选？现货与质控要点
已经有0人回复
-
SEL120-34A HCl：CDK8/CDK19高选择性激酶抑制剂
已经有0人回复
-
氢氧化钾活化操作
已经有18人回复
-
散金
已经有68人回复
-
生物质超级电容器
已经有0人回复
找到一些相关的精华帖子，希望有用哦~
如何把一个氨基酸序列转化成FASTA格式
已经有7人回复
- 从uniprot中下载的氨基酸序列的fasta格式中的PE,SV分别指什么
已经有4人回复
- 在一个TXT文本中，以fasta格式保存着1000多个氨基酸序列，每个序列的长度如何得出
已经有5人回复
- 请问如何将这些序列转化为FASTA格式？
已经有8人回复
- 如何处理大量的fasta格式的rna序列
已经有16人回复
- 菜鸟求助——着急求助解答：在NCBI上为什么下载不了FASTA格式的基因序列
已经有9人回复
- 【求助/交流】有登陆号，如何批量从Genbank中导出fasta格式的核酸序列
已经有6人回复
- 如何将序列转换为FASTA格式？
已经有11人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
