Uploaded image for project: 'IGB'
  1. IGB
  2. IGBF-3701

Muday time course: Check first sequences in renamed fastq VS. what is in SRA

    Details

    • Type: Task
    • Status: Closed (View Workflow)
    • Priority: Minor
    • Resolution: Done
    • Affects Version/s: None
    • Fix Version/s: None
    • Labels:
      None

      Description

      To ensure that the sequences are matching what we have versus what the SRA has, I came up with a simple check. We look for the first 10 sequences of a given SRA run across all the sequence files.

      For speed, I head the fastq.gz file to get only the first 5 sequences.
      Use zcat to keep files zipped.
      Grep for the sequence of interest.

      for f in *gz; do echo $f; zcat < $f | head -n 20 | grep "^CTGGCTTTTC" ; done

      Probability dictates that we should find just 1 result.

        Attachments

          Issue Links

            Activity

            robofjoy Robert Reid created issue -
            robofjoy Robert Reid made changes -
            Field Original Value New Value
            Epic Link IGBF-2993 [ 21429 ]
            robofjoy Robert Reid made changes -
            Link This issue blocks IGBF-3683 [ IGBF-3683 ]
            robofjoy Robert Reid made changes -
            Status To-Do [ 10305 ] In Progress [ 3 ]
            robofjoy Robert Reid made changes -
            Status In Progress [ 3 ] Needs 1st Level Review [ 10005 ]
            robofjoy Robert Reid made changes -
            Status Needs 1st Level Review [ 10005 ] First Level Review in Progress [ 10301 ]
            robofjoy Robert Reid made changes -
            Status First Level Review in Progress [ 10301 ] Ready for Pull Request [ 10304 ]
            robofjoy Robert Reid made changes -
            Status Ready for Pull Request [ 10304 ] Pull Request Submitted [ 10101 ]
            robofjoy Robert Reid made changes -
            Status Pull Request Submitted [ 10101 ] Reviewing Pull Request [ 10303 ]
            robofjoy Robert Reid made changes -
            Status Reviewing Pull Request [ 10303 ] Merged Needs Testing [ 10002 ]
            robofjoy Robert Reid made changes -
            Status Merged Needs Testing [ 10002 ] Post-merge Testing In Progress [ 10003 ]
            robofjoy Robert Reid made changes -
            Resolution Done [ 10000 ]
            Status Post-merge Testing In Progress [ 10003 ] Closed [ 6 ]

              People

              • Assignee:
                robofjoy Robert Reid
                Reporter:
                robofjoy Robert Reid
              • Votes:
                0 Vote for this issue
                Watchers:
                Start watching this issue

                Dates

                • Created:
                  Updated:
                  Resolved: