SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
fastq.gz stats READ-COUNT BASE-COUNT jgibbons1 Bioinformatics 9 10-30-2013 05:24 AM
How do you get the gene length that htseq-count is using? jebe RNA Sequencing 0 03-07-2013 08:37 AM
Tophat: options -N --read-edit-dist --read-gap-length Pradhaun Bioinformatics 0 01-04-2013 07:58 AM
picard error: Mismatch between read length and quals length writing read shawpa Bioinformatics 0 08-20-2012 05:52 AM
multiBamCov or htseq-count to count read per feature ? NicoBxl Bioinformatics 1 07-03-2012 02:05 AM

Reply
 
Thread Tools
Old 01-27-2014, 08:55 PM   #1
EpiBrass
Member
 
Location: Australia

Join Date: Nov 2013
Posts: 16
Default Count C versus read length

Hi,

I've had a good look through the bioinformatics thread but as far as I can see nobody has asked this question yet. I am using bisulphite sequencing in plants (CpG, CHG & CHH contexts).

I need to compare different sequencing runs post bisulphite conversion to determine the reproducibility between runs and library preps of the same sample (technical reps). It has been suggested that one way to approach this is to compare the number of C's in a read against the read length post trimming. This can then be graphed and in theory the longer each read is the more C's will be present. If the technical reps are not biased then the trend should be the same for all reps.

Does anyone have any suggestions of how to output from .fastq or .fasta the number Cs in a read and the read length into a tab-delimited file so it can then be graphed?

Many thanks,

Justin
EpiBrass is offline   Reply With Quote
Old 01-28-2014, 12:58 AM   #2
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480
Default

It sounds like you're describing an "M-bias" graph, an example of which is below (this one has bias):


If you create alignments with either bismark or bison you can make those sorts of graphs easily. With bison, it's the "bison_mbias -pdf something.bam" command, which requires that you have R and the ggplot2 package installed (otherwise, drop the "-pdf" part and it'll just output a tab-separated file that you can use for graphing). Similarly, bismark will produce them but I don't recall off-hand what the command is. I think bismark uses perl GD::Graph.

There's also a recent paper in Bioinformatics on BSeqQC or something like that that can also produce M-bias graphs.

I should not that these graphs aren't looking so much at reproducibility as variance in library setup (bisulfite conversion efficiency, etc.). If you really want to look at reproducibility, it would make more sense to directly compare the methylation calls and calculate how many of them are significantly different (it should by approximately whatever alpha level you use, so ~5% of sites for p<0.05 threshold).
dpryan is offline   Reply With Quote
Reply

Tags
base count, bisulphite sequencing, ngs, read length, wgbs

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 11:02 PM.


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