cgeb2001's DokuWiki!

User Tools

Site Tools

Trace: • bioinformatics_essay11 • bioinformatics_essay • others
gene_prediction_curation_with_igv

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revisionPrevious revision
Next revision		Previous revision
gene_prediction_curation_with_igv [2023/12/08 11:24] – [Strand-specific RNAseq reads] 134.190.232.149gene_prediction_curation_with_igv [2025/08/07 16:01] (current) – 134.190.145.228
Line 49: Line 49:
  
 You can play with the window size ''-W'' and overlap ''-s'' parameters to fine tune it to your own liking. You can either load in this track via ''File -> Load From File ..'' or through the JSON file (see below) You can play with the window size ''-W'' and overlap ''-s'' parameters to fine tune it to your own liking. You can either load in this track via ''File -> Load From File ..'' or through the JSON file (see below)
 +
 ==== Loading into IGV ==== ==== Loading into IGV ====
  
Line 228: Line 229:
  
 For some frustratingly unknow reason, IGV only displays the three forward frames, OR the three reverse frames. It is, as far as I know, not possible to display all six frames simultaneously. To swap between reverse and forward frames, click on small arrow of the //Sequence// track. To see the sequence track, you need to be sufficiently zoomed in. For some frustratingly unknow reason, IGV only displays the three forward frames, OR the three reverse frames. It is, as far as I know, not possible to display all six frames simultaneously. To swap between reverse and forward frames, click on small arrow of the //Sequence// track. To see the sequence track, you need to be sufficiently zoomed in.
 +
 ==== Curating gene models ==== ==== Curating gene models ====
  
Line 354: Line 356:
   * Does the translate CDS start with an M and end with a * (STOP)?   * Does the translate CDS start with an M and end with a * (STOP)?
   * Are there any premature * (STOP) in the CDS?   * Are there any premature * (STOP) in the CDS?
 +  * If Blastocystis genome, and there is no STOP, does it end with a TGTTTGTT motif?
  
 <code> <code>
Line 369: Line 372:
 </code> </code>
  
-NOTE: Currently, the checks performed in the script assume that you haven't annotated any UTRs in your GFF3. If you have UTR (''three_prime_UTR'' or ''five_prime_UTR'') features, for example the start coordinate of the first CDS will not match the start coordinate of the gene. It will thus throw an error, even though your GFF3 is perfectly fine. Something I'll have to fix in the future.+<del>NOTE: Currently, the checks performed in the script assume that you haven't annotated any UTRs in your GFF3. If you have UTR (''three_prime_UTR'' or ''five_prime_UTR'') features, for example the start coordinate of the first CDS will not match the start coordinate of the gene. It will thus throw an error, even though your GFF3 is perfectly fine. Something I'll have to fix in the future.</del>
  
-TODO: also make it so that it collects all errors are reports them at the end, rather than quitting at the first error.+UPDATE (Jan 2025): It should now account for UTRs
  
 ==== Re-streamlining gene names ==== ==== Re-streamlining gene names ====
gene_prediction_curation_with_igv.1702049069.txt.gz · Last modified: by 134.190.232.149