SortMeRNA 去除rRNA

单端测序文件

for i in *IP*fastq; do

 ~/Software/sortmerna-2.1b/sortmerna --ref /home/zhou/Software/sortmerna-2.1b/index/human_rRNA.fa,/home/zhou/Software/sortmerna-2.1b/index/human_rRNA_db: --reads $i --num_alignments 1 --fastx --aligned reads_rRNA --other ${i}_rmrRNA.fq --log -a 8 -m 24000 -v;

 done


双端测序

合并read files

/path_to_sortmerna/scripts/merge-paired-reads.sh reads_1.fq reads_2.fq reads_interleaved.fq

去除rRNA

 ~/Software/sortmerna-2.1b/sortmerna --ref /home/zhou/Software/sortmerna-2.1b/index/human_rRNA.fa,/home/zhou/Software/sortmerna-2.1b/index/human_rRNA_db: --reads $i --num_alignments 1 --fastx --aligned reads_rRNA --other ${i}_rmrRNA.fq --log -a 8 -m 24000 --paired_in -v;

分解merged的read flie

/path_to_sortmerna/scripts/unmerge-paired-reads.sh reads_non_rRNA.fq reads_non_rRNA_1.fq reads_non_rRNA_2.fq

你可能感兴趣的:(SortMeRNA 去除rRNA)