tDRnamer 1.0.0 (Dec 2021)
User: pchan
Run time: Fri Jan 28 01:25:04 2022
Command: tDRnamer --mode seq --seq /soe/pchan/home/data/tRNA/trfdb/fasta/human_trf.fa --db /www/lib/tDRnamer/hg38/hg38 --source euk --output trfdb_hg38 --minread 10 --minlen 13 --maxlen 60 --cores 4 --force

Status: Checking input options
Status: Checking sequences in /soe/pchan/home/data/tRNA/trfdb/fasta/human_trf.fa
Total number of sequences: 157
Missing sequence name: 0
Sequence name with special characters: 0
Below minimum sequence length (13 bp): 0
Above maximum sequence length (60 bp): 0
Invalid sequence: 0
Total number of sequences after filtering: 157

Status: Reformat sequences in trfdb_hg38-filtered-seq.fa
Status: Search sequences against reference database
Status: Align sequences with bowtie2
157 reads; of these:
  157 (100.00%) were unpaired; of these:
    2 (1.27%) aligned 0 times
    32 (20.38%) aligned exactly 1 time
    123 (78.34%) aligned >1 times
98.73% overall alignment rate

Status: Filtering and identifying tDRs
Status: 155 sequences are identified as tDRs.
Status: Grouping tDRs
Status: Grouping tDRs
Status: Naming tDRs and summarizing results
Status: Process completed.