bedtools getfasta -s -fi $GENOME -bed genes.bed -fo genes.fa -name 3.提取基因组中的一段序列 比如我只想提取一号染色体上的10000-50000这一段序列。 首先建立索引。 samtools faidx IRGSP-1.0_genome.fasta 然后根据染色体信息和物理位置直接提取。这里注意,如何要和目的4连用的话,要修改>后面的值与gff...
么有gff文件是无法提取基因序列的，gff文件记录了基因的位置信息，没有位置信息无法提取序列。
这不是很简单的问题吗？如果你能回答上面的问题的话