The storage and data transferring of large genome data are becoming important concerns for biomedical researchers. We present a novel multi-reference based genome compression method with a hierachical structure. Our approach works for the de facto standard alignment format (i.e., BAM) compression that is the pressing need at present. We align new sequences to a reference sequence using SOAP3, a GPU-based aligning software, and summarize mapping properties and information for exact mapped reads. To increase the exact aligning rate, we also realign the approximately mapped and unmapped reads by changing the reference sequence or shortening the read length. ...