
################################################################################
Started ShapeMapper v2.1.5 at 2021-03-26 13:25:11
Running from directory: /pine/scr/p/s/psirving/20210326-rnasep-acid-resequenced
args:  --target rnasep.fa --name API --modified --folder Sample_8 --untreated --folder Sample_10 --amplicon --output-parsed-mutations --overwrite --per-read-histograms
Created pipeline at 2021-03-26 13:25:11
Running PrimerLocator at 2021-03-26 13:25:11 . . .
. . . done at 2021-03-26 13:25:11
Running FastaFormatChecker at 2021-03-26 13:25:11 . . .
. . . done at 2021-03-26 13:25:11
Running BowtieIndexBuilder at 2021-03-26 13:25:11 . . .
. . . done at 2021-03-26 13:25:11
Running process group 4 at 2021-03-26 13:25:11 . . .
  Including these components:
    ProgressMonitor . . . started at 2021-03-26 13:25:11
    QualityTrimmer1 . . . started at 2021-03-26 13:25:11
    QualityTrimmer2 . . . started at 2021-03-26 13:25:11
    Interleaver . . . started at 2021-03-26 13:25:11
    Merger . . . started at 2021-03-26 13:25:11
    Tab6Interleaver . . . started at 2021-03-26 13:25:11
    BowtieAligner . . . started at 2021-03-26 13:25:11
    MutationParser_Modified . . . started at 2021-03-26 13:25:11
    SplitToFile1 . . . started at 2021-03-26 13:25:11
    MutationCounter_Modified . . . started at 2021-03-26 13:25:11
  /`````````````````````````````````````````````````````````````````````````````
  |Merger (sample: Modified) output message: 
  |----------------------------------------- 
  | 
  | java -Djava.library.path=/nas/longleaf/home/psirving/shapemapper-2.1.5/internals/thirdparty/miniconda/opt/bbmap-37.78/jni/ -ea -Xmx1000m -Xmx1000m -cp /nas/longleaf/home/psirving/shapemapper-2.1.5/internals/thirdparty/miniconda/opt/bbmap-37.78/current/ jgi.BBMerge vstrict=t in=stdin out=stdout outu=stdout interleaved=t usejni=t t=4
  | Executing jgi.BBMerge [vstrict=t, in=stdin, out=stdout, outu=stdout, interleaved=t, usejni=t, t=4]
  | Version 37.78 [vstrict=t, in=stdin, out=stdout, outu=stdout, interleaved=t, usejni=t, t=4]
  | 
  | Revised arguments: [maxbad=4, margin=3, minqo=8, qualiters=2, ratiomode=t, flatmode=f, minentropy=52, minoverlap=12, minoverlap0=4, maxratio=0.05, ratiomargin=12, ratiooffset=0.5, ratiominoverlapreduction=4, efilter=2, pfilter=0.008, minsecondratio=0.16, minapproxoverlap=22, in=stdin, out=stdout, outu=stdout, interleaved=t, usejni=t, t=4]
  | 
  | Set INTERLEAVED to true
  | Set threads to 4
  | Writing mergable reads merged.
  | Unspecified format for output stdout; defaulting to fastq.
  | Unspecified format for output stdout; defaulting to fastq.
  | Started output threads.
  | Unspecified format for input stdin; defaulting to fastq.
  | Total time: 1349.255 seconds.
  | 
  | Pairs:               	982276
  | Joined:              	839397   	85.454%
  | Ambiguous:           	139712   	14.223%
  | No Solution:         	3167       	0.322%
  | Too Short:           	0       	0.000%
  | 
  | Avg Insert:          	329.8
  | Standard Deviation:  	30.4
  | Mode:                	335
  | 
  | Insert range:        	38 - 542
  | 90th percentile:     	336
  | 75th percentile:     	335
  | 50th percentile:     	335
  | 25th percentile:     	333
  | 10th percentile:     	330
  \_____________________________________________________________________________
  /`````````````````````````````````````````````````````````````````````````````
  |BowtieAligner (sample: Modified) output message: 
  |------------------------------------------------ 
  | 
  | 982276 reads; of these:
  |   142879 (14.55%) were paired; of these:
  |     54691 (38.28%) aligned concordantly 0 times
  |     88136 (61.69%) aligned concordantly exactly 1 time
  |     52 (0.04%) aligned concordantly >1 times
  |     ----
  |     54691 pairs aligned concordantly 0 times; of these:
  |       76 (0.14%) aligned discordantly 1 time
  |     ----
  |     54615 pairs aligned 0 times concordantly or discordantly; of these:
  |       109230 mates make up the pairs; of these:
  |         71863 (65.79%) aligned 0 times
  |         37356 (34.20%) aligned exactly 1 time
  |         11 (0.01%) aligned >1 times
  |   839397 (85.45%) were unpaired; of these:
  |     29057 (3.46%) aligned 0 times
  |     809984 (96.50%) aligned exactly 1 time
  |     356 (0.04%) aligned >1 times
  | 91.03% overall alignment rate
  \_____________________________________________________________________________
  /`````````````````````````````````````````````````````````````````````````````
  |MutationCounter_Modified (RNA: rnasep) (sample: Modified) output message: 
  |------------------------------------------------------------------------- 
  | 
  | Read lengths
  | --------------------
  | bin_range	frequency
  | [0,49]	0.077227
  | [50,99]	0.050866
  | [100,149]	0.039843
  | [150,199]	0.034627
  | [200,249]	0.030445
  | [250,299]	0.025707
  | [300,349]	0.741284
  | [350,399]	0.000000
  | [400,449]	0.000000
  | [450,499]	0.000000
  | [500,549]	0.000000
  | [550,599]	0.000000
  | [600,649]	0.000000
  | [650,699]	0.000000
  | [700,749]	0.000000
  | [750,799]	0.000000
  | [800,849]	0.000000
  | [850,899]	0.000000
  | [900,949]	0.000000
  | [950,999]	0.000000
  | >=1000	0.000000
  | --------------------
  | 
  | Mutations per read
  | --------------------
  | bin_left	frequency
  | 0	0.190040
  | 1	0.204742
  | 2	0.200250
  | 3	0.160165
  | 4	0.105834
  | 5	0.061685
  | 6	0.033721
  | 7	0.018693
  | 8	0.010520
  | 9	0.006109
  | 10	0.003573
  | 11	0.002001
  | 12	0.001111
  | 13	0.000629
  | 14	0.000373
  | 15	0.000216
  | 16	0.000129
  | 17	0.000083
  | 18	0.000050
  | 19	0.000026
  | 20	0.000049
  | --------------------
  \_____________________________________________________________________________
. . . done at 2021-03-26 13:47:42
Running process group 5 at 2021-03-26 13:47:42 . . .
  Including these components:
    ProgressMonitor . . . started at 2021-03-26 13:47:42
    QualityTrimmer1 . . . started at 2021-03-26 13:47:42
    QualityTrimmer2 . . . started at 2021-03-26 13:47:42
    Interleaver . . . started at 2021-03-26 13:47:42
    Merger . . . started at 2021-03-26 13:47:42
    Tab6Interleaver . . . started at 2021-03-26 13:47:42
    BowtieAligner . . . started at 2021-03-26 13:47:42
    MutationParser_Untreated . . . started at 2021-03-26 13:47:42
    SplitToFile1 . . . started at 2021-03-26 13:47:42
    MutationCounter_Untreated . . . started at 2021-03-26 13:47:42
  /`````````````````````````````````````````````````````````````````````````````
  |Merger (sample: Untreated) output message: 
  |------------------------------------------ 
  | 
  | java -Djava.library.path=/nas/longleaf/home/psirving/shapemapper-2.1.5/internals/thirdparty/miniconda/opt/bbmap-37.78/jni/ -ea -Xmx1000m -Xmx1000m -cp /nas/longleaf/home/psirving/shapemapper-2.1.5/internals/thirdparty/miniconda/opt/bbmap-37.78/current/ jgi.BBMerge vstrict=t in=stdin out=stdout outu=stdout interleaved=t usejni=t t=4
  | Executing jgi.BBMerge [vstrict=t, in=stdin, out=stdout, outu=stdout, interleaved=t, usejni=t, t=4]
  | Version 37.78 [vstrict=t, in=stdin, out=stdout, outu=stdout, interleaved=t, usejni=t, t=4]
  | 
  | Revised arguments: [maxbad=4, margin=3, minqo=8, qualiters=2, ratiomode=t, flatmode=f, minentropy=52, minoverlap=12, minoverlap0=4, maxratio=0.05, ratiomargin=12, ratiooffset=0.5, ratiominoverlapreduction=4, efilter=2, pfilter=0.008, minsecondratio=0.16, minapproxoverlap=22, in=stdin, out=stdout, outu=stdout, interleaved=t, usejni=t, t=4]
  | 
  | Set INTERLEAVED to true
  | Set threads to 4
  | Writing mergable reads merged.
  | Unspecified format for output stdout; defaulting to fastq.
  | Unspecified format for output stdout; defaulting to fastq.
  | Started output threads.
  | Unspecified format for input stdin; defaulting to fastq.
  | Total time: 1415.967 seconds.
  | 
  | Pairs:               	960965
  | Joined:              	828734   	86.240%
  | Ambiguous:           	129236   	13.449%
  | No Solution:         	2995       	0.312%
  | Too Short:           	0       	0.000%
  | 
  | Avg Insert:          	324.8
  | Standard Deviation:  	51.2
  | Mode:                	335
  | 
  | Insert range:        	35 - 542
  | 90th percentile:     	336
  | 75th percentile:     	335
  | 50th percentile:     	335
  | 25th percentile:     	335
  | 10th percentile:     	330
  \_____________________________________________________________________________
  /`````````````````````````````````````````````````````````````````````````````
  |BowtieAligner (sample: Untreated) output message: 
  |------------------------------------------------- 
  | 
  | 960965 reads; of these:
  |   132231 (13.76%) were paired; of these:
  |     57231 (43.28%) aligned concordantly 0 times
  |     74972 (56.70%) aligned concordantly exactly 1 time
  |     28 (0.02%) aligned concordantly >1 times
  |     ----
  |     57231 pairs aligned concordantly 0 times; of these:
  |       206 (0.36%) aligned discordantly 1 time
  |     ----
  |     57025 pairs aligned 0 times concordantly or discordantly; of these:
  |       114050 mates make up the pairs; of these:
  |         76067 (66.70%) aligned 0 times
  |         37979 (33.30%) aligned exactly 1 time
  |         4 (0.00%) aligned >1 times
  |   828734 (86.24%) were unpaired; of these:
  |     39366 (4.75%) aligned 0 times
  |     789219 (95.23%) aligned exactly 1 time
  |     149 (0.02%) aligned >1 times
  | 89.44% overall alignment rate
  \_____________________________________________________________________________
  /`````````````````````````````````````````````````````````````````````````````
  |MutationCounter_Untreated (RNA: rnasep) (sample: Untreated) output message: 
  |--------------------------------------------------------------------------- 
  | 
  | Read lengths
  | --------------------
  | bin_range	frequency
  | [0,49]	0.088702
  | [50,99]	0.045644
  | [100,149]	0.033110
  | [150,199]	0.031036
  | [200,249]	0.028484
  | [250,299]	0.027207
  | [300,349]	0.745816
  | [350,399]	0.000000
  | [400,449]	0.000000
  | [450,499]	0.000000
  | [500,549]	0.000000
  | [550,599]	0.000000
  | [600,649]	0.000000
  | [650,699]	0.000000
  | [700,749]	0.000000
  | [750,799]	0.000000
  | [800,849]	0.000000
  | [850,899]	0.000000
  | [900,949]	0.000000
  | [950,999]	0.000000
  | >=1000	0.000000
  | --------------------
  | 
  | Mutations per read
  | --------------------
  | bin_left	frequency
  | 0	0.449604
  | 1	0.277715
  | 2	0.139302
  | 3	0.077432
  | 4	0.030522
  | 5	0.013959
  | 6	0.006524
  | 7	0.002851
  | 8	0.001226
  | 9	0.000544
  | 10	0.000189
  | 11	0.000077
  | 12	0.000026
  | 13	0.000011
  | 14	0.000006
  | 15	0.000002
  | 16	0.000002
  | 17	0.000003
  | 18	0.000002
  | 19	0.000002
  | 20	0.000001
  | --------------------
  \_____________________________________________________________________________
