Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • hstehr
    replied
    -s option for mpileup does exist in samtools 0.1.17

    In samtools 0.1.17 the mpileup command does indeed have the -s option, even though it is not mentioned in the manpage. It is listed under output options when typing 'samtools mpileup'. In version 0.1.16, the -s option is not available for mpileup, but only for pileup.

    I hope this answers the original question.

    However, I do not understand the encoding of the mapping qualities in the -s output. It does not seem to be in the usual Sanger/Phred encoding. Can anyone point me to where this is documented?

    Leave a comment:


  • earonesty
    replied
    Originally posted by nilshomer View Post
    Check out the paper: http://bioinformatics.oxfordjournals...s.btr076.short

    Anyhow, it uses the base qualities, not the mapping qualities. It resolves and weights the local ambiguities. Smith Waterman can be modeled as an HMM, with the BAQ running the forward/backward algorithms to compute the posterior probability of a query base and target base aligning based on all possible local alignments. This posterior is then converted to a base quality.
    Alignment quality is, effectively, summarized in the "mapping quality" MAQ statistic. BAQ is a "per base alignment quality",rather than a single statistic. The HMM model in BAQ is used to adjust base qualities , but the reason for adjustment is the alignment quality of each base.

    The resulting pileup has lowered qualities on those bases whose alignment was ambiguous.

    Using both MAQ and BAQ is, thus, a bit overlapped.

    That's all I was trying to point out.

    Leave a comment:


  • ninjapanda57
    replied
    Solution

    If you look at the MQ part of the info section in the VCF format that does equal mapping quality. The way that I extract that to work with the bed file is in perl

    all you do is

    call the line
    my @snp_array = $snp =~ /^(\S+)\s+(\S+)\s+(\S+)\s+(\S+)\s+\S+\s+(\S+)\s+/;
    my $mpileup_check = $snp_array[2];
    my ($chr, $pos, $id, $ref, $var, $cns_qual, $rd_depth, $map_qual, $snp_qual, @extra);
    if ($mpileup_check eq '.') {
    ($chr, $pos, $id, $ref, $var, $cns_qual, @extra) =split("\t", $snp);
    $extra[1] =~ /DP(\d+)/;
    $rd_depth = $1;
    $extra[1] =~ /MQ(\d+)/;
    $map_qual = $1;
    } else
    that should work, that's just in perl though.

    Hope this helps.

    Leave a comment:


  • nilshomer
    replied
    Check out the paper: http://bioinformatics.oxfordjournals...s.btr076.short

    Anyhow, it uses the base qualities, not the mapping qualities. It resolves and weights the local ambiguities. Smith Waterman can be modeled as an HMM, with the BAQ running the forward/backward algorithms to compute the posterior probability of a query base and target base aligning based on all possible local alignments. This posterior is then converted to a base quality.

    Leave a comment:


  • earonesty
    replied
    I thought the algorithm adjusted base qualitieis... based on mapping qualities. IE... it does it so you don't have to.

    Leave a comment:


  • nilshomer
    replied
    Originally posted by earonesty View Post
    mpileup is *supposed* to handle all the problems with mapping quality internally with the BAQ algorithm.
    BAQ adjusts the base qualities not the mapping qualities. The mapping qualities are set by the aligner/mapper.

    Leave a comment:


  • earonesty
    replied
    mpileup, pileup, Bio:B::Sam

    mpileup is *supposed* to handle all the problems with mapping quality internally with the BAQ algorithm. I think the resulting entries in the pileup have been quality-adjusted so that the phred scores should reflect the BAQ-adjusted values, right? Thus mpileup's output should be "usable without needing the mapping quality".

    Am I reading the mpileup help page incorrectly? It skips bases with BAQ quality < 13, etc. The final pileup does not reflect the original reads exactly, but rather a BAQ-corrected version.



    I've generally found good results, better than trying to handle mapping-quality issues on my own ... but it seems odd to handle all that, and then *remove* the option that would make verification easier.

    Bio:B::Sam's pileup() call is a great tool for walking a pileup, while having access to ALL possible information, and it's way faster than running and then parsing the results of samtools pileup.
    Last edited by earonesty; 06-29-2011, 08:26 AM. Reason: add a title

    Leave a comment:


  • FCAlive
    replied
    Bumpity bump bump bump bump

    Bumpity bump bump bump bump

    Leave a comment:


  • kweber2
    replied
    And I still have the same question too, if anyone has any answers...

    Leave a comment:


  • FCAlive
    replied
    I have the same question

    I have the same question.

    bump

    Leave a comment:


  • kweber2
    started a topic Getting map qualities with samtools mpileup?

    Getting map qualities with samtools mpileup?

    I've been using the -s option of samtools "pileup" to output the mapping qualities for each read at each position as a last column in the pileup file.

    I know that the "pileup" command is deprecated, and I need to switch to the "mpileup" command, but "mpileup" doesn't seem to have an option for including the map qualities. Does anyone have any advice? I know that the map quality is still reported once at the beginning of each read, but it seems like it will be tricky to map that quality to the correct base in each subsequent position. I'm hoping there's an easier alternative...

    Thanks for any help!

    Kris Weber
    Research Assistant
    Noble Lab, Department of Genome Sciences
    University of Washington

Latest Articles

Collapse

  • seqadmin
    Exploring the Dynamics of the Tumor Microenvironment
    by seqadmin




    The complexity of cancer is clearly demonstrated in the diverse ecosystem of the tumor microenvironment (TME). The TME is made up of numerous cell types and its development begins with the changes that happen during oncogenesis. “Genomic mutations, copy number changes, epigenetic alterations, and alternative gene expression occur to varying degrees within the affected tumor cells,” explained Andrea O’Hara, Ph.D., Strategic Technical Specialist at Azenta. “As...
    07-08-2024, 03:19 PM
  • seqadmin
    Exploring Human Diversity Through Large-Scale Omics
    by seqadmin


    In 2003, researchers from the Human Genome Project (HGP) announced the most comprehensive genome to date1. Although the genome wasn’t fully completed until nearly 20 years later2, numerous large-scale projects, such as the International HapMap Project and 1000 Genomes Project, continued the HGP's work, capturing extensive variation and genomic diversity within humans. Recently, newer initiatives have significantly increased in scale and expanded beyond genomics, offering a more detailed...
    06-25-2024, 06:43 AM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, 07-19-2024, 07:20 AM
0 responses
26 views
0 likes
Last Post seqadmin  
Started by seqadmin, 07-16-2024, 05:49 AM
0 responses
41 views
0 likes
Last Post seqadmin  
Started by seqadmin, 07-15-2024, 06:53 AM
0 responses
46 views
0 likes
Last Post seqadmin  
Started by seqadmin, 07-10-2024, 07:30 AM
0 responses
42 views
0 likes
Last Post seqadmin  
Working...
X