SelfBlastFilter.py

Created on Jan 21, 2011this is designed to run in a folder with a single fasta file and blast txt file that has been run blasting that fasta file against itself, with the number of descriptions > the number of times a given ID is hit.
ex. script:nohup pb blastn -dbnuc InputFasta_length_added.fas -query InputFasta_length_added.fas -out InputFasta_length_added75.txt -num_descriptions 75 -evalue 0.00001 -outfmt 10&
note that the fasta headers for the file being used CCDS are unique
Note that the input fasta file needs to be screened to make sure the is only onegene ID for each entry. I do this by adding their length to their end, but there are stillsome replicates, but these are OK since they are the same.