获取包含在bam、sam或fastq中的fast5文件的路径。

fast5seek的Python项目详细描述


#快速5搜索
[![pypi status](https://img.shields.io/pypi/v/fast5seek.svg)(https://pypi.python.org/pypi/fast5seek)
[![构建状态](https://travis-ci.org/mbhall88/fast5seek.svg?branch=master)(https://travis ci.org/mbhall88/fast5seek)
[![Github许可证](https://img.shields.io/github/license/mbhall88/fast5seek.svg)](https://github.com/mbhall88/fast5seek/blob/master/license)
[![推特关注](https://img.shields.io/twitter/follow/mbhall88.svg?style=social&logo=twitter&label=follow)(https://twitter.com/mbhall88)


此程序接受一个(或多个)fast5文件目录以及任意数量的fastq、sam或bam文件。输出是所有fast5文件的完整路径,fastq、bam或sam文件也在提供的fast5目录中。





`` bash
pip3 install fast5seek
``````


`` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` ` `/>
` in.fastq in.bam in.sam`文件,并从每个记录中提取读取id。然后,它会遍历`/path/to/fast5s'下的所有
fast5文件,并检查它们的读取ID是否在`<;in.fastq in.bam in.sam>;中的读取ID集中。如果是,则文件的
路径将写入其自己的行in`out.txt`.


如果未给出输出(`-o`),则将输出写入“stdout”。










\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\>从给定目录输出所有fast5文件的路径包含在fastq或bam/sam文件中。

-帮助显示此帮助消息并退出
-I FAST5YDIR [ FAST5XDIR…],FAST5XIDIR FAST5YDI[FAST5YDIR…] BR/>您想查询的FACT5文件目录。程序将递归地遍历子目录。
-r reference[引用…],--reference reference[引用…]
fastq或bam/sam文件。
-o output,--output output
要写入fast5路径的文件名。如果未输入
,它将写入到stdout的路径。
-m,--mapped only extract read id for mapped reads in bam/sam
文件。
--日志级别{0,1,2,3,4,5}
日志级别。0表示无,5表示调试。
默认值为4,它将报告信息、警告、错误、
和关键信息。
--无进度条不显示进度条。
````

不需要合并bam fastq sam文件。

fast5seek-i/myfast5/dir/1//other/fast5/dir/2/-r reads.sorted.bam reads2.bam



如果通配符位于同一目录中,则可以使用
其他程序,您可以执行类似于

fast5seek-i/path/to/fast5s/-r in.fastq xargs cp-t subset-dir/

打开
`fast5`文件,我们建议您首先将
`fast5seek`的输出保存到一个文件中保存,然后进行如下分析:


mkdir subset_dir-dir/
fast5seek-i/path/t o/fasts5s/-r in.fastq-o mapped_reads.txt
cat mapped_reads.txt xargs cp-t subset_subset_dir/




是的如果程序有任何问题,请[在上面打开一个问题](https://github.com/mbhall88/fast5seek/issues)。


欢迎加入QQ群-->: 979659372 Python中文网_新手群

推荐PyPI第三方库


热门话题
带有servlet、jsp和数据库的应用程序的java结构   java Admob本机Express广告   java无法在会话超时时获取会话范围的bean   java的“xmlns=”http://www.intellij.com/uidesigner/form/“`由Intellij以.xml格式标记为红色   迭代器在Java中重写Iterables<Obj>的正确方法   java如何刷新我的默认表模型   无效URL上的java PrettyFaces重定向循环   java如何将JSON数组传递给SpringMVC?   java如何从setOnClickListener获取数据并将其传递给倒计时程序   从java运行grep命令后bash获取错误   java[Ljava.lang.String;@5d79a22d结果   java错误p:ajax更新   请在这段代码中解释java instanceof关键字