. . . done at 2021-03-26 14:11:19
Running CalcProfile at 2021-03-26 14:11:19 . . .
. . . done at 2021-03-26 14:11:19
Running NormProfile at 2021-03-26 14:11:19 . . .
. . . done at 2021-03-26 14:11:20
Running TabToShape at 2021-03-26 14:11:20 . . .
. . . done at 2021-03-26 14:11:20
Running RenderFigures at 2021-03-26 14:11:20 . . .
/```````````````````````````````````````````````````````````````````````````````
|RenderFigures (RNA: rnasep) output message: 
|------------------------------------------- 
| 
| Quality control checks:
| 
| (See README.md for explanation)
| 
| Read depth check:
| 100.0% (294/294) nucleotides meet the minimum read depth of 5000
| PASS
| 
| Mutation rate check:
| 97.3% (286/294) nucleotides have positive mutation rates
| above background
| High background check:
| 0.0% (0/294) nucleotides have high background
| mutation rates.
| PASS: Not too many nucleotides with high background mutation rates.
| 
| Number highly reactive check:
| 22.8% (67/294) nucleotides show high apparent reactivity.
| PASS: The expected number of nucleotides or more are highly
|       reactive.
\_______________________________________________________________________________
. . . done at 2021-03-26 14:11:24
Running RenderMappedDepths at 2021-03-26 14:11:24 . . .
/```````````````````````````````````````````````````````````````````````````````
|RenderMappedDepths (RNA: rnasep) output message: 
|------------------------------------------------ 
| 
| Approximate maximum read depth for reads mapping near each amplicon primer pair
| (See *_mapped_depths.pdf for plots)
| (See *_per-amplicon_abundance.txt for exact numbers)
| 
| Modified:
|   pair forward                       left    right   reverse                       left    right   depth  
|   ---- ----------------------------- ------- ------- ----------------------------- ------- ------- ------ 
|   1    GGCCTTCGGGCCAA                0       13      GAACCGGACCGAAGCC              309     324     874k   
| 
| Untreated:
|   pair forward                       left    right   reverse                       left    right   depth  
|   ---- ----------------------------- ------- ------- ----------------------------- ------- ------- ------ 
|   1    GGCCTTCGGGCCAA                0       13      GAACCGGACCGAAGCC              309     324     831k   
| 
\_______________________________________________________________________________
. . . done at 2021-03-26 14:11:26
ShapeMapper run successfully completed at 2021-03-26 14:11:26.
