This paper is a continuation of recent attempts to understand, via mathematical modeling, the dynamics of marine bacteriophage infections. Previous authors have proposed systems of ordinary differential delay equations with delay dependent coefficients. In this paper we continue these studies in two respects. First, we show that the dynamics is sensitive to the phage mortality function, and in particular to the parameter we use to measure the density dependent phage mortality rate. Second, we incorporate spatial effects by deriving, in one spatial dimension, a delay reaction- diffusion model in which the delay term is rigorously derived by solving a von Foerster equation. Using this model, we formally compute the speed at which the viral infection spreads through the domain and investigate how this speed depends on the system parameters. Numerical simulations suggest that the minimum speed according to linear theory is the asymptotic speed of propagation. − KS(t)P (t), dI dt = −µiI(t )+ KS(t)P (t) − e −µiT KS(t − T )P (t − T ), dP dt = β − µpP (t) − KS(t)P (t )+ be −µiT KS(t − T )P (t − T ).