Advertisement
Guest User

Untitled

a guest
Jul 27th, 2017
52
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.41 KB | None | 0 0
  1. continue
  2. num_snps_skipped += 1
  3. samp_id = sline[id_index]
  4. ref_allele = sline[ref_allele_index]
  5. tum_allele = sline[tum_allele_index]
  6. snp = ref_allele + tum_allele
  7. if not snp in transitions:
  8. snp = nucleotide_complement[ref_allele] + nucleotide_complement[tum_allele]
  9. ref_trinuc = sline[ref_tri_index]
  10. if ref_trinuc == "NA":
  11. print "Warning: Reference allele not available on "+
  12. "line %d; skipping line"%line_number
  13. continue
  14. if not ref_trinuc[1] == snp[0]:
  15. print "Warning: Reference allele does not match reference "+
  16. "trinucleotide; skipping line %d"%line_number**
  17. continue
  18. snp_with_ctx = ref_trinuc[0] + snp + ref_trinuc[2]
  19. if not samp_id in signatures:
  20. signatures[samp_id] = [0 for i in substitution_order]
  21. if snp_with_ctx not in substitution_order:
  22. print "Warning: substitution on line " +
  23. "%d is %s, not "%(line_number,snp_with_ctx) +
  24. "found among possible substitutions. Skipping line."
  25.  
  26. Traceback (most recent call last):
  27. File "main.py", line 59, in <module>
  28. signatures = signature.make(args.in_file, substitution_order=stratton['substitution_order'], out_path = args.spectrum_output)
  29. File "/home/ateeqanees/Mutation/centos/mutation-signatures-master/signature.py", line 73, in make
  30. if not ref_trinuc[1] == snp[0]:
  31. IndexError: string index out of range
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement