deletions | additions       diff --git a/code/clean.py b/code/clean.py  index 7b8c058..9f76681 100644  --- a/code/clean.py  +++ b/code/clean.py 
...
epithelial = epithelial.filter(filter_fold_change, fc = fc_cutoff).saveas()  proximal_tube = proximal_tube.filter(filter_fold_change, fc = fc_cutoff).saveas()  kidney = kidney.filter(filter_fold_change, fc = fc_cutoff).saveas()  # Identify only those sites that are peaks in all three tissue types  combined = pybedtools.BedTool().multi_intersect(i pybedtools.BedTool().multi_intersect(  i  = [epithelial.fn, proximal_tube.fn, kidney.fn]) union = combined.filter(lambda x: int(x[3]) == 3).saveas()  union.cut(range(3)).saveas(data + "/sites-union.bed")  diff --git a/version-your-code.tex b/version-your-code.tex  index 329b501..e99e801 100644  --- a/version-your-code.tex  +++ b/version-your-code.tex 
...
epithelial = epithelial.filter(filter_fold_change, fc = fc_cutoff).saveas()  proximal_tube = proximal_tube.filter(filter_fold_change, fc = fc_cutoff).saveas()  kidney = kidney.filter(filter_fold_change, fc = fc_cutoff).saveas()  # Identify only those sites that are peaks in all three tissue types  combined = pybedtools.BedTool().multi_intersect(i pybedtools.BedTool().multi_intersect(  i  = [epithelial.fn, proximal_tube.fn, kidney.fn]) union = combined.filter(lambda x: int(x[3]) == 3).saveas()  union.cut(range(3)).saveas(data + "/sites-union.bed")  \end{verbatim} 
...
epithelial = epithelial.filter(filter_fold_change, fc = fc_cutoff).saveas()  proximal_tube = proximal_tube.filter(filter_fold_change, fc = fc_cutoff).saveas()  kidney = kidney.filter(filter_fold_change, fc = fc_cutoff).saveas()  # Identify only those sites that are peaks in all three tissue types  combined = pybedtools.BedTool().multi_intersect(i pybedtools.BedTool().multi_intersect(  i  = [epithelial.fn, proximal_tube.fn, kidney.fn]) union = combined.filter(lambda x: int(x[3]) == 3).saveas()  union.cut(range(3)).saveas(data + "/sites-union.bed")  \end{verbatim}