Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	3,983,006	0	A	C	34.5%	94.9 / 10.0	59	T627P (ACC→CCC)	apc3	Acetophenone carboxylase gamma subunit
Reads supporting (aligned to +/- strand):  ref base A (6/32);  new base C (20/0);  total (26/33)
Fisher's exact test for biased strand distribution p-value = 1.25e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACCCTGTACTCGTACAGCGAACCCGCCAGCCCGGTGGAGCTGGTGAACCTGGAGTGCTCGGTGATCGGCCACCTGCCACGCCCGCCGCAACCCGAGCTGCAGGGGCCGGCCGAGCCCCCCGCCCCCACCGCGCAATCCGTCCGCCCGATGCTGTTCAGTGCCGACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACAC  >  ACB122/3982880‑3983133                                                                                                                               |                                                                                                                                aCCCTGTACTCGTACAGCGAACCCGCCAGCCCGGTGGAGCTGGTGAACCTTGAGTGCTCGGTGATCGGCCACCTGCCACGCCCGCCGCAACCCGAGCTGCAGGGGCCGGCCGAGCCCCCCGCCCCCACCGCGCAAtc                                                                                                                       <  8:204723/137‑1 (MQ=255) aCCCTGTACTCGTACAGCGAACCCGCCAGCCCGGTGGAGCTGGTGAACCTGGAGTGCTCGGTGATCGGCCACCTGCCACGCCCGCCGCAACCCGAGCTGCAGGGGCCGGCCGAGCCCCCCGCCCCCACCGCGCAAtc                                                                                                                       <  7:236276/137‑1 (MQ=255)      gTACTCGTACAGCGAACCCGCCAGCCCGGTGGAGCTGGTGAACCTGGAGTGCTCGGTGATCGGCCACCTGCCACGCCCGCCGCAACCCGAGCTGCAGGGGCCGGCCGAGCCCCCCGCCCCCACCGCGCAAtccgtcc                                                                                                                  <  8:120092/137‑1 (MQ=255)        aCTCGTACAGCGAACCCGCCAGCCCGGTGGAGCTGGTGAACCTGGAGTGCTCGGTGATCGGCCACCTGCCATGCCCGCCGCAACCCGAGCTGCAGGGGCCGGCCGAGCCCCCCGCCCCCCCCCCGCAAACCGGCCGc                                                                                                                >  3:205079/1‑137 (MQ=255)        aCTCGTACAGCGAACCCGCCAGCCCGGTGGAGCTGGTGAACCTGGAGTGCTCGGTGATCGGCCACCTGCCACGCCCGCCGCAACCCGAGCTGCAGGGGCCGGCCGCGCCCCCCGCCCCCCCCCCGCCATCCGTCCGc                                                                                                                >  4:285369/1‑137 (MQ=255)        aCTCGTACAGCGAACCCGCCAGCCCGGTGGAGCTGGTGAACCTGGAGTGCTCGGTGATCGGCCACCTGCCACGCCCGCCGCAACCCGAGCTGCAGGGGCCGGCCGCGCCCCCCCCCCCCCCCGCCCAAaacgtcccc                                                                                                                >  4:141849/1‑135 (MQ=255)        aCTCGTACAGCGAACCCGCCAGCCCGGTGGAGCTGGTGAACCTGGAGTGCTCGGTGATCGGCCACCTGCCACGCCCGCCGCAACCCGAGCTGCAGGGGCCGGCCGCGCCCCCCCCCCCCACCCCCCCAACCGTCCGc                                                                                                                >  2:96978/1‑137 (MQ=255)          tCGTACAGCGAACCCGCCAGCCCGGTGGAGCTGGTGAACCTGGAGTGCTCGGTGATCGGCCACCTGCCACGCCCGCCGCAACCCGAGCTGCAGGGGCCGGCCGAGCCCCCCGCCCCCCCCGCGCAATCCGTCCGccc                                                                                                              >  7:89208/1‑137 (MQ=255)          tCGTACAGCGAACCCGCCAGCCCGGTGGAGCTGGTGAACCTGGAGTGCTCGGTGATCGGCCACCTGCCACGCCCGCCGCAACCCGAGCTGCAGGGGCCGGCCGAGCCCCCCGCCCCCCCCCCGCAAACCGGCCGccc                                                                                                              >  3:198542/1‑137 (MQ=255)           cGTACAGCGACCCCGCCAGCCCGGTGGAGCTGCTGAACCTGCAGTGATCGGTGATCCGCCACCTGCCACGCCCCCCGCAACCCGAGCAGCACGCGCCGGCCCAGCCCCCCCCCCAACCCGCCCAATCCGACCGccac                                                                                                             >  2:263714/1‑135 (MQ=255)           cGTACAGCGAACCCGCCAGCCCGGTGGAGCTGGTGAACCTGGAGTGCTCGGTGATCGGCCACCTGCCACGCCCGCCGCAACCCGAGCTGCAGGGGCCGGCCGCGCCCCCCGCCCCCCCCGCGCAATCCGGCCGCCCg                                                                                                             >  8:157358/1‑137 (MQ=255)           cGTACAGCGAACCCGCCAGCCCGGTGGAGCTGGTGAACCTGGAGTGCTCGGTGATCGGCCACCTGCCACGCCCGCCGCAACCCGAGCTGCAGGGGCCGGCCGAGCCCCCCGCCCCCACCGCGCAATCCGTCCGCCCg                                                                                                             <  4:26852/137‑1 (MQ=255)           cGTACAGCGAACCCGCCAGCCCGGTGGAGCTGGTGAACCTGGAGTGCGCGGTGATCGGCCACCTGCCACGCCCGCCGCAACCCGAGCTGCAGGGGCCGGCCGAGCCCCCCGCCCCCACCGCGCAATCCGTCCGCCCg                                                                                                             <  7:171726/137‑1 (MQ=255)           cGTACAGCGAACCCGCCAGCCCGGTGGAGCGGGTGAACCTGGAGTGCTCGGTGATCGGCCACCTGCCACGCCCGCCGCAACCCTAGCTGCAGGGGCCGGCCGAGCCCCCCGCCCCCACCTCGCAATCCGTCCGCCCg                                                                                                             <  6:308721/137‑1 (MQ=255)               cAGCGAACCCGCCAGCCCGGTGGAGCTGGTGAACCTGGAGTGCTCGGTGATCGGCCACCTGCGACGCCCGCCGCAACCCGAGCTGCAGGGGCCGGCCGAGGCCCCCGCCCCCCCCGCGCAATCCGACCGCCCgcggc                                                                                                         >  1:303821/1‑133 (MQ=255)               cAGCGAACCCGCCAGCCCGGTGGAGCTGGTGAACCTGGAGTGCTCGGTGATCGGCCACCTGCCACGCCCGCCGCAACCCGAGCTGCAGGGGCCGGCCGAGCCCCCCGCCCCCCCCGCCCAATCCGGGCGcccccagg                                                                                                         >  8:18651/1‑132 (MQ=255)               cAGCGAACCCGCCAGCCCGGTGGAGCTGGTGAACCTGGAGTGCTCGGTGATCGGCCACCTGCCACGCCCGCCGCAACCCGAGCTGCAGGGGCCGGCCGAGCCCCCCGCCCCCACCGCGCAATCCGTCCGCCCGATGc                                                                                                         >  4:204491/1‑137 (MQ=255)               cAGCGAACCCGCCAGCCCGGTGGAGCTGGTGAACCTGGAGTGCTCGGTGATCGGCCACCTGCCACGCCCGCCGCAACCCGAGCTGCAGGGGCCGGCCGAGCCCCCCCCCCCCCCCCCGCCATCCGGCCGcccccagg                                                                                                         >  4:320753/1‑132 (MQ=255)                aGCGAACCCGCCAGCCCGGTGGAGCTGGTGAACCTGGAGTGCTCGGTGATCGGCCACCTGCCACGCCCGCCGCAACCCGAGCTGCAGGGGCCGGCCGAGCCCCCCGCCCCCCCCGCGCCAACCGTCCGCCCGATGcg                                                                                                        >  4:136592/1‑136 (MQ=255)                aGCGAACCCGCCAGCCCGGTGGAGCTGGTGAACCTGGAGAGCTCGGTGATCGGCCACCTGCCACGCCCGCCGCCACCCGAGCTGCAGGGGCCGGCCGAGCCCCCCGCCCCCCCCGCGCCATCCGTCCGCCCgctgct                                                                                                        >  5:232969/1‑137 (MQ=255)                 gCGAACCCGCCAGCCCGGTGGAGCTGGTGAACCTGGAGTGCTCGGTGATCGGCCACCTGCCACGCCCGCCGCAACCCGAGCTGCAGGGGCCGGCCGAGCCCCCCGCCCCCACCGCGCAATCCGTCCGCCCGATGCTg                                                                                                       <  2:229187/137‑1 (MQ=255)                 cctaaCCCGGCAGCCCGGTGGAGCTGTTTAACCTGGATTTCTCGGGGATCGGCAAGCGGCCACGCCCGGCGCAACCGGAGGTTCAGTGGCCGTCCTAGCCCCCGGCCCCCACCGCGCAATCCTTCCGCCCGATGCTg                                                                                                       <  6:186698/134‑1 (MQ=255)                    aaCCCGCCAGCCCGGTGGAGCTGGTGAACCTGGAGTGCTCGGTGATCGGCCACCTGCCACGCCCGCCGCAACCCGAGCTGCAGGGGCCGGCCGCGCCCCCCGCCCCCCCCGCGCCATCCGTCCGCCCGATGCTGttg                                                                                                    >  7:86418/1‑136 (MQ=255)                    aaCCCGCCAGCCCGGTGGAGCTGGTGAACCTGGAGTGCTCGGTGATCGGCCACCTGCCACGCCCGCCGCAACCCGAGCTGCAGGGGCCGGCCGCGCCCCCCGCCCCCCCCGCGCCATCCGTCCGCCCGATGCGGttg                                                                                                    >  1:151271/1‑136 (MQ=255)                    aaCCCGCCAGCCCGGTGGAGCTGGTGAACCTGGAGTGCTCGGTGATCGGCCACCTGCCACGCCCGCCGCAACCCGAGCTGCAGGGGCCGGCCGAGCCCCCCGCCCCCACCGCGCAATCCGTCCGCCCGATGCTGtta                                                                                                    >  8:98221/1‑136 (MQ=255)                    aaCCCGCCAGCCCGGTGGAGCTGGTGAACCTGGAGTGCTCGGTGATCGGCCACCTGCCACGCCCGCCGCAACCCGAGCTGCAGGGGCCGGCCGAGCCCCCCGCCCCCACCGCGCAATCCGTCCGCCCGAAGCTGtta                                                                                                    >  8:297665/1‑136 (MQ=255)                           cAGCCCGGTGGAGCTGGTGAACCTGGAGTGCTCGGTGATCGGCCACCTGCCACGCCCGCCGCAACCCGAGCTGCAGGGGCCGGCCGAGCCCCCCGCCCCCACCGCGCAATCCGTCCGCCCGATGCTGTTCAGTGCCg                                                                                             <  5:341638/137‑1 (MQ=255)                              cccGTTGGAGCTGGTGACCCTGTAGTGCTCGGTGAGCGGCCACCTGCCACGCCCGCCGCAACCCGAGCTGCAGGGGCCGGCCGAGCCCCCCGCCCCCACCGCGCAATCCGTCCGCCCGATGCTTTTCAGTGCCGACg                                                                                          <  1:263714/137‑1 (MQ=255)                                   tGGAGCTGGTGAACCTGGAGTGCTCGGTGATCGGCCACCTGCCACGCCCGCCGCAACCCGAGCTGCAGGGGCCGGCCGAGCCCCCCGCCCCCACCGCGCAATCCGTCCGCCCGATGCTGTTCAGTGCCGACGGTGaa                                                                                     <  6:171055/137‑1 (MQ=255)                                   tGGAGCTGGTGAACCTGGAGTGCTCGGTGAGCGGCCACCTGCCACGCCCGCCGCAACCCGAGCTGCAGGGGCCGGCCGAGCCCCCCGCCCCCACCGCGCAATCCGTCCGCCCGATGCTGTTCAGTGCCGACGGTGaa                                                                                     <  6:271941/137‑1 (MQ=255)                                   tGGAGCTGGTGAACCTGGAGTGCTCGGGGATCGGCCACCTGCCACGCCCGCCGCAACCCGAGCTGCAGGGGCCGGCCGAGCCCCCCGCCCCCACCGCGCAATCCGTCCGCCCGATGCTGTTCAGTGCCGACGGTGaa                                                                                     <  4:140197/137‑1 (MQ=255)                                     gAGCTGGTGAACCTGGAGTGCTCGGTGATCGGCCACCTGCCACGCCCGCCGCAACCCGAGCTGCAGGGGCCGGCCGAGCCCCCCGCCCCCACCGCGCAATCCGTCCGCCCGATGCTGTTCAGTGCCGACGGTGAATg                                                                                   <  4:317343/137‑1 (MQ=255)                                          ggTGAACCTGGAGTGCTCGGTGATCGGCCACCTGCCACGCCCGCCGCAACCCGAGCTGCAGGGGCCGGCCGAGCCCCCCGCCCCCACCGCGCAATCCGTCCGCCCGATGCTGTTCAGTGCCGACGGTGAATGGCAGc                                                                              <  1:38556/137‑1 (MQ=255)                                               aCCTGGAGTGCGCGGTGATCGGCCACCTGCCACGCCCGCCGCAACCCGAGCTGCAGGGGCCGGCCGAGCCCCCCGCCCCCACCGCGCAATCCGTCCGCCCGATGCTGTTCAGTGCCGACGGTGAATGGCAGCCcaca                                                                         <  3:193554/137‑1 (MQ=255)                                                      gTGCTCGGTGATCGGCCCCCTGCCACGCCCGCCGCAACCCGAGCTGCAGGGGCCGGCCGAGCCCCCCGCCCCCACCGCGCAATCCGTCCGCCCGCTGCTGTTCAGGGCCGGCGGGGAATGGCAGCCCACCCCGgtgg                                                                  >  8:117265/1‑136 (MQ=255)                                                            ggTGATCGGCCACCTGCCACGCCCGCCGCAACCCGAGCTGCAGGGGCCGGCCGCGGCCCCCCCCCCCCCCGCGCAATCCGGCCCCCCGCTGCAGTTCAGTGCCGCCGGGGAATGGCAGCCCCCACCGGTGTTCAAcc                                                            >  8:227735/1‑136 (MQ=255)                                                            ggTGATCGGCCACCTGCCACGCCCGCCGCAACCCGAGCTGCAGGGGCCGGCCGCGCCCCCCGCCCCCCCCCCCCAATACCTCCGCCCGAGGCTGTTCAGTGCCCACCGGGAATGGCAGCCCCCACCGCGGTTCAAcc                                                            >  7:22484/1‑136 (MQ=255)                                                            ggTGATCGGCCACCTGCCACGCCCGCCGCAACCCGAGCTGCAGGGGCCGGCCGAGCCCGCCGCCCCCCCCCCGCCATCCGTCCGCCCCCTGCTGTTAAGTGCCCGCCGTGAATGGCAGCCCACACCGGGGGTCAAcc                                                            >  7:258709/1‑136 (MQ=255)                                                            ggTGATCGGCCACCTGCCACGCCCGCCGCAACCCGAGCTGCAGGGGCCGGCCGAGCCCCCCCCCCCCACCGCGCAATCCGTCCGCCCGCTGCTGTTCAGTGCCGACGGTGAATGGCAGCCCACCCCGGGGTTCAACg                                                            >  8:65415/1‑137 (MQ=255)                                                              tGATCGGCCACCTGCCACGCCCGCCGCAACCCGAGCTGCAGGGGCCGGCCGAGCCCCCCGCCCCCACCGCGCAATCCGTCCGCCCGATGCTGTTCAGTGCCGACGGTGAATGGCAGCCCACACCGGTGTTCAACGGc                                                          <  4:191310/137‑1 (MQ=255)                                                              tGATCGGCCACCTGCCACGCCCGCCGCAACCCGAGCTGCAGGGGCCGGCCGAGCCCCCCGCCCCCACCGCGCAATCCGTCCGCCCGATGCTGTTCAGTGCCGACGGTGAATGGCAGCCCACACCGGTGTTCAACGGc                                                          <  6:223303/137‑1 (MQ=255)                                                              tGATCGGCCACCTGCCACGCCCGCCGCAACCCGAGCTGCAGGGGCCGGCCGAGCCCCCCGCCCCCACCGCGCAATCCGTCCGCCCGATGCTGTTCAGTGCCGACGGTGAATGGCAGCCCACACCGGTGTTCAACGGc                                                          <  1:96978/137‑1 (MQ=255)                                                                      cACCTGCCACGCCCGCCGCAACCCGAGCTGCAGGGGCCGGCCGCGCCCCCCGCCCCCCCCGCGCAATCCGTCCGCCCGCAGGAGTTCAGTGCCGACGGGGGATGGCAGCCCCCACCCGTGTTCAACGGCGaccaccc                                                  >  7:332586/1‑136 (MQ=255)                                                                      cACCTGCCACGCCCGCCGCAACCCGAGCTGCAGGGGCCGGCCGAGCCCCCCGCCCCCCCCGCGCAAACCGTCCGCCCGGGGCTGGTCAGTGCCGCCGGGGAATGGCCGCCCACACCGGGGTTCAACGGCAacccccg                                                  >  4:105110/1‑136 (MQ=255)                                                                         cTGCCACGCCCGCCGCAACCCGAGCTGCAGGGGCCGGCCGAGCCCCCCGCCCCCACCGCGCAATCCGTCCGCCCGATGCTGTTCAGTGCCGACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGtt                                               <  1:129085/137‑1 (MQ=255)                                                                           gCCACGCCCGCCGCAACCCGAGCTGCAGGGGCCGGCCGAGCCCCCCGCCCCCCCCGCGCAAtccgtcc                                                                                                                  >  6:98639/1‑68 (MQ=255)                                                                           gCCACGCCCGCCGCAACCCGAGCTGCAGGGGCCGGCCGAGCCCCCCGCCCCCACCGCGCAAtccgtcc                                                                                                                  <  5:98639/68‑1 (MQ=255)                                                                                 cccgccGCAACCCGAGCTGCAGGGGACGGCCGAGCCCCCCGCCCCCACCGCGCAATCCGGCCGCCCGATGCTGTTCAGTGCCGACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGc                                       <  6:214588/137‑1 (MQ=255)                                                                                  ccgccgCAACCCTAGCTGCAGGGGCCGGCCGAGCCCCCCGCCCCCACCGCGCAATCCGTCCGCCCGATGCTGTTCAGTGCCGACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCa                                      <  8:86418/137‑1 (MQ=255)                                                                                        cAACCCGAGCTGCAGGTGCCGGCCGAGCCCCCCGCCCCCACCGCGCAATCCGTCCGCCCGATGCTGTTCAGTGCCGACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGt                                <  3:80217/137‑1 (MQ=255)                                                                                        cAACCCGAGCTGCAGGGGCCGGCCGAGCCCCCCGCCCCCACCGCGCAATCCGTCCGCCCGATGCTGTTCAGTGCCGACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGt                                <  2:297387/137‑1 (MQ=255)                                                                                                      ggggCCGGCCGAGCCCCCCGCCCCCACCGCGCAATCCGTCCGCCCGATGCTGTTCAGTGCCGACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCtgtg                  <  2:76025/137‑1 (MQ=255)                                                                                                              ccgAGCCCCCCGCCCCCACCGCGCAATCCGTCCGCCCGATGCTGTTCAGTGCCGACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAg          <  3:204491/137‑1 (MQ=255)                                                                                                                gccccccccGCCACCACCTCGCAATCCGTCCGCCCGATGCTGTTCAGTGCCGACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGa        <  4:205079/134‑1 (MQ=255)                                                                                                                 aGCCCCCCGCCCCCTCCGCGCAATCCGTCCGCCCGATGCTGTTCAGTGCCGACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGaa       <  3:269105/137‑1 (MQ=255)                                                                                                                 aGCCCCCCGCCCCCACCGCGCAATCCGTCCGCCCGATGCTGTTCAGTGCCGACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGaa       <  2:151271/137‑1 (MQ=255)                                                                                                                 aGCCCCCCGCCCCCACCGCGCAATCCGTCCGCCCGATGCTGTTCAGTGCCGACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGaa       <  5:321436/137‑1 (MQ=255)                                                                                                                      cccGCCCCCACCGCGCAATCCGTCCGCCCGATGCTGTTCAGTGCCGACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGacac  <  2:303821/137‑1 (MQ=255)                                                                                                                      cccGCCCCCACCGCGCAATCCGTCCGCCCGATGCTGTTCAGTGCCGACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGacac  <  7:227735/137‑1 (MQ=255)                                                                                                                               |                                                                                                                                ACCCTGTACTCGTACAGCGAACCCGCCAGCCCGGTGGAGCTGGTGAACCTGGAGTGCTCGGTGATCGGCCACCTGCCACGCCCGCCGCAACCCGAGCTGCAGGGGCCGGCCGAGCCCCCCGCCCCCACCGCGCAATCCGTCCGCCCGATGCTGTTCAGTGCCGACGGTGAATGGCAGCCCACACCGGTGTTCAACGGCAACCACCTGTTGCCTGGGCAGACCGTGCATGGCCCCTGTGTGATCGAGGAAGACAC  >  ACB122/3982880‑3983133

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