To improve on interpretation of alignments, the “by alignment score” (byAS) method was developed, based on the hypothesis that among multiple alignments from each sequencing record to the mixed genome index, the one with the greatest alignment score is the best match. These scores are used to split the original sequencing files into two sets of species-specific RNAseq reads for re-alignment with only the appropriate genome.