需求
已知某基因组序列，染色体或scaffold ID顺序不定，想要对其按数字排序。
image.png
原顺序：
image.png
想要的排序结果：
image.png
实现
使用bioawk，没有的话conda直接安装。
bioawk -c fastx '{print}' old.genome.fa | \
sort -k1,1V | awk '{print ">"$1;print $2}' >new.genome.fa
已知某基因组序列，染色体或scaffold ID顺序不定，想要对其按数字排序。
image.png
原顺序：
image.png
想要的排序结果：
image.png
使用bioawk，没有的话conda直接安装。
bioawk -c fastx '{print}' old.genome.fa | \
sort -k1,1V | awk '{print ">"$1;print $2}' >new.genome.fa
需求 已知某基因组序列，染色体或scaffold ID顺序不定，想要对其按数字排序。 原顺序： 想要的排序结果： ...
最近在做水稻基因组注释，对于注释结果进行整理过滤，还在探索中。 evm注释结果 将结果转化成蛋白文件： 提取长度大...
awk 命令使用记录 例1, 修改fasta序列ID 若后续只需要用到序列，不需要ID的信息，就可以把ID信息完...
数据： 研究物种的基因组文件：Protein.fasta 近源种的基因家族蛋白序列 :protein.fasta ...
使用HiCPro绘制全基因组互作矩阵为了图片的美观可以将挂载上的染色体先提取出来：genome.chr.fasta...
fasta文件格式：fasta 文件为一个 ID 对应一个 序列,可以是 转录本序列, 蛋白序列 ''' V350...
很简单的一个shell脚本，从UCSC里面单独下载X,Y染色体的fasta序列，写脚本从Y染色体序列里面模拟双端测...
很简单的一个shell脚本，从UCSC里面单独下载X,Y染色体的fasta序列，写脚本从Y染色体序列里面模拟双端测...
组装结果统计 下面对组装得到fasta 格式基因组序列进行长度、N50 等统计。 使用assembly-stats...
MCScanX 使用方法： 首先需要得到物种的蛋白序列以及基因组注释文件：fasta 与 gff 然后通过本地bl...
本文标题：使用bioawk对基因组fasta序列ID（染色体/scaffo
本文链接：https://www.haomeiwen.com/subject/uaefrltx.html
网友评论