We first confirmed that our model was able to capture the essential tumor population dynamics by showing a qualitative equivalence between the patient’s clinical course and our model simulation of similar tumor subpopulations consisting of 94% EGFRL858R, 6% BRAF V600E and assuming the existence of a very low initial frequency of 0.01% MET amplification of EGFRL858R, BRAFV600E and EGFRT790M in the presence of 1 μM erlotinib (Fig. 2A,B). This evidence concerns the gene BRAF and neoplasm.