diff --git a/diversity.cabal b/diversity.cabal
--- a/diversity.cabal
+++ b/diversity.cabal
@@ -2,7 +2,7 @@
 -- documentation, see http://haskell.org/cabal/users-guide/
 
 name:                diversity
-version:             0.3.5.1
+version:             0.4.0.0
 synopsis:            Return the diversity at each position for all sequences in a fasta file
 description:         Find the diversity of a collection of entities, mainly for use with fasta sequences. Produces a binary which works on fasta files to find the diversity of any order and rarefaction curves for a sliding window across all positions in the sequences. To analyze just a collection of entities, just use the whole sequences and list flag.
 homepage:            https://github.com/GregorySchwartz/diversity
diff --git a/src/src-exec/Main.hs b/src/src-exec/Main.hs
--- a/src/src-exec/Main.hs
+++ b/src/src-exec/Main.hs
@@ -78,7 +78,8 @@
          <> help "The start point and interval of subsamples in the\
                  \ rarefaction curve. For instance, '1 1' would be 1, 2, 3, ...\
                  \ '2 6' would be 2, 8, 14, ... Note: input is a string so\
-                 \ use quotations around the entry" )
+                 \ use quotations around the entry and it always has the\
+                 \ number of subsamples overall as the last point" )
       <*> switch
           ( long "fast-bin"
          <> short 'f'
diff --git a/src/src-lib/Math/Diversity/Diversity.hs b/src/src-lib/Math/Diversity/Diversity.hs
--- a/src/src-lib/Math/Diversity/Diversity.hs
+++ b/src/src-lib/Math/Diversity/Diversity.hs
@@ -54,7 +54,7 @@
                  -> [a]
                  -> [(Int, Double)]
 rarefactionCurve fastBin start interval xs =
-    map rarefact [start,(start + interval)..n_total]
+    map rarefact $ [start,(start + interval)..(n_total - 1)] ++ [n_total]
   where
     rarefact n
         | n == 0       = (fromIntegral n, 0)
@@ -85,7 +85,7 @@
                        -> [(a, b)]
                        -> [(Int, Double)]
 rarefactionSampleCurve fastBin start interval ls =
-    map rarefact [start,(start + interval)..t_total]
+    map rarefact $ [start,(start + interval)..(t_total - 1)] ++ [t_total]
   where
     rarefact t
         | t == 0       = (t, 0)