I have downloaded all fasta files for human chromosomes and merged them with cat. In order to use the "region" option in samtools view command, do I have to specify the exact header of the fasta file?
so for example if I want to analyse only a specific region in chromosome 12:
$ samtools view gi|89161190|ref|NC_000012.10|NC_000012 Homo sapiens chromosome 12, reference assembly, complete sequence:1000000-2000000
I guess it's better then to modify the header to something like ">chr12".
Any suggestions?
Thanks
so for example if I want to analyse only a specific region in chromosome 12:
$ samtools view gi|89161190|ref|NC_000012.10|NC_000012 Homo sapiens chromosome 12, reference assembly, complete sequence:1000000-2000000
I guess it's better then to modify the header to something like ">chr12".
Any suggestions?
Thanks
Comment