diff --git a/private/find_template_cell b/private/find_template_cell deleted file mode 100755 index 2349d215e9302d9bb859f146e3f73e55201b74e4..0000000000000000000000000000000000000000 Binary files a/private/find_template_cell and /dev/null differ diff --git a/private/genotype_mtx b/private/genotype_mtx deleted file mode 100755 index c535f52659c7803e1eb943eef8f990666301d1b4..0000000000000000000000000000000000000000 Binary files a/private/genotype_mtx and /dev/null differ diff --git a/private/genotype_mtx.nim b/private/genotype_mtx.nim index 66feb5819811ea34cbe1d6bd61512cc5f893315f..a935cb1e645cb349f6c4017effb65e942c98089f 100755 --- a/private/genotype_mtx.nim +++ b/private/genotype_mtx.nim @@ -8,8 +8,7 @@ import strutils import hts import utils import streams - -import fragments +import findGtNodes ## these outputs will be per chromosome @@ -57,8 +56,8 @@ proc getGtMtx(ibam:Bam, ivcf:VCF, barcodeTable:TableRef, outGtMtx:FileStream, ou outGtMtx.writeLine(join([$writtenVarintIndex, $cellIndex, $geno],sep = " ")) outTotalCountMtx.writeLine(join([$writtenVarintIndex, $cellIndex, $cellDP],sep = " ")) outAltCountMtx.writeLine(join([$writtenVarintIndex, $cellIndex, $calt],sep = " ")) - outSnpAnnot.writeLine(join([$rec.POS, $rec_ref, $rec_alt], sep="\t") ) nnsize.inc + outSnpAnnot.writeLine(join([$rec.POS, $rec_ref, $rec_alt], sep="\t") ) writtenVarintIndex.inc ## write to matrices @@ -101,13 +100,14 @@ while hts_getline(hf, cint(10), addr kstr) > 0: discard hf.hts_close() let s_Chrs = @["chr1"] +var chrom = "chr1" let out_prefix = "test_data/getMtx/" echo "generating gtMtx" try: - outGtMtx = openFileStream(out_prefix & "_gtMtx.mtx", fmWrite) - outSnpAnnot = openFileStream(out_prefix & "_snpAnnot.txt", fmWrite) - outTotalCountMtx = openFileStream(out_prefix & "_totalMtx.mtx", fmWrite) - outAltCountMtx = openFileStream(out_prefix & "_AltMtx.mtx", fmWrite) + outGtMtx = openFileStream(out_prefix & chrom & "_gtMtx.mtx", fmWrite) + outSnpAnnot = openFileStream(out_prefix & chrom & "_snpAnnot.txt", fmWrite) + outTotalCountMtx = openFileStream(out_prefix & chrom & "_totalMtx.mtx", fmWrite) + outAltCountMtx = openFileStream(out_prefix & chrom & "_AltMtx.mtx", fmWrite) except: stderr.write getCurrentExceptionMsg() diff --git a/private/infer_missing_snps b/private/infer_missing_snps deleted file mode 100755 index 85df612e74316fd0810c1ef6ddc70cf358b34890..0000000000000000000000000000000000000000 Binary files a/private/infer_missing_snps and /dev/null differ