CHECK report for Rhisat2 on merida1
This page was generated on 2020-04-15 12:48:37 -0400 (Wed, 15 Apr 2020).
Rhisat2 1.2.0 Charlotte Soneson
Snapshot Date: 2020-04-14 16:46:13 -0400 (Tue, 14 Apr 2020) |
URL: https://git.bioconductor.org/packages/Rhisat2 |
Branch: RELEASE_3_10 |
Last Commit: 458e19e |
Last Changed Date: 2019-10-29 13:11:28 -0400 (Tue, 29 Oct 2019) |
| malbec1 | Linux (Ubuntu 18.04.4 LTS) / x86_64 | OK | OK | OK | | |
tokay1 | Windows Server 2012 R2 Standard / x64 | OK | OK | OK | OK | |
merida1 | OS X 10.11.6 El Capitan / x86_64 | OK | OK | [ OK ] | OK | |
Summary
Command output
Installation output
Rhisat2.Rcheck/00install.out
Tests output
Rhisat2.Rcheck/tests/testthat.Rout
R version 3.6.3 (2020-02-29) -- "Holding the Windsock"
Copyright (C) 2020 The R Foundation for Statistical Computing
Platform: x86_64-apple-darwin15.6.0 (64-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> library(testthat)
> library(Rhisat2)
>
> test_check("Rhisat2")
*** Warning ***
'hisat2-build-s' was run directly. It is recommended that you run the wrapper script 'hisat2-build' instead.
[1] "HISAT2 version 2.1.0 by Daehwan Kim (infphilo@gmail.com, http://www.ccb.jhu.edu/people/infphilo)"
[2] "Usage: hisat2-build [options]* <reference_in> <ht2_index_base>"
[3] " reference_in comma-separated list of files with ref sequences"
[4] " hisat2_index_base write ht2 data to files with this dir/basename"
[5] "Options:"
[6] " -c reference sequences given on cmd line (as"
[7] " <reference_in>)"
[8] " -a/--noauto disable automatic -p/--bmax/--dcv memory-fitting"
[9] " -p number of threads"
[10] " --bmax <int> max bucket sz for blockwise suffix-array builder"
[11] " --bmaxdivn <int> max bucket sz as divisor of ref len (default: 4)"
[12] " --dcv <int> diff-cover period for blockwise (default: 1024)"
[13] " --nodc disable diff-cover (algorithm becomes quadratic)"
[14] " -r/--noref don't build .3/.4.ht2 (packed reference) portion"
[15] " -3/--justref just build .3/.4.ht2 (packed reference) portion"
[16] " -o/--offrate <int> SA is sampled every 2^offRate BWT chars (default: 5)"
[17] " -t/--ftabchars <int> # of chars consumed in initial lookup (default: 10)"
[18] " --localoffrate <int> SA (local) is sampled every 2^offRate BWT chars (default: 3)"
[19] " --localftabchars <int> # of chars consumed in initial lookup in a local index (default: 6)"
[20] " --snp <path> SNP file name"
[21] " --haplotype <path> haplotype file name"
[22] " --ss <path> Splice site file name"
[23] " --exon <path> Exon file name"
[24] " --seed <int> seed for random number generator"
[25] " -q/--quiet verbose output (for debugging)"
[26] " -h/--help print detailed description of tool and its options"
[27] " --usage print this usage message"
[28] " --version print version information and quit"
*** Warning ***
'hisat2-align' was run directly. It is recommended that you run the wrapper script 'hisat2' instead.
[1] "HISAT2 version 2.1.0 by Daehwan Kim (infphilo@gmail.com, www.ccb.jhu.edu/people/infphilo)"
[2] "Usage: "
[3] " hisat2-align [options]* -x <ht2-idx> {-1 <m1> -2 <m2> | -U <r>} [-S <sam>]"
[4] ""
[5] " <ht2-idx> Index filename prefix (minus trailing .X.ht2)."
[6] " <m1> Files with #1 mates, paired with files in <m2>."
[7] " <m2> Files with #2 mates, paired with files in <m1>."
[8] " <r> Files with unpaired reads."
[9] " <sam> File for SAM output (default: stdout)"
[10] ""
[11] " <m1>, <m2>, <r> can be comma-separated lists (no whitespace) and can be"
[12] " specified many times. E.g. '-U file1.fq,file2.fq -U file3.fq'."
[13] ""
[14] "Options (defaults in parentheses):"
[15] ""
[16] " Input:"
[17] " -q query input files are FASTQ .fq/.fastq (default)"
[18] " --qseq query input files are in Illumina's qseq format"
[19] " -f query input files are (multi-)FASTA .fa/.mfa"
[20] " -r query input files are raw one-sequence-per-line"
[21] " -c <m1>, <m2>, <r> are sequences themselves, not files"
[22] " -s/--skip <int> skip the first <int> reads/pairs in the input (none)"
[23] " -u/--upto <int> stop after first <int> reads/pairs (no limit)"
[24] " -5/--trim5 <int> trim <int> bases from 5'/left end of reads (0)"
[25] " -3/--trim3 <int> trim <int> bases from 3'/right end of reads (0)"
[26] " --phred33 qualities are Phred+33 (default)"
[27] " --phred64 qualities are Phred+64"
[28] " --int-quals qualities encoded as space-delimited integers"
[29] ""
[30] " Alignment:"
[31] " --n-ceil <func> func for max # non-A/C/G/Ts permitted in aln (L,0,0.15)"
[32] " --ignore-quals treat all quality values as 30 on Phred scale (off)"
[33] " --nofw do not align forward (original) version of read (off)"
[34] " --norc do not align reverse-complement version of read (off)"
[35] ""
[36] " Spliced Alignment:"
[37] " --pen-cansplice <int> penalty for a canonical splice site (0)"
[38] " --pen-noncansplice <int> penalty for a non-canonical splice site (12)"
[39] " --pen-canintronlen <func> penalty for long introns (G,-8,1) with canonical splice sites"
[40] " --pen-noncanintronlen <func> penalty for long introns (G,-8,1) with noncanonical splice sites"
[41] " --min-intronlen <int> minimum intron length (20)"
[42] " --max-intronlen <int> maximum intron length (500000)"
[43] " --known-splicesite-infile <path> provide a list of known splice sites"
[44] " --novel-splicesite-outfile <path> report a list of splice sites"
[45] " --novel-splicesite-infile <path> provide a list of novel splice sites"
[46] " --no-temp-splicesite disable the use of splice sites found"
[47] " --no-spliced-alignment disable spliced alignment"
[48] " --rna-strandness <string> specify strand-specific information (unstranded)"
[49] " --tmo reports only those alignments within known transcriptome"
[50] " --dta reports alignments tailored for transcript assemblers"
[51] " --dta-cufflinks reports alignments tailored specifically for cufflinks"
[52] " --avoid-pseudogene tries to avoid aligning reads to pseudogenes (experimental option)\023"
[53] " --no-templatelen-adjustment disables template length adjustment for RNA-seq reads"
[54] ""
[55] " Scoring:"
[56] " --mp <int>,<int> max and min penalties for mismatch; lower qual = lower penalty <6,2>"
[57] " --sp <int>,<int> max and min penalties for soft-clipping; lower qual = lower penalty <2,1>"
[58] " --no-softclip no soft-clipping"
[59] " --np <int> penalty for non-A/C/G/Ts in read/ref (1)"
[60] " --rdg <int>,<int> read gap open, extend penalties (5,3)"
[61] " --rfg <int>,<int> reference gap open, extend penalties (5,3)"
[62] " --score-min <func> min acceptable alignment score w/r/t read length"
[63] " (L,0.0,-0.2)"
[64] ""
[65] " Reporting:"
[66] " -k <int> (default: 5) report up to <int> alns per read"
[67] ""
[68] " Paired-end:"
[69] " -I/--minins <int> minimum fragment length (0), only valid with --no-spliced-alignment"
[70] " -X/--maxins <int> maximum fragment length (500), only valid with --no-spliced-alignment"
[71] " --fr/--rf/--ff -1, -2 mates align fw/rev, rev/fw, fw/fw (--fr)"
[72] " --no-mixed suppress unpaired alignments for paired reads"
[73] " --no-discordant suppress discordant alignments for paired reads"
[74] ""
[75] " Output:"
[76] " -t/--time print wall-clock time taken by search phases"
[77] " --summary-file print alignment summary to this file."
[78] " --new-summary print alignment summary in a new style, which is more machine-friendly."
[79] " --quiet print nothing to stderr except serious errors"
[80] " --met-file <path> send metrics to file at <path> (off)"
[81] " --met-stderr send metrics to stderr (off)"
[82] " --met <int> report internal counters & metrics every <int> secs (1)"
[83] " --no-head supppress header lines, i.e. lines starting with @"
[84] " --no-sq supppress @SQ header lines"
[85] " --rg-id <text> set read group id, reflected in @RG line and RG:Z: opt field"
[86] " --rg <text> add <text> (\"lab:value\") to @RG line of SAM header."
[87] " Note: @RG line only printed when --rg-id is set."
[88] " --omit-sec-seq put '*' in SEQ and QUAL fields for secondary alignments."
[89] ""
[90] " Performance:"
[91] " -o/--offrate <int> override offrate of index; must be >= index's offrate"
[92] " -p/--threads <int> number of alignment threads to launch (1)"
[93] " --reorder force SAM output order to match order of input reads"
[94] " --mm use memory-mapped I/O for index; many 'hisat2's can share"
[95] ""
[96] " Other:"
[97] " --qc-filter filter out reads that are bad according to QSEQ filter"
[98] " --seed <int> seed for random number generator (0)"
[99] " --non-deterministic seed rand. gen. arbitrarily instead of using read attributes"
[100] " --remove-chrname remove 'chr' from reference names in alignment"
[101] " --add-chrname add 'chr' to reference names in alignment "
[102] " --version print version information and quit"
[103] " -h/--help print this usage message"
[1] "/Library/Frameworks/R.framework/Versions/3.6/Resources/library/Rhisat2/hisat2 version 2.1.0"
[2] "64-bit"
[3] "Built on merida1.bioconductor.org"
[4] "Tue Apr 14 19:26:49 EDT 2020"
[5] "Compiler: InstalledDir: /Library/Developer/CommandLineTools/usr/bin"
[6] "Options: -O3 -m64 -msse2 -funroll-loops -g3 -DPOPCNT_CAPABILITY"
[7] "Sizeof {int, long, long long, void*, size_t, off_t}: {4, 8, 8, 8, 8, 8}"
Settings:
Output files: "/tmp/RtmpOs19rw/hisat2indexdir/index.*.ht2"
Line rate: 6 (line is 64 bytes)
Lines per side: 1 (side is 64 bytes)
Offset rate: 4 (one in 16)
FTable chars: 10
Strings: unpacked
Local offset rate: 3 (one in 8)
Local fTable chars: 6
Local sequence length: 57344
Local sequence overlap between two consecutive indexes: 1024
Endianness: little
Actual local endianness: little
Sanity checking: disabled
Assertions: disabled
Random seed: 0
Sizeofs: void*:8, int:4, long:8, size_t:8
Input files DNA, FASTA:
/Library/Frameworks/R.framework/Versions/3.6/Resources/library/Rhisat2/extdata/refs/chr1.fa
Reading reference sizes
Time reading reference sizes: 00:00:00
Calculating joined length
Writing header
Reserving space for joined string
Joining reference sequences
Time to join reference sequences: 00:00:00
Time to read SNPs and splice sites: 00:00:00
Using parameters --bmax 18650 --dcv 1024
Doing ahead-of-time memory usage test
Passed! Constructing with these parameters: --bmax 18650 --dcv 1024
Constructing suffix-array element generator
Converting suffix-array elements to index image
Allocating ftab, absorbFtab
Entering GFM loop
Exited GFM loop
fchr[A]: 0
fchr[C]: 26955
fchr[G]: 49554
fchr[T]: 72459
fchr[$]: 99465
Exiting GFM::buildToDisk()
Returning from initFromVector
Wrote 4227718 bytes to primary GFM file: /tmp/RtmpOs19rw/hisat2indexdir/index.1.ht2
Wrote 24872 bytes to secondary GFM file: /tmp/RtmpOs19rw/hisat2indexdir/index.2.ht2
Re-opening _in1 and _in2 as input streams
Returning from GFM constructor
Returning from initFromVector
Wrote 45363 bytes to primary GFM file: /tmp/RtmpOs19rw/hisat2indexdir/index.5.ht2
Wrote 25130 bytes to secondary GFM file: /tmp/RtmpOs19rw/hisat2indexdir/index.6.ht2
Re-opening _in5 and _in5 as input streams
Returning from HierEbwt constructor
Headers:
len: 99465
gbwtLen: 99466
nodes: 99466
sz: 24867
gbwtSz: 24867
lineRate: 6
offRate: 4
offMask: 0xfffffff0
ftabChars: 10
eftabLen: 0
eftabSz: 0
ftabLen: 1048577
ftabSz: 4194308
offsLen: 6217
offsSz: 24868
lineSz: 64
sideSz: 64
sideGbwtSz: 48
sideGbwtLen: 192
numSides: 519
numLines: 519
gbwtTotLen: 33216
gbwtTotSz: 33216
reverse: 0
linearFM: Yes
Total time for call to driver() for forward index: 00:00:00
Settings:
Output files: "/tmp/RtmpOs19rw/hisat2indexdir/index.*.ht2"
Line rate: 6 (line is 64 bytes)
Lines per side: 1 (side is 64 bytes)
Offset rate: 4 (one in 16)
FTable chars: 10
Strings: unpacked
Local offset rate: 3 (one in 8)
Local fTable chars: 6
Local sequence length: 57344
Local sequence overlap between two consecutive indexes: 1024
Endianness: little
Actual local endianness: little
Sanity checking: disabled
Assertions: disabled
Random seed: 0
Sizeofs: void*:8, int:4, long:8, size_t:8
Input files DNA, FASTA:
/Library/Frameworks/R.framework/Versions/3.6/Resources/library/Rhisat2/extdata/refs/chr1.fa
Reading reference sizes
Time reading reference sizes: 00:00:00
Calculating joined length
Writing header
Reserving space for joined string
Joining reference sequences
Time to join reference sequences: 00:00:00
Time to read SNPs and splice sites: 00:00:00
Using parameters --bmax 18650 --dcv 1024
Doing ahead-of-time memory usage test
Passed! Constructing with these parameters: --bmax 18650 --dcv 1024
Constructing suffix-array element generator
Converting suffix-array elements to index image
Allocating ftab, absorbFtab
Entering GFM loop
Exited GFM loop
fchr[A]: 0
fchr[C]: 26955
fchr[G]: 49554
fchr[T]: 72459
fchr[$]: 99465
Exiting GFM::buildToDisk()
Returning from initFromVector
Wrote 4227718 bytes to primary GFM file: /tmp/RtmpOs19rw/hisat2indexdir/index.1.ht2
Wrote 24872 bytes to secondary GFM file: /tmp/RtmpOs19rw/hisat2indexdir/index.2.ht2
Re-opening _in1 and _in2 as input streams
Returning from GFM constructor
Returning from initFromVector
Wrote 45363 bytes to primary GFM file: /tmp/RtmpOs19rw/hisat2indexdir/index.5.ht2
Wrote 25130 bytes to secondary GFM file: /tmp/RtmpOs19rw/hisat2indexdir/index.6.ht2
Re-opening _in5 and _in5 as input streams
Returning from HierEbwt constructor
Headers:
len: 99465
gbwtLen: 99466
nodes: 99466
sz: 24867
gbwtSz: 24867
lineRate: 6
offRate: 4
offMask: 0xfffffff0
ftabChars: 10
eftabLen: 0
eftabSz: 0
ftabLen: 1048577
ftabSz: 4194308
offsLen: 6217
offsSz: 24868
lineSz: 64
sideSz: 64
sideGbwtSz: 48
sideGbwtLen: 192
numSides: 519
numLines: 519
gbwtTotLen: 33216
gbwtTotSz: 33216
reverse: 0
linearFM: Yes
Total time for call to driver() for forward index: 00:00:00
4 reads; of these:
4 (100.00%) were unpaired; of these:
2 (50.00%) aligned 0 times
2 (50.00%) aligned exactly 1 time
0 (0.00%) aligned >1 times
50.00% overall alignment rate
══ testthat results ═══════════════════════════════════════════════════════════
[ OK: 32 | SKIPPED: 0 | WARNINGS: 1 | FAILED: 0 ]
>
> proc.time()
user system elapsed
21.380 1.016 22.522
Example timings
Rhisat2.Rcheck/Rhisat2-Ex.timings