seqIDs <-c(1:nrow(muleys))

muleys <- cbind(muleys,seqIDs)