This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Assembly based structural variants call from bacterial genome

Hi,

I am trying to call structural variants from assembled bacterial genome. The genome is assembled using Nanopore long reads and is polished with Illumina reads using Pilon. The coverage of nanopore reads is ~80X. I wonder if assembly based structural variants call would be better than mapping based approaches. Does anyone have suggestion based on their experiences?

Thanks

nanopore structural variants assembly
• 1,601 views • link
written by rthapa • rthapa

Assuming your assembly is correct and complete then whole genome alignment methods will be more accurate than read mapping methods. I have tried both on S. cerevisiae genomes and found this.

You can try my tool MUM&Co or Syri for a full range of SV types.

• 0 views • link
written by samuel.a.odonnell samuel.a.odonnell

Thank you. Is it possible to call SV by comparing multiple genomes with reference genome using MUM&Co? The script can take only single query genome at a time.

bash mumandco.sh -r ./yeast.tidy.fa -q ./yeast_tidy_DEL100.fa -g 12500000 -o DEL100_test
• 0 views • link
written by rthapa • rthapa

I've replied on the github. If you are looking to generate a multiple whole genome alignment and call SVs then you may want to look at the vg toolkit for genome graph variant calling

• 0 views • link
written by samuel.a.odonnell samuel.a.odonnell

0 answers

No answers yet.

Log in to answer this question.

More posts like this