亚洲av成人无遮挡网站在线观看,少妇性bbb搡bbb爽爽爽,亚洲av日韩精品久久久久久,兔费看少妇性l交大片免费,无码少妇一区二区三区

  免費注冊 查看新帖 |

Chinaunix

  平臺 論壇 博客 文庫
12下一頁
最近訪問板塊 發(fā)新帖
查看: 3536 | 回復: 16
打印 上一主題 下一主題

求助小腳本編寫 [復制鏈接]

論壇徽章:
1
15-16賽季CBA聯(lián)賽之天津
日期:2017-03-21 18:26:14
跳轉到指定樓層
1 [收藏(0)] [報告]
發(fā)表于 2017-02-26 09:40 |只看該作者 |倒序瀏覽
腳本讀入兩個文件 lst1和lst2
lst1:
RRA45973,RRA45973,CACTTCGA,,
RRA45974,RRA45974,CAGCGTTA,,
RRA45975,RRA45975,CATACCAA,,
RRA45976,RRA45976,CCAGTTCA,,
RRL06701,RRL06701,ATCCTGTA,,
需求把后面的8個堿基加兩個逗號刪除,堿基有的時候是6個。
lst2:
ATTACTCG        TATAGCCT
ATTACTCG        ATAGAGGC
ATTACTCG        CCTATCCT
ATTACTCG        GGCTCTGA
ATTACTCG        AGGCGAAG
要求把中間的空格換成逗號隔開。
寫一個腳本把做完修改的兩個文件對應行粘到一起,結果就像
RRA45973,RRA45973,ATTACTCG,TATAGCCT
.
.
.
請教各位大神應該怎么寫,小弟初學,自己寫的運行起來有問題。

論壇徽章:
307
程序設計版塊每周發(fā)帖之星
日期:2016-04-08 00:41:33操作系統(tǒng)版塊每日發(fā)帖之星
日期:2015-09-02 06:20:00每日論壇發(fā)貼之星
日期:2015-09-02 06:20:00程序設計版塊每日發(fā)帖之星
日期:2015-09-04 06:20:00每日論壇發(fā)貼之星
日期:2015-09-04 06:20:00每周論壇發(fā)貼之星
日期:2015-09-06 22:22:00程序設計版塊每日發(fā)帖之星
日期:2015-09-09 06:20:00程序設計版塊每日發(fā)帖之星
日期:2015-09-19 06:20:00程序設計版塊每日發(fā)帖之星
日期:2015-09-20 06:20:00每日論壇發(fā)貼之星
日期:2015-09-20 06:20:00程序設計版塊每日發(fā)帖之星
日期:2015-09-22 06:20:00程序設計版塊每日發(fā)帖之星
日期:2015-09-24 06:20:00
2 [報告]
發(fā)表于 2017-02-26 13:43 |只看該作者
能否將你的要求說的詳細一點,對生物信息的知識比較匱乏。另外,將希望得到的結果提出來。

論壇徽章:
1
15-16賽季CBA聯(lián)賽之天津
日期:2017-03-21 18:26:14
3 [報告]
發(fā)表于 2017-02-26 14:02 |只看該作者
回復 2# sunzhiguolu
就是說lst1那個文件每行的最后兩個逗號加上【ATGC】組成的8個(有的時候是六個)字符我要把它刪了,替換成list2中的相應行,但是list2中的每行要把中間的空格換成逗號,最終的結果就是我貼子上已經(jīng)粘的那個結果,只不過我只是貼了一行,已知lst1和lst2這兩個文件行數(shù)是一樣的,有多少行就輸出多少個合并后的行。

論壇徽章:
1
15-16賽季CBA聯(lián)賽之天津
日期:2017-03-21 18:26:14
4 [報告]
發(fā)表于 2017-02-26 14:10 |只看該作者
回復 2# sunzhiguolu
lst1每行改完后的結果就像這樣:RRA45973,RRA45973,
lst2每行改完后的結果就像這樣:ATTACTCG,TATAGCCT
最終希望的結果是把這兩行粘成一行
RRA45973,RRA45973,ATTACTCG,TATAGCCT
lst1中有多少行就打印出多少行。

