I just created an accidental test case. I had a nt file mixed in with my protein files. When diaomond tried to map it-- it didn't recognize it as a nt file and tried to align it without success. I think we should build in an assertion that a file passed is in fact what it claims to be (protein or nucleotide). Could potentially do this with BioPython... or maybe write our own.
I just created an accidental test case. I had a nt file mixed in with my protein files. When diaomond tried to map it-- it didn't recognize it as a nt file and tried to align it without success. I think we should build in an assertion that a file passed is in fact what it claims to be (protein or nucleotide). Could potentially do this with BioPython... or maybe write our own.