huget 0.1.0

gProfiler for human genomics
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55

use crate::readstruct::ReadSNP;
use reqwest::blocking::get;
use scraper::{Html, Selector};
use std::error::Error;
use std::fs::File;
use std::io::{BufRead, BufReader};

/*
Gaurav Sablok
codeprog@icloud.com
*/

#[tokio::main]
pub async fn mapid(pathstring: &str, mapsnpid: String) -> Result<String, Box<dyn Error>> {
    let filesnpopen = File::open(pathstring).expect("file not present");
    let filesnpread = BufReader::new(filesnpopen);

    let mut opensnpstruct: Vec<ReadSNP> = Vec::new();

    for i in filesnpread.lines() {
        let line = i.expect("line not present");
        if !line.starts_with("#") {
            let linestore = line
                .split("\t")
                .map(|x| x.to_string())
                .collect::<Vec<String>>();
            opensnpstruct.push(ReadSNP {
                rsid: linestore[0].clone(),
                chromosome: linestore[1].clone(),
                position: linestore[2].clone(),
                genotype: linestore[3].clone(),
            })
        }
    }

    let identifiedsnp = mapsnpid;
    let formatstring = format!("{}/{}", "https://www.ncbi.nlm.nih.gov/snp/", identifiedsnp);
    let response = get(formatstring).expect("string not found");
    let document = Html::parse_document(&response.text().expect("message not present"));
    let snpselect = Selector::parse(".sect_heading").expect("method failed");
    for element in document.select(&snpselect) {
        let vector_1 = element
            .text()
            .collect::<Vec<_>>()
            .join("-")
            .replace(" ", "")
            .replace("\n", "")
            .replace("Genomecontext", "")
            .replace("Selectflanklength", "")
            .replace(":", "");
        println!("{:?}", vector_1);
    }

    Ok("path to the mapping string has been written".to_string())
}