get odd specutils.analysis.snr result for sdss galaxy spectrum

I just want to calculate snr of an SDSS galaxy spectrum. Here's my main code.

`
from specutils import Spectrum1D
from specutils.analysis import snr

sdss_spec = Spectrum1D.read('./spec-0266-51630-0036.fits')
snr(sdss_spec)
`
It turned out that the result is "Quantity 49.439533 1e-51 erg3 / (Angstrom3 cm6 s3)>", I'm not sure that's a normal result. Or did I miss something?

Thanks a lot.

Reproduce the report with Spectrum1D.read('./spec-0266-51630-0036.fits') and specutils.analysis.snr on the supplied SDSS spectrum. Trace the snr entry point and inspect how the returned quantity and units are calculated; done means determining whether this result is expected or identifying the behavior that needs correction.

Written by the indexing model from the issue text.