使用tophat2和bowtie1寻找环形RNA时报错 - 做专业的生信问答网站

登录

该问题已被锁定！

2: 关注

3036: 浏览

使用tophat2和bowtie1寻找环形RNA时报错

使用tophat2和bowtie1寻找环形RNA时报错Mapping left_kept_reads to genome xxxx with Bowtie [FAILED]，原始代码

For single-read sequencing

TopHat2/TopHat-Fusion

tophat2 -a 6 --microexon-search -m 2 -p 10 -G knownGene.gtf -o tophat hg19_bowtie2_index RNA_seq.fastq
bamToFastq -i tophat/unmapped.bam -fq tophat/unmapped.fastq
tophat2 -o tophat_fusion -p 15 --fusion-search --keep-fasta-order --bowtie1 --no-coverage-search hg19_bowtie1_index tophat/unmapped.fastq

好问题 0 评论收藏举报

查看全部 2 个回答

孟浩巍超级管理员用户来自于: 北京市
2023-05-23 09:39

目前只看报错不能判断问题。

chatGPT的回答有一定的参考性，我有几方面的建议。

1. 考虑是不是index有问题，就是用测试数据是否能跑通？

2. 是不是能复现这个问题，就是每次运行都报这个错误，能否定位到具体是哪个reads引发的报错？

能复现错误，之后才好进一步排查。

赞同 0 0评论

关于作者

: m_79ebf5bf0169deebb7c741e4d3cb8be7 注册会员
这家伙很懒，还没有设置简介

0: 回答

0: 文章

1: 问题

问题动态

发布时间: 2023-05-22 21:30

更新时间: 2023-05-23 09:39

关注人数: 2 人关注

相关问题

cellranger使用问题: 3353 浏览 2 关注 1 回答 0 评论

虚拟机中使用GEC进行GWAS阈值矫正: 2657 浏览 2 关注 1 回答 0 评论

怎么寻找某个基因的不同可变剪切形式？: 2291 浏览 4 关注 2 回答 0 评论

ROSE算法寻找SEs: 2079 浏览 2 关注 1 回答 0 评论

使用Trimmomatic切除接头后做FASTQC: 2500 浏览 2 关注 1 回答 0 评论

mhcflurry的使用: 2405 浏览 2 关注 1 回答 0 评论

trimmomatric使用报错: 2810 浏览 2 关注 2 回答 0 评论

使用seurat包，导出特定cluster的细胞-基因counts矩阵。初始数据命名为pbmc: 2960 浏览 2 关注 1 回答 0 评论

ChIP-seq的bigwig文件可以使用rpkm>1过滤掉一些弱的信号吗？: 2714 浏览 2 关注 2 回答 0 评论

如何使用wget命令下载“清华云盘”内的文件: 5157 浏览 2 关注 2 回答 0 评论

推荐内容

问 bowtie2使用报错，求助啦: 2650 浏览 3 关注 3 回答 0 评论

问蛋白保守序列分析: 1706 浏览 2 关注 1 回答 0 评论

问 grep 能不能指定匹配第多少列？: 2314 浏览 3 关注 2 回答 1 评论

问动态库存在但调用报错问题: 2158 浏览 2 关注 1 回答 0 评论

问链特异性文库（mRNA/lncRNA/circRNA）如何将RNA类型分开？: 2183 浏览 2 关注 1 回答 0 评论

问群体结构矫正: 2241 浏览 2 关注 1 回答 0 评论

问用R画RNA_seq中变化显著的heatmap聚类热图: 2978 浏览 4 关注 3 回答 0 评论

问 computeMatrix画图问题: 3375 浏览 1 关注 2 回答 0 评论

问信号通路特征基因集的查询？用于作ssGSEA分析。: 2135 浏览 2 关注 1 回答 0 评论

All Rights Reserved Powered BY WeCenter V4.1.0 © 2025