Please Help Interpret Picard Bamindexstats Results
1
0
Entering edit mode
10.8 years ago
newDNASeqer ▴ 760

I have difficulty understanding the statistics generated by PiCard BamIndexStats, part of the my stats results is:

chr7_gl000195_random length= 182896 Aligned= 3135 Unaligned= 0 chr8 length= 146364022 Aligned= 1627297 Unaligned= 0 chr8_gl000196_random length= 38914 Aligned= 869 Unaligned= 0

Does this result mean that, for chr 7 and 8, 100% of the reads are aligned to the chromosomes regions? I used bwa mem to do the alignment, is it sort of too good to be true?

picard statistics bwa • 2.3k views
ADD COMMENT
0
Entering edit mode

Did you run with the LENIENT flag? I ask because Picard will discard unmapped reads if they have a MAPQ of zero.

ADD REPLY
0
Entering edit mode

Other things that would be helpful: can you tell us your command line for Picard, and perhaps include a few lines from the BAM?

ADD REPLY
0
Entering edit mode
10.8 years ago
always_learning ★ 1.1k

You may have bam file and corresponding index files for that cases here !! So as I think It will always show 100% Match here !! I may be wrong but would like to see other response !! Meanwhile what about other Chromosomes also !!

ADD COMMENT

Login before adding your answer.

Traffic: 1622 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6