mageri icon indicating copy to clipboard operation
mageri copied to clipboard

Error while running with UMI

Open VNagesh-Bio opened this issue 7 years ago • 6 comments

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

VNagesh-Bio avatar Aug 11 '17 16:08 VNagesh-Bio