Giter VIP home page Giter VIP logo

Comments (2)

heidihyang avatar heidihyang commented on September 15, 2024

Hi, I am no longer getting this error but the genotype folder is empty despite the program saying it's finished.

Below is my joblog:

(teflon_env) [hyangg@n1826 TEFLON]$ python ./teflon_genotype.py \

-wd $HOME/project-vlsork/TEFLON/
-d $HOME/project-vlsork/TEFLON/qlob.prep_TF/
-s $HOME/project-vlsork/TEFLON/samples.txt
-dt "pooled"
Lower-bound coverage threshold filters corresponding to samples ['ANG.5', 'WOO.4.w6'] is [1, 1]
NOTE: all sites with adjusted read counts > upper-bound coverage threshold will be marked -9
Upper-bound coverage threshold filters corresponding to samples ['ANG.5', 'WOO.4.w6'] is [82, 366]
NOTE: all sites with adjusted read counts > upper-bound coverage threshold will be marked -9
cdm: gunzip -c /u/home/h/hyangg/project-vlsork/TEFLON/qlob.prep_TF/qlob.pseudo2ref.pickle.gz > /u/home/h/hyangg/project-vlsork/TEFLON/qlob.prep_TF/qlob.pseudo2ref.pickle.gz.tmp
loading pickle: /u/home/h/hyangg/project-vlsork/TEFLON/qlob.prep_TF/qlob.pseudo2ref.pickle.gz.tmp
NOTE: this step can be time and memory intensive for large reference genomes
pickle loaded!
Converting coordinates from pseudospace to reference-based coordinates...
Traceback (most recent call last):
File "./teflon_genotype.py", line 122, in
main()
File "./teflon_genotype.py", line 116, in main
pt.pt_portal(countDir,genoDir,samples, posMap, stats, p2rC, l_thresh, h_thresh)
File "/u/project/vlsork/hyangg/TEFLON/teflon_scripts/genotyper_poolType.py", line 51, in pt_portal
p2rC.pseudo2refConvert_portal(outFILE1,posMap,outFILE2)
File "/u/project/vlsork/hyangg/TEFLON/teflon_scripts/pseudo2refConvert.py", line 50, in pseudo2refConvert_portal
ls[2]=pseudoMap[chrom][int(ls[2])]
IndexError: list index out of range

from teflon.

heidihyang avatar heidihyang commented on September 15, 2024

I solved the issue - RepeatMasker annotates for satellite repeats and low-complexity repeats in addition to TEs, but these were not in my TE hierarchy text. You can probably find a way to include them in the workflow but once I removed those from my TE reference bed file and re-ran the reference prep step everything worked smoothly.

from teflon.

Related Issues (16)

Recommend Projects

  • React photo React

    A declarative, efficient, and flexible JavaScript library for building user interfaces.

  • Vue.js photo Vue.js

    🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.

  • Typescript photo Typescript

    TypeScript is a superset of JavaScript that compiles to clean JavaScript output.

  • TensorFlow photo TensorFlow

    An Open Source Machine Learning Framework for Everyone

  • Django photo Django

    The Web framework for perfectionists with deadlines.

  • D3 photo D3

    Bring data to life with SVG, Canvas and HTML. 📊📈🎉

Recommend Topics

  • javascript

    JavaScript (JS) is a lightweight interpreted programming language with first-class functions.

  • web

    Some thing interesting about web. New door for the world.

  • server

    A server is a program made to process requests and deliver data to clients.

  • Machine learning

    Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.

  • Game

    Some thing interesting about game, make everyone happy.

Recommend Org

  • Facebook photo Facebook

    We are working to build community through open source technology. NB: members must have two-factor auth.

  • Microsoft photo Microsoft

    Open source projects and samples from Microsoft.

  • Google photo Google

    Google ❤️ Open Source for everyone.

  • D3 photo D3

    Data-Driven Documents codes.