SEQanswers

Go Back   SEQanswers > Sequencing Technologies/Companies > SOLiD



Similar Threads
Thread Thread Starter Forum Replies Last Post
TopHat -paired end vs single end reads adarshjose RNA Sequencing 10 06-12-2012 07:15 PM
Tophat paired end alignment help: comma or blank to separate two pairs? mediator Bioinformatics 3 01-31-2012 09:44 AM
How long should paired-end alignment run? agc Bioinformatics 11 09-07-2011 01:31 AM
MAQ paired-end alignment parameters baohua100 Bioinformatics 0 08-19-2009 01:28 PM
how paired end alignment works? totalnew Bioinformatics 8 04-27-2009 01:46 PM

Reply
 
Thread Tools
Old 03-11-2011, 03:04 AM   #1
adeslat
Junior Member
 
Location: washington area

Join Date: Mar 2011
Posts: 7
Default Tophat paired end alignment.

Dear forum,

I have RNASeq SOLID data that is paired end. I want to use Tophat to perform the alignment. Does anyone out there know how to
ensure that "even if an end of a pair could not be aligned, it should still have an entry"

Aaron Quinlan requested that I ensure that this is so -- because I went to use BEDtools with the -pe option and it failed saying my data was not sorted -- he suggested that it may be because I did not ensure that the pair every seq had an entry even if it was not aligned....

Best,
Anne
adeslat is offline   Reply With Quote
Old 03-11-2011, 06:46 AM   #2
Jon_Keats
Senior Member
 
Location: Phoenix, AZ

Join Date: Mar 2010
Posts: 279
Default

Check the sam header SO identifier. In the current release the sort nomenclature is not matching the SAM specification and will cause some applications the through errors. Either resort using picard, not samtools, or just use sed to modify the SO header to conform to specification
Jon_Keats is offline   Reply With Quote
Old 03-11-2011, 06:51 AM   #3
adeslat
Junior Member
 
Location: washington area

Join Date: Mar 2011
Posts: 7
Smile

Quote:
Originally Posted by Jon_Keats View Post
Check the sam header SO identifier. In the current release the sort nomenclature is not matching the SAM specification and will cause some applications the through errors. Either resort using picard, not samtools, or just use sed to modify the SO header to conform to specification
Jon,
Of all the people to reply -- thanks -- I've been reading your various how-tos -- they are simply wonderful. Thank you -- I will do as you suggest and report back -- picard works better?

Best,
Anne
adeslat is offline   Reply With Quote
Old 03-11-2011, 06:53 AM   #4
Jon_Keats
Senior Member
 
Location: Phoenix, AZ

Join Date: Mar 2010
Posts: 279
Default

In the case of this error that I discovered trying to get picard to remove the duplicates I found that samtools did not correct the SO header while picard did, not sure why that occurs as a bwa generated bam file when sorted with samtools gets the right header.
Jon_Keats is offline   Reply With Quote
Old 03-20-2011, 10:41 PM   #5
adeslat
Junior Member
 
Location: washington area

Join Date: Mar 2011
Posts: 7
Default elsewhere documented errors with samtools

Hi Jon,

I was looking at "Scripture" and found looking through their course material at:

http://molecularevolution.org/resour...pture_activity

the following:

"There is currently an incompatibility between samtools sort and Picard (which is used by Scripture); therefore you need to sort the BAM file using Picard and not samtools:"

By the way -- the instructions on this site are spot on, clean and complete.

The reference is the same Nature Biotechnology issue where the Steven Salzberg and Cole Trapnell so nicely introduced the "Tuxedo" suite in the letter "Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation"
Cole Trapnell, Brian A Williams,Geo Pertea, Ali Mortazavi, Gordon Kwan, Marijke J van Baren, Steven L Salzberg, Barbara J Wold & Lior Pachter
Nature Biotechnology 28, 511–515 (2010) doi:10.1038/nbt.1621

Scripture is discussed in:
"Ab initio reconstruction of cell type–specific transcriptomes in mouse reveals the conserved multi-exonic structure of lincRNAs" - pp503 - 510
Mitchell Guttman, Manuel Garber, Joshua Z Levin, Julie Donaghey, James Robinson, Xian Adiconis, Lin Fan, Magdalena J Koziol, Andreas Gnirke, Chad Nusbaum, John L Rinn, Eric S Lander & Aviv Regev
doi:10.1038/nbt.1633

Both are excellent articles....

Anne
adeslat is offline   Reply With Quote
Old 08-07-2011, 04:05 AM   #6
Joker!sAce
Member
 
Location: Denmark

Join Date: Feb 2011
Posts: 21
Default

Hi,

I looked at this tutorial and was looking for the "scripture_alpha2.jar" file. I cannot find that file anywhere. Any insights or help would be appreciated.

Thanks!
Joker!sAce is offline   Reply With Quote
Reply

Tags
paired-end, rna-seq, solid, tophat rna-seq

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 04:23 PM.


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