| Previous changeset 2:d6de2d1f4af9 (2013-03-03) Next changeset 4:fafa105e5f58 (2013-03-19) |
|
Commit message:
Reference skips have a base quality score |
|
modified:
pileup_to_vcf.py |
| b |
| diff -r d6de2d1f4af9 -r aa76c8dd97e6 pileup_to_vcf.py --- a/pileup_to_vcf.py Sun Mar 03 20:38:43 2013 -0600 +++ b/pileup_to_vcf.py Tue Mar 19 12:49:45 2013 -0500 |
| [ |
| @@ -158,6 +158,7 @@ if debug: print >> sys.stderr, "%3d\t%s\t%3d\t%s %3d" % (bi,base,qi,quals[qi],ord(quals[qi]) - 33) bi += 1 if base in '<>' : #reference skip (between paired reads) + qi += 1 pass; elif base in '.,' : # match reference on forward/reverse strand mc += 1 |