Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Slacanch
    replied
    well i don't really care about the sequence, it's for profile comparison.
    as long as the positions are right, i'm ok with it.

    maybe i shouldn't have called it pileup, but coverage. just a measure of how many reads overlap each position, in order to be able to compare profiles at certain loci (such as transcriptional readthrough for example).

    for future reference, i discovered a script that does more or less this in the "misc" folder of samtools. it's called interpolate_sam.pl.

    i'll check it more in detail and report back.

    Leave a comment:


  • dpryan
    replied
    What would a pileup of the uncovered region between mates even look like. They couldn't contribute any sequence, so what would you show for them? The only conceivable reason I could think of to do this is in peak calling for chip-seq or something like that, but then you don't care about the actual sequence so there'd be no point in going to full pileup route.

    Leave a comment:


  • Slacanch
    started a topic generating pileup of paired end fragments

    generating pileup of paired end fragments

    Dear all,
    i'm working with paired end rnaseq data. i wanted to generate a pileup with samtools mpileup but noticed that only the mapping reads get counted, and not the unmapped part between two mates, like so:

    ______READ1--------READ2________ <- mapping
    ______1 1 1 1 0 0 0 1 1 1 1________ <- pileup

    is generating pileups using the whole fragment (read + mate + part in between them) frowned upon for some reason? would it be correct to do it? if yes are there any tools?

    alternatively, is there a tool to merge mates and generate the complete fragment in a sam file?

    Thank you!
    Last edited by Slacanch; 01-20-2015, 08:42 AM.

Latest Articles

Collapse

  • seqadmin
    Recent Advances in Sequencing Technologies
    by seqadmin



    Innovations in next-generation sequencing technologies and techniques are driving more precise and comprehensive exploration of complex biological systems. Current advancements include improved accessibility for long-read sequencing and significant progress in single-cell and 3D genomics. This article explores some of the most impactful developments in the field over the past year.

    Long-Read Sequencing
    Long-read sequencing has seen remarkable advancements,...
    12-02-2024, 01:49 PM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, Today, 08:24 AM
0 responses
8 views
0 likes
Last Post seqadmin  
Started by seqadmin, Yesterday, 07:41 AM
0 responses
7 views
0 likes
Last Post seqadmin  
Started by seqadmin, 12-11-2024, 07:45 AM
0 responses
12 views
0 likes
Last Post seqadmin  
Started by seqadmin, 12-10-2024, 07:59 AM
0 responses
13 views
0 likes
Last Post seqadmin  
Working...
X