Comments (2)
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.
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)
- teflon.v0.4.py IndexError: list index out of range HOT 3
- Killed HOT 1
- Advice on parsing results and getting allele frequencies HOT 2
- Index BAM with BWA or Samtools? HOT 8
- The genotype folder empty. HOT 4
- Interpreting the breakpoint coordinates HOT 2
- malformed mega_clustered.bed HOT 6
- No reference genome in prefix.prep_MP, but multiple in prep_RM HOT 3
- error after collapsing samples HOT 3
- teflon_genotype.py fails
- IUPAC code in TE database
- Error when I try to run sample_pipeline.sh
- Standard deviation warning HOT 1
- New error
- IndexError in pseudoGenerate.py due to Incorrect BED file format
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
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.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from teflon.