Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	4,135,607	A→C	37.0%	Q32P (CAG→CCG)	tgnC →	(Z)‑2‑((N‑methylformamido)methylene)‑5‑hydroxybu tyrolactone dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	4,135,607	0	A	C	37.0%	31.6 / 7.1	46	Q32P (CAG→CCG)	tgnC	(Z)‑2‑((N‑methylformamido)methylene)‑5‑hydroxybu tyrolactone dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (10/17);  new base C (9/9);  total (19/27)
Fisher's exact test for biased strand distribution p-value = 5.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGCCCGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAGCTGCTGCGTCGTCTG  >  ACB122/4135473‑4135740                                                                                                                                       |                                                                                                                                      acgCCCGCACCCTCCCCGAATCCGCCACGGGGATTGTTGCATGACCCGCGGGCGTTGCCAGATGTGCATCGACGGCCAAGGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCa                                                                                                                                       <  4:462334/133‑1 (MQ=255) aaGCCCGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCAGCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCa                                                                                                                                       <  7:389927/135‑1 (MQ=255)   gCCCGCACCCTCACCGAATCCGCCAAGTGGATTGTTGCATGACGCTCGTTCGTTTCCAGATGTGCATCCACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCAACCGCCCAgg                                                                                                                                     <  6:129235/135‑1 (MQ=255)    cccGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCAGGCCACCGCCCAGGc                                                                                                                                    >  2:166346/1‑135 (MQ=255)      cGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGCCAGCCTCACCCCGGCCCCCGCCCCGGccg                                                                                                                                  >  1:366792/1‑134 (MQ=255)      cGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCt                                                                                                                                  >  2:345488/1‑135 (MQ=255)           ccTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGgc                                                                                                                              >  4:109068/1‑134 (MQ=255)              cACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCACGCCGGGgcgcc                                                                                                                           >  8:112346/1‑133 (MQ=255)                ccGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCACGCCGGGGCGCCAcg                                                                                                                        >  5:165529/1‑134 (MQ=255)                  gAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCCGCCGGGGCGCaccagc                                                                                                                      >  3:328606/1‑130 (MQ=255)                   aaTCGGCCAAGGGGATTGTTGCATGACCGTGGTTGTTTTCCAGATGTGAAACGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGCGGGGGCGCAACTGcc                                                                                                                     <  8:246271/135‑1 (MQ=255)                     tCCGCCACGGGGTTTTTTGAATGACCGTCGTGCGTTTCCAGTTGTGCAGCGACGGCCAATGGGGCGATCCCCAGAGCGCCAAGACCTTCGACAGCCACACCCCGCCCCCCGCCCAGGCCGGGGCGCAACTGCCCg                                                                                                                   <  4:229113/135‑1 (MQ=255)                     gccGCCCAGGGGATTGTTGCATGACCCTCGTGCGTTGCCAGATGGGCATCGACGGCCAATGGCGCGAGGCCCAGAGCGCCAAGACCTTCGACAGCCTCACCCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCg                                                                                                                   <  3:324778/134‑1 (MQ=255)                         ccAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCAGCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGc                                                                                                               <  4:36011/135‑1 (MQ=255)                          cAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGCCTGGGCGCAACTGCCCgccgcc                                                                                                              >  3:310725/1‑135 (MQ=255)                            cggggATTGTTGCTGGCCCCTCGTTCGTTTCCAGATGTGCATCGCCGCCCACTGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGa                                                                                                            <  8:327020/134‑1 (MQ=255)                                 aTTGTTGCATGACCCTCGTTCTTTTCCTCATGTGCATCTACGGCCTAGGGGGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTAAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAg                                                                                                       <  4:283858/135‑1 (MQ=255)                                    gttgCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGCCGGGGCGCCACTGCCCGCCGCCGCTGAAGCCg                                                                                                    >  5:152978/1‑135 (MQ=255)                                     ttgCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCa                                                                                                                                       >  8:183678/1‑99 (MQ=255)                                     ttgCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCa                                                                                                                                       <  7:183678/99‑1 (MQ=255)                                        cATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGAtgt                                                                                                >  4:40792/1‑135 (MQ=255)                                           gACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCACCCCGGCCACCGCCCCCGCCAGGGCGCAACAGCCCGACGCCGAGGAACCCGCAGTGGa                                                                                             >  8:290341/1‑135 (MQ=255)                                             cccTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAg                                                                                            <  8:450377/134‑1 (MQ=255)                                                        tgcccgATGTGCAGCGCCGGCCACTGGCGCGAGGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAg                                                                                <  8:159193/130‑1 (MQ=255)                                                          ccctttgtcCATCGCCGGCCAATGGCGCGAGGCCAAGAGCGCCAAGACCTTCGACCGCCCCACCCCGCCCACCGCCCCGGCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAgg                                                                               <  3:109068/129‑1 (MQ=255)                                                           cAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCCCCCAGGCCAGGGCGCAACTGCCCCCCGCCGAATAAGACGAATTGGAGCTGGCCGGTAAGgcc                                                                             >  7:408462/1‑135 (MQ=255)                                                                   cATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGCCGGGGCGGCCCTGCCCGCCGCCGCGGGAGCCGCGGTGGGGCTGGCCGGTGAGGCCGccccccg                                                                     >  7:271907/1‑131 (MQ=255)                                                                    aTCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGCCTGGGCGCCACGGCCCGCCGCCGAGGCAGCCGCTGTGGAGCTGGCCGGTCAGGCCGCCCCGCGt                                                                    >  1:86476/1‑135 (MQ=255)                                                                    aTCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGCCGCCGCTGAAGCCGATGTGGAGCTGGCCGTTCACGCCGCCCAGCGt                                                                    >  7:155546/1‑135 (MQ=255)                                                                      cgccgcccACTGGCGCGATGCCCCGAGCGGCACCCCCTCCGACAGCCTCACCCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGCCGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGc                                                                  <  5:238837/129‑1 (MQ=255)                                                                        aCGCCCACTGGCGCGAGGCCCCGAGCGGAACGCCCTTCGACCGCCTCACCCCGCCCCCCGCCCCGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGc                                                                  <  2:172302/133‑1 (MQ=255)                                                                         cGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCCCCGCCCCGGCCTGGGCGCCACTGCCCGCCGCCGCTGAAGCCGATGTGGAGCTGGCCGGTCAGGCCGCCCACCGGGccgt                                                               >  2:203105/1‑133 (MQ=255)                                                                              aaTGGCGCGAGGCCCAGAGCGGCAAGACCTGCGACAGCCTCAACCCGGCCCCCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACa                                                          <  6:464864/135‑1 (MQ=255)                                                                                   cgcgATGCCCAGAGCGGCAAGACCTTCGACAGCCTCCACCCGGCCACCGCCCCGGCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAg                                                     <  7:328316/135‑1 (MQ=255)                                                                                        ggCCCAGCGCGGCAAGACCTTCGACAGCCGCACCCCGGCCCCCCCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTg                                                <  6:192886/134‑1 (MQ=255)                                                                                            cAGCGCCGCAACACCTTCGACAGCCTCACCCCGCCCCCCGCCCCGGCGGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGgcgc                                            <  8:60432/135‑1 (MQ=255)                                                                                            cAGAGCGGCAAGACCTGCGACAGCCTCACCCCGGCCACCGCCCCGGCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGGGCCTTCGACAGCAAGGCGTGgcgc                                            <  6:170843/135‑1 (MQ=255)                                                                                            cAGAGCGGCAAGACCTCCGACAGCCTCACCCCGGCCACCGCCCTGGCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGgcgc                                            <  6:163055/135‑1 (MQ=255)                                                                                                 cGGCAAGACCTGCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCAt                                       <  4:328606/135‑1 (MQ=255)                                                                                                      agcccTTCGACAGCCCCACCCCGGCCCCCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCg                                  <  2:151982/132‑1 (MQ=255)                                                                                                           ttCGACCGCCTCACCCCGCCCACCCCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTAcc                             <  1:82660/135‑1 (MQ=255)                                                                                                           ttCGACAGCCGCACCCCGGCCACCGCCCCGGCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTAcc                             <  2:40996/135‑1 (MQ=255)                                                                                                                  gCCTCACCCCGGCCCCCGCCCCGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCcgcg                      <  5:287772/135‑1 (MQ=255)                                                                                                                  gCCTCACCCCGGCCCCCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCcgcg                      <  1:148946/135‑1 (MQ=255)                                                                                                                  cccccAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCcgcg                      <  2:150318/131‑1 (MQ=255)                                                                                                                      cAACCCGGCCACCGCCCCGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGCTGTGGAGCTGGCCGTTGAGGCCCCCCCGCGTGCCGTCGCCCGCAAGGCGTGGCGCAGCATCCCCGCGACCGCCCCCGGCaa                  >  8:189562/1‑135 (MQ=255)                                                                                                                      cAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGCCGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCaa                  >  4:321367/1‑135 (MQ=255)                                                                                                                                aCCGCCCAGGCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGCCAGCAAGGCGTGGCGCAGCATCCCCGCTACCGCCCGCCGCCAGCTgcggcgg        >  4:41711/1‑134 (MQ=255)                                                                                                                                      cAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAGCTGCTGCGTCGTCTg  <  7:343202/135‑1 (MQ=255)                                                                                                                                      cAGGCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAGCTGCTGCGTCGTCTg  <  5:70470/135‑1 (MQ=255)                                                                                                                                       |                                                                                                                                      AAGCCCGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAGCTGCTGCGTCGTCTG  >  ACB122/4135473‑4135740

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: ‑