首页 > 其他分享 >seqkit 软件根据序列ID删除指定的序列

seqkit 软件根据序列ID删除指定的序列

时间:2023-10-14 09:33:16浏览次数:31  
标签:test1 ## pc1 ID fa base 序列 root seqkit

 

001、 单个删除

(base) [root@pc1 test1]# ls
a.fa
(base) [root@pc1 test1]# cat a.fa                         ## 测试文件
>chr1
tttcccggg
>chr2
tttgggjjj
cccjjjjjj
>chr3
ccc
>chr4
aaaaatt
(base) [root@pc1 test1]# seqkit grep -v -p "chr1" a.fa      ## 删除chr1
>chr2
tttgggjjjcccjjjjjj
>chr3
ccc
>chr4
aaaaatt

 

(base) [root@pc1 test1]# ls
a.fa
(base) [root@pc1 test1]# cat a.fa     ## 测试fasta
>chr1
tttcccggg
>chr2
tttgggjjj
cccjjjjjj
>chr3
ccc
>chr4
aaaaatt
(base) [root@pc1 test1]# seqkit grep -v -p "chr2" a.fa     ## 删除chr2
>chr1
tttcccggg
>chr3
ccc
>chr4
aaaaatt

 

(base) [root@pc1 test1]# ls
a.fa
(base) [root@pc1 test1]# cat a.fa                      ## 测试fasta
>chr1
tttcccggg
>chr2
tttgggjjj
cccjjjjjj
>chr3
ccc
>chr4
aaaaatt
(base) [root@pc1 test1]# seqkit grep -v -p "chr4" a.fa    ## 删除chr4
>chr1
tttcccggg
>chr2
tttgggjjjcccjjjjjj
>chr3
ccc

 

002、批量删除

(base) [root@pc1 test1]# ls
a.fa  rm.list
(base) [root@pc1 test1]# cat a.fa         ## 测试fasta
>chr1
tttcccggg
>chr2
tttgggjjj
cccjjjjjj
>chr3
ccc
>chr4
aaaaatt
(base) [root@pc1 test1]# cat rm.list     ## 删除列表
chr2
chr4
(base) [root@pc1 test1]# seqkit grep -v -f rm.list a.fa     ## 删除命令
[INFO] 2 patterns loaded from file
>chr1
tttcccggg
>chr3
ccc

 。

 

标签:test1,##,pc1,ID,fa,base,序列,root,seqkit
From: https://www.cnblogs.com/liujiaxin2018/p/17763700.html

相关文章

  • python实现fasta文件碱基序列每行按照指定数目输出
     001、(base)[root@pc1test1]#lsa.fatest.py(base)[root@pc1test1]#cata.fa##测试fasta>chr1tttcccggg>chr2tttgggjjjcccjjjjjj>chr3ccc>chr4aaaaatt(base)[root@pc1test1]#cattest.py##程序#!/usr/bin/envpython3#......
  • 在Python中使用LSTM和PyTorch进行时间序列预测|附代码数据
    全文链接:http://tecdat.cn/?p=8145最近我们被客户要求撰写关于LSTM的研究报告,包括一些图形和统计输出。顾名思义,时间序列数据是一种随时间变化的数据类型。例如,24小时内的温度,一个月内各种产品的价格,一年中特定公司的股票价格诸如长期短期记忆网络(LSTM)之类的高级深度学习模型能......
  • python 实现统计fasta文件每一条序列的长度
     001、a、[root@pc1test1]#lsa.fatest.py[root@pc1test1]#cata.fa##测试fasta>chr1tttcccggg>chr2tttgggccc>chr3cccttt>chr4aaaaattt[root@pc1test1]#cattest.py##统计每条序列的长度#!/usr/bin/envpython3#-*-coding:......
  • python 中序列ID从fasta文件中批量提取序列数据
     001、[root@pc1test1]#lsa.fachr.listtest.py[root@pc1test1]#cata.fa##测试fasta文件>chr1tttcccggg>chr2tttgggccc>chr3cccttt>chr4aaaaattt[root@pc1test1]#catchr.list##序列IDchr2chr4 [root@pc1......
  • IdeaMaven换源(保姆级别演示)
    (1)在Idea中找到settings这个按钮,然后点击进去,如下图所示。(2)然后在里面找到Buuild,Execution,Deployment,按钮(3)点击里面的BuildTools找到Maven,将Usersettingsfile这个框后面的Override勾选上(4)点击OK(5)再在这个Usersettingsfile里面的路径.m下现新建先新建一个文本文档,命名为setti......
  • Codeforces Global Round 11 A. Avoiding Zero
    给一个大小为\(n\)的数组\(a_1,a_2,\cdots,a_n\)。你需要构造一个大小为\(n\)的数组\(b\)且满足以下条件:数组\(b\)是数组\(a\)的冲排列对于\(\forallk=1,2,\cdots,n\),\(\sum_{i=1}^{k}b_i\neq0\)。输出任意一组构造,或者回答不可能。若\(\sum_{i......
  • IDEA_多窗口_二叉树目录
    IDEAIDEA打开两个项目File——>Open/OpenRecent——>选择项目是替换目前正打开的项目窗口-ThisWindow/保留目前已打开的项目,重新打开一个新的窗口-NewWindowIDEA文件夹分支显示多个空文件夹创建时,内无文件的目录会叠加一起,点击设置按钮、TreeAppearance......
  • Exception in thread "main" java.security.InvalidKeyException: Wrong key size问题
    问题描述在Java里面使用DES加密算法,然后就爆出这个错误:问题解决换用了另外一种加密解密的函数:SecretKeySpec;即将原来的这种:换成了这种:我是觉得使用DES加密算法时,它一直显示key的字节长度不对,就想着换一种表述方式,又看到了别的友友的经验分享,就换成这样试了试(直接放进mai......
  • Android 的ViewBinding实现的原理
    AndroidViewBinding是一种用于替代传统的findViewById和findViewById的视图绑定方法。它允许你以类型安全的方式访问应用布局中的视图元素,而无需手动查找它们。ViewBinding的实现原理如下:布局文件解析:在编译期间,AndroidGradle插件会扫描项目中的布局文件(XML文件),并为每个......
  • 日本开始大规模补助台积电、Rapidus、英特尔等巨头 | 百能云芯
    日本经济新闻报导,台积电考虑在熊本县菊阳町附近兴建日本二厂,这一计划在日本政府看来十分重要,以至于他们准备提供高达约9000亿日元的补助资金。在这一巨额资金的支持下,台积电计划在熊本县菊阳町附近建设日本二厂,预计从2027年开始生产6纳米芯片。台积电是半导体制造领域......