Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
Local realignment around indels (454 technology) desmo Bioinformatics 2 05-17-2012 06:31 AM
GATK: VCF file for Local realignment around indels jorge Bioinformatics 2 10-11-2011 12:15 AM
Samtools mpileup creates extra large file after local realignment Hkins552 Bioinformatics 7 07-21-2011 02:20 PM
GATK recalibration score and local realignment: what is the first step? m_elena_bioinfo Bioinformatics 0 01-24-2011 03:59 AM
Local realignment using GATK and smra seq_GA Bioinformatics 28 01-17-2011 09:10 AM

Thread Tools
Old 07-07-2011, 03:46 AM   #1
Junior Member
Location: UK

Join Date: Jun 2011
Posts: 2
Question Interpreting Local Realignment results from GATK

Hi there. I need some help in interpreting the result of Local Realignment (around indels) using GATK please. I am new to using the package and trying to understand it. I have run the local realignment around indels according to the GATK manual to reduce the rate of mismatched alignment and false positive SNPs. After running the suspicious-interval-finding step, the output file gives the following 3 messages:
-2757553 reads were filtered out during traversal out of 73748806 total (3.74%)
- 639389 reads (0.87% of total) failing ZeroMappingQualityReadFilter
- 2118164 reads (2.87% of total) failing UnmappedReadFilter
Now here are the problems/questions:
Q1. Word counting (wc l) of my input.bam file shows that I have 96,378,717 number of reads. Why the programme reports only 73,748,806 reads in the first message?
Q2. Similarly, the number of reads with zero MAPQ in the input.bam file is 25,809,346 whereas the output says 639,389 in the 2nd message. Similarly my count of the number of unmapped read is 23,055,530 which is much higher compared to that reported in the 3rd messages (2,118,164). Why these anomalies?
Q4. After realignment I get the message:
0 reads were filtered out during traversal out of 12,369,995 total (0.00%)
Word count ( wc l) of realigned bam file shows that there are 12,369,991 number of reads, which is almost same as reported. But I expected to see at least 70,991,253 (73,748,806 −2,757,553) number of reads in the realigned bam file considering the read numbers reported in first message are correct. Such a low number of reads in the realigned bam file compared to input bam file is a real concern. Can anyone suggest what is going on?

alma is offline   Reply With Quote

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 05:23 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO