Here is an answer to billstevens (many questions got intertwined):
It seems that you've already launched the tophat command (unsuccessfully) before, which led to the creation of an annotation.gtf.gff file from the annotation.gtf file. If you remove this file, the command should run smoothly.
Best,
Julien
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
You can also just build them yourself:
gtf_to_fasta /home/chris/databases/refAnnotation.gtf /media/chris/human/Sequence/hg19.fa /home/chris/databases/refAnnotationout
bowtie2-build /home/chris/databases/refAnnotationout /home/chris/databases/refAnnotationbt
Leave a comment:
-
Originally posted by Dario1984 View PostThe GTF file is fine. Yes, read the man page for ln. The target and the link parameters are separated by a space.
Code:$ man ln
Leave a comment:
-
Bowtie error on tophat
Originally posted by Dario1984 View PostYou're not allowed to put filenames in the PATH variable. You put the path to the directory which contains those programs.
I think there is a problem with the version of tophat for OSX seems that doesnt recognise bowtie2, I used this command:
rsb-33-073:tophat-2.0.7.OSX_x86_64 jaime$ ./tophat2 /Users/jaime/Desktop/bowtie2_indexes /Users/jaime/Desktop/RNAseq/inoculated_libraries
[2013-02-20 09:37:27] Beginning TopHat run (v2.0.7)
-----------------------------------------------
[2013-02-20 09:37:27] Checking for Bowtie
Bowtie 2 not found, checking for older version..
Error: Bowtie not found on this system.
I put the bowtie2 in the tophat directory and try to run from there but bowtie is not recognised by tophat
thanks
Leave a comment:
-
The GTF file is fine. Yes, read the man page for ln. The target and the link parameters are separated by a space.
Code:$ man ln
Leave a comment:
-
Thanks Dario. Did I get the soft-linking of the genes.gtf file? and is there a space between genome* and .(full stop sign). Sorry I'm a newbie...
Leave a comment:
-
Originally posted by Ajayi Oyeyemi View Post[2013-02-18 16:06:49] Checking for Bowtie index files
Error: Could not find Bowtie 2 index files (genome.*.bt2)
Instead of
Code:$ ln -s ./Bos_taurus/NCBI/Btau_4.6.1/Sequence/Bowtie2Index/genome.*.
Code:$ ln -s ./Bos_taurus/NCBI/Btau_4.6.1/Sequence/Bowtie2Index/genome* .
Leave a comment:
-
Originally posted by Dario1984 View PostThat is not true. I downloaded bosTau7.fa and had a look at the beginning of the file with R.
Code:> readLines("F:\\bosTau7.fa", n=10) [1] ">chr1" [2] "taccccactcacacttatggatagatcaactaaacagaaaattaacaagg" [3] "aaaaaaaaaAAAAAAAAAACATGTATCTATAAAGCTCACTAAAGCAAAGC" [4] "ACAGTAAAATGAGGTATGGCAATAAAAGGCGTTAAATTTATCaaaaaaaa" [5] "aaaaaaaaagtgctgaaaaattaaaaaaaaaaaaaagaaaaaaaTATAGA" [6] "GATAAAATTAAATTTTTGTATTTATTCTTGCCGTTGAAAATTATGGGCAA" [7] "CTTTTCATTTTTTCAAATTTAATATTGTGTTACTCAAAAATACCCTAAAG" [8] "TTGATCTCTGTTTAAATCTAGAACACTTTTTGTTAAGTTTATGTCGATAT" [9] "GTCTTAGCTCTTTTTTTTAATTAGACATGATATTTTCTTTCATCACATTT" [10] "TCAATATGTAGGAGAGCTAGTGATATTTGTGTAATAATTTTGTACTCGGC"
hanks Dario, let me be much more clearer. I downloaded the cow genome from the igenome in cufflinks website and I unpacked the tar.gz file. I created a link using the ln -s function to the gtf file and the genome. See below the result I got.
[ooa4@cbsum1c1b009 ooa4]$ ln -s ./Bos_taurus/NCBI/Btau_4.6.1/Sequence/Bowtie2Index/genome.*.
[ooa4@cbsum1c1b009 ooa4]$ ln -s ./Bos_taurus/NCBI/Btau_4.6.1/Annotation/Genes/genes.gtf
[ooa4@cbsum1c1b009 ooa4]$ ls -al
total 43380296
drwxrwxr-x 9 ooa4 ooa4 4096 Feb 18 16:05 .
drwxrwxrwx. 3 root root 4096 Feb 7 12:46 ..
drwxrwxr-x 3 ooa4 ooa4 4096 Feb 18 15:57 Bos_taurus
-rw-rw-r-- 1 ooa4 ooa4 27915273767 Feb 13 15:47 Bos_taurus_NCBI_Btau_4.6.1.tar.gz
-rw-rw-r-- 1 ooa4 ooa4 32225697 Feb 18 15:46 bovineee
lrwxrwxrwx 1 ooa4 ooa4 55 Feb 18 16:05 genes.gtf -> ./Bos_taurus/NCBI/Btau_4.6.1/Annotation/Genes/genes.gtf
lrwxrwxrwx 1 ooa4 ooa4 60 Feb 18 16:04 genome.*. -> ./Bos_taurus/NCBI/Btau_4.6.1/Sequence/Bowtie2Index/genome.*.
drwxrwxr-x 3 ooa4 ooa4 4096 Feb 12 16:31 index
drwxrwxr-x 3 ooa4 ooa4 4096 Feb 16 08:44 output1
drwxrwxr-x 3 ooa4 ooa4 4096 Feb 16 08:45 output1B
-rwxrwxr-x 1 ooa4 ooa4 5790 Mar 20 2012 README.txt
-rw-rw-r-- 1 ooa4 ooa4 6683898300 Feb 8 15:00 SRR594497.fastq
-rw-rw-r-- 1 ooa4 ooa4 995599407 Feb 8 11:48 SRR594497.sra
-rw-rw-r-- 1 ooa4 ooa4 7338744848 Feb 8 15:07 SRR594499.fastq
-rw-rw-r-- 1 ooa4 ooa4 1436256918 Feb 8 11:47 SRR594499.sra
drwxrwxr-x 4 ooa4 ooa4 4096 Feb 16 10:09 top_gtf
drwxr-xr-x 4 root root 4096 Feb 13 09:13 tophat_output2
drwxrwxr-x 3 ooa4 ooa4 4096 Feb 16 10:02 unique_output
-rw-rw-r-- 1 ooa4 ooa4 19337728 Feb 13 15:47 wget-log
[ooa4@cbsum1c1b009 ooa4]$ tophat2 -p 8 -G genes.gtf -o SRR594497_tophat genome SRR594497.fastq
[2013-02-18 16:06:49] Beginning TopHat run (v2.0.7)
-----------------------------------------------
[2013-02-18 16:06:49] Checking for Bowtie
Bowtie version: 2.0.6.0
[2013-02-18 16:06:49] Checking for Samtools
Samtools version: 0.1.18.0
[2013-02-18 16:06:49] Checking for Bowtie index files
Error: Could not find Bowtie 2 index files (genome.*.bt2)
I'm using the trapnell protocol and I'm following it to the letter.
Any advice?
Leave a comment:
-
You need to have the path to the bowtie2 executable file in your PATH environment variable.
Leave a comment:
-
Bowtie error on tophat
Hi everyone,
Im new in RNA seq data analysis Im just learning how to install tophat, bowtie and so on, doing the troubleshooting of tophat recommended by their website I obtain this:
2013-02-18 16:10:51] Checking for Bowtie
Bowtie 2 not found, checking for older version..
Error: Bowtie not found on this system.
the line that I used is this:
./tophat -r 20 test_ref reads_1.fq reads_2.fq
Prevously I Installed bowtie2 and guessed only to change bowtie files to tophat directory was enough but doesnt work so what should I do??
Thanks
Leave a comment:
-
That is not true. I downloaded bosTau7.fa and had a look at the beginning of the file with R.
Code:> readLines("F:\\bosTau7.fa", n=10) [1] ">chr1" [2] "taccccactcacacttatggatagatcaactaaacagaaaattaacaagg" [3] "aaaaaaaaaAAAAAAAAAACATGTATCTATAAAGCTCACTAAAGCAAAGC" [4] "ACAGTAAAATGAGGTATGGCAATAAAAGGCGTTAAATTTATCaaaaaaaa" [5] "aaaaaaaaagtgctgaaaaattaaaaaaaaaaaaaagaaaaaaaTATAGA" [6] "GATAAAATTAAATTTTTGTATTTATTCTTGCCGTTGAAAATTATGGGCAA" [7] "CTTTTCATTTTTTCAAATTTAATATTGTGTTACTCAAAAATACCCTAAAG" [8] "TTGATCTCTGTTTAAATCTAGAACACTTTTTGTTAAGTTTATGTCGATAT" [9] "GTCTTAGCTCTTTTTTTTAATTAGACATGATATTTTCTTTCATCACATTT" [10] "TCAATATGTAGGAGAGCTAGTGATATTTGTGTAATAATTTTGTACTCGGC"
Leave a comment:
-
index and gtf conflict.
Originally posted by Dario1984 View PostSee the pattern in the error message. It ends with bt2. They don't have the cow index on their website, so you will have to make it.
I downloaded the bovine genome from this site http://hgdownload.soe.ucsc.edu/golde...sTau7/bigZips/ and I used the bowtie2-build to construct my index files. However, my gtf file was also downloaded from this site http://genome.ucsc.edu/cgi-bin/hgTables ( I used the latest version BaylorBtau_4.6.1/bosTau7). I faced a lot of problems during alignment with tophat and someone mentioned that chromosome numbering differs( for genome is 1 and for the gtf file chr 1). I don't know how to fix this. Can someone point out where I can get a gtf file that is compatible the site I got my genome from?
Please I really need your help as I'm new to rna-seq data analysis.
Thanks.
Leave a comment:
-
Error: Could not find Bowtie 2 index files (/home/workdir/ooa4/rnaseq/.*.bt2)
Leave a comment:
Latest Articles
Collapse
-
by seqadmin
The sequencing world is rapidly changing due to declining costs, enhanced accuracies, and the advent of newer, cutting-edge instruments. Equally important to these developments are improvements in sequencing analysis, a process that converts vast amounts of raw data into a comprehensible and meaningful form. This complex task requires expertise and the right analysis tools. In this article, we highlight the progress and innovation in sequencing analysis by reviewing several of the...-
Channel: Articles
05-06-2024, 07:48 AM -
-
by seqadmin
The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist...-
Channel: Articles
04-22-2024, 07:01 AM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, 05-14-2024, 07:03 AM
|
0 responses
15 views
0 likes
|
Last Post
by seqadmin
05-14-2024, 07:03 AM
|
||
Started by seqadmin, 05-10-2024, 06:35 AM
|
0 responses
37 views
0 likes
|
Last Post
by seqadmin
05-10-2024, 06:35 AM
|
||
Started by seqadmin, 05-09-2024, 02:46 PM
|
0 responses
45 views
0 likes
|
Last Post
by seqadmin
05-09-2024, 02:46 PM
|
||
Started by seqadmin, 05-07-2024, 06:57 AM
|
0 responses
39 views
0 likes
|
Last Post
by seqadmin
05-07-2024, 06:57 AM
|
Leave a comment: