在开始之前，先得做一些基础配置。
安装Rust工具链
curl --proto '=https' --tlsv1.2 -sSf https://sh.rustup.rs | sh
我使用VScode进行开发，并安装了如下插件
- rust-analyzer
- rust syntax
- crates
- better toml
- rust test lens
- Tabnine
接着，尝试使用Rust构建第一个程序。我的目标是读取Fasta，统计其中的A, T, C, G各有多少个。
完成这个需求，有两个思路，第一个是调用成熟的库，第二个是自己造轮子写一个函数。
对于思路1，我使用关键词，Rust和Bio搜索到一个库叫做，rust-bio。 https://rust-bio.github.io/
接下来，使用cargo new创建项目
cargo new fasta_count
第二步: 添加依赖。 编辑 Cargo.toml , 添加依赖信息
[dependencies]
bio = "0.37.1"
第三步: 学习rust-bio文档，查找相关函数，然后在 src/main.rs中编写代码。
use bio::io::fasta;
use std::io;
fn main() {
let reader = fasta::Reader::new(io::stdin());
let mut nb_a = 0;
let mut nb_t = 0;
let mut nb_c = 0;
let mut nb_g = 0;
for result in reader.records(){
let record = result.expect("Error during fasta record parsing");
for &base in record.seq() { //return sequence
if base == b'a' || base == b'A' {
nb_a += 1;
} else if base == b'c' || base == b'C'{
nb_c += 1;
} else if base == b'g' || base == b'G'{
nb_g += 1;
} else if base == b't' || base == b'T'{
nb_t += 1;
}
}
}
println!("A:{}, C:{}, G:{}, T:{}", nb_a, nb_c, nb_g, nb_t);
}
我使用标准库io::stdin从标准输入中读取数据，然后利用fasta::Reader的new出一个新的reader实例。对该实例的records对象进行遍历，得到fasta的中记录，record。最后对记录中序列进行遍历，得到不同碱基的数目。
第四步：编译代码。cargo会先将依赖环境下载到当前项目下，然后在进行编译。
cargo build
编译结果存放在 target的debug目录下，名为fasta_count。我们运行程序，得到结果。
cat TAIR10.fa | target/debug/fasta_count
# A:38223602, C:21551439, G:21528650, T:38177852
第一种思路的代码，是我看了极客时间「陈天·Rust编程第一课」中第一个示例"使用HTTP获取HTML然后保存成Markdown"后自己尝试的第一个和生信相关程序。对于第二种思路，我个人还是决定后续学习更多Rust基础再来尝试，因为我目前都不会用基础库来读取文件。
目前我的初体验是，Rust的语法更接近C/C++，它的编译器非常强大，能够很好的指出的我代码中的错误，便于检索debug。
