Check that BAM is indexed
flashton2003 opened this issue · comments
flashton2003 commented
Hello,
Just a small suggestion, but currently, if the BAM file is not indexed, clairvoyante produces 100 lines of error traceback, with this comment at the top.
[main_samview] random alignment retrieval only works for indexed BAM or CRAM files.
No read has been process, either the genome region you specified has no read cover, or please check the correctness of your BAM input (/home/ubuntu/tm_data/replimune/variant_calling/2020.04.02/RH018A_NGS_Sequence_Data/RH018A_NGS_Sequence_Data.primary_aln.bam).
Could I suggest that you add an explicit check for the index and return a more easily comprehensible error if it isn't there.
Best,
Phil
Ruibang Luo commented
Improved in Clair.