论文
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
欢迎大家关注我的公众号
小明的数据分析笔记本
小明的数据分析笔记本 公众号 主要分享：1、R语言和python做数据分析和数据可视化的简单小例子；2、园艺植物相关转录组学、基因组学、群体遗传学文献阅读笔记；3、生物信息学入门学习资料及自己的学习笔记！