Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • #16
    Originally posted by SylvainL View Post
    Looking at your reads, it seems you had a problem when you extracted them from the CRAm files. They look exactly similar with a shift of 3-4 bp...
    I peeked into the original cram file using the following command and the result is the same Fastq. So i don't know if the Cram decompression to Fastq has caused the problem:

    fazal@fazal-Precision-T1700:/media/fazal/backup/BCL11A/Cram$ java -jar /opt/cramtools-3.0.jar fastq -I 18418_2#1.cram | head -20
    @HS32_18418:2:2307:11553:47098#1/1
    CCTAACCCTAACCCTAACCCTAACCCTAACCCTAACCCTAACCCTAACCCTAACCCTAACCCTAACCCTAACCCT
    +
    BBBBBBF/FFFFFFFFFFFFFFFFFFFFFFFFFFFBBFFFB</<FFFF//FFFFFFFBBF<FFFFFFBF/B<FF/
    @HS32_18418:2:2307:11553:47098#1/2
    TAACCCTACCCCTAACCCTAACCCTAACCCTAACCCTAACCCTAACCCTAACCCTTACCCTAACCCTAACCCAAC
    +
    </FBB<///F<BF/B/F<FB<<FB<B<FFF/BFFBFBBB<FFB<FFFFBFFFFFF<FFFF/FFBFFF<FFBBBBB
    @HS32_18418:2:1201:20716:93279#1/2
    ACCCTAACCCTAACCCTAACCCTAACCCTAACCCTAACCCTAACCCTAACCCTAACCCTAACCCTAACCCTAACC
    +
    BBBBBFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFBFFFFF//FFB
    @HS32_18418:2:1201:20716:93279#1/1
    AACCCTAACCCTAACCCTAACCCTAACCCTAACCCTAACCCTAACCCTAACCCTAACCCTAACCCTAACCCAAAC
    +
    /<FFF</<FBF<</FFFBF<FFFFBBFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFBBBBB
    @HS32_18418:2:1102:8324:84406#1/1
    CTAACCCTAACCCTAACCCTAACCCTAACCCTAACCCTAACCCTAACCCTAACCCTAACCCTAACCCAAACCCTA
    +
    BBBB</<<FFFFFBBF<<BFFBB/F/<<<<F/FFFFFBBB<F//B/F/<FBFFB//</BFFBB/</////</7/<

    Comment


    • #17
      Originally posted by HESmith View Post
      These sequences (CCCTAA) are from the telomere repeat. CRAMs are sorted files so, depending on the reference used for alignment, it's not unusual to see all of these reads at the beginning of the file.

      What IS unusual is that both reads are from the same strand, when one should be the reverse complement of the other (e.g., TTAGGG). Not sure how that happened, but it's likely to be the reason why realignment failed.
      Hi HESmith,
      Bear with me for the silly and naive question i am about to ask: In PE chipseq, pair of reads should be completely reverse complementary to each other or there is like a region in which the reads are reverse complementary?

      Comment


      • #18
        Originally posted by fh331 View Post
        Hi HESmith,
        Bear with me for the silly and naive question i am about to ask: In PE chipseq, pair of reads should be completely reverse complementary to each other or there is like a region in which the reads are reverse complementary?
        They will only be completely rev-comp of each other, IF the insert size=number of F/R sequencing cycles (not likely for every fragment in your library).

        Comment


        • #19
          Originally posted by SylvainL View Post
          Are you sure the reads in file 1 and file 2 are in the same order? Just print the first 5 reads of each file...

          You can also try to map only one file (not considered as paired-end then) to see the percentage of mapped reads...
          Hi SylvainL,
          I tried to align single file from the PE. Here's the output report and worked fine i think:

          fazal@fazal-Precision-T1700:/media/fazal/backup/BCL11A/FastQ_Files$ bowtie -m 1 -S /media/fazal/backup/BCL11A/Bowtie_Indices/Bowtie1_Index/human_g1k_v37.fasta 18418_2-1_1.fastq > /media/fazal/backup/BCL11A/Sam_from_FastQ/18418_2-1_1.sam
          # reads processed: 40095362
          # reads with at least one reported alignment: 33362937 (83.21%)
          # reads that failed to align: 2434332 (6.07%)
          # reads with alignments suppressed due to -m: 4298093 (10.72%)
          Reported 33362937 alignments to 1 output stream(s)

          No I don't know what's going wrong with aligning the two PEs together?

          I am kind of getting frustrated

          Comment


          • #20
            Like @HESmith said, you may have a problem with the orientation of your read pairs.

            The default orientation for valid alignments in bowtie1 is --fr. If your read pairs got converted somehow to ff, then none of the alignments would be valid.

            Comment


            • #21
              Probably the wrong orientation came from the exportation of the reads from the CRAM files. As fanli said, change the setting in bowtie (or reverse complement the reads 2)

              Comment

              Latest Articles

              Collapse

              • seqadmin
                Techniques and Challenges in Conservation Genomics
                by seqadmin



                The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

                Avian Conservation
                Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
                03-08-2024, 10:41 AM
              • seqadmin
                The Impact of AI in Genomic Medicine
                by seqadmin



                Artificial intelligence (AI) has evolved from a futuristic vision to a mainstream technology, highlighted by the introduction of tools like OpenAI's ChatGPT and Google's Gemini. In recent years, AI has become increasingly integrated into the field of genomics. This integration has enabled new scientific discoveries while simultaneously raising important ethical questions1. Interviews with two researchers at the center of this intersection provide insightful perspectives into...
                02-26-2024, 02:07 PM

              ad_right_rmr

              Collapse

              News

              Collapse

              Topics Statistics Last Post
              Started by seqadmin, 03-14-2024, 06:13 AM
              0 responses
              34 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 03-08-2024, 08:03 AM
              0 responses
              72 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 03-07-2024, 08:13 AM
              0 responses
              81 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 03-06-2024, 09:51 AM
              0 responses
              68 views
              0 likes
              Last Post seqadmin  
              Working...
              X