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

Do not die in --cycleOnly mode when the lane is not paired end, just count all single end basequals into the first column and leave the second column filled with 0s 

git-svn-id: file:///humgen/gsa-scr1/gsa-engineering/svn_contents/trunk@2778 348d0f76-0448-11de-a6fe-93d51630548a
This commit is contained in:
asivache 2010-02-04 19:48:12 +00:00
parent a3cd56897d
commit 152f65b362
1 changed files with 5 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
java/src/org/broadinstitute/sting/oneoffprojects/walkers/QualityScoreByStrandWalker.java
View File

@ -82,7 +82,11 @@ public class QualityScoreByStrandWalker extends LocusWalker<StrandedCounts,Stran
if ( ! useCycle ) {
counts.updateLocus( (int) e.getQual(), ! e.getRead().getReadNegativeStrandFlag() );
} else {
counts.updateLocus( (int) e.getQual(), ! e.getRead().getFirstOfPairFlag() );
if ( e.getRead().getReadPairedFlag() ) {
counts.updateLocus( (int) e.getQual(), e.getRead().getFirstOfPairFlag() );
} else {
counts.updateLocus( (int) e.getQual(), true );
}
}
}