-# ../samtools pileup ex1.bam|perl -ape '$_=$F[4];s/\^.//g;$x=tr/A-Z/X/;$y=tr/a-z/x/;s/(\d+)(X)/$x-=$1,$2/eg;s/(\d+)(x)/$y-=$1,$2/eg;$_=join("\t",@F[0,1],$x,$y)."\n"'
\ No newline at end of file
+# ../samtools pileup ex1.bam|perl -ape '$_=$F[4];s/(\d+)(??{".{$1}"})|\^.//g;@_=(tr/A-Z//,tr/a-z//);$_=join("\t",@F[0,1],@_)."\n"'
+
+# ../samtools pileup -cf ex1.fa ex1.bam|perl -ape '$_=$F[8];s/\^.//g;s/(\d+)(??{".{$1}"})|\^.//g;@_=(tr/A-Za-z//,tr/,.//);$_=join("\t",@F[0,1],@_)."\n"'