Hello all,
I am trying to identify a human functional homolog of a specific bacteria. I assembled and scaffolded the bacterial genome using SPAdes. I am now trying to align the assembled scaffolds onto a human-derived Lactobacillus species genome available on NCBI to identify a human-derived species who is a functional homolog of the assembled genome.
My question is, is there an efficient way which I can find a published genome which would be the closest to my assembly for the purpose of alignment?
I tried running a BLAST, however, I was told the query was too long. The total length I'm working with is 2.1 million bp.
Thank you in advance!
•
written
by
Kacper •
1 answer
Find all lactobaccilus genomes put all genomes into a single file, run minimap2 on your contigs, and align against the combined genome.
Finally summarize the BAM with idxstats that will tell you which genome matches the most often.
Some from command line but you can also download the genomes from the NCBI website, see the FAQ
Log in to answer this question.
-
Species identifcation using 16s rRNA
written by A_heath •
-
how to annotate draft genome
written by zhangdengwei
-
chimeric assembly of a mitogenome
written by vicmansj86 •
-
Genome annotation using COG
written by Paul
-
Genome Annotation of a new bacterial strain
written by Paul
-
Submission of a draft genome with or without repetitive contigs?
written by minions-b •
-
Problem with WGS
written by aggsandy8 •
-
characterizing a plasmid genome
written by SJ Basu •
-
Augustus: Choosing Species
written by emmajosenl •
-
SPAdes insert size
written by treitlis •
