Listing fasta MSAs with certain taxa/header counts
0
0
Entering edit mode
4.0 years ago
herber4 ▴ 30

Hello all,

I am working on a phylogenetic analysis with 23 bacterial genomes.

I generated the gene clusters with panaroo and am now making gene trees for individual clusters.

The issue is that I am working with the single copy orthologues and have a directory full of .aln.fas multiple sequence alignments. The directory has thousands of alignments in it, and i only need the alignments with a gene sequence (header) for each taxa (n=23).

is there a one liner capable of grabbing all of the .aln.fas files with only 23 headers and piping them in to a list?

Therefore allowing me to easily put them into a new directory and generate the gene trees.

best,

AAH

gene genome next-gen alignment fasta • 393 views
ADD COMMENT

Login before adding your answer.

Traffic: 1182 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6