Dear all I have a draft genome sequence having 18 contigs which I have concatinated and found operons using softberry module fgenes. Likewise, I have one reference (complete genome), of similar bacteria and i have identified operons in this genome too. Now the question is, if I have some genes such as essential genes and I know the position (coordinates) of those genes in both draft and comeplete genome, then I want to compare the coempelete operon of particular gene in both reference and draft genome to know whether it occurs at the correct operon or not. I have to perform this analysis for arounf 70 genes. Can i use multigene blast? but how? Any help will really be appreciated.