論壇徽章:
28
15-16賽季CBA聯(lián)賽之八一
日期:2016-02-22 19:10:4215-16賽季CBA聯(lián)賽之深圳
日期:2016-12-01 10:34:0415-16賽季CBA聯(lián)賽之新疆
日期:2016-12-07 10:24:2915-16賽季CBA聯(lián)賽之同曦
日期:2016-12-15 12:06:43CU十四周年紀念徽章
日期:2016-12-18 13:03:4415-16賽季CBA聯(lián)賽之吉林
日期:2017-01-03 15:52:2515-16賽季CBA聯(lián)賽之遼寧
日期:2017-01-04 14:58:2415-16賽季CBA聯(lián)賽之遼寧
日期:2017-01-15 09:42:512016科比退役紀念章
日期:2017-02-06 17:21:50黑曼巴
日期:2017-02-10 15:46:1215-16賽季CBA聯(lián)賽之上海
日期:2017-03-18 10:14:5415-16賽季CBA聯(lián)賽之青島
日期:2017-03-18 22:00:44
5 [報告]
發(fā)表于 2017-02-26 14:25 |只看該作者

  1. paste -d '' <(sed 's/\w*,*$//' lst1) <(sed -r 's/\s+/,/' lst2)
復制代碼

?

論壇徽章:
1
15-16賽季CBA聯(lián)賽之天津
日期:2017-03-21 18:26:14
6 [報告]
發(fā)表于 2017-02-26 14:35 |只看該作者
回復 5# moperyblue
因為這是這兩個文件的一部分,如果只是用paste命令用正則的話會把一些不需要修改的地方也改了,我粘出來的只是這兩個文件需要修改的部分。所以我還是希望能有一個只是針對上面我給出的兩個文件的小腳本。

論壇徽章:
28
15-16賽季CBA聯(lián)賽之八一
日期:2016-02-22 19:10:4215-16賽季CBA聯(lián)賽之深圳
日期:2016-12-01 10:34:0415-16賽季CBA聯(lián)賽之新疆
日期:2016-12-07 10:24:2915-16賽季CBA聯(lián)賽之同曦
日期:2016-12-15 12:06:43CU十四周年紀念徽章
日期:2016-12-18 13:03:4415-16賽季CBA聯(lián)賽之吉林
日期:2017-01-03 15:52:2515-16賽季CBA聯(lián)賽之遼寧
日期:2017-01-04 14:58:2415-16賽季CBA聯(lián)賽之遼寧
日期:2017-01-15 09:42:512016科比退役紀念章
日期:2017-02-06 17:21:50黑曼巴
日期:2017-02-10 15:46:1215-16賽季CBA聯(lián)賽之上海
日期:2017-03-18 10:14:5415-16賽季CBA聯(lián)賽之青島
日期:2017-03-18 22:00:44
7 [報告]
發(fā)表于 2017-02-26 17:44 |只看該作者
回復 6# kevin_furant

可以只針對符合條件的記錄修改后 再 連接起來

論壇徽章:
0
8 [報告]
發(fā)表于 2017-02-26 18:04 |只看該作者
感觉怎么像index

論壇徽章:
1
15-16賽季CBA聯(lián)賽之天津
日期:2017-03-21 18:26:14
9 [報告]
發(fā)表于 2017-02-26 18:09 |只看該作者
回復 8# 華小飛_Perl
其實就是index,拆分雙端index的時候有些東西需要改,顯得麻煩所以就想搞一搞。樓上的方法其實就能實現(xiàn),不過我正在學perl,我想看一下別人怎么寫的  為什么我自己寫的不能達到效果。

論壇徽章:
0
10 [報告]
發(fā)表于 2017-02-26 19:04 |只看該作者
回復 9# kevin_furant

樓上大神們的建議很不錯,雙端拆分的話應該是miseq / nextseq / miniseq平臺的下機數(shù)據(jù)咯~
您需要登錄后才可以回帖 登錄 | 注冊

本版積分規(guī)則 發(fā)表回復

  

北京盛拓優(yōu)訊信息技術有限公司. 版權所有 京ICP備16024965號-6 北京市公安局海淀分局網(wǎng)監(jiān)中心備案編號:11010802020122 niuxiaotong@pcpop.com 17352615567
未成年舉報專區(qū)
中國互聯(lián)網(wǎng)協(xié)會會員  聯(lián)系我們:huangweiwei@itpub.net
感謝所有關心和支持過ChinaUnix的朋友們 轉載本站內(nèi)容請注明原作者名及出處

清除 Cookies - ChinaUnix - Archiver - WAP - TOP