| Previous changeset 64:c6dd3215ebe0 (2017-12-06) Next changeset 66:43a1aa648537 (2017-12-07) |
|
Commit message:
Uploaded |
|
modified:
new_imgt.r |
| b |
| diff -r c6dd3215ebe0 -r dc06e94bc1e3 new_imgt.r --- a/new_imgt.r Wed Dec 06 08:52:00 2017 -0500 +++ b/new_imgt.r Wed Dec 06 09:55:37 2017 -0500 |
| [ |
| @@ -4,7 +4,12 @@ merged.file = args[2] gene = args[3] -merged = read.table(merged.file, header=T, sep="\t", fill=T, stringsAsFactors=F, comment.char="") +merged = read.table(merged.file, header=T, sep="\t", fill=T, stringsAsFactors=F, comment.char="", quote="") + +if(!("Sequence.ID" %in% names(merged))){ #change-o db + print("Change-O DB changing 'SEQUENCE_ID' to 'Sequence.ID'") + names(merged)[which(names[merged] == "SEQUENCE_ID")] = "Sequence.ID" +} if(gene != "-"){ merged = merged[grepl(paste("^", gene, sep=""), merged$best_match),] |