Hi, Sej,

I found that I have to use sam files with only mapped reads as input in the Paralyzer to make it work. I need to do

samtools view -b -F 4 input.bam > output_mapped.bam

to generate bam files with only mapped reads and then convert the bam to sam file by

samtools view -h output_mapped.bam > output_mapped.sam

to generate sam file that is working with Paralyzer.

Thanks for your input.

Xiao

strong textHi, ATpoint,

Thanks a lot for your reply. I am sorry for my late reply to your question (I am new here and I did not get email notification somehow for messages). I share my command line which leads to error as below:

(base) bieniaszs-ipro:hg19_NL4_3 bieniaszlab$ PARAlyzer 128G hnRNPU_NL43.ini
Running PARalyzer v1.5
Parsing SAM file(s)...Done
Creating Read Groups & Clusters...Exception in thread "main" java.lang.Exception: * is not found in the .2bit file
at MyTwoBitParser.loadChromosome(MyTwoBitParser.java:58)
at MyTwoBitParser.getSequence(MyTwoBitParser.java:118)
at PARalyze.main(PARalyze.java:204)

My 2bit file was generated by Fatotwobit tool of human hg19, I did not use filter file and my .ini file is as below:

BANDWIDTH=3
CONVERSION=T>C
MINIMUM_READ_COUNT_PER_GROUP=5
MINIMUM_READ_COUNT_PER_CLUSTER=5
MINIMUM_READ_COUNT_FOR_KDE=5
MINIMUM_CLUSTER_SIZE=11
MINIMUM_CONVERSION_LOCATIONS_FOR_CLUSTER=1
MINIMUM_CONVERSION_COUNT_FOR_CLUSTER=1
MINIMUM_READ_COUNT_FOR_CLUSTER_INCLUSION=5
MINIMUM_READ_LENGTH=1
MAXIMUM_NUMBER_OF_NON_CONVERSION_MISMATCHES=0

EXTEND_BY_READ

#ADDITIONAL_NUCLEOTIDES_BEYOND_SIGNAL=20

SAM_FILE=/Users/bieniaszlab/miniconda3/data/hnRNP_U/hg19_NL4_3/hnRNPU_NL43.sort.sam=COLLAPSED
GENOME_2BIT_FILE=/Users/bieniaszlab/miniconda3/data/hnRNP_U/hg19_NL4_3/NL43.2bit

OUTPUT_DISTRIBUTIONS_FILE=/Users/bieniaszlab/miniconda3/data/hnRNP_U/hg19_NL4_3/hnRNPU_NL43_distribution.csv
OUTPUT_GROUPS_FILE=/Users/bieniaszlab/miniconda3/data/hnRNP_U/hg19_NL4_3/hnRNPU_NL43_groups.csv
OUTPUT_CLUSTERS_FILE=/Users/bieniaszlab/miniconda3/data/hnRNP_U/hg19_NL4_3/hnRNPU_NL43_clusters.csv

Regards,
Xiao

Hi, ATpoint,

I really appreciate your stimulating input! I found that I have to use sam files with only mapped reads as input in the Paralyzer to make it work. I need to do samtools view -b -F 4 input.bam > output_mapped.bam to generate bam files with only mapped reads and then convert the bam to sam file by samtools view -h output_mapped.bam > output_mapped.sam to generate sam file that is working with Paralyzer. As you mentioned that in SAM/BAM files, typically * is used as chromosome name to indicate unmapped reads. This is likely giving me a problem in my Paralyzer run. My Paralyzer runs well with sam files with only mapped reads!

Hi, Seokju,

Sorry for the late reply but I just saw your message today. I did not include =COLLAPSED in my PARalyzer ini file. I use fastx_toolkit to collapse reads. Here is one of my ini file for PARalyzer below:

BANDWIDTH=3
CONVERSION=T>C
MINIMUM_READ_COUNT_PER_GROUP=10
MINIMUM_READ_COUNT_PER_CLUSTER=5
MINIMUM_READ_COUNT_FOR_KDE=5
MINIMUM_CLUSTER_SIZE=15
MINIMUM_CONVERSION_LOCATIONS_FOR_CLUSTER=2
MINIMUM_CONVERSION_COUNT_FOR_CLUSTER=2
MINIMUM_READ_COUNT_FOR_CLUSTER_INCLUSION=5
MINIMUM_READ_LENGTH=1
MAXIMUM_NUMBER_OF_NON_CONVERSION_MISMATCHES=1

EXTEND_BY_READ

GENOME_2BIT_FILE=/Users/bieniaszlab/Documents/Xiao/CLIP/Index/hg19/hg19.2bit
SAM_FILE=/Users/bieniaszlab/Documents/Xiao/CLIP/CLIP17_IRCLIP4/Nextseq_Midoutput/CLIP17_synNCp15_NC43FSFS_hg19_m2.sorted.sam
OUTPUT_DISTRIBUTIONS_FILE=/Users/bieniaszlab/Documents/Xiao/CLIP/CLIP17_IRCLIP4/CLIP17_synNCp15_NC43FSFS_hg19_m2_distri.csv
OUTPUT_GROUPS_FILE=/Users/bieniaszlab/Documents/Xiao/CLIP/CLIP17_IRCLIP4/Nextseq_Midoutput/CLIP17_synNCp15_NC43FSFS_hg19_m2_group.csv
OUTPUT_CLUSTERS_FILE=/Users/bieniaszlab/Documents/Xiao/CLIP/CLIP17_IRCLIP4/Nextseq_Midoutput/CLIP17_synNCp15_NC43FSFS_hg19_m2_clusters.csv