Index with samtools f
Tested on small dataset.