bwa计算wig文件

bwa会输出mapping不上的reads例如:app

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

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

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

因此会导求出来的wig文件有 * 这个做为scaffold,致使转化bw的时候出问题ping