Chapitre 16. FASTX grep : Création d'un programme utilitaire pour sélectionner des séquences
Un collègue m'a demandé un jour de trouver toutes les séquences d'ARN dans un fichier FASTQ dont la description ou le nom contenait la chaîne LSU (pour long subunit RNA). Bien qu'il soit possible de résoudre ce problème pour les fichiers FASTQ en utilisant le programme grep 1 pour trouver toutes les lignes d'un fichier correspondant à un certain motif, écrire une solution en Python te permet de créer un programme qui pourrait être étendu pour gérer d'autres formats, comme FASTA, ainsi que pour sélectionner des enregistrements en fonction d'autres critères, comme la longueur ou le contenu GC. En outre, tu peux ajouter des options pour modifier le format de la séquence de sortie et introduire des commodités pour l'utilisateur, comme deviner le format du fichier d'entrée en se basant sur l'extension du fichier.
Dans ce chapitre, tu apprendras :
-
A propos de la structure d'un fichier FASTQ
-
Comment effectuer une correspondance d'expression régulière insensible à la casse ?
-
A propos des idées DWIM (Do What I Mean) et DRY (Don't Repeat Yourself) dans le code
-
Comment utiliser les opérations
andetorpour réduire les valeurs booléennes et les bits ?
Recherche de lignes dans un fichier à l'aide de grep
Le programme grep peut trouver toutes les lignes d'un ...
Become an O’Reilly member and get unlimited access to this title plus top books and audiobooks from O’Reilly and nearly 200 top publishers, thousands of courses curated by job role, 150+ live events each month,
and much more.
Read now
Unlock full access