Сборка de novo

Команды

Команда Значение
velveth work_directory 31 -fastq -short chr2t.fastq Выделение k-меров.
velvetg work_directory Сборка контигов.
makeblastdb -in /nfs/srv/databases/ngs/Human/chr2.fasta -dbtype nucl -out chr2 Создание базы данных из файла хромосомы.
blastn -query contigs.fa -db chr12 -outfmt 6 -out blast.out Поиск по базе.
work_directory - это папка, в которую velvetg кладёт информацию о созданных контигах.
N50: 237
Ссылка на таблицу со всеми находками Ссылка на таблицу с находками, откартированными ровно один раз

Контиги, картировавшиеся ровно один раз

Их много - 255, многие из них перекрываются. Тому, что они не соединились, могут быть объяснением две причины:
1) Длина перекрывания меньше длины k-мера
2) Ошибки секвенирования различных чтений, переходящие в ошибки в контигах