mageri
mageri copied to clipboard
Error while running with UMI
Hey Mike,
I am getting the following exception when running MAGERI:
[Fri Aug 11 12:22:47 EDT 2017 +00m00s] [my_project] Started analysis. [Fri Aug 11 12:22:47 EDT 2017 +00m00s] [my_project] Pre-processing sample group my_sample. [Fri Aug 11 12:22:47 EDT 2017 +00m00s] [Indexer] Building UMI index, 0 reads processed, 0.0% extracted.. Exception in thread "main" java.lang.RuntimeException: Error while parsing quality at com.milaboratory.core.sequencing.io.fastq.SFastqReader.parse(SFastqReader.java:258) at com.milaboratory.core.sequencing.io.fastq.PFastqReader.take(PFastqReader.java:213) at com.antigenomics.mageri.core.input.PMigReader$PairedReaderWrapper.take(PMigReader.java:182) at com.antigenomics.mageri.core.input.PMigReader$PairedReaderWrapper.take(PMigReader.java:166) at cc.redberry.pipe.blocks.O2ITransmitter.run(O2ITransmitter.java:66) at java.lang.Thread.run(Thread.java:744) Caused by: com.milaboratory.core.sequence.quality.WrongQualityStringException: [-1] at com.milaboratory.core.sequence.quality.SequenceQualityPhred.parse(SequenceQualityPhred.java:130) at com.milaboratory.core.sequencing.io.fastq.SFastqReader.parse(SFastqReader.java:256) ... 5 more
I am running this on paired end reads.
My Read 1 looks like this:
@NB501788:53:H3HYMBGX3:1:11101:14803:1046 1:N:0:TACAGGTC+NATGCTGG UMI:GATCAC:CCCFFD
CTCCTNAGATACTGTTATCGTGCAGCGCNNNNNNNNNNNNNNNNNNNNNNNNNNTTAAAGAAATATGCA
+
AAAAA#EEEEEEEEEEEEEEEEEAEEEE##########################AEEEEEEEEEEEEEE
@NB501788:53:H3HYMBGX3:1:11101:11212:1046 1:N:0:TACAGGTC+NATGCTGG UMI:TGGAAT:CCCFFD
GCCGTNATGCAGTAGCAGCGAGGCATTCNNNNNNNNNNNNNNNNNNNNNNNNNNGGCTACTTCTTATACT
+
AAA/A#EEE/EEEEEEEE/EEEEEEA/E##########################EEEEEEEAAEEEEEEE
And Read 2 looks similar but with 2:N:0:....
Any help in this matter will be greatly appreciated.
Thanks, Vaishnavi