Does anyone know what this output in tophat mean? How much of a problem is it?
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
I am curious too...
I am curious as well as to what this truly means and if it should be a concern. I encountered it on my latest runs in which I am using the Ensembl version of the human genome assembly. An example of the warning displayed is:
[Sun Aug 08 23:20:23 2010] Beginning TopHat run (v1.0.14)
-----------------------------------------------
[Sun Aug 08 23:20:23 2010] Preparing output location brain/
[Sun Aug 08 23:20:23 2010] Checking for Bowtie index files
[Sun Aug 08 23:20:23 2010] Checking for reference FASTA file
[Sun Aug 08 23:20:23 2010] Checking for Bowtie
Bowtie version: 0.12.5.0
[Sun Aug 08 23:20:23 2010] Checking reads
seed length: 36bp
format: fastq
quality scale: phred33 (default)
[Sun Aug 08 23:23:07 2010] Mapping reads against Homo_sapiens.GRCh37.59 with Bowtie
[Sun Aug 08 23:26:43 2010] Joining segment hits
[Sun Aug 08 23:27:09 2010] Mapping reads against Homo_sapiens.GRCh37.59 with Bowtie
[Sun Aug 08 23:31:02 2010] Joining segment hits
[Sun Aug 08 23:31:26 2010] Searching for junctions via segment mapping
[Sun Aug 08 23:46:58 2010] Retrieving sequences for splices
[Sun Aug 08 23:52:48 2010] Indexing splices
Warning: Encountered reference sequence with only gaps
Warning: Encountered reference sequence with only gaps
Warning: Encountered reference sequence with only gaps
Warning: Encountered reference sequence with only gaps
Warning: Encountered reference sequence with only gaps
Warning: Encountered reference sequence with only gaps
Warning: Encountered reference sequence with only gaps
Warning: Encountered reference sequence with only gaps
Warning: Encountered reference sequence with only gaps
Warning: Encountered reference sequence with only gaps
Warning: Encountered reference sequence with only gaps
Warning: Encountered reference sequence with only gaps
Warning: Encountered reference sequence with only gaps
Warning: Encountered reference sequence with only gaps
Warning: Encountered reference sequence with only gaps
Warning: Encountered reference sequence with only gaps
Warning: Encountered reference sequence with only gaps
Warning: Encountered reference sequence with only gaps
Warning: Encountered reference sequence with only gaps
Warning: Encountered reference sequence with only gaps
Warning: Encountered reference sequence with only gaps
Warning: Encountered reference sequence with only gaps
Warning: Encountered reference sequence with only gaps
Warning: Encountered reference sequence with only gaps
Warning: Encountered reference sequence with only gaps
Warning: Encountered reference sequence with only gaps
[Sun Aug 08 23:54:58 2010] Mapping reads against segment_juncs with Bowtie
[Sun Aug 08 23:58:35 2010] Joining segment hits
[Sun Aug 08 23:59:00 2010] Mapping reads against segment_juncs with Bowtie
[Mon Aug 09 00:02:49 2010] Joining segment hits
[Mon Aug 09 00:03:13 2010] Reporting output tracks
-----------------------------------------------
Run complete [00:45:40 elapsed]
-steve
-
Originally posted by digitonin View PostI never did find an answer to this. I didn't see any problems associated it with it either though.
First of all, sorry for bumping such a old post, but this problem still persists...
Does anyone found a reason why does this happens?
I have obtained this same warning message when using the current Drosophila melanogaster reference (r6.03; BDGP5.78)
Comment
-
Hi,
I am using bowtie 1.1.2 and I get the same message -
encountered reference sequence with gaps
My genome is not repmasked and has the following type of header
>ta_IWGSC_CSSassembly_1as_v2_44039
Does such a warning interfere with the bowtie-build in anyway?
Comment
-
This will happen whenever you have cases like this
Code:>chr1 >chr2 ACAGCTACT
Code:>chr3 NNNNNNN >chr4 ACGTAGCTGACT
Comment
Latest Articles
Collapse
-
by seqadmin
Like all molecular biology applications, next-generation sequencing (NGS) workflows require diligent quality control (QC) measures to ensure accurate and reproducible results. Proper QC begins at nucleic acid extraction and continues all the way through to data analysis. This article outlines the key QC steps in an NGS workflow, along with the commonly used tools and techniques.
Nucleic Acid Quality Control
Preparing for NGS starts with isolating the...-
Channel: Articles
02-10-2025, 01:58 PM -
-
by seqadmin
In recent years, precision medicine has become a major focus for researchers and healthcare professionals. This approach offers personalized treatment and wellness plans by utilizing insights from each person's unique biology and lifestyle to deliver more effective care. Its advancement relies on innovative technologies that enable a deeper understanding of individual variability. In a joint documentary with our colleagues at Biocompare, we examined the foundational principles of precision...-
Channel: Articles
01-27-2025, 07:46 AM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Genetic Mapping of Plasmodium knowlesi Identifies Essential Genes and Drug Resistance Mechanisms
by seqadmin
Started by seqadmin, 02-07-2025, 09:30 AM
|
0 responses
67 views
0 likes
|
Last Post
by seqadmin
02-07-2025, 09:30 AM
|
||
Started by seqadmin, 02-05-2025, 10:34 AM
|
0 responses
104 views
0 likes
|
Last Post
by seqadmin
02-05-2025, 10:34 AM
|
||
Started by seqadmin, 02-03-2025, 09:07 AM
|
0 responses
83 views
0 likes
|
Last Post
by seqadmin
02-03-2025, 09:07 AM
|
||
Started by seqadmin, 01-31-2025, 08:31 AM
|
0 responses
45 views
0 likes
|
Last Post
by seqadmin
01-31-2025, 08:31 AM
|
Comment