Unconfigured Ad

Collapse
This topic is closed.
X
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • GenoMax
    Senior Member
    • Feb 2008
    • 7142

    @Phil: Is a better fix coming along for this? That kind of seems like a hacky solution.

    Comment

    • ewels
      Phil Ewels
      • Mar 2011
      • 32

      Originally posted by GenoMax View Post
      @Phil: Is a better fix coming along for this? That kind of seems like a hacky solution.
      No idea You'll have to ask Simon about that..

      Comment

      • fnn4
        Junior Member
        • Jul 2014
        • 8

        Hello,

        Does any one know where I can find average Q score on miseq raw reads fastQC report? I know fastQC gives per sequence Q scores. But it doesn't tell the average Q scores of the raw reads though. Thank you very much in advance.

        Comment

        • ewels
          Phil Ewels
          • Mar 2011
          • 32

          Hi @fnn4,

          Do you mean a single mean quality score for the entire FastQ file? The Per Sequence Quality Scores plot shows the distribution of average Q scores across all reads in the file. I don't think that FastQC reports an average for the entire file, but it should be simple enough to do using the raw data reported in fastqc_data.txt - just sum the (Q score * read count) from the per-sequence section and divide by the sum of the read count.

          What's your reason for wanting this statistic, out of interest?

          Phil

          Comment

          • fnn4
            Junior Member
            • Jul 2014
            • 8

            Hi Phil,

            Thank you for your reply. Yes, that is what I was asking for. When we submit samples, we have to pass a certain threshold for average Q score.
            Originally posted by ewels View Post
            Hi @fnn4,

            Do you mean a single mean quality score for the entire FastQ file? The Per Sequence Quality Scores plot shows the distribution of average Q scores across all reads in the file. I don't think that FastQC reports an average for the entire file, but it should be simple enough to do using the raw data reported in fastqc_data.txt - just sum the (Q score * read count) from the per-sequence section and divide by the sum of the read count.

            What's your reason for wanting this statistic, out of interest?

            Phil

            Comment

            • Guillefriis86
              Junior Member
              • Jan 2014
              • 2

              FastQC weird Kmer count in GBS data

              Hi @simonandrews and everybody elese, I'm using FastQC to check the quality of my GBS reads, focusing in th R1 sequences, and I'm getting weird results in the Kmer content. After trimming and removing adapters and barcodes, I got a first Kmer profile showing a higher than expected count for several 6-mers, like this:
              Click image for larger version

Name:	kmer_profiles.png
Views:	1
Size:	21.4 KB
ID:	307201

              I thought it may be the adapters of the other side of the read for the sequencing of the R2, so I trimmed all the reads longer than 130 pb to 130 pb. And strangely, I got new Kmer peaks, not shown before:

              Click image for larger version

Name:	kmer_profiles_crop.png
Views:	1
Size:	32.9 KB
ID:	307202

              Does anybody know why and what are these new peaks, or what it's the origin of them? May be related with FastQC algorithms? Somebody found similar problems? Thanks a lot guys!!

              Comment

              • hiretabletsae
                Junior Member
                • Feb 2019
                • 1

                I most want to say that you are doing a great job.
                Last edited by hiretabletsae; 02-06-2019, 10:07 PM.
                https://www.hiretablets.ae/touch-screen-monitor-rental

                Comment

                Latest Articles

                Collapse

                ad_right_rmr

                Collapse

                News

                Collapse

                Topics Statistics Last Post
                Started by SEQadmin2, Yesterday, 10:09 AM
                0 responses
                10 views
                0 reactions
                Last Post SEQadmin2  
                Started by SEQadmin2, 06-04-2026, 08:59 AM
                0 responses
                18 views
                0 reactions
                Last Post SEQadmin2  
                Started by SEQadmin2, 06-02-2026, 12:03 PM
                0 responses
                26 views
                0 reactions
                Last Post SEQadmin2  
                Started by SEQadmin2, 06-02-2026, 11:40 AM
                0 responses
                21 views
                0 reactions
                Last Post SEQadmin2  
                Working...