三代测序PacBioONT reads过滤和修剪-Chooper

chopper简介

chopper是NanoFilt和NanoLyse的Rust语言版本,适用于长reads测序(如PacBio和纳米孔测序ONT)的过滤和修剪fastq文件。

chopper相对于python编写的NanoFilt和NanoLyse,运行输出相同结果的时间更短,且NanoFilt和NanoLyse现已停止更新,选择使用chopper会有更高的效率。

github地址: https://github.com/wdecoster/chopper

chopper

chopper安装

# conda安装
conda install -c bioconda chopper -y# linux发行版zip包安装
mkdir chopper && cd chopper 
wget https://github.com/wdecoster/chopper/releases/download/v0.8.0/chopper-linux.zip
unzip chopper-linux.zip
# 添加执行权限
chmod +x chopper
# 加入环境变量
echo "export PATH=/software/chopper/:$PATH" >> ~/.bashrc
source ~/.bashrcchopper  -h

运行示例

# 对于fastq.gz文件作为输入+管道+过滤质量低于10、长度低于1000bp的序列,gzip压缩为.gz文件
gunzip -c sample.fastq.gz | chopper -q 10 -l 1000 \
| gzip > sample.filtered.fastq.gz# fastq文件作为输入
chopper -q 10 -l 500 -i sample.fastq > sample.filtered.fastq.gz# fastq.gz文件作为输入,并对输出文件进行压缩
chopper -q 10 -l 500 -i sample.fastq.gz \
| gzip > sample.filtered.fastq.gz

参数说明:

-l ,–minlength :过滤掉小于指定长度的序列
-q , --quality :过滤掉低于指定质量序列
-i, input: 输入文件,默认从stdin读取
–maxlength :过滤掉超过此长度的序列, 默认2147483647
–threads:线程数,默认4
–contam: 参考fasta文件检查潜在的污染,默认None
–minGC:过滤掉低于指定GC含量此百分比的序列,默认0
–maxGC:过滤掉大于指定GC含量此百分比的序列, 默认1.0
–headcrop:从头部切掉n bp,默认0
–tailcrop:尾部切掉n bp,默认0

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/865271.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

思维,CF 739A - Alyona and mex

一、题目 1、题目描述 2、输入输出 2.1输入 2.2输出 3、原题链接 739A - Alyona and mex 二、解题报告 1、思路分析 我们考虑区间mex运算的值最大也就是区间长度,所以我们最大值的上界就是所有区间中的最小长度,假如记为mi 我们一定可以构造出答案…

开启IT世界的第一步:高考新生的暑期学习指南

目录 前言 了解IT领域 学习编程语言 实践项目 学习资源 阅读专业书籍 培养良好的学习习惯 结语 最后 - 投票 前言 七月的钟声敲响,各省的高考分数已揭晓,意味着一段紧张而又充满奋斗的旅程画上了句号。然而,高考的结束并不意味…

Java 泛型:上界通配符和下界通配符的用途和限制

Java 泛型&#xff1a;上界通配符和下界通配符的用途和限制 Java泛型中的上界通配符&#xff08;<? extends T>&#xff09;和下界通配符&#xff08;<? super T>&#xff09;是两个强大的工具&#xff0c;它们在不同的场景中各自发挥重要作用。理解它们的用途和…

zabbix 配置钉钉告警

1.申请一个钉钉企业版 2.群内申请一个机器人 下载电脑版钉钉&#xff0c;登录后&#xff0c;在要接收群消息的群里&#xff0c;点击右上角设置图标&#xff0c;下滑找到机器人&#xff0c;添加一个机器人&#xff0c;保存机器人的webhook地址 保存这里的加签字符串 保存这里的…

通信软件开发之业务知识:PON口割接什么意思?

一 PON口割接&#xff08;原创总结&#xff09; 在通信领域&#xff0c;PON口割接指的是对无源光网络&#xff08;Passive Optical Network&#xff0c;PON&#xff09;端口进行的切换或调整操作。简单来说&#xff0c;就是对光纤网络中的某个端口进行重新连接或重新分配&…

153. 寻找旋转排序数组中的最小值(中等)

153. 寻找旋转排序数组中的最小值 1. 题目描述2.详细题解3.代码实现3.1 Python3.2 Java 1. 题目描述 题目中转&#xff1a;153. 寻找旋转排序数组中的最小值 2.详细题解 如果不考虑 O ( l o g n ) O(log n) O(logn)的时间复杂度&#xff0c;直接 O ( n ) O(n) O(n)时间复杂…

【Hugging Face全面拥抱LangChain:全新官方合作包】

文末有福利&#xff01; ❝ 最近Hugging Face官宣发布langchain_huggingface&#xff0c;这是一个由 Hugging Face 和 LangChain 共同维护的 LangChain 合作伙伴包。这个新的 Python 包旨在将 Hugging Face 最新功能引入 LangChain 并保持同步。 通过Hugging Face官方包的加持&…

【Linux 命令】文件比较 diff

diff 命令是 Unix 和类 Unix 系统&#xff08;如 Linux 和 macOS&#xff09;中用于比较文件内容差异的一个非常有用的命令行工具。它可以逐行比较两个文件的内容&#xff0c;并输出它们之间的差异。这些差异通常以行为单位显示&#xff0c;并且会标记出哪些行是唯一的、添加的…

