当前位置: 首页 > news >正文

【minimap2】一定要注意组合参数

当我需要minimap2在输出sam文件中包含secondary alignment时,我认为默认的输出开关就应该是开着的,因此没有设置--secondary=yes,使用以下参数:

minimap2 -ax sr -t $threads ${INDEX} ${fastq_dictory}/${prefix}_qc2_1.fq.gz ${fastq_dictory}/${prefix}_qc2_2.fq.gz > align/${prefix}.${align_software}.${parameters}.sam

怎么看也没有secondary alignment,我以为是没有sort到一起,但是考虑到两点:

  1. 整体name-sort会导致占用内存很大,因为你总要从头看到尾才能知道有没有漏掉的同名reads

  2. minimap2在比对时默认就应该是按照name-sort做完一个做下一个,不存在在很远的地方重新再比一次的必要性

因此并没有直接使用samtools sort -n,那么我就设置上--secondary=yes,此时输出文件大了很多倍说明此时secondary alignment有被输出出来,经过less检查也确认了这一点,我起初还以为是-p的原因,直到打开man页面看到了-sr后面醒目的粗体括号我才想起来原来这里面把原设定改了,非常无语于自己无能的记忆力。

出名还是要趁早啊,老了记性就差了,差很多。

http://www.gsyq.cn/news/57599.html

相关文章:

  • 3-数据库
  • 1-计算机网络
  • html空间如何添加图片
  • Trick——字符串
  • BLOG-1-电梯调度算法
  • 线圈公司质量排名,线圈生产厂家性价比推荐,变压器公司行业排名,变压器生产厂家客户推荐。滤波器公司技术排名,开关变压器公司性价比榜单
  • 2022年春季研究资助计划征集技术提案
  • 2025年集成房屋设计公司十大排名,岗亭加工厂家十大排行榜,专业岗亭定制工厂怎么选?彩钢移动厕所厂家推荐。
  • 图像的颜色模式
  • noip13
  • host with linux
  • 深入解析:算法基础篇(8)贪心算法
  • 20251122
  • AI填补核聚变传感器数据空白技术解析
  • 电梯调度程序分析
  • 帮同学签了个到,我发现竟然能盗光他所有账号
  • hive sql开发难不难
  • KingbaseES电科金仓数据库SQL调优 - 实践
  • 深入解析:职场发展—如何避雷垃圾公司
  • ssh隧道
  • onenote复制纯文本的方法
  • 测试经验分享,登录功能+购物车+限时秒杀(测试点)汇总 - 详解
  • SGD 、Momentum 、Adam
  • P22_损失函数与反向传播
  • 团队作业3-需求改进系统设计
  • C#扩展成员全面解析:从方法到属性的演进
  • 白银市一对一培训机构推荐,2026年最新课外辅导全面测评口碑排名榜
  • 天水市一对一培训机构推荐,2026最新课外辅导机构口碑深度测评排名榜
  • Spring BeanFactoryPostProcessor 接口
  • 详细介绍:第七篇:匹配篇 | 怎么像做产品一样,为每个岗位“定制”你的简历?