Unconfigured Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • hujiacun
    Member
    • Mar 2015
    • 17

    how can I retrieve one read sequence?

    Hi, I use Integrative genomic Viewer to find one Read which mate is mapped=no.
    This no mapped pair-end read may be a result from structural variation. So, if I know one read ID and its one end sequence, How can I get its pair-end sequence from raw Data.

    Thanks!
  • GenoMax
    Senior Member
    • Feb 2008
    • 7142

    #2
    If your sequence files are not compressed:

    Code:
    $ cat file.fastq | grep -A 3 "read_ID_you_want"
    If they are compressed:

    Code:
    $ zcat file.fastq.gz | grep -A 3 "read_ID_you_want"
    If you only want sequence (not full fastq record) then replace 3 with 1.

    Comment

    Latest Articles

    Collapse

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by SEQadmin2, Today, 11:58 AM
    0 responses
    6 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 06-05-2026, 10:09 AM
    0 responses
    21 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 06-04-2026, 08:59 AM
    0 responses
    34 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 06-02-2026, 12:03 PM
    0 responses
    55 views
    0 reactions
    Last Post SEQadmin2  
    Working...