Question: Fasta file substitution
1
max1415r • 10 wrote:
I have a batch protein FASTA file. I wan to create output file with position for each amino acid with possible mutation (remaining 19 aa) in tab-delimited file.
eg. sequence
>sp|Q6NUK1|SCMC1_HUMAN Calcium-binding mitochondrial carrier protein
MLRWLRDFVLPTAACQDAEQPTRYETLFQALDRNGDGVVDIGELQEGLRNLGIPLGQDAE
>sp|Q6KCM7|SCMC2_HUMAN Calcium-binding mitochondrial carrier protein
MLCLCLYVPVIGEAQTEFQYFESKGLPAELKSIFKLSVFIPSQEFSTYRQWKQKIVQAGD
Output file; (Protein ID, position, amino acid, substitution(19). It will be given for all protein(around 4000)
Q6NUK1 | 1 | M | A |
Q6NUK1 | 1 | M | R |
Q6NUK1 | 1 | M | N |
Q6NUK1 | 1 | M | D |
Q6NUK1 | 1 | M | C |
Q6NUK1 | 1 | M | Q |
Q6NUK1 | 1 | M | E |
Q6NUK1 | 1 | M | G |
Q6NUK1 | 1 | M | H |
Q6NUK1 | 1 | M | I |
Q6NUK1 | 1 | M | L |
Q6NUK1 | 1 | M | K |
Q6NUK1 | 1 | M | F |
Q6NUK1 | 1 | M | P |
Q6NUK1 | 1 | M | S |
Q6NUK1 | 1 | M | T |
Q6NUK1 | 1 | M | W |
Q6NUK1 | 1 | M | Y |
Q6NUK1 | 1 | M | V |