If the well-annotated species and yours are highly related, then you'll likely have luck just aligning your reads against it and using it's reference annotation to call variants in promoter regions. Otherwise, you'd be better served by assembling your species and trying to align the contigs (you'd then find differences in the results). Since these are subspecies we're talking about, I assume the former possibility will work well and prove easier than the latter.