zzh
/
gatk-3.8
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Fix minor mapping quality bug. 

git-svn-id: file:///humgen/gsa-scr1/gsa-engineering/svn_contents/trunk@1973 348d0f76-0448-11de-a6fe-93d51630548a
This commit is contained in:
hanna 2009-11-04 14:33:23 +00:00
parent 21b1ece45d
commit 43c3ee61d5
2 changed files with 5 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

3
c/bwa/bwa_gateway.cpp
View File

@ -24,6 +24,9 @@ BWA::BWA(const char* ann_filename,
bwts[1] = bwt_restore_bwt(reverse_bwt_filename);
bwt_restore_sa(reverse_sa_filename, bwts[1]);
load_default_options();
// initialize the bwase subsystem
bwase_initialize();
}
BWA::~BWA() {

4
java/src/org/broadinstitute/sting/alignment/bwa/AlignmentValidationWalker.java
View File

@ -55,7 +55,7 @@ public class AlignmentValidationWalker extends ReadWalker<Integer,Integer> {
@Override
public boolean filter(char[] ref, SAMRecord read) {
return true;
//return read.getReadName().contains("SL-XBC:1:76:604:397#0");
//return read.getReadName().contains("SL-XBC:1:61:1719:1512#0");
}
/**
@ -84,7 +84,7 @@ public class AlignmentValidationWalker extends ReadWalker<Integer,Integer> {
matches &= (alignment.getCigar().equals(read.getCigar()));
int mapQDelta = Math.abs(alignment.getMappingQuality()-read.getMappingQuality());
maxMapQDelta = Math.max(mapQDelta,maxMapQDelta);
//matches &= (alignment.getMappingQuality() == read.getMappingQuality());
matches &= (alignment.getMappingQuality() == read.getMappingQuality());
if(matches) break;
}
}