BOL

perl -e '$/="\n>"; while (<>) { s/>//g; my ($id, $seq) = split /\n/; print ">$_" if length $seq; }' < all_p_ctg.fasta > all_p_ctg_corrected.fasta #Convert #Fasta #Zero

#Zero coverage regions ~/Tools/bedtools2/bin/bedtools genomecov -ibam PvsH_aln_sorted.bam -bga | awk '$4==0' > ZERO_coverage_regions.txt