diff --git a/micall/monitor/coverage_plots.R b/micall/monitor/coverage_plots.R index da4feb12a..fa24c17f1 100755 --- a/micall/monitor/coverage_plots.R +++ b/micall/monitor/coverage_plots.R @@ -225,7 +225,7 @@ data <- read.csv(file=input.csv, header=TRUE, sep=',') #seed,region,q-cutoff,query.aa.pos,refseq.aa.pos,A,C,D,E,F,G,H,I,K,L,M,N,P,Q,R,S,T,V,W,Y,* #HCV1A-H77-core-seed,HCV1A-H77-core,0,0,1,0,0,0,0,0,0,0,3,0,0,1542,0,0,0,0,0,0,0,0,0,0 -alphabet <- c('A', 'C', 'D', 'E', 'F', 'G', 'H', 'I', 'K', 'L', 'M', 'N', 'P', 'Q', 'R', 'S', 'T', 'V', 'W', 'Y') # X. represents stop codon '*' +alphabet <- c('A', 'C', 'D', 'E', 'F', 'G', 'H', 'I', 'K', 'L', 'M', 'N', 'P', 'Q', 'R', 'S', 'T', 'V', 'W', 'Y', 'X.') # X. represents stop codon '*' data$coverage <- apply(data[ , which(is.element(names(data), alphabet))], 1, sum)