#!/bin/bash

locus=$(zcat < GCF_017654675.1_Xenopus_laevis_v10.1_genomic.gbff.gz | grep -n LOCUS | tail -n 1 | sed 's/:.*//')


mtDNA=$(zcat < GCF_017654675.1_Xenopus_laevis_v10.1_genomic.gbff.gz | sed -n "${locus}",'$p')


echo "${mtDNA}" | sed -n '/FEATURES/,/ORIGIN/ p' | head -n $(($( echo "${mtDNA}" | sed -n '/FEATURES/,/ORIGIN/ p' | wc -l )-1)) | tail -n +2 | 
sed -n '/[[:digit:]]\.\.[[:digit:]]/p' | sed '/^[[:space:]]*\//d'| sed 's/[[:space:]]*//' | sed 's/[[:space:]]/\t/g' | cut -f 1 | sort | uniq -c