Testing out Pilon

Objective: Test out using Pilon on the MiSeq and PGM data Approach: Running Pilon using bash commands, evaluate revise reference relative to HGAP and GenBank LT2

Testing single paried end dataset



In [10]:

    
%%bash
java -Xmx4G -jar ../utilities/pilon-1.10.jar \
                --genome ../data/RM8375/ref/CFSAN008157.HGAP.fasta \
                --frags ../analysis/bioinf/sequence_purity/mapping/SRR1555296.bam \
                --changes --vcf --tracks \
                --fix "all" --debug #note --fix "all" default









    



Pilon version 1.10 Tue Nov 18 20:14:04 2014 -0500
Genome: ../data/RM8375/ref/CFSAN008157.HGAP.fasta
OpenJDK 64-Bit Server VM warning: INFO: os::commit_memory(0x00000007b2500000, 20447232, 0) failed; error='Cannot allocate memory' (errno=12)
#
# There is insufficient memory for the Java Runtime Environment to continue.
# Native memory allocation (malloc) failed to allocate 20447232 bytes for committing reserved memory.
# An error report file with more information is saved as:
# /notebooks/dev/hs_err_pid246.log

Error with head, needed to rerun bwa with revised params file

Error with memory allocation, not sure if too many applicaitons running on laptop, or docker memory restrictions

re-running after closing windows and restarting finder (hung up on moving files to trash)
Need to look into changing amount of memory allocated to ipython notebook



In [ ]:



In [12]:

    
%%bash 
python ../analysis/bioinf/sequence_purity/run_bwa_mem_pe.py \
        ../analysis/bioinf/sequence_purity/bwa_mem_pipeline_params.txt

Testing pilon on PGM data



In [5]:

    
%%bash
java -Xmx8G -jar ../utilities/pilon-1.10.jar \
                --genome ../data/RM8375/ref/CFSAN008157.HGAP.fasta \
                --unpaired ../data/RM8375/PGM/bam/IonXpress_001_R_2014_03_23_18_22_09_user_SN2-17-8375_Orthogonal_Measurement_1_Run_2_PacBioRef2.bam \
                --changes --vcf --tracks \
                --fix "all" --debug #note --fix "all" default









    



Process is terminated.