24小时热门版块排行榜    

查看: 2260  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 正版小辉辉 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 大佬帮着分析今年的评审意见?该课题连续申请了三年了,明年要不要换个方向。 100+3 cai-11138 2026-08-31 11/550 2026-09-06 20:49 by archvillain
[文学芳草园] 两块石头 +7 阿美_Lml888 2026-09-03 10/500 2026-09-06 20:08 by 阿美_Lml888
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 j5mmowWXNWxv 2026-09-05 4/200 2026-09-06 19:21 by eObJPa4gWmp6
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 jPDp0sc2B7zQ 2026-09-05 6/300 2026-09-06 18:21 by eObJPa4gWmp6
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 jPDp0sc2B7zQ 2026-09-05 4/200 2026-09-06 18:09 by eObJPa4gWmp6
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 6/300 2026-09-06 18:09 by eObJPa4gWmp6
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-04 8/400 2026-09-06 17:09 by eObJPa4gWmp6
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 j5mmowWXNWxv 2026-09-05 4/200 2026-09-06 12:49 by eObJPa4gWmp6
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 7/350 2026-09-06 11:20 by eObJPa4gWmp6
[基金申请] 科研人应该花精力去思考如何解决问题,而不是去凝练问题 +10 瞬息宇宙 2026-09-01 19/950 2026-09-06 11:15 by jiangxuan2004
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 8/400 2026-09-06 11:09 by eObJPa4gWmp6
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 jPDp0sc2B7zQ 2026-09-04 8/400 2026-09-06 11:00 by eObJPa4gWmp6
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 7ZYIEW2YWn72 2026-09-04 9/450 2026-09-06 06:00 by E1iiBLMU2XnD
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +8 CRJ0Aq2FtVH0 2026-09-03 15/750 2026-09-06 05:48 by E1iiBLMU2XnD
[基金申请] 面上没中,邀请各位路过的虫友分析一下分数 +12 jackleilei 2026-09-03 13/650 2026-09-05 18:10 by fireguard
[文学芳草园] 初秋的晨风 +4 阿美_Lml888 2026-09-04 6/300 2026-09-05 13:51 by 阿美_Lml888
[基金申请] 学科评审组评审是指会评吗? +6 瞬息宇宙 2026-08-31 9/450 2026-09-05 13:11 by haizai99
[硕博家园] Ei源刊怎么投 +4 Fengshun9711 2026-09-04 4/200 2026-09-04 19:30 by 研途知予
[硕博家园] 哈尔滨工业大学韩晓军教授课题组招收2027年硕士推免生及博士研究生 +3 灯灯灯灯DDDD 2026-09-03 3/150 2026-09-03 21:22 by 科研皮皮猪
[基金申请] 要骂人了,新模版改版就是要淡化问题凝练这种虚的东西,结果有个评委还在说凝练得不够 +9 瞬息宇宙 2026-08-31 17/850 2026-09-02 14:04 by anjeeshine
信息提示
请填处理意见