This is a test version of Biostars. For the public version, visit https://www.biostars.org.
After Novo assembly

After assembling (novo assembler ) a sequence of Tuberculosis I found that the number of base pair is greater than that number of base pair is the reference strain Is this logical?

assembly

Is this logical?

Yes.

And computational biology is a quantitative science. Please tell us exactly how much greater your assembly is and which reference genome you have used.

i used Mycobacterium tuberculosis H37Rv this reference has 4.4Mb and my sequence has 7.3 Mb?

Have you tried to compare the two to see how the assemblies are different? Use Mauve to compare.
Did you have an excess of sequence (> 100x gross coverage) that went into this assembly?

Yes, you should map your contigs to the refseq and then identify contigs which do NOT map. Blast these contigs to identify their origin. Furthermore, your may plot GC-content of your contigs versus coverage. Do you see more than one cluster in the scatter plot?

1 answer

There could be any number of reasons resulting in an increase in the number of bases - duplication at various levels, insertions, etc. Why is comparing base count any kind of metric?

Why is comparing base count any kind of metric?

Tuberculosis genomes are extremely well conserved, much more than any other bacterial species.

Ah, I see. Thank you :)

Log in to answer this question.