SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Yes .. BBMap can do that! GenoMax Bioinformatics 213 10-17-2018 12:06 AM
BBmap dedupe help JamesSeward Bioinformatics 7 07-15-2016 11:20 PM
Introducing CalcTrueQuality, a quality-score recalibrator Brian Bushnell Bioinformatics 15 11-10-2015 05:15 PM
bbmap aligner cmccabe Bioinformatics 22 11-03-2015 01:03 PM
BBMap for BitSeq dietmar13 Bioinformatics 1 04-30-2015 09:40 AM

Reply
 
Thread Tools
Old 05-08-2018, 01:57 PM   #1
Coyk
Member
 
Location: West Virginia

Join Date: Jun 2011
Posts: 14
Default calctruequality in bbmap

I'm trying to recalibrate Q scores of a NextSeq run using MiSeq contigs assembled with Tadpole. The commands I used to map the reads to the reference were as follows:

>bbmap.sh in=concatABC.fastq.gz outm=mapped.sam ref=./Lpe09_06TdpAssemblies/contigs09_06.fa ignorequality maxindel=100 minratio=0.4 ambig=toss qahist=qahist_raw.txt qhist=qhist_raw.txt mhist=mhist_raw.txt

The command I used to get the calibration matrices was:
>calctruequality.sh in=mapped.sam

I get the following output
java -ea -Xmx57992m -Xms57992m -cp /home/me/bbmap/current/ jgi.CalcTrueQuality in=mapped.sam
Executing jgi.CalcTrueQuality [in=mapped.sam]

Exception in thread "Thread-2" Exception in thread "Thread-1" java.lang.AssertionError: TODO: Encountered a read with 'M' in cigar string but no MD tag and no ScafMap loaded.
at stream.SamLine.toShortMatch(SamLine.java:1212)
at stream.SamLine.toRead (SamLine.java:2015)
at stream.SamLine.toRead (SamLine.java:1875)
at stream.SamReadStreamer$ProcessThread.makeReads (SamReadStreamer.java:206)
at stream.SamReadStreamer$ProcessThread.run (SamReadStreamer.java:135)
java.lang.AssertionError: TODO: Encountered a read with 'M' in cigar string but no MD tag and no ScafMap loaded.
at stream.SamLine.toShortMatch(SamLine.java:1212)
at stream.SamLine.toRead (SamLine.java:2015)
at stream.SamLine.toRead (SamLine.java:1875)
at stream.SamReadStreamer$ProcessThread.makeReads (SamReadStreamer.java:206)
at stream.SamReadStreamer$ProcessThread.run (SamReadStreamer.java:135)
Exception in thread "Thread-3" java.lang.AssertionError: TODO: Encountered a read with 'M' in cigar string but no MD tag and no ScafMap loaded.
at stream.SamLine.toShortMatch(SamLine.java:1212)
at stream.SamLine.toRead (SamLine.java:2015)
at stream.SamLine.toRead (SamLine.java:1875)
at stream.SamReadStreamer$ProcessThread.makeReads (SamReadStreamer.java:206)
at stream.SamReadStreamer$ProcessThread.run (SamReadStreamer.java:135)

The program never quit so I don't know if it had stopped or was hanging but I let it sit there for over an hour. I didn't see /ref/qual/ files being generated. I'm not sure how to proceed. Since the mapped sam file was generated with bbmap using contigs from Tadpole, I'm not sure why there are Ms in the cigar strings. The infile I used for bbmap was made by first concatenating all Read1 from 3 libraries then all Read2 from 3 libraries then those two files concatenated.
Coyk is offline   Reply With Quote
Reply

Tags
bbmap, bbtools, calctruequality, nextseq 500

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 06:26 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2018, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO