The current GATK version is 3.7-0
Examples: Monday, today, last week, Mar 26, 3/26/04

Howdy, Stranger!

It looks like you're new here. If you want to get involved, click one of these buttons!

Get notifications!


You can opt in to receive email notifications, for example when your questions get answered or when there are new announcements, by following the instructions given here.

Did you remember to?


1. Search using the upper-right search box, e.g. using the error message.
2. Try the latest version of tools.
3. Include tool and Java versions.
4. Tell us whether you are following GATK Best Practices.
5. Include relevant details, e.g. platform, DNA- or RNA-Seq, WES (+capture kit) or WGS (PCR-free or PCR+), paired- or single-end, read length, expected average coverage, somatic data, etc.
6. For tool errors, include the error stacktrace as well as the exact command.
7. For format issues, include the result of running ValidateSamFile for BAMs or ValidateVariants for VCFs.
8. For weird results, include an illustrative example, e.g. attach IGV screenshots according to Article#5484.
9. For a seeming variant that is uncalled, include results of following Article#1235.

Did we ask for a bug report?


Then follow instructions in Article#1894.

Formatting tip!


Wrap blocks of code, error messages and BAM/VCF snippets--especially content with hashes (#)--with lines with three backticks ( ``` ) each to make a code block as demonstrated here.

Jump to another community
Picard 2.9.4 is now available. Download and read release notes here.
GATK 3.7 is here! Be sure to read the Version Highlights and optionally the full Release Notes.

Pileup format

arshiarshi Member
edited August 2012 in Ask the GATK team

I am running Pileup with the verbose option. I have two questions regarding it.
(1) Why are all the value in the mapping quality column 0 ?
(2)There is another column, not mentioned in the description of pileup, separated by '@'. What does this column mean ?

11 86988 A A D 0 C37@931@1036@0
11 86989 G G D 0 C37@932@1036@0
11 86990 T T B 0 C37@933@1036@0
11 86991 G G D 0 C37@934@1036@0
11 86992 A A B 0 C37@935@1036@0
11 86993 C C C 0 C37@936@1036@0
11 86994 C CCC D=A 0 C37@937@1036@0,38@0@100@0,39@0@100@0

Thanks,
Arshi

Tagged:

Best Answer

  • ebanksebanks Broad InstituteMember, Broadie, Dev
    Accepted Answer

    The mapping quality isn't emitted by default, so you can't be seeing them at all with that command line. I think perhaps you are seeing 0s because there are no RODs (e.g. VCFs) being input.

Answers

  • ebanksebanks Broad InstituteMember, Broadie, Dev

    Thanks for reporting this. I'm just about to add documentation for the verbose output. Here's what it will say:
    In addition to the standard pileup output, adds 'verbose' output too. The verbose output contains the number of spanning deletions, and for each read in the pileup it has the read name, offset in the base string, read length, and read mapping quality. These per read items are delimited with an '@' character.

  • Thanks !
    Could you also help me with the issue where I am getting all 0 mapping qualities ?. I have checked my data in IGV, and very few of my reads have 0 mapping qualities.
    I get the correct phred quality scores, though.
    This is how I am running Pileup,
    java -Xmx8g -path/toGATK.jar/ \
    -T Pileup \
    -R path/toGATK/resources/hg19.fa \
    -I a.bam \
    -o a.pileup

    Thanks

  • ebanksebanks Broad InstituteMember, Broadie, Dev
    Accepted Answer

    The mapping quality isn't emitted by default, so you can't be seeing them at all with that command line. I think perhaps you are seeing 0s because there are no RODs (e.g. VCFs) being input.

  • Thanks a lot Eric !. The INDEL is a great option in pileup. I am also trying to get all the INDELS and SNPs through UnifiedGenotyper (-glm BOTH). Is there a way that GATK can output the number of Indels at each position. Similar to a pileup format ?. I am interested in both known and predictive INDELS and their count.
    Perhaps I can use the .vcf file from UnifiedGenotyper ?.

    Thanks,
    Arshi

  • Just to be a littel clear, I tried the --metadata option in Pileup and used the 1000G_indel.vcf file as RODs.

  • ebanksebanks Broad InstituteMember, Broadie, Dev

    Hmm, no I don't think you can do what you want with the GATK right now.

  • Ok. Thanks for your quick reply !

Sign In or Register to comment.