Question: What seems to work best to get consensus sequence after mapping: variant calling or de-novo assembly?
1
gravatar for Kevin D
3.6 years ago by
Kevin D30
INRA France
Kevin D30 wrote:

Hi everyone!

I have a question for the community and I wish I have your opinion and some piece of advise about it. I have a genomic sequence from a plant species on which I mapped Illumina paired reads (2*150 bp). Those reads come from another species from the same genus. Both species are closely-related and I expect some SNP and indels. I perfomed the mapping with BWA. Now I got the bam file. I wonder if I should do variant calling using the combined samtools mpileup/bcftools tools or extract the mapped reads and try a de-novo assembly with those reads using SOAP-denovo2 for instance... In one hand, SOAP will take the insert size into account but it is harder to set. On the other hand, samtools mpileup looks easier to handle but may be less accurate (?) I don't really know...

At the end, I would like to get the genomic sequence of the new species in fasta format.

NB : I prefer not to perform a whole genome assembly since the reads quantity may not be sufficient to get an interesting coverage.

Thanks in advance,

Kevin

ADD COMMENTlink written 3.6 years ago by Kevin D30
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1944 users visited in the last hour