rasusa 4.1.0

Randomly subsample reads or alignments
1
2
3
4
5
6
7
8
9
10
11
12
13
14

use noodles_util::alignment;

#[test]
fn test_get_len() -> Result<(), Box<dyn std::error::Error>> {
    let mut reader =
        alignment::io::reader::Builder::default().build_from_path("tests/cases/test.bam")?;
    let header = reader.read_header()?;
    for result in reader.records(&header) {
        let record = result?;
        let l = record.sequence().len();
        println!("len: {}", l);
    }
    Ok(())
}