The current GATK version is 3.8-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!

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.

#### ☞ Got a problem?

1. Search using the upper-right search box, e.g. using the error message.
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.

#### ☞ 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.

GATK version 4.beta.3 (i.e. the third beta release) is out. See the GATK4 beta page for download and details.

# Pileup format

Member
edited August 2012

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:

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 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.

• Member

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

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.

• Member

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

• Member

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