三代测序数据比对工具lra
论文
lra: A long read aligner for sequences and contigs
https://journals.plos.org/ploscompbiol/article?id=10.1371/journal.pcbi.1009078
https://github.com/ChaissonLab/LRA
安装直接使用conda
conda install -c bioconda lra
参考序列构建索引
lra index -ONT pome_mt_potential.fa
比对
lra align -ONT pome_mt_potential.fa candidate_mt_longreads.fasta -t 8 -p s > output.sam
samtools sam转bam
samtools view -S -b -o output.bam output.sam
bam sort
samtools sort output.bam -o output.sorted.bam
提取map上的reads
参考 https://gist.github.com/davfre/8596159
samtools view -F 0x04 -b output.sorted.bam > output.aligned.bam
本文参与 腾讯云自媒体同步曝光计划，分享自微信公众号。
原始发表：2022-03-23，如有侵权请联系 cloudcommunity@tencent.com 删除
评论
登录后参与评论
推荐阅读
目录
腾讯云开发者
