title: Fitting distributions, goodness of fit, p-value. Is it possible to do this with Scipy (Python)?
body: <p>INTRODUCTION: I'm a bioinformatician. In my analysis which I perform on all human genes (about 20 000) I search for a particular short sequence motif to check how many times this motif occurs in each gene. </p>
  
  <p>Genes are 'written' in a linear sequence in four letters (A,T,G,C). For example: CGTAGGGGGTTTAC... This is the four-letter alphabet of genetic code which is like the secret language of each cell, it’s how DNA actually stores information.</p>
  
  <p>I suspect that frequent repetitions of a particular short motif sequence (AGTGGAC) in some genes are crucial in a specific biochemical process in the cell. Since the motif itself is very short it is difficult with computational tools to distinguish between true functional examples in genes and those that look similar by chance. To avoid this problem I get sequences of all genes and concatenated into a single string and shuffled. The length of each of the original genes was stored. Then for each of the original sequence lengths, a random sequence was constructed by repeatedly picking A or T or G or C at random from the concatenated sequence and transferring it to the random sequence. In this way, the resulting set of randomized sequences has the same length distribution, as well as the overall A,T,G,C composition. Then I search for the motif in these randomized sequences. I performed this procedure 1000 times and averaged the results.</p>
  
  <ul>
  <li>15000 genes that do not contain a given motif</li>
  <li>5000 genes that contain 1 motif</li>
  <li>3000 genes that contain 2 motifs</li>
  <li>1000 genes that contain 3 motifs</li>
  <li>...</li>
  <li>1 gene that contain 6 motifs</li>
  </ul>
  
  <p>So even after 1000 times randomization of true genetic code, there aren't any genes which have more than 6 motifs. But in the true genetic code, there are a few genes which contain more then 20 occurrences of the motif, which suggest that these repetition might be functional and it's unlikely to find them in such an abundance by pure chance.</p>
  
  <p>PROBLEM:</p>
  
  <p>I would like to know the probability of finding a gene with let's say 20 occurrences of the motif in my distribution. So I want to know the probability to find such a gene by chance. I would like to implement this in Python, but I don't know how.</p>
  
  <p>Can I do such an analysis in Python?</p>
  
  <p>Any help would be appreciated. </p>