Function test_create_plot_data

src/plot.rs:1014–1060

()

    #[test]
    fn test_create_plot_data() {
        let region = Region {
            target: "chr1".to_string(),
            start: 0,
            end: 20,
        };
        let (reads, reference, total_reads, coverage, subsampled_reads) = create_plot_data(
            "tests/sample_1/reads.bam",
            "tests/sample_1/reference.fa",
            &region,
            100,
            None,
            0.0,
        )
        .unwrap();
        let expected_reference = Reference {
            start: 0,
            reference: "TTGCCGGGGTGGGGAGAGAG".to_string(),
        };
        let expected_read = Read {
            name: "sim_Som1-5-2_chr1_1_1acd6f".to_string(),
            cigar: PlotCigar::from_str("16=|iAA|80=|1T|1=").unwrap(),
            position: 4,
            flags: 99,
            mapq: 30,
            row: Some(1),
            end_position: 106,
            mpos: 789264,
            aux: AuxRecord(HashMap::new()),
            raw_cigar: "16M2I82M".to_string(),
        };

        let expected_reads = vec![EncodedRead::from_reads(vec![expected_read])];
        let expected_coverage = Coverage {
            start: 0,
            matches: "0|0|0|0|1|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0".to_string(),
            a: "".to_string(),
            t: "".to_string(),
            c: "".to_string(),
            g: "".to_string(),
        };
        assert_eq!(reference, expected_reference);
        assert_eq!(coverage, expected_coverage);
        assert_eq!(reads, expected_reads);
        assert_eq!(total_reads, 1);
        assert_eq!(subsampled_reads, 1);
    }

    #[test]
    fn test_create_plot_data_2() {