bwa計算wig文件

bwa會輸出mapping不上的reads例如：

A00268:187:H5HVYDSXX:2:1101:22824:1532  77      *       0       0       *       *       0       0       TGGGCCTCGGCCAGCGGGTTGATCAGCACCTGCTCGTTGCGCAGGTCCATCTTGCGGTTGACCTTGGGCGGCCTTTTGGCGACGAAGACGCCGTTCTCGACATAGAGGTTGCCACGCTGCAGCTGCCGGCCATCGAACCAGCGGGCGTTC  FF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF  AS:i:0  XS:i:0

77 = 64+8+4+1(first segement+next segement in the template unmapped +segement unmapped+paired)

A00268:187:H5HVYDSXX:2:1101:22824:1532  141     *       0       0       *       *       0       0       CCGCGACTTGCCCGCTGAAGCCTCCCGCCCTTGCGATTTCTGATCCCCCTCCTGCTGCTCGCCAGCGGCCTGGCCCAGGCACAGGCGACCGGGCAGCCCTACCAGGCCAGCGCCCCGCTGTCCTGGGAACTGAAGAACGCCCGCTGGTTC  FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF  AS:i:0  XS:i:0

因此會導求出來的wig文件有 * 這個做爲scaffold，致使轉化bw的時候出問題