Entering edit mode
3.0 years ago
Fizzah
▴
30
Hello there, I am trying to generate histogram for paired end reads via bbmap tool program. However I am getting this output:
root@cemb-ls1-pgs:/mnt/cemb/S017679/raw# /home/fizzah/bbmap/reformat.sh BU21R1.fastq qchist=qchist.txt
java -ea -Xms300m -cp /home/fizzah/bbmap/current/ jgi.ReformatReads BU21R1.fastq qchist=qchist.txt
Executing jgi.ReformatReads [BU21R1.fastq, qchist=qchist.txt]
Set qcount histogram output to qchist.txt
No output stream specified. To write to stdout, please specify 'out=stdout.fq' or similar.
Input is being processed as unpaired
Input: 37926174 reads 5726852274 bases
Output: 37926174 reads (100.00%) 5726852274 bases (100.00%)
Time: 50.960 seconds.
Reads Processed: 37926k 744.23k reads/sec
Bases Processed: 5726m 112.38m bases/sec
when I specify stream output, it says:
File qchist.txt exists and overwrite=false
Please guide about its solution
BBMap does not produce graphical output. What you will get is a text based histogram.
In general for all BBMap programs you have to explicitly add option
ow=t
to overwrite existing files. Be aware that this option instantly wipes the contents of existing files so if you are not careful you can lose data.Ok. So can you name some tool that can give graphical outcome based histogram for paired end reads?