def precursor_sequence(loci, reference):
region = "<STR_LIT>" % (loci[<NUM_LIT:1>], loci[<NUM_LIT:2>], loci[<NUM_LIT:3>], loci[<NUM_LIT:4>])<EOL>precursor = pybedtools.BedTool(str(region), from_string=True).sequence(fi=reference, s=True)<EOL>return open(precursor.seqfn).read().split("<STR_LIT:\n>")[<NUM_LIT:1>]<EOL>
Get sequence from genome