Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

different result with mito annotation #77

Open
tanpham15 opened this issue Feb 19, 2024 · 0 comments
Open

different result with mito annotation #77

tanpham15 opened this issue Feb 19, 2024 · 0 comments

Comments

@tanpham15
Copy link

Dear Marcela Uliano-Silva,

After running mitohifi (in both v3.2 and 3.2.1), I used final mitogenome with mitos (on https://usegalaxy.eu/) to get annotation in bed.

Here is the list parameter I used in:

  • mitohifi: genetic code (-o) 5
  • mitos (usegalaxy) genetic code 5, Reference data: Refseq39+MiTFi tRNA model
  • mitos2 (usegalaxy) genetic code 5, Reference data: Refseq63 Metazoa/ Reference data: Refseq89 Metazoa

Finally, comparation with mitos I saw a slightly different between results of mitohifi and mitos:

1, Mitohifi is better when found all 37 genes, in mitos2 (with both refseq63 and refseq89), they missed "trnH" - translational exceptions: cox1,-1,3375,4910 (CGA, None), with mitos (refseq39) they also found all 37 genes
2, In mitohifi (png), it indicate "tRNA-Ser2", but in all three results with mitos, they found "trnS1(act)".

Could you please let me know how can I get a report for gene length (in bed file)? and can mitohifi make a circle for final mitogenome annotation?

Thank you very much in advance!

Best regards

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant