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

登录

该问题已被锁定！

2: 关注

1157: 浏览

使用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 人关注

相关问题

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

关于tophat2 + cufflinks 得到的差异基因和featureCounts差异基因数目相差太大原因: 950 浏览 2 关注 2 回答 0 评论

使用MAS5.0处理过的芯片数据，能不能直接取对数后用limma包找差异表达基因，为什么不能？: 867 浏览 2 关注 1 回答 0 评论

使用bwa index 时遇到问题: 853 浏览 2 关注 1 回答 2 评论

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

meme寻找大鼠的motif: 925 浏览 2 关注 2 回答 0 评论

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

ASR祖先序列重建，最后一步使用PAML时出现一些问题: 1311 浏览 1 关注 3 回答 0 评论

无对照的转录组数据如何寻找差异表达基因: 1120 浏览 2 关注 2 回答 0 评论

不同样品中寻找特异的OTU: 1014 浏览 2 关注 2 回答 0 评论

推荐内容

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

问群体内同源基因的所有变异: 716 浏览 2 关注 1 回答 0 评论

问 rna-seq数据校正: 828 浏览 2 关注 2 回答 0 评论

问 IOBR包输入基因表达矩阵要求: 754 浏览 2 关注 1 回答 0 评论

问 L1插入位点的百分比密度分布: 770 浏览 2 关注 1 回答 0 评论

问安装R包sf报错: 1287 浏览 2 关注 2 回答 0 评论

问 how to speed up following code: 819 浏览 2 关注 1 回答 0 评论

问 sc-ATAC数据质控: 1130 浏览 2 关注 1 回答 0 评论

问 atac重复样品可视化: 872 浏览 2 关注 1 回答 0 评论

问利用GATK的HaplotypeCaller生成GVCF一直报错怎么回事: 1112 浏览 2 关注 2 回答 0 评论

All Rights Reserved Powered BY WeCenter V4.1.0 © 2024