Hi Yilong,
Thanks so much for your reply.
I tried to run 'fakemut' from maq using this command.
Code:
maq fakemut chr22.fa > out.fakeref.fasta 2> out.fake.snp
Do you know what does each column means from "out.fake.snp"?
I can guess, but I am not quite sure.
I can't seem find any explanation about it from the web?
Code:
chr22 16051056 G c 99
chr22 16054178 - T 99
chr22 16055216 C G 99
chr22 16056830 A C 99
chr22 16056930 A G 99
chr22 16057491 G A 99
chr22 16058778 C t 99