diff --git a/src/intersect_SNP.R b/src/intersect_SNP.R
index 6cd39fb2904864a5a62398fd7720b76f4854071c..b2f7df28490762189716463a36e18774e9f873b6 100755
--- a/src/intersect_SNP.R
+++ b/src/intersect_SNP.R
@@ -17,7 +17,7 @@ snp <- snp_b %>%
   filter(cords %in% only_b)
 
 fastafile <- read.fasta(file = args[3],
-                        as.string = TRUE) 
+                        as.string = TRUE)
 
 snp$seq_list <- snp %>%
   apply(1, FUN = function(x, fastafile, POS, CHROM){
@@ -38,5 +38,6 @@ snp$seq_list <- snp %>%
   )
 
 snp %>%
+  arrange(desc(tumor_sample.AF)) %>%
   write_csv(paste0("only_", args[2]))