| 查看: 837 | 回复: 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代码可以直接使用？ |
回复此楼» 猜你喜欢
生物质炭与活化剂氢氧化钾混合注意事项
已经有0人回复
-
科研人必读：新一代 KRAS G12C 抑制剂 Olomorasib 的研究进展
已经有0人回复
-
化学工程及工业化学论文润色/翻译怎么收费?
已经有242人回复
-
BT-11（Omilancor）研究化合物全解析：从分子结构到实验应用要点
已经有0人回复
-
富马酸贝达喹啉：二芳基喹啉类工具化合物的化学特征
已经有0人回复
-
BT投稿with editor状态快两周了，正常吗
已经有25人回复
-
Silevertinib（BDTX-1535）结构拆解：炔基双环片段为何成为合成瓶颈
已经有0人回复
-
RO4929097 解析：双苯并氮杂䓬骨架与 Notch 通路实验设计要点
已经有0人回复
-
分子量 1157 的 PROTAC 怎么配不析出：SGK3-PROTAC1 的溶解、分装与图谱复核
已经有0人回复
-
logP 6.1 的羧酸怎么配母液：GSK 650394 的溶解窗口与稀释顺序
已经有0人回复
-
从"水不溶"到 71 mg/mL：雷特格韦的 pH 依赖溶解与螯合色谱边界
已经有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人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
