Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,167,839	C→G	13.8%	L20V (CTC→GTC)	eda →	KHG/KDPG aldolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,167,839	0	C	G	13.8%	81.6 / 7.3	29	L20V (CTC→GTC)	eda	KHG/KDPG aldolase
Reads supporting (aligned to +/- strand):  ref base C (17/8);  new base G (2/2);  total (19/10)
Fisher's exact test for biased strand distribution p-value = 5.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.38e-01

CGCCGCGCTGGCGGGCAACGACCTTACTGAAATGCCGATTCGCGCTTTTCTTCACGACCCCCTGGACATCTACTGGTGCCCATGAGCCAAGGATCCACAGTCATGACCACCCTTGAACGCCCACAGCCCAAGCTCTCGATGGCTGATAAAGCCGCGCGGATCGACGCCATTTGCAAAAAGGCGCGCATCCTGCCGGTAATCACCATTGCTCGTGAGGAAGACATCCTGCCGCTGGCTGACGCCCTGGCCGCTGGCGGTATCCGT  >  NC_002947/1167707‑1167970                                                                                                                                     |                                                                                                                                    cgccgcgcTGGCGGGCAACGACCTTACTGAAATGCCGATTCGCGCTTTTCTTCACGACCCCCTGGACATCTACTGGTGCCCATGAGCCGAGGATCCACAGTCATGACCACCCTTGAACGCCCACAGCCCAAGctct                                                                                                                                  >  1:437399/1‑136 (MQ=255) cgccgcgcTGGCGGGCAACGACCTTACTGAAATGCCGATTCGCGCTTTTCTTCACGACCCCCTGGACATCTACTGGTGCCCATGAGCCAAGGATCCACAGTCATGACCACCCTTGAACGCCCACAGCCCAAGctct                                                                                                                                  <  1:507834/136‑1 (MQ=255)  gccgcgcTGGCGGGCAACGACCTTACTGAAATGCCGATTCGCGCTTTTCTTCACGACCCCCTGGACATCTACTGGTGCCCATGAGCCAAGGATCCACAGTCATGACCACCCTTGAACGCCCACAGCCCAAGGtctc                                                                                                                                 >  2:507093/1‑136 (MQ=255)  gccgcgcTGGCGGGCAACGACCTTACTGAAATGCCGATTCGCGCTTTTCTTCACGACCCCCTGGACATCTACTGGTGCCCATGAGCCAAGGATCCACAGTCAGGACCACCCTTGAACGCCCACAGCCCAAGGtctc                                                                                                                                 >  2:504387/1‑136 (MQ=255)   ccgcgcTGGCGGGCAACGACCTTACTGAAATGCCGATTCGCGCTTTTCTTCACGACCCCCTGGACATCTACTGGTG‑CCATGAGCCAAGGATCCACAGTCATGACCACCCTTGAACGCCCACAGCCCAAGCTCTCGa                                                                                                                               <  1:182572/136‑1 (MQ=255)                    gACCTTACTGAAATGCCGATTCGCGCTTTTCTTCACGACCCCCTGGACATCTACTGGTGCCCATGAGCCAAGGATCCACAGTCATGACCACCCTTGAACGCCCACAGCCCAAGCTCTCGATGGCTGATAAAGCcgc                                                                                                               >  2:202567/1‑136 (MQ=255)                                   ccGATTCGCGCTTTTCTTCACGACCCCCTGGACATCTACTGGTGCCCATGAGCCAAGGATCCACAGTCATGACCACCCTTGAACGCCCACAGCCCAAGCTCTCGATGGCTGATAAAGCCGCGCGGATCGACGCCAt                                                                                                >  2:445717/1‑136 (MQ=255)                                        tCGCGCTTTTCTTCACGACCCCCTGGACATCTACTGGTGCCCATGAGCCAAGGATCCACAGTCATGACCACCCTTGAACGCCCACAGCCCAAGCTCTCGATGGCTGATAAAGCCGCGCGGATCGACGCCATTTGCa                                                                                           >  2:256875/1‑136 (MQ=255)                                                            cccTGGACATCTACTGGTGCCCATGAGCCAAGGATCCACAGTCATGACCACCCTTGAACGCCCACAGCCCAAGCTCTCGATGGCTGATAAAGCCGCGCGGATCGACGCCATTTGCAAAAAGGCGCGCATCCTGCCg                                                                       <  2:421755/136‑1 (MQ=255)                                                                       tACTGGTGCCCATGAGCCAAGGATCCACAGTCATGACCACCCTTGAACGCCCACAGCCCAAGCTCTCGATGGCTGATAAAGCCGCGCGGATCGACGCCATTTGCAAAAAGGCGCGCATCCTGCCGGTAATCACCAt                                                            <  1:259362/136‑1 (MQ=255)                                                                       gaCTGGTGCCCATGAGCCAAGGATCCACAGTCATGACCACCCTTGAACGCCCACAGCCCAAGCTCTCGATGGCTGATAAAGCCGCGCGGATCGACGCCATTTGCAAAAAGGCGCGCATCCTGCCGGTAATCACCAt                                                            <  1:227965/135‑1 (MQ=255)                                                                         cTGGTGCCCATGAGCCAAGGATCCACAGTCATGACCACCCTTGAACGCCCACAGCCCAAGCTCTCGATGGCTGATAAAGCCGCGCGGATCGACGCCATTTGCAAAAAGGCGCGCATCCTGCCGGTAATCACCATTg                                                          >  1:204138/1‑136 (MQ=255)                                                                           ggTGCCCATGAGCCAAGGATCCACAGTCATGACCACCCTTGAACGCCCACAGCCCAAGCTCTCGATGGCTGATAAAGCCGCGCGGATCGACGCCATTTGCAAAAAGGCGCGCATCCTGCCGGTAATCACCATTGCt                                                        >  1:5362/1‑136 (MQ=255)                                                                              gCCCATGAGCCAAGGATCCACATTCATGACCACCCTTGAACGCCCACAGCCCAAGCTCTCGATGGCTGATAAAGCCGCGCGGATCGACGCCATTTTCAAAAAGGCGCGCATCCTGCCGGTAATACCCATTGCTCGt                                                     >  2:211909/1‑136 (MQ=255)                                                                              gCCCATGAGCCAAGGATCCACAGTCATGACCACCCTTGAACGCCCACAGCCCAAGCTCTCGATGGCTGATAAAGCCGCGCGGATCGACGCCATTTGCAAAAAGGCGCGCATCCTGCCGGTAATCACCATTGCTCGt                                                     >  2:388880/1‑136 (MQ=255)                                                                               cccATGAGCCAAGGATCCACAGTCATGACCACCCTTGAACGCCCACAGCCCAAGCTCTCGATGGCTGATAAAGCCGCGCGGATCGACGCCATTTGCAAAAAGGCGCGCATCCTGCCGGTAATCACCATTGCTCgcg                                                    >  2:480909/1‑134 (MQ=255)                                                                               cccATGAGCCAAGGATCCACAGTCATGACCACCCTTGAACGCCCACAGCCCAAGCTCTCGATGGCTGATAAAGCCGCGCGGATCGACGCCATTTGCAAAAAGGCGCGCATCCTGCCGGTAATCACCATTGCTCGTg                                                    >  1:388723/1‑136 (MQ=255)                                                                                  aTGAGCCAAGGATCCACAGTCATGACCACCCTTGAACGCCCACAGCCCAAGCTCTCGATGGCTGATAAAGCCGCGCGGATCGACGCCATTTGCAAAAAGGCGCGCATCCTGCCGGTAATCACCATTGCTCGTGAgg                                                 >  1:82453/1‑136 (MQ=255)                                                                                   tGAGCCAAGGATCCACAGTCATGACCACCCTTGAACGCCCACAGCCCAAGCTCTCGATGGCTGATAAAGCCGCGCGGATCGACGCCATTTGCAAAAAGGCGCGCATCCTGCCGGTAATCACCATTGCTCGTGAGGa                                                >  1:385415/1‑136 (MQ=255)                                                                                                   aGTCATGACCACCCTTGAACGCCCACAGCCCAAGCTCTCGATGGCTGATAAAGCCGCGCGGATCGACGCCATTTGCAAAAAGGCGCGCATCCTGCCGGTAATCACCATTGCTCGTGAGGAAGACATCCTGCCgctg                                >  2:350466/1‑136 (MQ=255)                                                                                                         gACCACCCTTGAACGCCCACAGCCCAAGCTCTCGATGGCTGATAAAGCCGCGCGGATCGACGCCATTTGCAAAAAGGCGCGCATCCTGCCGGTAATCACCATTGCTCGTGAGGAAGACATCCTGCCGCTGGCTGAc                          <  2:190826/136‑1 (MQ=255)                                                                                                          accaccCTTGAACGCCCACAGCCCAAGCTCTCGATGGCTGATAAAGCCGCGCGGATCGACGCCATTTGCAAAAAGGCGCGCATCCTGCCGGTAATCACCATTGCTCGTGAGGAAGACATCCTGCCGCTGGCTGACg                         >  1:560337/1‑136 (MQ=255)                                                                                                                cTTGAACGCCCACAGCCCAAGCTCTCGATGGCTGATAAAGCCGCGCGGATCGACGCCATTTGCAAAAAGGCGCGCATCCTGCCGGTAATCACCATTGCTCGTGAGGAAGACATCCTGCCGCTGGCTGACGCCCTgg                   >  2:377969/1‑136 (MQ=255)                                                                                                                  tGAACGCCCACAGCCCAAGCTCTCGATGGCTGATAAAGCCGCGCGGATCGACGCCATTTGCAAAAAGGCGCGCATCCTGCCGGTAATCACCATTGCTCGTGAGGAAGACATCCTGCCGCTGGCTGACGCCCTGGcc                 >  1:309151/1‑136 (MQ=255)                                                                                                                    aaCGCCCACAGCCCAAGCTCTCGATGGCTGATAAAGCCGCGCGGATCGACGCCATTTGCAAAAAGGCGCGCATCCTGCCGGTAATCACCATTGCTCGTGAGGAAGACATCCTGCCGCTGGCTGACGCCCTGGCCGc               >  1:376349/1‑136 (MQ=255)                                                                                                                       gCCCACAGCCCAAGGTCTCGATGGCTGATAAAGCCGCGCGGATCGACGCCATTTGCAAAAAGGCGCGCATCCTGCCGGTAATCACCATTGCTCGTGAGGAAGACATCCTGCCGCTGGCTGACGCCCTGGCCGCTgg            <  1:504387/136‑1 (MQ=255)                                                                                                                       gCCCACAGCCCAAGGTCTCGATGGCTGATAAAGCCGCGCGGATCGACGCCATTTGCAAAAAGGCGCGCATCCTGCCGGTAATCACCATTGCTCGTGAGGAAGACATCCTGCCGCTGGCTGACGCCCTGGCCGCTgg            <  1:507093/136‑1 (MQ=255)                                                                                                                       gCCCACAGCCCAAGCTCTCGATGGCTGATAAAGCCGCGCGGATCGACGCCATTTGCAAAAAGGCGCGCATCCTGCCGGTAATCACCATTGCTCGTGAGGAAGACATCCTGCCGCTGGCTGACGCCCTGGCCGCTgg            <  1:350466/136‑1 (MQ=255)                                                                                                                            caGCCCAAGCTCTCGATGGCTGATAAAGCCGCGCGGATCGACGCCATTTGCAAAAAGGCGCGCATCCTGCCGGTAATCACCATTGCTCGTGAGGAAGACATCCTGCCGCTGGCTGACGCCCTGGCCGCTGGCGGTa       >  2:519058/1‑136 (MQ=255)                                                                                                                             aGCCCAAGCTCTCGATGGCTGATAAAGCCGCGCGGATCGACGCCATTTGCAAAAAGGCGCGCATCCTGCCGGTAATCACCATTGCTCGTGAGGAAGACATCCTGCCGCTGGCTGACGCCCTGGCCGCTGGCGGTAt      <  1:213660/136‑1 (MQ=255)                                                                                                                                 cAAGCTCTCGATGGCTGATAAAGCCGCGCGGATCGACGCCATTTGCAAAAAGGCGCGCATCCTGCCGGTAATCACCATTGCTCGTGAGGAAGACATCCTGCCGCTGGCTGACGCCCTGGCCGCTGGCGGTATCCGt  <  1:256875/136‑1 (MQ=255)                                                                                                                                     |                                                                                                                                    CGCCGCGCTGGCGGGCAACGACCTTACTGAAATGCCGATTCGCGCTTTTCTTCACGACCCCCTGGACATCTACTGGTGCCCATGAGCCAAGGATCCACAGTCATGACCACCCTTGAACGCCCACAGCCCAAGCTCTCGATGGCTGATAAAGCCGCGCGGATCGACGCCATTTGCAAAAAGGCGCGCATCCTGCCGGTAATCACCATTGCTCGTGAGGAAGACATCCTGCCGCTGGCTGACGCCCTGGCCGCTGGCGGTATCCGT  >  NC_002947/1167707‑1167970

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