我這里使用python和biopython堡称，關(guān)于python和Biopython的安裝可以參考我的B站視頻
image.png
輸入文件就是葉綠體基因組的fasta文件瞎抛，經(jīng)過注釋你已經(jīng)知道LSC/IR/SSC的位置坐標，比如自己的葉綠體基因組序列
1——80000是大單拷貝區(qū)
80001——100000是反向重復(fù)區(qū)1
100001——130000是小單拷貝區(qū)
首先是輸出葉綠體基因組的總長度
image.png
輸出總的GC含量
image.png
大單拷貝區(qū)的GC含量
image.png
反向重復(fù)區(qū)
image.png
小單拷貝區(qū)
image.png
這里需要注意的是位置坐標是起始位置-1：終止位置
歡迎大家關(guān)注我的公眾號
小明的數(shù)據(jù)分析筆記本
小明的數(shù)據(jù)分析筆記本 公眾號 主要分享：1却紧、R語言和python做數(shù)據(jù)分析和數(shù)據(jù)可視化的簡單小例子桐臊；2钞艇、園藝植物相關(guān)轉(zhuǎn)錄組學、基因組學豪硅、群體遺傳學文獻閱讀筆記哩照；3、生物信息學入門學習資料及自己的學習筆記懒浮！
