Soy nuevo en esto, pero sé que para muchas secuencias es más fácil usar Perl para editar las cabeceras FASTA para evitar errores en análisis posteriores. Mi problema es no saber cómo meter una línea de comandos. Tengo archivos de entrada:
>gi|343457071|gb|JF424046.1| Streptomyces abikoensis strain KCTC 9662 recombinase A (recA) gene, partial cds
TCCACCGGGTCGACCGCTCTCGACGTCGCGCTCGGTGTCGGCGGCCTGCCGCGCGGCCGCGTGGTGGAGA
TCTACGGACCGGAGTCCTCCGGTAAGACGACCCTGACGTTGCACGCCGTGGCCAACGCCCAGCGGGCCGG
CGGCACCGTCGCCTTCGTGGACGCCGAGCACGCCCTCGACCCCGAGTACGCCAGAAAGCTCGGCGTCGAC
ATCGACAACCTGATCCTTTCCCAGCCGGACAACGGCGAGCAGGCTCTCGAGATCGTCGACATGCTGGTCC
>gi|78173104|gb|DQ234054.1| Streptomyces argillaceus strain ATCC 12956 RecA (recA) gene, complete cds
CCATGGCAGGCACCGACCGCGAGAAGGCCCTGGACGCCGCACTCGCACAGATTGAACGGCAATTCGGCAA
GGGCGCGGTCATGCGCATGGGCGACCGCTCGAAGGAGCCCATCGAGGTCATCCCGACCGGGTCGACCGCG
CTCGACGTGGCCCTCGGCGTCGGCGGTCTGCCGCGCGGCCGTGTCATCGAGGTCTACGGACCCGAGTCCT
CCGGCAAGACGACCCTGACCCTGCACGCGGTGGCGAACGCCCAGAGGGCCGGCGGCCAGGTGGCGTTCGT
GGACGCCGAGCACGCCCTCGACCCCGAGTACGCGCAGAAGCTCGGCGTGGACATCGACAACCTGATCCTG
TCCCAGCCGGACAACGGCGAGCAGGCCCTGGAGATCGTGGACATGCTCGTCCGCTCCGGGGCCCTCGACC
TGATCGTCATCGACTCCGTCGCCGCGCTCGTCCCGCGTGCGGAGATCGAGGGCGAGATGGGCGACAGCCA ...