METODE

    A partir de la nostra seqüència en format FASTA , primer cal passar-la a format tabulat per tal de calcular el número de parell de bases
    mitjançant la comanda:
    awk '{print f "%s", $0}' AC090953.fa > AC090953.tbl

     Un cop s'obté la seqüència  en format tbl  es podem realitzar les comandes per determinar la longitutd de la seqüència i el contingut de g+c

    - Longitud de la seqüència: awk '{print length {$2}}' AC090953.tbl
    - Contingut de G+C: awk '{print $2}' AC090953.tbl | fold -1 | sort | uniq -c | gawk '{print $2, $1/ número de bases de la seqüència }'
 

    Enrera