![](https://img-blog.csdnimg.cn/20201014180756928.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
生信
文章平均质量分 50
weixin_46530443
这个作者很懒,什么都没留下…
展开
-
使用maker进行基因组注释
所遇到的100个问题。(我使用过的maker主要是老师和师姐装的,如果自己装,最少还要遇到一百个问题)1.maker要在/dev/shm下运行,有点类似计算机的运行内存;2.如果提示某个perl脚本环境有问题,依照提示修改脚本环境变量相应行即可;3.运行new_species.pl遇到报错:设置图中的AUGUSTUS_CONFIG_PATH即可;4.注释出的雌株雄株基因数目相差过大,合并第二轮augustus和snap的模型后进行最后一步maker注释;5.使用maker2zff 遇到报错:原创 2021-12-09 10:37:04 · 2437 阅读 · 0 评论 -
bamtools拆分bam文件
bamtools拆分bam文件sRNA使用shortstack比对到基因组上的比对文件是合并的bam文件,需要按照样品将其拆分,使用bamtools进行拆分:bamtools split -in merged_alignments.bam -tag RG运行时间较久,有更好的方式再来更新。...原创 2021-12-07 14:19:19 · 1243 阅读 · 1 评论 -
linux下处理fastq格式文件
1.统计每条序列长度:awk '{if(NR%4==2){print}}' s1_1.fq |awk '{print length($1)}'原创 2021-09-15 09:45:59 · 506 阅读 · 0 评论 -
计算反转录转座子插入时间三:MEGA批量化处理
得到了成对的LTRs后,可根据两条LTR序列的不同,根据K-2模型计算K——每个核苷酸位点的平均替代数,使用MEGA进行计算。由于文件众多,需要多次处理,使用MEGA-cc进行处理。安装方法参考https://blog.csdn.net/g_r_c/article/details/22602071安装完成后,打开MEGA,点击右下角PROTOTYPE,进入模拟模式,选择序列类型(我使用Nucleotide non-coding)ALIGN——MUSCLE,使用默认参数点击确定生成比对配置文件*.mao。原创 2021-09-14 18:46:25 · 1668 阅读 · 6 评论 -
计算反转录转座子插入时间二:提取成对LTRs序列
上一篇讲了反转录转座子插入时间的计算原理,这篇主要是计算过程。 首先使用ltr_finder和ltr_harvest对基因组中的反转录转座子进行预测,之后用LTR_retriever整合预测结果,流程主要参考徐洲更先生的博文,输出结果如下:图源:https://www.jianshu.com/p/f962d5c40fdf其中*.fa.pass.list.gff3内包含有所有完整的反转录转座子,内容如下:每个block表示一个反转录转座子,第一行为整个LTR-RT的区域,4、5行分别...原创 2021-09-13 21:42:20 · 3595 阅读 · 1 评论 -
计算反转录转座子插入时间一:计算原理
转座子(Transposable elements, TEs)又称为转座原件,是一类特殊的DNA序列,可以通过转录或逆转录,在内切酶的作用下从染色体一个位置跳跃到另一个位置。根据转座方式不同,转座子可以分为三类[1]:I型转座子又称为反转录转座子(retrotransposon),在转座时,先以DNA为模板转录成一段mRNA,再以其为模板反转录成cDNA,在整合酶的作用下整合到基因组的新位置上。II型转座子又称转座子(transposon),在转座酶的作用下,从原来的位置解离下来,再重新整合到染色体上原创 2021-09-12 18:07:11 · 3198 阅读 · 1 评论