Entering edit mode
2.2 years ago
waik
•
0
I commanded
$BWA2 mem -t 1 -R "@RG\tID:sample\tSM:sample\tPL:platform" ../resource/reference/hg38.chr21.fa sample_1.fastq.gz sample_2.fastq.gz > sample.mapped.sam
but there was error
bash: mem: command not found
I tried to find solutions. but no clues.
Make sure the
$BWA2
bash variable is set withecho $BWA2
.and use simple quotes for option -R , otherwise, tabs will be expanded in the SAM header which can break some programs.