d嘤嘤不想求异或喵(牛客周赛49)

题意&#xff1a; 嘤嘤有两个整数 l,r&#xff0c;她想知道区间 [l,r] 所有整数的异或和是多少. 分析&#xff1a; 样例1只有一个数输出1 样例2 1^201^10113 样例3 1^2^301^10^1111^11000 #include<bits/stdc.h> using namespace std; typedef long long ll; ll f(l…

14-15 为什么我们现在对阅读如此难以接受

写出来感觉很奇怪&#xff0c;但最近我感觉自己失去了阅读能力。长篇文本对我来说尤其具有挑战性。句子很难读完。更别提章节了。章节有很多段落&#xff0c;而段落又由许多句子组成。 啊。 即使在极少数情况下&#xff0c;我读完了一章&#xff0c;下一页上已经有另一章等着…

Unity3D滚动列表的编码与实现详解

Unity3D作为一款强大的游戏开发引擎&#xff0c;提供了丰富的UI系统&#xff08;UGUI&#xff09;来支持开发者构建各种用户交互界面。滚动列表作为常见的UI元素之一&#xff0c;在游戏、应用程序中广泛使用&#xff0c;如角色列表、物品背包、设置菜单等。本文将详细介绍Unity…

Asterinas 驱动框架

Asterinas 驱动框架 总线初始化 在整个OS初始化的时候&#xff1a; Rust pub fn init() { arch::before_all_init(); logger::init(); #[cfg(feature “intel_tdx”)] let td_info init_tdx().unwrap(); #[cfg(feature “intel_tdx”)] early_println!( “td gpaw: {}, td a…

有sdwan可以不用专线吗?sdwan和专线的区别优势

SD-WAN&#xff08;Software-Defined Wide Area Network&#xff0c;软件定义广域网&#xff09;确实可以在很大程度上替代传统的专线连接&#xff0c;尤其是在追求成本效益和网络灵活性的场景下。SD-WAN的核心优势在于其智能化和自动化的能力&#xff0c;这使得它能够优化数据…

华清远见人工智能课程:项目优势助力,学习更高效!

在人工智能飞速发展的今天&#xff0c;学习人工智能成为新的高薪赛道。我们都知道人工智能的学习离不开项目练手&#xff0c;只有通过实际项目的操作&#xff0c;才能真正掌握人工智能的核心技能。但遗憾的是&#xff0c;很多人工智能课程只注重理论知识的传授&#xff0c;缺乏…

中介子方程四十

XXFXXuXXWXXuXXdXXrXXαXXuXpXXKXηXiXXnXXyXηXuXXrXXaXnXXαXXnXaXXrXXuXηXyXXnXXiXηXKXXpXuXXαXXrXXdXXuXWXπXXWXeXyXeXbXπXpXXNXXqXeXXrXXαXXuXpXXKXηXiXXnXXyXηXuXXrXXaXnXXαXXnXaXXrXXuXηXyXXnXXiXηXKXXpXuXXαXXrXXeXqXXNXXpXπXbXeXyXeXWXXπXWXuXXdXXrXXα…

本地文件上传Github的方法

本文仅用于个人回忆本地文件上传GitHub方法&#xff0c;使用HTTP上传而不是SSH&#xff0c;其余文件已配置完成&#xff0c;如果你和我一样&#xff0c;那么请往下看&#xff0c;这条博客不说废话。 Step1 选择需要上传的文件&#xff0c;右键选择Git Bash Here Step2 创建一…

GAMES104:04游戏引擎中的渲染系统1:游戏渲染基础-学习笔记

文章目录 概览&#xff1a;游戏引擎中的渲染系统四个课时概览 一&#xff0c;渲染管线流程二&#xff0c;了解GPUSIMD 和 SIMTGPU 架构CPU到GPU的数据传输GPU性能限制 三&#xff0c;可见性Renderable可渲染对象提高渲染效率Visibility Culling 可见性裁剪 四&#xff0c;纹理压…

电路模型和电路定律

电路---为了某种需要由某些电工设备或元件按一定方式组合起来的电流的通路 实际电路的两个作用 1.电能的传输&#xff0c;分配和转换 2.传递和处理信号 电路中的几个基本概念 激励---电源或信号源的电压或电流&#xff0c;也称为输入 响应---由激励在电路各部分产生的电流…

怎么参与场外期权?

今天期权懂带你了解怎么参与场外期权&#xff1f; 目前个人投资者暂时还不能直接参与场外个股期权&#xff0c;因为场外个股期权现在只能机构来进行交易。 所以个人投资者目前只能通过机构通道来进行操作&#xff0c;类似期权懂&#xff0c;找到期权懂经理&#xff0c;然后通…

手机删除的照片怎么找回?3个急救指南,让你重新拥有

一不小心手滑&#xff0c;手机里的几百张珍贵照片就消失了&#xff01;是不是心如刀割&#xff0c;感觉错过了几个亿&#xff1f;手机删除的照片怎么找回&#xff1f;别担心&#xff0c;小菜一碟&#xff0c;恢复照片就像打游戏一样&#xff0c;只要掌握正确的“攻略”&#xf…