版塊導航: 正在加載中...

登錄注冊

應《網(wǎng)絡安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網(wǎng)跟帖服務。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

emanlee

木蟲 (小有名氣)

應助: 28 (小學生)
金幣: 3521.6
散金: 100
紅花: 2
帖子: 116
在線: 125.6小時
蟲號: 1466309
注冊: 2011-10-29
性別: GG
專業(yè): 計算機應用技術

[求助] 從兩個不同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比對？
是否有現(xiàn)成的perl或者python，或者C代碼可以直接使用？

回復此樓

» 猜你喜歡

不合理蛙科研實驗中的趣事：實驗器材的 “烏龍” 已經有0人回復
不合理蛙科研實驗中的趣事：和實驗材料的 “斗智斗勇” 已經有0人回復
化學工程及工業(yè)化學論文潤色/翻譯怎么收費? 已經有115人回復
蛋白質檢測：精準分析，解鎖生物分子的密碼已經有0人回復
不合理蛙科研實驗之小鼠實驗：嚴謹設計，解析生命機制的重要載體已經有0人回復
不合理蛙科研實驗之重金屬檢測：精準篩查，守護健康與環(huán)境的防線已經有0人回復
不合理蛙科研實驗之“蛙測重金屬我背鍋三千” 已經有0人回復
不合理蛙科研實驗之“鼠逃三次我發(fā)三篇SCI” 已經有0人回復
納米粒度分析已經有3人回復
林科院林化所招收材料與化工專業(yè)碩士，坐標南京已經有0人回復

1樓 2013-08-25 08:33:22

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

智能機器人

Robot (super robot)

我們都愛小木蟲

找到一些相關的精華帖子，希望有用哦~

如何把一個氨基酸序列轉化成FASTA格式已經有7人回復
從uniprot中下載的氨基酸序列的fasta格式中的PE,SV分別指什么已經有4人回復
在一個TXT文本中，以fasta格式保存著1000多個氨基酸序列，每個序列的長度如何得出已經有5人回復
請問如何將這些序列轉化為FASTA格式？已經有8人回復
如何處理大量的fasta格式的rna序列已經有16人回復
菜鳥求助——著急求助解答：在NCBI上為什么下載不了FASTA格式的基因序列已經有9人回復
【求助/交流】有登陸號，如何批量從Genbank中導出fasta格式的核酸序列已經有6人回復
如何將序列轉換為FASTA格式？已經有11人回復

點擊這里搜索更多相關資源

科研從小木蟲開始，人人為我，我為人人

相關版塊跳轉我要訂閱樓主 emanlee 的主題更新