SEQanswers

Go Back   SEQanswers > Applications Forums > RNA Sequencing



Similar Threads
Thread Thread Starter Forum Replies Last Post
heterozygosity rate joxcargator73 Bioinformatics 2 10-08-2014 04:24 PM
Optimise TopHat2 speed and results sindrle RNA Sequencing 8 08-18-2013 07:45 AM
Tophat2 normal alignment rate JJP RNA Sequencing 7 08-12-2013 09:16 AM
Tophat2 Alignment Rate SHeaph Bioinformatics 5 11-21-2012 01:05 AM
error rate der_eiskern Illumina/Solexa 0 12-11-2009 02:51 PM

Reply
 
Thread Tools
Old 11-14-2014, 05:49 AM   #1
duer777
Junior Member
 
Location: Belgium

Join Date: Nov 2014
Posts: 2
Post rate my tophat2 results :-)

Hi,
New in RNA-seq and without any help in my lab, i would like to ask you your advice about my mapping results. (In my opinion, it's ok, but i will be very happy if experimented poeple can confirm it or not).
here are align summary tophat results and accepted_hits.bam analysed by samtool flagstat.

Left reads:
Input : 16641634
Mapped : 15985940 (96.1% of input)
of these: 4863655 (30.4%) have multiple alignments (7966 have >20)
Right reads:
Input : 16641634
Mapped : 15922959 (95.7% of input)
of these: 4834402 (30.4%) have multiple alignments (7972 have >20)
95.9% overall read mapping rate.

Aligned pairs: 15508899
of these: 4652618 (30.0%) have multiple alignments
87521 ( 0.6%) are discordant alignments
92.7% concordant pair alignment rate.

=============

47082928 + 0 in total (QC-passed reads + QC-failed reads)
0 + 0 duplicates
47082928 + 0 mapped (100.00%:-nan%)
47082928 + 0 paired in sequencing
23593110 + 0 read1
23489818 + 0 read2
44621168 + 0 properly paired (94.77%:-nan%)
45651258 + 0 with itself and mate mapped
1431670 + 0 singletons (3.04%:-nan%)
2790766 + 0 with mate mapped to a different chr
86220 + 0 with mate mapped to a different chr (mapQ>=5)

Thanks for your advice.
duer777 is offline   Reply With Quote
Reply

Tags
rna-seq advice, tophat2 bowtie2

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 09:59 PM.


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