Go Back   SEQanswers > Applications Forums > RNA Sequencing

Similar Threads
Thread Thread Starter Forum Replies Last Post
Read counts from SAM file mapped to de novo assembled transcripts using HTSeq-count alan_sm RNA Sequencing 2 06-12-2015 08:54 PM
fastq.gz stats READ-COUNT BASE-COUNT jgibbons1 Bioinformatics 9 10-30-2013 05:24 AM
Finding Neighboring Genes Issues! jplotnik Bioinformatics 1 10-09-2012 05:18 AM
multiBamCov or htseq-count to count read per feature ? NicoBxl Bioinformatics 1 07-03-2012 02:05 AM
small rna seq with total lane read count different between conditions anle Bioinformatics 0 01-25-2012 01:25 PM

Thread Tools
Old 11-10-2013, 07:28 AM   #1
Location: Santa Fe, NM

Join Date: Nov 2013
Posts: 20
Default six huge read count outliers on neighboring transcripts on one lane - theories?

This is a purely academic question since these particular transcripts are not of particular interest to us (and they're only bad on one out of 71 lanes), but I'm curious as to what the theory may be for what happened to these particular reads.

The setup is that we've run 71 Arabidopsis samples, one per lane, at Otogenetics, and they sent the data over to DNANexus for read mapping against TAIR9. I've downloaded the mapped read counts and used them in DESeq, edgeR, etc. The data is really good, quite consistent across the lanes for almost all transcripts.

Except this one lane, lane #30, has six huge outliers in neighboring transcripts on chromosome 4:

at4g12470 799 bp from Chr4:7,401,109..7,401,907
at4g12480 833 bp from Chr4:7,406,105..7,406,937
at4g12490 786 bp from Chr4:7,409,621..7,410,406
at4g12500 778 bp from Chr4:7,414,150..7,414,927
at4g12510 568 bp from Chr4:7,417,236..7,417,803
at4g12520 683 bp from Chr4:7,421,056..7,421,738

I've attached a plot below (SVG) that shows these six plus the neighboring two transcripts that look perfectly normal. If this were a microarray assay I'd say there was a scratch on the plate; but with RNA-Seq, any simple explanation?

Sam Hokin
Computational Scientist, Carnegie and NCGR

Last edited by samhokin; 11-10-2013 at 07:33 AM.
samhokin is offline   Reply With Quote
Old 11-12-2013, 08:39 AM   #2
Junior Member
Location: Atlanta, GA

Join Date: Nov 2011
Posts: 2
Default Data outlyer

This is Natalie from Otogenetics. I would like to discuss your data with you to see if there is anything on our end that could help. Unfortunately I cannot tell which sample has your outlying data. Could you contact us at or give us a call? I would be happy to discuss the sample offline.
phdveg is offline   Reply With Quote


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 05:17 PM.

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