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 &gt; the number of times a given ID is hit.<br>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&amp;<br>note that the fasta headers for the file being used CCDS are unique <br>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.