Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • Missing FASTQ files in European 1KGP trio?

    I think that some FASTQ files containing exome reads from the European trio in the 1KGP project are missing.

    Read groups in an exome file (NA12878.chrom22.ILLUMINA.bwa.CEU.exon_targetted.20100311.bam):

    Code:
    daggy@daggy-Studio-XPS-8100:~/Downloads$ samtools view -H NA12878/exon/NA12878.chrom22.ILLUMINA.bwa.CEU.exon_targetted.20100311.bam
    ...
    @RG     ID:SRR013620    PL:ILLUMINA     LB:Solexa-5324  PI:0    DS:SRP000033    SM:NA12878      CN:BI
    @RG     ID:SRR013890    PL:ILLUMINA     LB:Solexa-7288  PI:0    DS:SRP000033    SM:NA12878      CN:BI
    @RG     ID:SRR014820    PL:ILLUMINA     LB:Solexa-10292 PI:300  DS:SRP000033    SM:NA12878      CN:BI
    @RG     ID:SRR015370    PL:ILLUMINA     LB:2773962043   PI:0    DS:SRP000033    SM:NA12878      CN:WUGSC
    @RG     ID:SRR016607    PL:ILLUMINA     LB:Solexa-11435 PI:201  DS:SRP000033    SM:NA12878      CN:BI
    @RG     ID:SRR017279    PL:ILLUMINA     LB:Solexa-7288  PI:0    DS:SRP000033    SM:NA12878      CN:BI
    If you take a look at the FASTQ raw reads, you will find all read groups but SRR013620 for NA12878. For NA12891, the read group SRR013888 is missing; and for NA12892, SRR013619. The African trio (NA19238-40) is OK.

    But the reads in read group SRR013620 are present in the BAM file:
    Code:
    daggy@daggy-Studio-XPS-8100:~/Downloads$ samtools view NA12878.chrom22.ILLUMINA.bwa.CEU.exon_targetted.20100311.bam | grep SRR013620 | wc
     142742 3118113 41310387
    After searching on Google, I found the read group SRR013620 at an EBI site.

    The files at the 1KGP FTP site are filtered ("filt" in the file names). Those at the EBI site do not have that substring suggesting that they are unfiltered. A comparison of the files SRR013890.filt.fastq.gz and SRR013890/SRR013890.fastq.gz shows that the 1KGP is smaller (294 MB vs 303 MB), confirming the hypothesis.

    I found in another thread the sequence index file of 1KGP that the md5sum is missing (described here), suggesting that the missing FASTQ files are not supposed to exist.

    Can anybody tell me what is going on here? Are some FASTQ files missing, or located somewhere else? Can I use the FASTQ files at the EBI site, and do they have lower quality if unfiltered?

Latest Articles

Collapse

  • seqadmin
    Essential Discoveries and Tools in Epitranscriptomics
    by seqadmin




    The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist...
    04-22-2024, 07:01 AM
  • seqadmin
    Current Approaches to Protein Sequencing
    by seqadmin


    Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...
    04-04-2024, 04:25 PM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, Today, 08:47 AM
0 responses
12 views
0 likes
Last Post seqadmin  
Started by seqadmin, 04-11-2024, 12:08 PM
0 responses
60 views
0 likes
Last Post seqadmin  
Started by seqadmin, 04-10-2024, 10:19 PM
0 responses
59 views
0 likes
Last Post seqadmin  
Started by seqadmin, 04-10-2024, 09:21 AM
0 responses
54 views
0 likes
Last Post seqadmin  
Working...
X