| 1835 | } |
| 1836 | |
| 1837 | void filterWithQualities(string & pileupShort,string & qualityS, int minimalQualityPerPosition) { |
| 1838 | string newString; |
| 1839 | for (unsigned int i = 0; i< pileupShort.length(); i++) { |
| 1840 | int q = int (qualityS[i]); |
| 1841 | if (q>=minimalQualityPerPosition) |
| 1842 | newString+=pileupShort[i]; |
| 1843 | } |
| 1844 | pileupShort=newString; |
| 1845 | return; |
| 1846 | |
| 1847 | } |
| 1848 | |
| 1849 | char complement(char n) { |
| 1850 | switch(n) { |