Skip to content

Commit a67e566

Browse files
committed
cleaning
1 parent b4b499d commit a67e566

File tree

1 file changed

+0
-35
lines changed

1 file changed

+0
-35
lines changed

sequence/bamreader.cpp

Lines changed: 0 additions & 35 deletions
Original file line numberDiff line numberDiff line change
@@ -51,40 +51,5 @@ bool BamReader::next()
5151
sam_close(fp_in);
5252
}
5353

54-
55-
56-
// LOOP SNIPPET
57-
// while(sam_read1(fp_in,bamHdr,aln) > 0){
58-
59-
// int32_t pos = aln->core.pos +1; //left most position of alignment in zero based coordianate (+1)
60-
// char *chr = bamHdr->target_name[aln->core.tid] ; //contig name (chromosome)
61-
// uint32_t len = aln->core.l_qseq; //length of the read.
62-
63-
// uint8_t *q = bam_get_seq(aln); //quality string
64-
// uint32_t q2 = aln->core.qual ; //mapping quality
65-
66-
// qDebug()<<chr;
67-
68-
69-
// char *qseq = (char *)malloc(len);
70-
71-
// for(int i=0; i< len ; i++){
72-
// qseq[i] = seq_nt16_str[bam_seqi(q,i)]; //gets nucleotide id and converts them into IUPAC id.
73-
// }
74-
75-
76-
77-
// qDebug()<<chr<<pos<<len<<qseq<<q<<q2;
78-
79-
// // if(strcmp(chrom, chr) == 0){
80-
81-
// // if(locus > pos+len){
82-
// // printf("%s\t%d\t%d\t%s\t%s\t%d\n",chr,pos,len,qseq,q,q2);
83-
// // }
84-
// // }
85-
// }
86-
87-
88-
8954
return false;
9055
}

0 commit comments

Comments
 (0)