1
$\begingroup$

I have MAG geneome sequence released in NCBI. In the released sequence there are 4915 contigs present. I want to get the full length genome sequence as FASTA file to proceed with further works like circular genome viewing and other molecular typing. Kindly let me know how to proceed in this regard. I am not much aware of bioinformatics so need even the basic information in this regard.

Thanks

$\endgroup$
3
  • $\begingroup$ Please edit the question to limit it to a specific problem with enough detail to identify an adequate answer. $\endgroup$
    – Community Bot
    Mar 1, 2023 at 13:15
  • $\begingroup$ Please explain what genome you are working with in particular and how you have assembled this to date. Circular genome viewing ... bacteria? $\endgroup$
    – M__
    Mar 1, 2023 at 13:19
  • $\begingroup$ Just to note @Laura is an expert on this subject - if its bacteria. This post might help bioinformatics.stackexchange.com/questions/20336/… $\endgroup$
    – M__
    Mar 2, 2023 at 13:34

1 Answer 1

1
$\begingroup$

Simple answer is, it might not be possible to do it.

The large number of contigs suggests to me that a short-read sequencing method (such as Illumina sequencing) was used. Alone, short-read sequencing technologies cannot resolve repetitive sequences that are longer than the length of the reads. This is why you have a large number of contigs.*

*You might already have done this, but if you did not, it might be worth checking that your MAG is not actually a mixture of different genomes. This can also somewhat increase the contig number and, more importantly, can skew your interpretation if you assume it's only one organism. There are various ways to test this, but the simplest I know is the program CheckM. To run it you will need a way to use Unix/Linux (could be on a virtual machine). If you are not familiar with Unix, bash and conda, I would heartily suggest software carpentry lessons such as The Unix Shell (https://swcarpentry.github.io/shell-novice/) and possibly this conda lesson although it is in beta and I haven't tested it https://carpentries-incubator.github.io/introduction-to-conda-for-data-scientists/

I'm not familiar with the term 'molecular typing', but it looks like you're interested in what organisms / genes there are in your sample? You probably might not need a full genome for that.

Small probability of helping - but it might also be good to have a look at the sizes of your contigs. You might luck out and have a few large contigs (hundred thousand - million basepairs) which could be compared to existing database genomes.

$\endgroup$

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.