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

为什么 cutadapt 不过滤任何读取?

如何解决为什么 cutadapt 不过滤任何读取?

编码新手;背景是微生物生物信息学。我正在尝试运行 cutadapt 来过滤我的碱基对读取,但它不会过滤任何内容并且文件最终为空。我已经尝试使用 -h for cutadapt 上的建议更改脚本数十次,但没有任何效果。我的读数是 250 个碱基对。最新的 cutadapt 在我的其他文件上运行良好,但不适用于这些 Illumina 1.9 MySeq 250 bp 读数。

在编辑和运行脚本的某个时刻,它说读取未正确配对,尽管它没有理由这样做。那么为什么它不过滤任何读取,为什么它说我的读取配对不正确?我该如何解决?提前致谢!

python /project/geobiology/sw_local/cutadapt/bin/cutadapt -m 100 -M 250 -q 10 --quality-base=64 -f fastq -o ${sample}_1_trimmed.fq.gz -p ${sample}_2_trimmed.fq.gz ${sample}_1.fastq ${sample}_2.fastq >> cutadapt_primer_trimming_stats.txt 2>&1

版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 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”。这是什么意思?