Entering edit mode
4.8 years ago
worarado.kan
▴
20
I try running mapping to reference using tophat -cluffink after I run bowtie2-build I get GDDH13_1-1_formatted.1.bt2 GDDH13_1-1_formatted.2.bt2 GDDH13_1-1_formatted.3.bt2 GDDH13_1-1_formatted.4.bt2 GDDH13_1-1_formatted.rev.1.bt2 GDDH13_1-1_formatted.rev.2.bt2
Next I have to run tophat2 but I have a problem, so what is the command i should use for running tophat2?
Please recommend me.
tophat-2.0.13.Linux_x86_64/tophat -G ref.gtf --output-dir out index_prefix input_reads.[fasta|fastq]] Try using this command
Thank you for your message
Checking for Bowtie index files (genome).. Error: Could not find Bowtie 2 index files (GDDH13_1-1_formatted.2.bt2.*.bt2l) [hscience@localhost Tophat2]$
it always showed error. I dont understand why can not find index
Are you running from the Tophot2 folder? It's better to run from the folder that contains the data so you don't need to add the path to the directory. Also you can use
ll
orls -l
in your folder to make sure the files are readable.http://linuxcommand.org/lc3_lts0090.php
These links might help:
http://www.sthda.com/english/wiki/tophat2-download-build-reference-genome-and-align-the-reads-to-the-reference-genome
https://github.com/Jeanielmj/bioinformatics-workshop/wiki/Read-Alignment-with-TopHat2
Thank you so much for information