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

Merge branch 'master' of ssh://nickel.broadinstitute.org/humgen/gsa-scr1/gsa-engineering/git/unstable

This commit is contained in:
Ryan Poplin 2011-11-10 13:39:24 -05:00
parent 26762d6c6f 39678b6a20
commit 07dbf0bd40
1 changed files with 11 additions and 8 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

19
public/java/src/org/broadinstitute/sting/gatk/filters/MalformedReadFilter.java
View File

@ -50,19 +50,20 @@ public class MalformedReadFilter extends ReadFilter {
public boolean filterOut(SAMRecord read) { public boolean filterOut(SAMRecord read) {
// slowly changing the behavior to blow up first and filtering out if a parameter is explicitly provided // slowly changing the behavior to blow up first and filtering out if a parameter is explicitly provided
if (!checkMismatchingBasesAndQuals(read)) {
if (!filterMismatchingBaseAndQuals)
throw new UserException.MalformedBAM(read, "BAM file has a read with mismatching number of bases and base qualities. Offender: " + read.getReadName() +" [" + read.getReadLength() + " bases] [" +read.getBaseQualities().length +"] quals");
else
return true;
}
return !checkInvalidAlignmentStart(read) || return !checkInvalidAlignmentStart(read) ||
!checkInvalidAlignmentEnd(read) || !checkInvalidAlignmentEnd(read) ||
!checkAlignmentDisagreesWithHeader(this.header,read) || !checkAlignmentDisagreesWithHeader(this.header,read) ||
!checkHasReadGroup(read) ||
!checkMismatchingBasesAndQuals(read, filterMismatchingBaseAndQuals) ||
!checkCigarDisagreesWithAlignment(read); !checkCigarDisagreesWithAlignment(read);
} }
private static boolean checkHasReadGroup(SAMRecord read) {
if ( read.getReadGroup() == null )
throw new UserException.ReadMissingReadGroup(read);
return true;
}
/** /**
* Check for the case in which the alignment start is inconsistent with the read unmapped flag. * Check for the case in which the alignment start is inconsistent with the read unmapped flag.
* @param read The read to validate. * @param read The read to validate.
@ -127,7 +128,9 @@ public class MalformedReadFilter extends ReadFilter {
* @param read the read to validate * @param read the read to validate
* @return true if they have the same number. False otherwise. * @return true if they have the same number. False otherwise.
*/ */
private static boolean checkMismatchingBasesAndQuals(SAMRecord read) { private static boolean checkMismatchingBasesAndQuals(SAMRecord read, boolean filterMismatchingBaseAndQuals) {
if (!filterMismatchingBaseAndQuals)
throw new UserException.MalformedBAM(read, "BAM file has a read with mismatching number of bases and base qualities. Offender: " + read.getReadName() +" [" + read.getReadLength() + " bases] [" +read.getBaseQualities().length +"] quals");
return (read.getReadLength() == read.getBaseQualities().length); return (read.getReadLength() == read.getBaseQualities().length);
} }
} }