quick_ed_validation.c

#include <stdio.h>
#include <string.h>

int main(int argc, char *argv[]) {
	char *a = strdup("TTATG-ACAAATCGTGAGAAGCG-TGGGCA-C-TCAGGATGCCACAGAAATCAGAGCATCTGTC-GGAAGAATGATTGGTCGAATCGGGAGATTCTACATCCAAATGTGCACTGAACTCAAACTCAGTGATTATGATGGACGACTAATCCAGAATAGCATAACAATAGAGGGGGATGGTGCTTTCTGC-TTTGATGAGAGAAGAAATAAATA--TAGAAG--CATCCGAGTGC-TGGAAGGTACCTAAGAACA-GGGAGGTACCATAT--AGAAGAGTAGACGGAAAGATGATATCAGAGAACTCATCCTTTATGACAAAGAAGAAATAAGG-GAGTTTGGCGCCAAGCAAACAATGGCGAAGATGCAACAGCAGGTCTTACTCATATCATGATTTGGCAATTCCAACCTGAATGATGCCACAT--CAGA--ACAAGACCGCTTGTTCGCACCGGAATCGATCCCAGACT-TG--CTCTAATGCAAGGTTCAACACTTCCCAGAAGGTCTGGTGCCGCAGGTGCTGCGTAGCTGAAAGGGAGTTGGAACAATAGCAATGGAGTTAATCATCAGAATGATCAAACGTGGAATCAATGACCG-AATTTCTGTA-TGGT--GAATGGACGAAGCCCCGACAAGGCTGCTTATGAAAGAAATGTGCAATATCCTCAAAGGGAAAATTTCAAACA--TGCCCAGAGGGCAATGATGGATCTCAAGTAAGAGAAAGTCGAAACCCAGG-AACGCTGAGATTGAAG-GC-C-TTTTCCTGGCACCGTTACCGCACTCATTCTGGCCCCGATCAGTTGCACATAAATCCTGCCTGCCTGCTTGCTTATCTGGGCTTGCA-TGAGGA-TGGGCATGACTTTGAAAGGGAAGGGTACTCACTGGTCTATGGATAGACCCATTCAAATTACTCCCAAACAGCCAAGGGTGTATTGCTGCTGATGAGACCAAACTTGAAAACCCAG--CACAAGAGTCA-ATTG-GTCACTGAAATGCCACTCTGCTGCATTTGAAGGATTTAAGAGTATCAAG-TTCAT-A-A-GGAAGAAAGTGAT-TCCAC-AGGAAAGCT-TCCACAAGA-GGGTCGAGATTGCTTCAAATGAGAATGTGGAGAATCAT-GGATCCCAATACCCTGGAACTGAGAAGCAGATACTGGTACATAAGGACCAGGAGTGGAGGAAATACCAATCAACAAAAAGGCATCCGCA-GCCAGATCGGCCCTGTGCAGCCTACATTGTCTCATATTCCAATAG-GATCT-CC-GTCGAAAAGAGCAACCGTTATGGCAGCATTCAGCGGGAACAATGAAGGACGGACATCCGACATGCGAACAGGACGGTTATAAGAAATGAT-GAAAGTGCAAAGCCAGAAGATTTGTCCTTCC--GGGCGGGGAGTGCTGTCGCGCTCAACCGA-AAGGCAACGAACCCGATCGTCTGCGCCTTTTTCGTGTATAATGAGTAATGAAGGGTCTTA-TTCTGTC--A-ACAATGCAGAGGAGTA");
	
	char* b = strdup("TTATGAACAGAT-G-GA-AACTGATGGGGATCGCCAGAATGCAACTGAGATTAGGGCATCCGTCGGGAAG-ATGATTGATGGAATTGGGAGATTCTACATCCAAATGTGCACTGAACTTAAACTCAGTGATCATGAAGGGCGATTGATCCAGAACAGCTTGACAATAGA-GAAAATGGTGCTCTCTGCTTTTGATGAAAGAAGGAATAAATACCTGGAAGAACACCCCAGCGCGGGGAAAGATCCCAAGAAAACTGGGGGACCCATATACAGGAGAGTAGATGGAAA-ATG-GATGAGGGAACTCGTTCTTTATGACAAAGAAGAAATAAGGCGAATCTGGCGCCAAGCCAACAATGGTGAGGATGCGACAGCTGGTCTAACTCACATAATGATCTGGC-ATTCCAATTTGAATGATGCAACATACCAGAGGACAAGAGCTCTTGTTCGAACTGGAATGGATCCCAGAATGTGCTCTCTGATGCAGGGCTCGACTCTCCCTAGAAGGTCCGGAGCTGCAGGTGCTGC--AG-TCAAA-GGAATCGGGACAATGGTGATGGA---ACTGATCAGAATGGTCAAACGGGGGATCAACGATCGAAATTTCTGGAGAGGTGAGAATGGGCGGA-AAACAAGAA-G-TGCTTATGAGAG-AATGTGCAACATTCTTAAA-GGAAAATTTCAAACAGCTGCACAAAGAGCAATGGTGGA--TCAAGTGAGAGAAAGTCGGAACCCAGGAAATGCTGAGATCGAAGATCTCATATTTTTGGCA-AGAT-CTGCATTGA-TATTGAGAGGGTCAGTTGCTCACAAATCTTGCCTACCTGCCTGTGCATATGGACCTGCAGTATCCAGTGGGTACGACTTCGAAAAAGAGGGATATTC-CTTG-GTGGGAATAGACCCTTTCAAACTACTTCAAAATAGCCAA--ATATACAGC--CTAATCAGACCTAAC--GAGAATCCAGCACACAAGAGTCAGCTGGTGTGGATGGCATGCCATTCTGCTGCATTTGAA-GATTTAAGATTGTTAAGCTTCATCAGAGGGAAAAAAGT-ATCTCCTCGGGGGAAACTGTCAACTAGAGGGGTACAAATTGCTTCAAATGAGAACATGGATAAT-ATGGGAT-CGAGCACTCTTGAACTGAGAAGCGGGTACTGGGCCATAAGGACCAGGAGTGGAGGAAACACTAATCAAC-AGAGGGCCTCCGCAGGCCA-AAC--CAGTGTGCAACCTACGTT-T-TC-T-GTACAA-AGAAACCTCCCATTTGAAAAG-TCAACCATCATGGCAGCATTCACTGGAAATACGGAGGGAAGAACTTCAGACATGAGGGCA-GA-AATCATAAG-AATGATGGAAGGTGCAAAACCAGAAGAAGTGTCATTCCGGGGGAGGGGAGTTTTCGAGCTCTCAGACGAGAAGGCAGCGAACCCGATCG--TGC-CCTCTTT--TG-AT-ATGAGTAATGAAGGATCTTATTTCT-TCGGAGACAATGCAGAAGAGTA");
	
	int len = strlen(a);
	
	int c = 0;
	int c2 = 0;
	for (int i =0; i < len; i++) {
		if (b[i] == '-') {
			c2++;
		}
		else if (a[i] != b[i]) {
			c++;
		}
	}
	
	printf("Distance: %i\n",c);
	printf("Insertions: %i\n",c2);
	printf("lenA: %i\n",strlen(a));
	printf("lenB: %i\n",strlen(b));
}