微信公众号搜"智元新知"关注
微信扫一扫可直接关注哦!

kneaddata 和 bowtie2 错误:跳过读取

如何解决kneaddata 和 bowtie2 错误:跳过读取

我在 bowtie2 中使用 kneaddata。 使用带有修剪选项的认参数 (trimmomatic_options = MINLEN:12 SLIDINGWINDOW:4:15) 我的日志文件显示以下类型的警告重复多次:

Warning: skipping read '####' because length (1) <= # seed mismatches (0)
Warning: skipping read '####' because it was < 2 characters long

日志文件的结尾有以下信息:

7995997 reads; of these:
  7995997 (100.00%) were unpaired; of these:
    7506309 (93.88%) aligned 0 times
    333734 (4.17%) aligned exactly 1 time
    155954 (1.95%) aligned >1 times
6.12% overall alignment rate
pair1_aligned : 216036
pair2_aligned : 216036
orphan1_unaligned : 47096
orphan2_unaligned : 14117
orphan2_aligned : 43260
pair2_unaligned : 3722548
pair1_unaligned : 3722548
orphan1_aligned : 14356

我有以下问题:

  1. 这是否意味着我的读取有问题?如果是这样,如何找出答案?
  2. 我应该避免使用trimmomatic吗?您认为这会导致读取在这里被过滤掉吗?
  3. 在接下来的步骤中我会丢失这些读数吗?

版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。

相关推荐


Selenium Web驱动程序和Java。元素在(x,y)点处不可单击。其他元素将获得点击?
Python-如何使用点“。” 访问字典成员?
Java 字符串是不可变的。到底是什么意思?
Java中的“ final”关键字如何工作?(我仍然可以修改对象。)
“loop:”在Java代码中。这是什么,为什么要编译?
java.lang.ClassNotFoundException:sun.jdbc.odbc.JdbcOdbcDriver发生异常。为什么?
这是用Java进行XML解析的最佳库。
Java的PriorityQueue的内置迭代器不会以任何特定顺序遍历数据结构。为什么?
如何在Java中聆听按键时移动图像。
Java“Program to an interface”。这是什么意思?