首页
学习
活动
专区
工具
TVP
发布
精选内容/技术社群/优惠产品,尽在小程序
立即前往

根据id快速提取fastq序列

根据fastq序列的id,从原始fastq中提取序列这个操作,应该是大家在处理序列文件的过程中经常遇到的。如果大家用过Biopython,应该知道Bio模块在做fastq这些文件的处理时非常方便。...还是举个例子比较好,我从比对筛选过滤之后的bam文件中提取了第一列序列名,保存为id.name文件,想根据这个id文件从原始的fastq文件(单端)raw.fastq中把序列提出来。...这里id.name中id数目42万左右,raw.fastq序列数1000万左右: $ wc -l id.name426648 id.name$ wc -l raw.fastq 41867248...include=t 这里很多参数的意义都很明了,include=t是提取id.name中的序列,include=f是提取非id.name中的序列,这里我们应该用t。...filterbyname.sh还可以操作双端fastq,具体的可以用-h查看参数,作者写的很清楚呢。

3.1K30
您找到你想要的搜索结果了吗?
是的
没有找到

Android根据资源名获取资源ID

接触过Android开发的同学们都知道在Android中访问程序资源基本都是通过资源ID来访问。这样开发起来很简单,并且可以不去考虑各种分辨率,语言等不同资源显式指定。...痛点 但是,有时候也会有一些问题,比如我们根据服务器端的值取图片,但是服务器端绝对不会返回给我们的是资源id,最多是一种和文件名相关联的值,操作资源少的时候,可以维护一个容器进行值与资源ID的映射,但是多的话...便捷的方法 在这种情况下,使用文件名来得到资源ID显得事半功倍。 通过调用Resources的getIdentifier可以很轻松地得到资源ID。...完整的资源名为package:type/entry,如果资源名这个参数有完整地指定,后面的defType和defPackage可以省略。...defType和defPackage省略时,需要将其设置成null 注意这个方法不提倡,因为直接通过资源ID访问资源会更加效率高 如果资源没有找到,返回0,在Android资源ID中0不是合法的资源ID

3.4K10

Js排序算法_js 排序算法

它的时间复杂度也是 O(nlogn),但它在时间复杂度为 O(nlogn) 级的几种排序算法中,大多数情况下效率更高,所以快速排序的应用非常广泛。...注意: 快速排序不一定是最快的排序方法,这取决于需要排序的数据结构、数据量。不过,大多数情况下,面试官和工作场所用它的概率也是相对较高的,所以我们应该花时间把它学透彻。...当左、右两个部分各数据排序完成后,整个数组的排序也就完成了。 接下来通过一个例子理解这些步骤。假设有一个含有未排序元素 [7, -2, 4, 1, 6, 5, 0, -4, 2] 的数组。...空间复杂度在快速排序中平均也是O(log2n))。 从空间性能上看,尽管快速排序只需要一个元素的辅助空间,但快速排序需要一个栈空间来实现递归。...最好的情况下,即快速排序的每一趟排序都将元素序列均匀地分割成长度相近的两个子表,所需栈的最大深度为log(n+1);但最坏的情况下,栈的最大深度为n。这样,快速排序的空间复杂度为O(log2n))。

25.2K20

ORACLE:根据id查询所有子孙数据,或者根据id查询所有父数据(start with connect by prior)

. ---> 最常见的例子就是省市区一体表,就是通过id、pid、level来进行控制,从而一张表来存储数据.我们进行拿数据的时候,不用再连表拿取,直接通过(start with connect by...INSERT INTO REGION VALUES ('11', '绍兴市', '3', '2'); INSERT INTO REGION VALUES ('12', '西湖区', '3', '3'); 三、根据...id查询所有的子数据 需求:我输入山东省的id,会把山东省及下面的市区都查询出来 select * from REGION start with id = 2 connect by prior id...= pid -- prior 右边是子级id,就往子级的方向查询 ORDER BY id; 结果展示 四、根据id查询所有的父数据 需求:我输入黄岛区的id,会把黄岛区及其所在的市省国查询出来 select...* from REGION start with id = 8 connect by prior pid = id -- prior 右边是父id,就往父级的方向查询 ORDER BY id; 结果展示

2.1K10

基因类型注释根据基因ID就好了

生物信息学数据库种类繁多,其中基因ID是很多人比较困惑的,尤其是很多产品居然还不是基因ID的问题,比如表达芯片是探针,所以我策划了一系列ID转换教程,见文末!...我的包里面有一个函数大家比较感兴趣,就是为什么可以根据基因ID拿到其染色体坐标呢?而且还可以得到其基因类型。...(IDs, ID_type) annoGene(IDs, ID_type,out_file ='tmp.html') annoGene(IDs, ID_type,out_file ='tmp.csv')...配合着详细的介绍: 第三个万能芯片探针ID注释平台R包 第二个万能芯片探针ID注释平台R包 第一个万能芯片探针ID注释平台R包 GEO数据库中国区镜像横空出世 因为这些包暂时托管在GitHub平台,但是非常多的朋友访问...需要注意的是,这个函数的type参数,其实是有3个选择,这里我演示的是选择soft这个来源的基因注释信息。 并不是所有的平台都是有soft注释,也不是所有的平台都被我的这个工具囊括哦。

1.5K30

nginx配置根据参数转发

因此,设计解决方案为:由程序将需要跳转的完整url作为参数pic_url的值,加入整个url内,所以需要在代理中实现中转,将请求转发给pic_url url: 原: https://n1-test.xxx.com...,比如arg_pic_url可以获取到连接中pic_url这个参数的值 3....$query_string 可以获取路径之后包含所有参数kv的string 4.但是如果要获取的参数本身也是一个url,且该url包含有自带的参数,例如在这个例子里,直接使用$pic_url只能捕获http...SOURCE=EMTM这一段值,它后方的&之后的参数是无法捕捉到。...http{}部分添加一行DNS解析即可,注意,要写在nginx配置的http{}内: resolver 8.8.8.8 ipv6=off; 再次尝试,图片可以正常访问,检查浏览器控制台可以看到各项url参数均正常携带

6.6K20
领券