Les différentes méthodes de séquençage à haut débit ne sont pas totalement infaillibles dans leur capacité à faire le dit séquençage! Elles sont toutes basées sur des processus qui à divers moments de la méthode, peuvent engendrer des problèmes comme des échantillons de qualité discutable, des erreurs de lecture, des séquences partielles dans l'ensemble d'un fichier (rappel: chaque fichier contient des millions de séquences), etc. Il nous faudra dans une étape subséquente retirer toute donnée de séquence qui pourra perturber les reste de notre analyse; pour établir les critères de filtration, il nous faut regarder le profil de nos données brutes. Quoi utiliser? Dans la communauté, dans le moment, les deux applications les plus utilisées sont FastQC et fastP, qui toutes les deux vont créer des rapports sous la forme de pages web qui permettront de vous donner un aperçu de vos données et qui vous permettront de déterminer les meilleurs critères de filtration dans l'étape suivante.

Pour une première fois, on a une situation où l'on doit faire à répétition la même opération sur tous les fichiers FASTQ que nous avons téléchargé. On pourrait faire 32 fois la même opération mais, soyons sérieux, c'est à ça que sert une grappe de calcul et Python Pourquoi ne pas utiliser Python pour écrire tous les fichiers bash à soumettre à SLURM? À la base, il faudrait encoder statiquement tous les paramètres à soumettre dans notre script, ce qui le rend peu flexible… Introduisons un autre outil: un fichier de configuration à lire par notre script. Comment écrire ce fichier de paramètres? En utilisant un fichier texte écrit dans un format structuré qui pourra être lu par le script pour construire une structure de données directement utilisable. Un tel format est le format YAML, très flexible et facilement interprété par Python.

• Plus à venir…