diff --git a/htsinfer/mapping.py b/htsinfer/mapping.py index e1cef46..52780f3 100644 --- a/htsinfer/mapping.py +++ b/htsinfer/mapping.py @@ -51,11 +51,7 @@ def __init__( self.star_dirs: List[Path] = [] def evaluate(self): - """Infer read orientation. - - Returns: - Orientation results object. - """ + """Align FASTQ files to reference transcripts with STAR.""" # get transcripts for current organims transcripts = self.subset_transcripts_by_organism() @@ -270,6 +266,10 @@ def prepare_star_alignment_commands( ) -> Dict[Path, List[str]]: """Prepare STAR alignment commands. + Input FASTQ files are assumed to be sorted according to reference names + or coordinates, the order of input reads is kept with the option + "PairedKeepInputOrder", no additional sorting of aligned reads is done. + Args: index_dir: Path to directory containing STAR index.