#bioinformatics
User can start by running ```1181102378_BLAST_project.py```.This program accepts a FASTA file from the user and does 3 things:
- run BLAST on the sequence,
- perform multiple sequence alignment on the top 10 homologues,
- extract the positions of non-conserved regions.
A FASTA file is provided for the users (HBA2.fasta) and has been written in the input message to prompt the user to enter the correct file name. If there are ambiguous codes in the sequence, like letters other than A, T, G and C, the program will terminate. Else, all 3 functions will run, and the required files will be generated in the output/ folder. Contained in the folder sample output/ are what the output should be.
The program output are:
qblast_blastp.xml: the top 50 homologues of the user input sequence generated from BLASThitlist.fasta: the top 10 homologues extracted from the xml filealigned.fasta: the output of multiple sequence alignmentnonconserved-regions.txt: the positions of nonconserved regions will be symbolised with a "," instead of a "*" underneath the sequences, the first sequence is the query sequence, followed by the rest of the subject sequences.
There are several files included in the folder that are meant for the multiple sequence alignment program. These include clustalo.exe and its dependent .dll files.