Commit 58576fb7 authored by peguerin's avatar peguerin
Browse files

results supernova run second try

parent aa392300
--------------------------------------------------------------------------------
SUMMARY
--------------------------------------------------------------------------------
- Sat Jul 13 16:08:21 2019
- [serran]
- software release = 2.1.1(6bb16452a)
- likely sequencers = HiSeq2000/2500
- assembly checksum = -3,183,485,510,844,980,441
--------------------------------------------------------------------------------
INPUT
- 250.02 M = READS = number of reads; ideal 800M-1200M for human
- 138.50 b = MEAN READ LEN = mean read length after trimming; ideal 140
- 50.40 x = RAW COV = raw coverage; ideal ~56
- 39.93 x = EFFECTIVE COV = effective read coverage; ideal ~42 for raw 56x
- 84.35 % = READ TWO Q30 = fraction of Q30 bases in read 2; ideal 75-85
- 352.00 b = MEDIAN INSERT = median insert size; ideal 350-400
- 87.74 % = PROPER PAIRS = fraction of proper read pairs; ideal >= 75
- 1.00 = BARCODE FRACTION = fraction of barcodes used; between 0 and 1
- 744.16 Mb = EST GENOME SIZE = estimated genome size
- 7.19 % = REPETITIVE FRAC = genome repetitivity index
- 0.08 % = HIGH AT FRACTION = high AT index
- 40.37 % = ASSEMBLY GC CONTENT = GC content of assembly
- 1.42 % = DINUCLEOTIDE FRACTION = dinucleotide content
- 14.63 Kb = MOLECULE LEN = weighted mean molecule size; ideal 50-100
- 23.94 = P10 = molecule count extending 10 kb on both sides
- 102.00 b = HETDIST = mean distance between heterozygous SNPs
- 5.55 % = UNBAR = fraction of reads that are not barcoded
- 194.00 = BARCODE N50 = N50 reads per barcode
- 7.57 % = DUPS = fraction of reads that are duplicates
- 54.85 % = PHASED = nonduplicate and phased reads; ideal 45-50
--------------------------------------------------------------------------------
OUTPUT
- 7.47 K = LONG SCAFFOLDS = number of scaffolds >= 10 kb
- 19.18 Kb = EDGE N50 = N50 edge size
- 65.41 Kb = CONTIG N50 = N50 contig size
- 259.87 Kb = PHASEBLOCK N50 = N50 phase block size
- 140.88 Kb = SCAFFOLD N50 = N50 scaffold size
- 15.39 % = MISSING 10KB = % of base assembly missing from scaffolds >= 10 kb
- 564.95 Mb = ASSEMBLY SIZE = assembly size (only scaffolds >= 10 kb)
--------------------------------------------------------------------------------
ALARMS
- The length-weighted mean molecule length is 14632.04 bases. The molecule
length estimation was successful, however, ideally we would expect a larger
value. Standard methods starting from blood can yield 100 kb or larger DNA, but
it can be difficult to obtain long DNA from other sample types. Short molecules
may reduce the scaffold and phase block N50 length, and could result in
misassemblies. We have observed assembly quality to improve with longer DNA.
--------------------------------------------------------------------------------
assembly_size,barcode_fraction,bases_per_read,bridge_1_50,bridge_50,bridge_model_badness_of_fit,checksum,contig_N50,dinucleotide_percent,dup_perc,edge_N50,effective_coverage,effective_coverage_median,est_genome_size,gc_percent,hetdist,high_AT_index,likely_sequencers,lw_mean_mol_len,m10,median_ins_sz,nreads,p10,phase_block_N50,placed_perc,proper_pairs_perc,q30_r2_perc,raw_coverage,read_rate_DF_1_threaded,read_rate_IO_10_threaded,read_rate_IO_1_threaded,repfrac,rpb_N50,scaffold_N50,scaffolds_10kb_plus,scaffolds_1kb_plus,unbar_perc,valid_bc_perc
564954197,1,138.5,518.12,0.001,1.05403,-3183485510844980441,65406,1.41542,7.56675,19184,39.9341,40.8665,7.44164e+08,40.3723,102,0.083919,"HiSeq2000/2500",14632,15.3899,352,250017372,23.9412,259871,54.8465,87.7445,84.3476,50.3956,106.177,884.351,939.628,7.18967,194,140884,7470,37832,5.54651,94.4535
This diff is collapsed.
......@@ -4,3 +4,4 @@ nohup /media/bigvol/peguerin/programmes/supernova-2.1.1/supernova mkoutput --sty
## second try
nohup ./supernova-2.1.1/supernova run --id=serran --fastqs=x_dir/ --localmem=470 --maxreads=250000000 &
nohup ./supernova-2.1.1/supernova mkoutput --style=pseudohap2 --asmdir=/media/bigvol/peguerin/serran/outs/assembly --outprefix=serranus &
\ No newline at end of file
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment