genomewalker
8/7/2015 - 11:49 AM

get_orf_coverage_from_bam-coAssembly.awk

#Read a hist file from bedtools coverage -hist and create a coverage per orf file.
#Coverage per each gene is calculated as: coverage_orfA = sum(depth_of_coverage * fraction_covered)
{
    if ($0 i ~ /^all/){
        next;
    }else{
        split($4,a,"_");
        b=$1"_"$6"_"$2+1"_"$3"_orf-"a[2]"\t"$13;
        c[b]=c[b] + ($11*$14)
    }
}END{
    for (var in c){
        print var"\t"c[var]
    }
}