SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Geneious and chimera melcar Bioinformatics 1 11-13-2014 05:19 PM
Time to run variant caller plugin for Ion Torrent g3n3 Ion Torrent 2 08-21-2013 03:05 PM
CLC Genome Finishing Plugin mcnelson.phd Bioinformatics 0 02-14-2013 01:13 PM
Assemble RNA-Seq data with TopHat in Geneious Server! Geneious Vendor Forum 0 09-19-2011 10:19 AM
Hello - Anyone ever used Geneious?! jeremyrualosr Bioinformatics 5 04-07-2011 08:57 PM

Reply
 
Thread Tools
Old 11-03-2015, 06:54 AM   #1
zhanghao
Member
 
Location: BEIJING

Join Date: Nov 2015
Posts: 17
Default Tophat plugin in geneious

I used tophat plugin in geneious to map my RNAseq data, but I found the result is much worse than using command line (tophat), many introns are not well mapped. You can see the picture of my setting in Geneious and in Assembly Report, I found the command used in geneious is "home/bac/geneious9.0data/plugins/com.biomatters.plugins.tophat.TophatPlugin/com/biomatters/plugins/tophat/linux64/bowtie2-build input.fa input

/home/bac/geneious9.0data/plugins/com.biomatters.plugins.tophat.TophatPlugin/com/biomatters/plugins/tophat/linux64/tophat --b2-sensitive -N 2 --read-gap-length 2 --read-edit-dist 2 -a 8 -m 0 -i 70 -I 500000 -F 0.15 -g 40 -p 30 --segment-length 25 --library-type fr-unstranded input forwardReads.fastq"

I checked this settings, and found 2 problems:

1. When I used command line version of tophat, I used the reference (PH1_genome.fasta) and its annotation (gene.gtf) seperately. First, index the reference file (fasta) and then run tophat with the option -G gene.gtf (the CDS annotation), for this parameter, topHat will first extract the transcript sequences and use Bowtie to align reads to this virtual transcriptome first. But in geneious, I import the reference (PH1_genome.fasta) first and then import gene.gtf to add on the reference, so in geneious I just used the combined file as reference. In this condition, can geneious extract the transcript sequences first? And I didn't find the command "-G gene.gtf" in Assembly Report in geneious, so is this the root cause of this problem? Maybe I should import the sequence (PH1_genome.fasta) and annotation (gene.gtf) seperately, and add command line "-G gene.gtf" in the option "Additional Command Line Parameters". Is it right???

2. From tophat manual, it said "When running TopHat with paired reads it is critical that the *_1 files an the *_2 files appear in separate comma-delimited lists", but in geneious, I select the 2 documents, but in the command in Assembly Report, I found it is just one document "forwardReads.fastq". This will influence the result?
Attached Images
File Type: png th.png (143.3 KB, 3 views)
zhanghao is offline   Reply With Quote
Reply

Tags
geneious, tophat

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 01:30 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2019, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO