24小时热门版块排行榜    

Znn3bq.jpeg
汕头大学海洋科学接受调剂
查看: 2073  |  回复: 17

正版小辉辉

铁虫 (初入文坛)

[求助] shell求助---删除特定行中的字符串

我有如下一个文件:
>S6_1
GTAGGGAATCTTCCGCAATGGACGAAAGTCTAACGGAGCAACAGAAACCCGGGTAGTCC
>S6_2
GTAGGGAATCTTCCGCAATGGACGAAAGTCTAGATACCCGTGTAGTCC
>S6_3
GTAGGGAATCTTCCGCAATGGACGAAAGTCTGTTAGATACCCTAGTAGTCC

我想删除每一行(碱基)中的最后17个字符(AGATACCCTAGTAGTCC),但是不知道用语言如何实现?
貌似用sed和awk比较容易。
回复此楼

» 猜你喜欢

哈哈
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

ljw4010

荣誉版主 (职业作家)

小木虫从头派教主

前提是每行AGATACCCTAGTAGTCC唯一

发自小木虫Android客户端

» 本帖已获得的红花(最新10朵)

不要说话
5楼2017-07-13 11:12:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

ljw4010

荣誉版主 (职业作家)

小木虫从头派教主

不要说话
2楼2017-07-12 22:53:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

正版小辉辉

铁虫 (初入文坛)

引用回帖:
2楼: Originally posted by ljw4010 at 2017-07-12 22:53:09
sed可以搞

我也明白,可是不会写代码

发自小木虫Android客户端
哈哈
3楼2017-07-12 22:56:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ljw4010

荣誉版主 (职业作家)

小木虫从头派教主

sed -i 's/AGATACCCTAGTAGTCC//g'   filename

发自小木虫Android客户端
不要说话
4楼2017-07-13 11:11:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

正版小辉辉

铁虫 (初入文坛)

送红花一朵
引用回帖:
5楼: Originally posted by ljw4010 at 2017-07-13 11:12:27
前提是每行AGATACCCTAGTAGTCC唯一

但是如果不一样呢?这样是不是就不好办了?
哈哈
6楼2017-07-13 11:28:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ljw4010

荣誉版主 (职业作家)

小木虫从头派教主

不一样就这么搞,分组正则,我没验过,你试下,sed -i 's/(.*)AGATACCCTAGTAGTCC$/\1/g'? ?filename

发自小木虫Android客户端
不要说话
7楼2017-07-13 11:39:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ljw4010

荣誉版主 (职业作家)

小木虫从头派教主

sed -i 's/(.*)AGATACCCTAGTAGTCC$/\1/g' filename

发自小木虫Android客户端
不要说话
8楼2017-07-13 11:40:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

正版小辉辉

铁虫 (初入文坛)

引用回帖:
8楼: Originally posted by ljw4010 at 2017-07-13 11:40:50
sed -i 's/(.*)AGATACCCTAGTAGTCC$/\1/g' filename

qiime@qiime-190-virtual-box:~$ sed -i 's/(.*)ATTAGATACCCTAGTAGTCC$/\1/g' /home/qiime/S1_S14.fasta
sed: -e expression #1, char 36: invalid reference \1 on `s' command's RHS
哈哈
9楼2017-07-13 12:19:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

正版小辉辉

铁虫 (初入文坛)

引用回帖:
9楼: Originally posted by 正版小辉辉 at 2017-07-13 12:19:34
qiime@qiime-190-virtual-box:~$ sed -i 's/(.*)ATTAGATACCCTAGTAGTCC$/\1/g' /home/qiime/S1_S14.fasta
sed: -e expression #1, char 36: invalid reference \1 on `s' command's RHS...

这样看起来,有错误出现
哈哈
10楼2017-07-13 12:20:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 正版小辉辉 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 考研求调剂 +12 子木呐 2026-04-12 13/650 2026-04-14 01:19 by 王珺璞
[考研] 调剂求收留 +31 果然有我 2026-04-10 32/1600 2026-04-13 21:35 by gretl
[考研] 一志愿沪9,326求生物学调剂 +9 刘墨墨 2026-04-13 9/450 2026-04-13 18:59 by lbsjt
[考研] 一志愿211 0703化学 346分求调剂 +26 土豆er? 2026-04-09 29/1450 2026-04-13 15:15 by 独醉梦孤城
[考研] 297工科,求调剂? +13 河南农业大学-能 2026-04-12 13/650 2026-04-13 14:12 by dingyanbo1
[考研] 本科南方医科大学 一志愿985 药学学硕284分 求调剂 +5 弱水听文 2026-04-09 5/250 2026-04-12 13:16 by shengxi123
[考研] 药学305求调剂 +8 玛卡巴卡boom 2026-04-10 8/400 2026-04-12 00:07 by zhouwenxian
[考研] 331求调剂 +5 王国帅 2026-04-11 5/250 2026-04-11 22:56 by 溪涧流水
[考研] 296求调剂 +14 汪!?! 2026-04-08 15/750 2026-04-11 20:28 by dongdian1
[考研] 11408。358求调剂 +3 TMYzds 2026-04-07 3/150 2026-04-11 17:10 by 氮气气气
[考研] 352 求调剂 +6 yzion 2026-04-11 8/400 2026-04-11 16:24 by 明月此时有
[考研] 281求调剂 +11 觉得好的吧 2026-04-10 11/550 2026-04-11 09:35 by 逆水乘风
[考研] 一志愿211,化学310分,本科重点双非,求调剂 +23 努力奋斗112 2026-04-08 23/1150 2026-04-10 23:29 by 314126402
[考研] 287求调剂 +15 Fnhc 2026-04-07 21/1050 2026-04-10 19:09 by chemisry
[考研] 初试261 +3 Asht少 2026-04-10 6/300 2026-04-10 16:38 by Asht少
[考研] 285求调剂 +9 AZMK 2026-04-07 11/550 2026-04-10 15:24 by AZMK
[考研] 材料专硕调剂 +16 哈哈哈吼吼吼哈 2026-04-07 17/850 2026-04-09 21:16 by wutongshun
[考研] 328求调剂 +17 lftmya 2026-04-07 18/900 2026-04-09 08:05 by 5268321
[考研] 电子信息346 +4 zuoshaodian 2026-04-08 4/200 2026-04-08 11:54 by zzucheup
[考研] 277、学硕,求调剂 数一104, +11 瓶子PZ 2026-04-07 12/600 2026-04-07 23:30 by 一只好果子?
信息提示
请填处理意见