2018mapping例

makeblastdb -in dm6.fa -dbtype nucl

blastn -outfmt 6 -db dm6.fa -query read.fa

#low complexity な配列もヒットさせたり、ナノポアのようにクオリティの低いリードもマッチさせたい場合
blastn -outfmt 6 -num_threads 4 -word_size 10 -dust no -db dm6.fa -query read.fa
lastdb dm6.fa dm6.fa

lastal -f BlastTab+ dm6.fa read.fa

#10回以上のリピートにもヒットさせたくて、初期ヒットの上限を1万個まで上げる場合
lastal -m 10000 -P 4 -f BlastTab+ dm6.fa read.fa
bwa index dm6.fa

bwa mem dm6.fa read.fq 
bowtie2-build dm6.fa dm6.fa

bowtie2 -x dm6.fa -U read.fq

#ローカルアライメントモードでヒットした箇所全てを出力する場合
bowtie2 -a --local -p 4 -x dm6.fa -U read.fq
bowtie2-build dm6.fa dm6.fa

tophat2 dm6.fa read.fq
hisat2-build dm6.fa dm6.fa

hisat2 -x dm6.fa -U read.fq
minimap2 -a -x map-ont dm6.fa read.fq
minialign -xont dm6.fa read.fq
  • 2018mapping例.txt
  • 最終更新: 2019/07/15 02:46
  • by 127.0.0.1