我需要从各个基因组序列数据中得出一个简短的串联重复序列(STR)。我正在尝试使用BioPython。当我搜索任何STR等位基因或序列时,我想知道该序列在FASTA文件中的位置(该序列所在的染色体编号,以及该序列开始和结束的碱基对位置),但是我不清楚该怎么做。使用BioPython。有人告诉我FASTA文件不包含染色体信息,那么我还能做什么?我也尝试过CRAM / BAM文件,但是BioPython不接受这些文件。我只想打印该FASTA文件的每个序列的位置。
我该怎么做?
我需要从各个基因组序列数据中得出一个简短的串联重复序列(STR)。我正在尝试使用BioPython。当我搜索任何STR等位基因或序列时,我想知道该序列在FASTA文件中的位置(该序列所在的染色体编号,以及该序列开始和结束的碱基对位置),但是我不清楚该怎么做。使用BioPython。有人告诉我FASTA文件不包含染色体信息,那么我还能做什么?我也尝试过CRAM / BAM文件,但是BioPython不接受这些文件。我只想打印该FASTA文件的每个序列的位置。
我该怎么做?