what is the best coverage that would be suitable to carryout denovo assembly (velvet) on illumina resequencing reads- would 20x be ok?
-
-
I've found Velvet to work best at around 40x, using 150bp reads. But it depends somewhat on the read length; with longer reads you may be able to get away with less coverage.Originally posted by mmmm View Post -
thank you- I have 250PE and 300PE reads (but of lower coverage 20x) so will try velvet and see
I was wondering is it ok to combine fastq files of 250PE and 300PE redas (of the same sample) using cat command, would this combination affect the analysis (either velvet or bwa??Comment
-
for velvet it might be better to leave the 2 sets of reads separate, in case the library prep gave a different insert length for the different runs. You can still use both sets of reads in the same run of velvet though, just flag one set as
'-shortPaired', and the second set as '-shortPaired2' when running velveth, but make sure you have velvet compiled with at least 'CATEGORIES=2'.Comment
-
thank you- will give it ago
so if I have 3 diff PE reads (for the same samples but from diff runs)- I will make CATEGORIES = 3
velveth output/ kmer -fastq -shortPaired -separate run1.read1.fastq run1.read2.fastq -shortPaired2 -sparate run2.read1.fastq run2.read2.fastq and so on?Comment
-
Researchers using sequencing and genomics tools often have to make trade-offs. They can choose between speed or scale, short reads or long-range information, or targeted panels or a view of the whole transcriptome. New technologies that have been released this year are built to address those tough choices.
We asked six companies the same four questions to learn about their latest products. The new technologies bring a lot to the table, including rethinking sequencing...
Comment