#Extract all the fasta header name with certain ids
kraken --db ../../../../DATABASE/minikraken_20171019_8GB.tgz out.fa
more out.fa_class.txt | grep "227859" | awk '{print $2}' > all_real_ids.txt
minimap2 -t 36 -k19 -w5 -A1 -B2 -O3,13 -E2,1 -s200 -z200 -N50 --min-occ-floor=100 finaal_output.fasta finaal_output.fasta > finaal_self_align.paf