# HG changeset patch
# User davidvanzessen
# Date 1513072400 18000
# Node ID 7ce82833977c80ab6274e6de167a3342e0727ec5
# Parent  a9d2ed661541b2e923a1cd6f4c207e7853dde88b
Uploaded

diff -r a9d2ed661541 -r 7ce82833977c RScript.r
--- a/RScript.r	Thu Oct 05 09:06:21 2017 -0400
+++ b/RScript.r	Tue Dec 12 04:53:20 2017 -0500
@@ -34,6 +34,7 @@
 }
 
 dat = dat[,c("Patient",  "Receptor", "Sample", "Cell_Count", "Clone_Molecule_Count_From_Spikes", "Log10_Frequency", "Total_Read_Count", "J_Segment_Major_Gene", "V_Segment_Major_Gene", "CDR3_Sense_Sequence", "Clone_Sequence")] 
+dat = dat[!(nchar(dat$Clone_Sequence) < 2),]
 
 dat$dsPerM = 0
 dat = dat[!is.na(dat$Patient),]