DESeq_out_liv_vs_skelm.txtという名前で格納されている。替えるとDEGが抽出できる。比較ではFDRが十分に小さいものが得られなかった。ここまででDEG抽出は終了である。すると以下の図が描ける。左が肝臓で、右が骨格筋のものである。発現量(リード数)を多い順に並び替えてプロットした。遺伝子ばかりが良く発現しているこの図は納得がいく。ひとつ、ハエの脂肪体のデータを載せる。章ではえゲノムのカバレッジを求めた。読み込ませてグラフを描くだけである。つけてRが読み込めるようにする。