由於語法渲染問題而影響閱讀體驗, 請移步博客閱讀~
本文GitPage地址

Processing Sanger Sequencing file (abi) wiht Python or R

In R, you can using sangerseqR to reading ab1 file.
An example can see: Karobben 2020

Python

Biopython

  1. from Bio import SeqIO
  2. handle = open("test.ab1", "rb")
  3. for record in SeqIO.parse(handle, "abi"):
  4. print(record)

Trace plot

Biopython: Click hear

Package sanger-sequencing

Link
Can’t find any examples…

Perl

ABIF

Convert abi to Scf

reference: stock overflow

  1. convert_trace -out_format scf < trace.ab1 > trace.scf

Convert ab1 to excel

Working Manual: ab1_organizer

  1. pip install ab1_organizer
  2. ab1_organizer.py -f ./path/to/zip/file -t ./path/to/order/table.xlsx

Base Calling

GitHub: tracy
Documentation
Paper

(PS: It seems like we need to download some files from google. So, as you can see, it’s not easy for the people in China mainland to install it.)
It is a great app for basecalling, alignment, assembly and deconvolution of sequencing chromatogram files.

Installation

  1. apt install \
  2. build-essential g++ \
  3. cmake \
  4. git-all \
  5. liblzma-dev \
  6. zlib1g-dev \
  7. libbz2-dev \
  8. liblzma-dev \
  9. libboost-date-time-dev \
  10. libboost-program-options-dev \
  11. libboost-system-dev \
  12. libboost-filesystem-dev \
  13. libboost-iostreams-dev
  14. git clone --recursive https://github.com/gear-genomics/tracy.git
  15. cd tracy/
  16. make all
  17. make install
  18. ./bin/tracy -h

Enjoy~

本文由Python腳本GitHub/語雀自動更新

由於語法渲染問題而影響閱讀體驗, 請移步博客閱讀~
本文GitPage地址

GitHub: Karobben
Blog:Karobben
BiliBili:史上最不正經的生物狗