Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB143	1,639,779	T→G	28.5%	N212H (AAC→CAC)	gcd ←	Quinoprotein glucose dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB143	1,639,779	0	T	G	28.5%	52.7 / 5.2	36	N212H (AAC→CAC)	gcd	Quinoprotein glucose dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (6/19);  new base G (9/2);  total (15/21)
Fisher's exact test for biased strand distribution p-value = 2.44e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGGTCGTAGCGCCAGATTTCTGCGCCAGTGTCCGGGTCCAGCGCCAGCAGGCGACTGTGCGCGGTGCAGGCGTAAAGCATGCCGTTGACCTTCAGCGGGGTGTTCTGGTTGGTCAGCTCCACCGGGTCGTTTTCGGTTGGCAGATCACCGGTGCGAATGCGCCAGGCTTCTTCCAGACGGTAGGCGTTCTGCGGGGTGATCTGGCGCAGCGGCGAGTAGCGGTCGCCATGCTCGGTACGGCCGTAGGCCTGCCATTCGC  >  ACB143/1639647‑1639907                                                                                                                                     |                                                                                                                                 gCGGGTCGTAGCGCCAGATTTCTGCGCCAGTGTCCGGGGCCAGCGCCAGCAGGCGACTGTGCGCGGTGCAGGCGTAAAGCATGCCGTTGACCTTCAGCGGGGTGTTCTGGTTGGTCAGCTCCACCGGGTCGTTTTc                                                                                                                               <  6:168242/136‑1 (MQ=255)  cGGGTCGTAGCGCCAGATTTCTGCGCCAGTGTCCGGGTCCAGCGCCAGCAGGCGACTGTGCGCGGTGCAGGCGTAAAGCATGCCGTTGACCTTCAGCGGGGTGTTCTGGTTGGTCAGCTCCACCGGGTCGTTTTc                                                                                                                               <  2:174146/135‑1 (MQ=255)     gTCGTAGCGCCAGATTTCTGCGCCAGTGTCCGGGTCCAGCGCCAGCAGGCGACTGTGCGCGGTGCAGGCGTAAAGCATGCCGTTGACCTTCAGCGGGGTGTTCTGGTTGGTCAGCTCCACCGGGTCGTGTTggggt                                                                                                                           >  8:9056/1‑134 (MQ=255)     gTCGTAGCGCCAGATTTCTGCGCCAGTGTCCGGGTCCAGCGCCAGCAGGCGACTGTGCGCGGTGCAGGCGTAAAGCATGCCGTTGACCTTCAGCGGGGTGTTCTGGTTGGTCAGCTCCACCGGGTCGTGTTCGGtt                                                                                                                           >  5:69598/1‑136 (MQ=255)            cgcCAGATTTCTGCGCCAGTGTCCGGGTCCAGCGCCAGCAGGCGACTGTGCGCGGTGCAGGCGTAAAGCATGCCGTTGACCTTCAGCGGGGTGTTCTGGTTGGTCAGCTCCACCGGGTCGTTTTCGGTTGGTAGAt                                                                                                                    <  3:256567/136‑1 (MQ=255)                      cTGCGCCAGTGTCCGGGTCCAGCGCCAGCAGGCGACTGTGCGCGGTGCAGGCGTAAAGCATGCCGTTGACCTTCAGCGGGGTGTTCTGGTTGGTCAGCTCCACCGGGTCGTTTTCGGTTGGCAGATCACCGGTGc                                                                                                           <  1:2096/135‑1 (MQ=255)                          gcCAGTGTCCGGGTCCAGCGCCAGCAGGCGACTGTGCGCGGTGCAGGCGTAAAGCATGCCGTTGACCTTCAGCGGGGTGTTCTGGTTGGTCAGCTCCACCGGGTCGTTTTCGGTTGGCAGATCACCGGTGCGAATg                                                                                                      <  7:27607/136‑1 (MQ=255)                            cAGTGTCCGGGTCCAGCGCCAGCAGGCGACTGTGCGCGGTGCAGGCGTAAAGCATGCCGTTGACCTTCAGCGGGGTGTTCTGGTTGGTCAGCTCCACCGGGTCGTTTTGGGGTGGGCGATCCCCCGTGCGAATgcg                                                                                                    >  4:40509/1‑136 (MQ=255)                                 tCCGGGTCCAGCGCCAGCAGGCGACTGTGCGCGGTGCAGGCGTAAAGCATGCCGTTGACCTTCAGCGGGGTGTTCTGGTTGGTCAGCTCCACCGGGTCGTTTTCGGTTGGCAGATCACCGGTGCGAATGCGCCAgg                                                                                               <  7:30640/136‑1 (MQ=255)                                     ggTCCAGCCCCAGCAGGCGACTGTGCGCGGTGCAGGCGTAAAGCATGCCGTTGACCTTCAGCGGGGTGTTCTGGTTGGTCAGCTCCACCGGGTCGTTTTCGGTTGGCAGATCACCGGTGCGAATGCGCCAGGcttc                                                                                           <  5:120041/136‑1 (MQ=255)                                                 gcagGCGACTGTGCGCGGTGCAGGCGTAAAGCATGCCGTTGACCTTCAGCGGGGTGTGCTTGTTGGTCAGCTCCACAGGGTCGTTTTCGGTTGGCAGATCACCGGTGCGAATGCGCCAGGCTTCTTCCAGACGGTa                                                                               <  3:232683/136‑1 (MQ=255)                                                               gcgGTGCAGGCGTAAAGCATGCCGTTGACCTTCAGCGGGGTGTTCTGGTTGGTCAGCTCCACCGGGTCGTGTTGGGTTGGCAGATCACCGGTGCGCATGCGCCACGCTTCTTACAGACGGTTGGCGTTGTGCgggg                                                                 >  6:240188/1‑136 (MQ=255)                                                                              aGCATGCCGTTGACCTTCAGCGGGGTGTTCTGGTTGGTCAGCTCCACCGGGTCGTTTTCGGTTGGCAGATAACCGGTGCGAATGCGCCAGGCTTCTTCCATACGGGAGGCGGTCTGCGGGGTGATCTGGCGCAgcg                                                  >  7:191010/1‑136 (MQ=255)                                                                                cATGCCGTTGACCTTCAGCGGGGTGTTCTGGTTGGTCAGCTCCACCGGGTCGTTTTGGGTTGGCAGATCACCGGTGCGAATTCGCCACGCGTATTCCAGCCGGTAGGGGTTCTGCGGGGGGATCAGGCGCAgcggc                                                >  5:118911/1‑136 (MQ=255)                                                                                cATGCCGTTGACCTTCAGCGGGGTGTTCTGGTTGGTCAGCTCCACCGGGTCGTTTTCGGTTTGCATAACAAAGGTGCCAAAGCGCCATGCATCTTCCACACCGTATGCCTTCTGCGGGGGGATATTGCGCAccggc                                                >  2:34213/1‑136 (MQ=255)                                                                                cATGCCGTTGACCTTCAGCGGGGTGTTCTGGTTGGTCAGCTCCACCGGGTCGTTTTCGGTTGGCAGATCACCCGTGCGAATTCGCCAGGCTTGTTACAGACAGTATGCGTTCTGCGGGGTGATCTGGCGCAgcggc                                                >  2:98325/1‑136 (MQ=255)                                                                                cATGCCGTTGACCTTCAGCGGGGTGTTCTGGTTGGTCAGCTCCACCGGGTCGTGTTGGGTTGGGAGATAACCGGTGCGAATGCGCCCGGCGTCTTCCAGCCGGTAGGCGTTGTGCGGGGTGATATGGCGCAgcggg                                                >  6:212146/1‑135 (MQ=255)                                                                                cATGCCGTTGACCTTCAGCGGGGTGTTCTGGTTGGTCAGCTCCACCGGGTCGTGTTGGGGTGGCAGATCACCCGTGCGAATTCGCCCGGCTTCTTGCAGACGGGAGGGGGTGGGCGGGGGGATCTGGCGCCccggc                                                >  6:251304/1‑136 (MQ=255)                                                                                     cGTTGACCTTCAGCGGGGTGTTCTGGTTGGTCAGCTCCACCGGGTCGTGTTGGGTTGGCAGATGACCCGTGCGAATTCGCCAGGCTTGTTAAAGACCGTAGGGGTTCTGCGGGGTGTTATTGCGCAGCGGCGAGtt                                           >  2:74402/1‑135 (MQ=255)                                                                                     cGTTGACCTTCAGCGGGGTGTTCTGGCTGGTCAGCTCCACCGGGTCGTTTTCGGTTGGCAGATCACCGGTGCGAATGCGCCAGGCTTCTTCCAGACGGTAGGCGTTCTGCGGGGTGATCTGGCGCAGCGGCGAGTa                                           <  2:203943/136‑1 (MQ=255)                                                                                      gTTGACCTTCAGCGGGGTGTTCTGGTTGGTCAGCTCCACCGGGTCGTGTTGGGTTGGCAGATCACCGGGGCGAATGCGCCACGCGTGTTCCACCCGGTATGCGTTCTGGGGGGTGGTCGGGCGCCGCGGCGCGTAg                                          >  5:223902/1‑136 (MQ=255)                                                                                           ccTTCAGCGGGGTGTTCTGGGTGGTCAGCTCCACCGGGTCGTTTTCGGTTGGCAGATCACCGGTGCGAATGCGCCAGGCTTCTTCCAGACGGTAGGCGTTCTGCGGGGTGATCTGGCGCAGCGGCGAGTAGCGGTc                                     <  8:173447/136‑1 (MQ=255)                                                                                             ttCAGCGGGGTGTTCTGGTTGGTCAGCTCCACCGGGTCGTTTTGGGTTTGGAGATAACCGGGGGGAATTCGCCCGGCTTCTTACAGCCCGTAGGCGGTTTTCGGGGGGGTCTGGGGCAGCGGCGCGTAGCGGTCGc                                   >  5:81132/1‑136 (MQ=255)                                                                                             ttCAGCGGGGTGTTCTGGTTGGTCAGCTCCACCGGGTCGTGTTAGGTTGGCAGATCACCCGTGCGAATGCGCCACGCTTCTTACAGACGGGAGGCGGTCTGGGGGGTGTTATGGCGCAGCGGCGAGTAGCGGTCGc                                   >  5:69119/1‑136 (MQ=255)                                                                                               ccgCGGGGTGGTCTGGTTGGGCAGCTCCACCGGGTCGGGTTCGGTTGGCAGATCACCGGTGCGCATGCGCCAGGCTTCTTCCAGACGGTAGGCGTTCTGCGGGGTGATCGGGCGCAGCGGCGAGTAGCGGTCGcc                                  <  3:219857/133‑1 (MQ=255)                                                                                               cAGCGGGGTGTTCTGGTTGGTCAGCTCCACCGGGTCGTTTTCGGTTGGCAGATCACCGGTGCGAATGCGCCAGGCTTCTTCCAGACGGTAGGCGTTCTGCGGGGTGATCTGGCGCAGCGGCGAGTAGCGGTCGcc                                  <  5:165079/135‑1 (MQ=255)                                                                                               cAGCGGGGTGTTCTGGTTGGTCAGCTCCACCGGGTCGTTTTCGGTTGGCAGATCACCGGTGCGAATGCGCCAGGCTTCTTCCAGACGGTAGGCGTTCTGCGGGGTGATCTGGCGCAGCGGCGAGTAGCGGTCGcc                                  <  5:154206/135‑1 (MQ=255)                                                                                                aGCGGGGTGTTCTGGTTGGTCAGCTCCACCGGGTCGTTTTCGGTTGGCAGATCACCGGTGCGAATGCGCCAGGCTTCTTCCAGACGGTAGGCGTTCTGCGGGGTGATCTGGCGCAGCGGCGAGTAGCGGTCGCCAt                                <  1:164403/136‑1 (MQ=255)                                                                                                     ggggtTCTGGTTGGTCAGCTCCACCGGGTCGTTTTCGGTTGGCAGATCACCGGTGCGAATGCGCCAGGCTTCTTCCAGACGGTAGGCGTTCTGCGGGGTGATCTGGCGCAGCGGCGAGTAGCGGTCGCCATGCTCg                           <  6:236198/133‑1 (MQ=255)                                                                                                     ggTGTTCTGGTTGGTCAGCTCCACCGGGTCGTTTTCGGTTGGCAGATCACCGGTGCGAATGCGCCAGGCTTCTTCCAGACGGTAGGCGTTCTGCGGGGTGATCTGGCGCAGCGGCGAGTAGCGGTCGCCATGCTCg                           <  6:175111/136‑1 (MQ=255)                                                                                                     ggTGTTCTGGTTGGTCAGCTCCACCGGGTCGTTTTCGGTTGGCAGATCACCGGTGCGAATGCGCCAGGCTTCTTCCAGACGGTAGGCGTTCTGCGGGGTGATCTGGCGCAGCGGCGAGTAGCGGTCGCCATGCTCg                           <  5:171162/136‑1 (MQ=255)                                                                                                     ggTGTTCTGGTTGGTCAGCTCCACCGGGTCGTTTTCGGTTGGCAGATCACCGGTGCGAATGCGCCAGGCTTCTTCCAGACGGTAGGCGTTCTGCGGGGTGATCTGGCGCAGCGGCGAGTAGCGGTCGCCATGCTCg                           <  4:68754/136‑1 (MQ=255)                                                                                                     ggTGTTCTGGGTGGTCAGCTCCACCGGGTCGTTATCGGTGGGCAGAGCCCCGGTGCGAATGCGCCAGGCTTCTTCCAGACGGTAGGCGTTCTGCGGGGTGATCTGGCGCAGCGGCGAGTAGCGGTCGCCATGCTCg                           <  2:127276/136‑1 (MQ=255)                                                                                                         ttCTGGTTGGTCAGCTCCACCGGGTCGTGTTGGGGTGGCAGATAACCGGTGCGCATGCGCCAGGCTTGTTGCAGACGGTAGGCGTTGTGCGGGGTGTTCTTGCGCAGCGGCGAGTAGCGGTCGCCATGCTCCGTAc                       >  2:38862/1‑136 (MQ=255)                                                                                                                tggtCAGCGCCACCGGGTCGTGTTCGGTTGGCAGATCACCGGTGCGAATGCGCCAGGCTTCTTCCAGACGGTAGGCGTTCTGCGGGGTGATCTGGCGCAGCGGCGAGTAGCGGTCGCCATGCTCGGTACGGCCGTa                <  6:155363/136‑1 (MQ=255)                                                                                                                              gggTCGTTTTCGGTTGGCAGATCACCGGTGCGAATGCGCCAGGCTTCTTCCAGACGGTAGGCGTTCTGCGGGGTGATCTGGCGCAGCGGCGAGTAGCGGTCGCCATGCTCGGTACGGCCGTAGGCCTGCCATTCGc  <  5:210468/136‑1 (MQ=255)                                                                                                                                     |                                                                                                                                 GCGGGTCGTAGCGCCAGATTTCTGCGCCAGTGTCCGGGTCCAGCGCCAGCAGGCGACTGTGCGCGGTGCAGGCGTAAAGCATGCCGTTGACCTTCAGCGGGGTGTTCTGGTTGGTCAGCTCCACCGGGTCGTTTTCGGTTGGCAGATCACCGGTGCGAATGCGCCAGGCTTCTTCCAGACGGTAGGCGTTCTGCGGGGTGATCTGGCGCAGCGGCGAGTAGCGGTCGCCATGCTCGGTACGGCCGTAGGCCTGCCATTCGC  >  ACB143/1639647‑1639907

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