I currently have where you can make a gene tree based off an HMM, and it gives that tree only if the marker is present in those genomes, which is nice for looking at distribution/evolution of that marker. I also do presence/absence of a suite of metabolic markers. I like Mike Lee's example of showing presence/absence with a highlighted tree of all genomes, so you can broadly see where the given marker is NOT located.
Steps:
Can pass an assembly accession file to ncbi-genome-download. Can also get metadata from here
Calling ORFs/annotations because just pulling down nucleotide genbank files
Search for the marker across proteins
Make ribosomal protein tree of all genomes in the set
Highlight with color the presence/absence, and also give an output of # among that clade (such as phyla)
Add a downloading step and ORF calling. Currently assume just pass the .faa files
Add step that for a given marker, makes the ribosomal tree of only those genomes that have the marker, can do this based on the list that comes out the .faa file, and then can compare phylogenies
I currently have where you can make a gene tree based off an HMM, and it gives that tree only if the marker is present in those genomes, which is nice for looking at distribution/evolution of that marker. I also do presence/absence of a suite of metabolic markers. I like Mike Lee's example of showing presence/absence with a highlighted tree of all genomes, so you can broadly see where the given marker is NOT located.
Steps: