When parsing depths from a BAM file, the reference sequence ID is retrieved from the BAM header info, however, that reference sequence may not always be an NCBI sequence. This is an issue with IRMA output where BAM files have reference sequence IDs that are fairly generic (e.g. "A_PB2"). In this case, it may be preferable to not try to retrieve any ref seq and omit it entirely from the output or spoof it based on a consensus sequence.
When parsing depths from a BAM file, the reference sequence ID is retrieved from the BAM header info, however, that reference sequence may not always be an NCBI sequence. This is an issue with IRMA output where BAM files have reference sequence IDs that are fairly generic (e.g. "A_PB2"). In this case, it may be preferable to not try to retrieve any ref seq and omit it entirely from the output or spoof it based on a consensus sequence.