def parse_node_coverage(line):
    # S       s34     CGTGACT LN:i:7  SN:Z:1  SO:i:122101     SR:i:0  dc:f:0
    # "nodeid","nodelen","chromo","pos","rrank",assemb
    """
    Parse the gaf alignment
    Input: line from gaf alignment
    Output: tuple of nodeid, nodelen, start_chromo, start_pos, coverage

    """
    line_comp = line.strip().split()
    nodeid = line_comp[1]
    nodelen = len(line_comp[2])
    start_chromo = line_comp[4].split(":")[2]
    start_pos = line_comp[5].split(":")[2]
    rrank = line_comp[-2].split(":")[2]
    coverage = line_comp[-1].split(":")[2]

    return nodeid, nodelen, start_chromo, start_pos, rrank, coverage