Open WeiWei1112 opened 6 years ago
You may want to check the assembly graph to see what is around that contigs and why there were not assembled.
You may want to check the assembly graph to see what is around that contigs and why there were not assembled.
Thanks for your reply! Do you know what does the "--careful" option do?
Sure. Per SPAdes manual (http://cab.spbu.ru/files/release3.13.0/manual.html):
--careful
Tries to reduce the number of mismatches and short indels. Also runs MismatchCorrector – a post processing tool, which uses BWA tool (comes with SPAdes). This option is recommended only for assembly of small genomes. We strongly recommend not to use it for large and medium-size eukaryotic genomes
确定。根据 SPAdes 手册 (http://cab.spbu.ru/files/release3.13.0/manual.html):
--careful Tries to reduce the number of mismatches and short indels. Also runs MismatchCorrector – a post processing tool, which uses BWA tool (comes with SPAdes). This option is recommended only for assembly of small genomes. We strongly recommend not to use it for large and medium-size eukaryotic genomes
确定。根据 SPAdes 手册 (http://cab.spbu.ru/files/release3.13.0/manual.html):
--careful Tries to reduce the number of mismatches and short indels. Also runs MismatchCorrector – a post processing tool, which uses BWA tool (comes with SPAdes). This option is recommended only for assembly of small genomes. We strongly recommend not to use it for large and medium-size eukaryotic genomes
Hello, if I'm assembling MDA data, is it better to use the --careful option ?
Hi, I want to use Spades to assemble a ~40kb plasmid sequenced with Illumina Miseq (2X250nt). When I use the default option:
spades.py -1 ../BWA/B13_vector_unmapped.R1.fastq -2 ../BWA/B13_vector_unmapped.R2.fastq --phred-offset 33 -o B13
I get many several small contigs like below:But if I add the "--careful" option, I can get
>B13_GGCTTAAG-TCGTGACC_1_length_34210_cov_184.621189
Apparently the result is much better with "careful" option. But the question is that do I throw away possible true variants if I just correct mismatches and indels? My sequence contains a gene cluster in which the genes might be similar but slightly different to each other. Do you have any suggestions for the option settings? Thank you so much!
Wei Wei