reverse - 使用 samtools faidx 提取反向 BLAST 匹配
问题描述
我正在使用 samtools faidx 提取与 BLAST 输出文件给出的范围匹配的序列(格式为表格 -outfmt 6)。不幸的是,有正向和反向的 BLAST 匹配。samtools faidx 可以毫无问题地处理前向范围。反向范围导致 samtools 给出错误消息。
语法是:
samtools faidx file.faa "$scaffold":"$start"-"$end"
在反向匹配的情况下(“$start”>“$end”),错误信息是:
>$Scaffold:$start-$end/rc
[faidx] Zero length sequence: $Scaffold:$start-$end
有谁知道 samtools 是否可以处理带有特定标志的反向输入?我还没有找到适合这种情况的脚本,是的,这在生物科学中应该很常见!
解决方案
我没有在 samtools 中看到云处理反向输入的任何特定“标签”。
尝试:
samtools faidx file.faa "$scaffold": "$end"-"$start"
它仍然是您想要的区域。如果您希望提取的序列反向(和互补)方向,请尝试一些在线工具。
推荐阅读
- ios - 在 Swift 中使用按钮添加变量和返回数字
- java - SSL VPN 连接到使用 Java 的 Fortinet 防火墙
- javascript - 如何在 JSON 对象中连接双引号字符串
- mobile - Arduino MKR NB 1500 - 信号丢失后重新连接
- python - 如何修复国际象棋词典验证器中的索引错误
- javascript - 如何在带有 next.js 的数组映射函数中使用 if 语句
- ios - 如何检测和获取 UITextField 内容 (Swift)
- cookies - 附加到“Set-Cookie”属性的错误⚠️
- node.js - SyntaxError:在 mocha 和 chai http 测试中位置 0 的 JSON 中出现意外的令牌 A
- javascript - 初学者 javascript 问题:箭头函数作为可选参数