def getForwardreadpos(targetPos,queryStart,hitStart,hitStrand):
    """
    This takes a nucleotide position which relates to the database hit
    (i.e the V region gene from IMGT) and identifies its position in the query sequence.
    It requires a target position (i.e. the first base of start codon, where the first base=1) and
    queryStart and hitStart data from the SearchIO object.
    It returns the Pythonesque position in the query.
    """
    if hitStrand == 1:
        offset =  hitStart - queryStart
    elif hitStrand == -1:
        offset =  - (hitStart + queryStart)
    else:
        raise ValueError('hitStrand must be -1 or 1')
    queryPos = targetPos + offset
    return(queryPos)