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
    Best Practices for Single-Cell Sequencing Analysis
    by seqadmin



    While isolating and preparing single cells for sequencing was historically the bottleneck, recent technological advancements have shifted the challenge to data analysis. This highlights the rapidly evolving nature of single-cell sequencing. The inherent complexity of single-cell analysis has intensified with the surge in data volume and the incorporation of diverse and more complex datasets. This article explores the challenges in analysis, examines common pitfalls, offers...
    Yesterday, 07:15 AM
  • seqadmin
    Latest Developments in Precision Medicine
    by seqadmin



    Technological advances have led to drastic improvements in the field of precision medicine, enabling more personalized approaches to treatment. This article explores four leading groups that are overcoming many of the challenges of genomic profiling and precision medicine through their innovative platforms and technologies.

    Somatic Genomics
    “We have such a tremendous amount of genetic diversity that exists within each of us, and not just between us as individuals,”...
    05-24-2024, 01:16 PM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, Today, 06:58 AM
0 responses
13 views
0 likes
Last Post seqadmin  
Started by seqadmin, Yesterday, 08:18 AM
0 responses
19 views
0 likes
Last Post seqadmin  
Started by seqadmin, Yesterday, 08:04 AM
0 responses
18 views
0 likes
Last Post seqadmin  
Started by seqadmin, 06-03-2024, 06:55 AM
0 responses
13 views
0 likes
Last Post seqadmin  
Working...
X