Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • dpryan
    replied
    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).

    Leave a comment:


  • EpiBrass
    started a topic Count C versus read length

    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

Latest Articles

Collapse

  • seqadmin
    Best Practices for Single-Cell Sequencing Analysis
    by seqadmin



    While isolating and preparing single cells for sequencing was historically the bottleneck, recent technological advancements have shifted the challenge to data analysis. This highlights the rapidly evolving nature of single-cell sequencing. The inherent complexity of single-cell analysis has intensified with the surge in data volume and the incorporation of diverse and more complex datasets. This article explores the challenges in analysis, examines common pitfalls, offers...
    06-06-2024, 07:15 AM
  • seqadmin
    Latest Developments in Precision Medicine
    by seqadmin



    Technological advances have led to drastic improvements in the field of precision medicine, enabling more personalized approaches to treatment. This article explores four leading groups that are overcoming many of the challenges of genomic profiling and precision medicine through their innovative platforms and technologies.

    Somatic Genomics
    “We have such a tremendous amount of genetic diversity that exists within each of us, and not just between us as individuals,”...
    05-24-2024, 01:16 PM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, Yesterday, 07:49 AM
0 responses
14 views
0 likes
Last Post seqadmin  
Started by seqadmin, 06-20-2024, 07:23 AM
0 responses
14 views
0 likes
Last Post seqadmin  
Started by seqadmin, 06-17-2024, 06:54 AM
0 responses
16 views
0 likes
Last Post seqadmin  
Started by seqadmin, 06-14-2024, 07:24 AM
0 responses
25 views
0 likes
Last Post seqadmin  
Working...
X