HMMRATAC
HMMRATAC copied to clipboard
genome.info is created with 0 bytes
samtools view -H ATACseq.sorted.bam| perl -ne 'if(/^@SQ.*?SN:(\w+)\s+LN:(\d+)/){print $1,"\t",$2,"\n"}' > genome.info results in an empty genome.info file
I am using UBUNTU 20.04.5
Does anyone have any suggestion?