Good morning to everyone, I'm using the boxplot function as suggested in the RnaSeqEdgeRQL to show a summary graph of the quality of my input .FASTQ files, which contain my RNASeq input reads (Go to the Paragraph 'Accuracy of base-calling'). As you can see in the provided example, the x-axis shows the bases positions of the analysed reads. So, in the suggested example, with reads up to 100 bp in length, the scales ranges from 1 to 100.
In my case, my reads are from Roche454 pyrosequencing, so the reads are longer (up to 400 bp, 150-200 on average).
I used the following code:
QS <- qualityScores("pathofmyfile/FD.fastq") boxplot(QS, ylab="Quality score", xlab="Base position", main="pathofmyfile/FD.fastq", cex=0.25, col="red")
To get the resulting plot:
But my reads are longer than 150 bp... why aren't all the positions shown? Am I missing something? Could it depend on the previous qualityScores function? I get this warning. I don't think it's related since it involves the y-axis, but I'm reporting it anyway.
qualityScores Rsubread 1.32.4 Scan the input file... Totally 105191 reads were scanned; the sampling interval is 10. Now extract read quality information... Warning: the Phred score offset (33) seems wrong : 0. Completed successfully. Quality scores for 10000 reads (equally spaced in the file) are returned. However, the Phred score offset (33) seemed to be wrong. The quality scores can be meaningless.
Any idea why my x-axis isn't complete? Thanks in advance!