In silico region of difference (RD) analysis of Mycobacterium tuberculosis complex from sequence reads