当前位置: 首页 > 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.jsqmd.com/news/47842/

相关文章:

  • 3-数据库
  • 4-java
  • 重构高阶智驾:天瞳威视以国产芯片,解锁Robotaxi平民化路径 - 实践
  • 1-计算机网络
  • 实用指南:MCU定点计算深度解析:原理、技巧与实现
  • 2-操作系统
  • html空间如何添加图片
  • html空间可以设置边框吗
  • 【SpringBoot从初学者到专家的成长23】利用SpringBoot构建高效的Web应用-拥抱你的第一个SpringBoot项目
  • PyCharm,Run Configurations,Python interpreter下拉框会显示哪些地方的python.exe
  • Deepseek大模型结合Chrome搜索爬取2025AI投资趋势数据 - 指南
  • 聚焦GeoAI!6本遥感地学好刊全解析,助你精准投稿
  • call 与 delegatecall - all-in
  • 20232417 2025-2026-1 《网络与系统攻防技术》实验六实验报告
  • Trick——字符串
  • 2022年春季研究资助计划征集技术提案
  • BLOG-1-电梯调度算法
  • 线圈公司质量排名,线圈生产厂家性价比推荐,变压器公司行业排名,变压器生产厂家客户推荐。滤波器公司技术排名,开关变压器公司性价比榜单
  • IDA之修复结构体
  • 2025年专业集成房屋安装公司排名,集成活动房十大靠谱厂家推荐榜,石家庄集成房屋公司十大排行:义广达彩钢口碑推荐榜,河北集成活动房厂家十大推荐榜
  • 真的好怪题解:P14314 [Aboi Round 2] Oneshot
  • ElasticSearch索引库操作 - 努力-
  • 2025年集成房屋设计公司十大排名,岗亭加工厂家十大排行榜,专业岗亭定制工厂怎么选?彩钢移动厕所厂家推荐。
  • 礼盒拖车公司推荐,礼盒拖车定制公司排行榜,礼盒拖车厂家口碑推荐,礼盒拖车生产厂家-航利通达
  • 360笔试
  • 图像的颜色模式
  • .NET+AI | MEAI | Function Caling 实操(4)
  • 高频变压器公司口碑榜单,电感公司技术排名,电感厂家交付效率排名,磁性元器件公司客户推荐,电感器公司产能排名,线圈公司行业排名-汉翔电子
  • MinIo介绍 - 努力-
  • BLOG1