24小时热门版块排行榜    

查看: 2387  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 正版小辉辉 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 XJXMuRjD7cK7 2026-11-13 4/200 2026-11-15 06:50 by xSrQNSSjcgoX
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 XJXMuRjD7cK7 2026-11-13 3/150 2026-11-15 05:33 by xSrQNSSjcgoX
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 uKOSCocCb6Fm 2026-11-12 7/350 2026-11-15 03:32 by xSrQNSSjcgoX
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 uKOSCocCb6Fm 2026-11-12 6/300 2026-11-15 02:59 by xSrQNSSjcgoX
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 uKOSCocCb6Fm 2026-11-12 10/500 2026-11-15 02:15 by xSrQNSSjcgoX
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 uKOSCocCb6Fm 2026-11-12 9/450 2026-11-15 02:04 by xSrQNSSjcgoX
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +7 0TEMgxdT1eIw 2026-11-09 10/500 2026-11-15 00:39 by Equinoxhua
[硕博家园] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 G2rhSvqvhNzl 2026-11-10 6/300 2026-11-15 00:27 by Equinoxhua
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 G2rhSvqvhNzl 2026-11-10 6/300 2026-11-15 00:24 by Equinoxhua
[考博] 27博士申请(市政/环境) +3 ladf 2026-11-12 3/150 2026-11-13 08:19 by Equinoxhua
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 TgvkbZZl9tqB 2026-11-11 4/200 2026-11-12 12:51 by ECLTXnJnRSXf
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 TgvkbZZl9tqB 2026-11-11 4/200 2026-11-12 12:48 by ECLTXnJnRSXf
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 TgvkbZZl9tqB 2026-11-11 6/300 2026-11-12 12:42 by ECLTXnJnRSXf
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 cdqGD6ftHmbl 2026-11-11 5/250 2026-11-12 12:28 by ECLTXnJnRSXf
[基金申请] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 Z7qfCsdIiD1T 2026-11-11 5/250 2026-11-12 11:44 by HajghKOMynmz
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 Z7qfCsdIiD1T 2026-11-10 6/300 2026-11-12 11:22 by HajghKOMynmz
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 C6CRn66RYYXS 2026-11-10 7/350 2026-11-12 10:16 by HajghKOMynmz
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 YHiw1gou1Eda 2026-11-10 6/300 2026-11-12 10:08 by HajghKOMynmz
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +7 FChKX8meu1ds 2026-11-09 9/450 2026-11-12 10:03 by HajghKOMynmz
[论文投稿] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 cdqGD6ftHmbl 2026-11-11 3/150 2026-11-12 07:53 by HajghKOMynmz
信息提示
请填处理意见