gatoravi

ggplot(t_m) + geom_bar(stat = "identity",  aes(x=V3.x, y=V4)) + geom_bar(stat = "identity",  aes(x=V5, y=V6), fill = "red") + facet_wrap(~ V17) +  labs(title = paste("BetaBinomial_pp = ", t_m$V2.y, " DNG_ST_pp = ", t_m$V3.y), x = t_m$V18, y = "Read Count")

gatoravi

#Avinash Ramu, WUSTL
#! /bin/bash

#The cutoff file contains the number of DeNovoGear calls for a particular posterior probability cutoff. Extract the same number of SamTools calls. If there is a tie in the last Likelihood Ratio value then pick all calls with that CLR value 

cutoff_f=1e-8Cutoffs
STop=STop_nonX_sorted

for cutoff in `cat $cutoff_f`
do