pub fn write_gbk_format_sequence(sequence: &str, file: &mut File) -> Result<()>
writes the DNA sequence in gbk format with numbering