Skip to content

Commit

Permalink
Update run-marker-blast-annotation-illumina.md
Browse files Browse the repository at this point in the history
  • Loading branch information
kbeutel authored Feb 29, 2024
1 parent 0ffa740 commit a8460b9
Showing 1 changed file with 5 additions and 3 deletions.
Original file line number Diff line number Diff line change
@@ -1,8 +1,10 @@
## Run Marker BLAST Annotation

[BLAST (Basic Local Alignment Search Tool)](https://blast.ncbi.nlm.nih.gov/Blast.cgi) (external link) determines if markers are useful by aligning the probe sequences to the genome. A marker is not useful if it aligns to multiple locations on a genome.
[BLAST (Basic Local Alignment Search Tool)](https://blast.ncbi.nlm.nih.gov/Blast.cgi) (external link) is used to determine if markers are problematic by aligning the probe sequences to the genome. A marker is likely problematic if it aligns to multiple locations on a genome. BLAST will run on the marker positions in [markerPositions.txt](../#/documentation/RunTheGenvisisWorkflow--create-marker-positions-illumina).

BLAST will run on the marker positions in [markerPositions.txt](../#/documentation/RunTheGenvisisWorkflow--create-marker-positions-illumina).
If you do not have a [probe sequence file](../#/documentation/GetStarted--find-required-raw-data-files), Genvisis cannot run the Run Marker BLAST Annotation step. Verify that Run Marker BLAST Annotation is unchecked in the Genvisis workflow. If Genvisis does not run the Run Marker BLAST Annotation step, it will automatically identify markers that are likely to be problematic based on mappability (which uses the reference genome instead of probe sequences to identify regions of the genome that appear multiple times and thus may be problematic).

BLAST is superior to mappability. If a probe set file is available, we recommend running Run Marker BLAST Annotation rather than relying on mappability.

#### 1: Parse Sample Files step [parse-samples] must have been run already or must be selected
This step uses the sampRAF files that were created in the [Parse Samples Files](../#/documentation/RunTheGenvisisWorkflow--parse-sample-files-illumina) step.
Expand All @@ -14,7 +16,7 @@ Markers are 50 bp in length
Checks if subunits of the marker have matches in the sample genome?

#### 3: Full path to an Illumina manifest file (e.g., HumanExome-12-v1-o-B.csv)
You must provide a manifest with probe sequences.
You must provide a manifest with probe sequences. Contact Illumina if you need assistance locating a manifest file that contains probe sequences.

#### 4: (optional) Marker positions file to override MarkerSet positions with (e.g. for a different genome build) (e.g. markerPositions.txt)
Use this option if you want to BLAST on a different genome build than that of your project. For example, if the project was created in hg37 and the positions in *markerPositions.txt* are hg37, you can provide an alternative markerPositions file with positions in hg38 for this part.
Expand Down

0 comments on commit a8460b9

Please sign in to comment.