| 查看: 796 | 回复: 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代码可以直接使用？ |
回复此楼» 猜你喜欢
留学--博士招生
已经有16人回复
-
化学工程，本211，求调剂，ab区不限
已经有4人回复
-
化学工程及工业化学论文润色/翻译怎么收费?
已经有208人回复
-
湖北师范大学招调剂啦
已经有8人回复
-
邵阳学院食品与化学工程学院硕士调剂
已经有0人回复
-
天津商业大学 生物与医药课题组招生
已经有0人回复
-
生物化工学硕招收调剂
已经有0人回复
-
广东石油化工学院2026年硕士研究生需要多名生物，制药工程，食品专业的调剂生
已经有0人回复
-
317分 一志愿江南大学 化学工程学硕 求调剂
已经有6人回复
-
化工申博
已经有3人回复
找到一些相关的精华帖子，希望有用哦~
如何把一个氨基酸序列转化成FASTA格式
已经有7人回复
- 从uniprot中下载的氨基酸序列的fasta格式中的PE,SV分别指什么
已经有4人回复
- 在一个TXT文本中，以fasta格式保存着1000多个氨基酸序列，每个序列的长度如何得出
已经有5人回复
- 请问如何将这些序列转化为FASTA格式？
已经有8人回复
- 如何处理大量的fasta格式的rna序列
已经有16人回复
- 菜鸟求助——着急求助解答：在NCBI上为什么下载不了FASTA格式的基因序列
已经有9人回复
- 【求助/交流】有登陆号，如何批量从Genbank中导出fasta格式的核酸序列
已经有6人回复
- 如何将序列转换为FASTA格式？
已经有11人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
