如何在bam文件中查找一段碱基序列？

已知一段碱基序列，想要在bam文件中查看包含这段序列的所有read（类似于在word文件中查找文本内容），应当如何操作？
目前尝试过的命令：

$ samtools view -h ref.c2.bam |less /AATGATACGGCGACCACCGAGATCT
$ /CTGCGACGGCTGCT: No such file or directory

报错，显示查不到相应内容，是因为bam文件不支持less的查找功能吗？有没有别的方法？

写回答
好问题 0 提建议
追加酬金
关注问题
分享
邀请回答
编辑收藏删除
收藏举报

报告相同问题？

关注问题

RNA-seq在subread比对时提示不能找到gene identifie 同时输出的txt文件格式混乱 linux
2022-12-24 22:13

回答 3 已采纳 1、警告信息 "Unable to find gene identifier for a feature" 表示 Subread 在处理 GTF 文件时无法找到 gene_id 属性，也就是说你使用的
怎么把linux下一个路径下的几十个数据文件分批运行？求教?
2016-10-02 14:20

回答 1 已采纳 http://blog.csdn.net/zhouzhaoxiong1227/article/details/46236117
Trooper of Bam
2017-08-05 09:09

回答 1 已采纳 https://www.baidu.com/link?url=Z1-DZbxkRHAWOCdIZgXpyRSGrF1es1wxYn7b1t0D1SW8ZWcVS5otZS4p2XVbnf_5Y9wqt
linux 查看vcf文件,Linux生信练习4-vcf
2021-05-10 05:06

刘万祥ExcelPro的博客原始数据准备cd biosoft/bowtie2/bowtie2-2.4.1-linux-x86_64/example/reads/../../bowtie2 -x ../index/lambda_virus -1 reads_1.fq -2 reads_2.fq | samtools sort -@ 5 -o tmp.bam -bcftools mpileup -f ../...
Troopers of Bam
2017-07-07 07:41

回答 1 已采纳 http://poj.org/problem?id=1531
尝试写入文件时索引超出范围
2018-02-14 04:22

回答 1 已采纳 The problem seems to be coming from these two lines inside your ingest function: n := strings.Tri
Golang从Postgresql函数获取值，该函数返回一个int但也接收一个字符串数组 postgresql
2018-08-13 18:02

回答 1 已采纳 As the error says: "message": "sql: Scan error on column index 0: converting driver.Value ty
Linux系统下RNA-seq分析（利用bowtie2比对）
2024-03-12 01:28

Luminosite Eternelle的博客 3.bowtie2-build 是用来建立索引的命令，/home/lumino/TEST1/TEST1.1/tair10.fa 是你的基因组 FASTA 文件的路径，tair10 是你希望为索引文件使用的前缀。综上所述，从提供的数据来看，这些数据的质量在映射率和读取...
几种人工神经网络之间的异同问题神经网络
2023-04-21 21:21

回答 1 已采纳基于Monster 组和GPT的调写： MLP 网（多层感知器网络）是一种前向反馈神经网络，由多个全连接的神经元层组成。每一层的神经元接收上一层的输出作为输入，并将它们线性组合并通过一个激活函数产生输
bash脚本，谁帮忙能解释一下第三行是什么意思？谢谢
2016-10-03 15:12

回答 3 已采纳 ``` 读取/路径/*.bam里的每个bam个文件，每个文件名作为i 对i按照/分段，取第8个，然后按照.分割，取第一个。如果你的路径是8层，那么这段代码提取的是文件的文件名（不包括后缀
更新Symfony2中的安全标识（ACL） php symfony
2015-01-06 13:25

回答 1 已采纳 Follow this Link: https://github.com/symfony/symfony/commit/8d39213f4cca19466f84a5656a199eee98602
如何对基因组序列进行注释
2018-09-06 22:48

weixin_33770878的博客基因组组装完成后，或者是完成了草图，就不可避免遇到一个问题，需要对基因组序列进行注释。注释之前首先得构建基因模型，有三种策略：从头注释(de novo prediction)：通过已有的概率模型来预测基因结构，在预测...
PHP只有Hit Counter？ php
2013-09-23 12:45

回答 3 已采纳 You can use the following as a basic hit counter: <?php $counter = file_get_contents('./fi
使用ArchR分析单细胞ATAC-seq数据(第一章)
2020-05-25 10:11

徐洲更hoptop的博客 "fragment"是ATAC-seq实验中的一个重要概念，它指的是通过Tn5转座酶对DNA分子进行酶切，然后经由双端测序得到的序列。 fragment产生过程我们会根据Tn5插入导致的偏移从read比对得到的位置推断出fragment的...
生物信息学软件_高通量测序技术|生信的发展，常用数据格式及分析软件
2020-12-03 21:17

weixin_39819152的博客 WeDiscover发起基因检测技术经典书籍共读活动，第一季共读李金明教授的《高通量测序技术》，今天解读第四章前两节关于生物信息学的发展，常用数据存储格式及分析软件。随着高通量测序技术的快速发展，其已由实验室...
【生物信息学】外泌体miRNA成长训练营(14天）
2022-07-15 19:50

道法自然功到自成的博客 Dicer酶的作用是在细胞质中，将Pre-miRNA剪切成大约含有22个核苷酸的双链miRNA结构，也就是miRNA和miRNA*的双链结构，Dicer酶是RNA酶3家族成员，在进化上高度保守，序列分析表明，各物种的Dicer都具有相似的结构域，...
单倍型定相软件Haplotype phasing
2020-12-20 10:15

wangchuang2017的博客 Haplotype phasing software Share: haplotype phasing software The Eagle software estimates haplotype phase either using a phased reference panel or within a genotyped cohort. Haplotype-specific ...
没有解决我的问题, 去提问

问题事件

关注

码龄粉丝数原力等级 --

被采纳

被点赞

采纳率
系统已结题 10月2日
关注

码龄粉丝数原力等级 --

被采纳

被点赞

采纳率
创建了问题 9月24日

悬赏问题

¥15 在matlab中Application Compiler后的软件无法打开
¥15 想问一下STM32创建工程模板时遇到得问题
¥15 Fiddler抓包443
¥20 Qt Quick Android 项目报错及显示问题
¥15 而且都没有 OpenCVConfig.cmake文件我是不是需要安装opencv，如何解决？
¥15 oracleBIEE analytics
¥15 H.264选择性加密例程
¥50 windows的SFTP服务器如何能批量同步用户信息？
¥15 centos7.9升级python3.0的问题
¥15 安装CentOS6时卡住

如何在bam文件中查找一段碱基序列？

0条回答 默认 最新

问题事件

悬赏问题

0条回答默认最新