Convertir FASTA en FASTQ

Comprendre les formats de fichiers FASTA et FASTQ

FASTA est un format basé sur du texte largement utilisé pour représenter des séquences de nucléotides ou de peptides. Chaque séquence dans un fichier FASTA commence par une ligne de description unique (qui commence par >), suivie par des lignes de données de séquence. FASTQ est un autre format utilisé en bioinformatique, qui étend le format FASTA en incluant des scores de qualité pour chaque nucléotide, ce qui le rend essentiel pour stocker les données brutes de séquençage provenant d'instruments de séquençage à haut débit.

Pourquoi convertir FASTA en FASTQ ?

La conversion de FASTA en FASTQ est nécessaire lorsque les outils en aval requièrent des scores de qualité, même si seuls des scores fictifs ou uniformes sont attribués. Cela est courant lors de la simulation de données de séquençage ou la préparation de données pour certains pipelines bioinformatiques.

Comment convertir FASTA en FASTQ

Étant donné que les fichiers FASTA ne contiennent pas d'informations sur la qualité, la conversion en FASTQ implique généralement l'attribution d'un score de qualité par défaut à chaque base. Ceci peut être réalisé à l'aide d'outils en ligne de commande ou de logiciels spécialisés.

Logiciels recommandés pour la conversion de FASTA en FASTQ

• Seqtk : Un outil léger et efficace en ligne de commande pour traiter des séquences. Pour convertir, utilisez :
  seqtk seq -F '#' input.fasta > output.fastq
  Cela attribue un score de qualité par défaut ('#', correspondant à un score Phred faible mais valide) à chaque base.
• BBMap (reformat.sh) : Un autre outil polyvalent. Utilisez :
  reformat.sh in=input.fasta out=output.fastq qout=30
  Cela attribue un score de qualité Phred de 30 à toutes les bases.
• Plateforme Galaxy : Pour ceux qui préfèrent une interface graphique, téléchargez votre fichier FASTA et utilisez l'outil FASTA-to-FASTQ, en spécifiant le score de qualité souhaité.

Exemple étape par étape avec Seqtk

1. Installez Seqtk depuis son dépôt GitHub ou via des gestionnaires de paquets.
2. Ouvrez votre terminal et exécutez :
   seqtk seq -F 'I' input.fasta > output.fastq
   Cela attribue un score de qualité élevé ('I', Phred 40) à chaque base.
3. Vérifiez output.fastq pour les séquences converties avec leurs scores de qualité.

Conclusion

La conversion de FASTA en FASTQ est simple avec des outils comme Seqtk ou BBMap. Rappelez-vous, les scores de qualité résultants sont artificiels sauf si vous disposez de données de qualité réelles.