KREPA6 functions in RNA editing catalytic complex structural organization and gRNA utilization in T. brucei

Custom software package PARERS v2.0 (Pipeline for Analysis of RNA Edited RNA Sequences), was written using Python 3.12. This pipeline pairs and merges read sequences using BBMerge. Merged reads are checked for matching forward and reverse primer sequences so that any reads with mismatched primers can be excluded from downstream analyses. It can also take advantage of Unique Molecular Identifier (UMI) barcodes for the purposes of deduplication, so that only one instance of each UMI pair can be analyzed. Non-T (A/C/G) sequences are aligned to the corresponding A/C/G maxicircle sequence and any sequences containing a non-T mismatch can be excluded from analysis. Sequence alignments are generated by PARERS v2.0, and data output includes bar graphs and bubble plots made with R version 4.4.0.定制化软件包PARERS v2.0（全称RNA编辑RNA序列分析流程，Pipeline for Analysis of RNA Edited RNA Sequences）基于Python 3.12编写完成。本流程采用BBMerge工具对测序读段序列进行配对与合并；随后对合并后的读段开展正反向引物序列匹配性校验，将引物序列不匹配的读段排除在下游分析之外。此外，该流程可借助唯一分子标识符（Unique Molecular Identifier, UMI）条形码实现去重，仅保留每个UMI对的唯一实例用于后续分析。针对非T（A/C/G）序列，流程会将其与对应的A/C/G大环（maxicircle）序列进行比对，将存在非T错配的序列排除出分析范围。PARERS v2.0可生成序列比对结果，最终输出的数据包含基于R 4.4.0版本绘制的柱状图与气泡图。