Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	6,014,847	A→C	27.2%	K133N (AAA→AAC)	puuB_8 →	Gamma‑glutamylputrescine oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	6,014,847	0	A	C	27.2%	90.4 / 5.0	56	K133N (AAA→AAC)	puuB_8	Gamma‑glutamylputrescine oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (28/11);  new base C (2/15);  total (30/26)
Fisher's exact test for biased strand distribution p-value = 3.72e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCATCGGAGAAGAAGGCGTGCGCACCATGAAACTGATGGGCCTGGAGGCCGTACAGATCGTGCGCGAGCGGGTAGAAAAGCATGCTATCGCCTGCGACCTGGCCTGGGGCTACTGCGACCTGGCCAACAAACCCGCCGAGCTGCTGGGCTTTGCCCAAGACGCCGAGGAGCTGCGCGGTCTGGGCTATGCACATGAGCTGCGGCTGGTCGGCAAGGACGATATCCACAGCGTGGTGGGTGCCGACTGCTATGTGGGGGGGC  >  ACB122/6014717‑6014977                                                                                                                                   |                                                                                                                                   tCATCGGAGAAGAAGGCGTGCGCACCATGAAACTGATGGGCCTGGAGGCCGTACAGATCGTGCGCGAGCGGGTAGAAAAGCATGCTATCGCCTGCGACCTGGCCTGGGGCTACTGCGACCTGGCCAACAAACccg                                                                                                                                <  5:136738/135‑1 (MQ=255) tCATCGGAGAAGAAGGCGTGCGCACCATGAAACTGATGGGCCTGGAGGCCGTACAGATCGTGCGCGAGCGGGTAGAAAAGCATGCTATCGCCTGCGACCTGGCCTGGGGCTACTGCGACCTGGCCAACAAACccg                                                                                                                                <  8:133169/135‑1 (MQ=255)        agaagaAGGCGTGCGCACCATGAAACTGATGGGCCTGGAGGCCGTACAGATCGTGCGCGAGCGGGTAGAAAAGCATGTTATCGCCTGCGACCTGGCCTGGGGCTACTGCGACCTGGCCAACAAACCCGCCGAgct                                                                                                                         >  5:13486/1‑135 (MQ=255)                       cACCATGAAACTGATGGGCCTGGAGGCCGTACAGATCGTGCGCGAGCGGGTAGAAAAGCATGCTATCGCCTGCGACCTGGCCTGGGGCTACTGCGACCTGGCCAACAACCCCGCCGAGCTGCTGGGCTTTGCCCa                                                                                                          <  8:381671/135‑1 (MQ=255)                       cACCATGAAACTGATGGGCCTGGAGGCCGTACAGATCGTGCGCGAGCGGGTAGAAAAGCATGCTATCGCCTGCGACCTGGCCTGGGGCTACTGCGACCTGGCCAACAAACCCGCCGAGCTGCTGGGCTTTGcc                                                                                                            >  7:326876/1‑133 (MQ=255)                          cATGAAACTGATGGGCCTGGAGGCCGTACAGATCGTGCGCGAGCGGGTAGAAAAGCATGCTATCGCCTGCGACCTGGCCTGGGGCTACTGCGACCTGGCCAACAAACCCGCCGAGCTGCTGGGCTTTGCCCAAGa                                                                                                       >  5:147416/1‑135 (MQ=255)                               acgtGAGGGCCCGGGAGGCCGAACACATCGTGCGCGAGGGGGTAGAAACCCATGCTATCGCCTGCGACCTGGCCTGGGGCAACTGCGACCTGGCCAACAACCCCGCCGAGCTGCGGGGCTTTGCCCAAGACGCCg                                                                                                  <  8:58786/132‑1 (MQ=255)                                     tGGGCCGGTCGGCCGTAAAGATCGGGCGGGAGCGGGAATAAACCCAGGCTATCCCCTGCCACCTGGCCTGGGGGAACTGCGACCGGGCCAACAACCCCGCCGCCCTGCTGGGGTTTGCCCAAGACGCCGAGGAGc                                                                                            <  6:262098/135‑1 (MQ=255)                                       ggCCTGGAGGCCGTACTGCTCGTGCGCGAGCGGGTAGAAAAGCATGCTATCGCGTGCGACCTGGCCTGGGGCTACTGCGACCTGGCCAACAACCCCGCCGAGCTGCTGGGCTTGGCCCAAGACGCCGAGGCGCTg                                                                                          >  5:56740/1‑135 (MQ=255)                                       ggCCTGGAGGCCGTACAGATCGTGCGCGAGCGGGTAGAAAAGCATGCTATCGCCTGCGACCTGGCCTGGGGCTACTGCGACCTGGCCAACAAACCCGCCGAGCTGCTGGGCTTTGCCCAAGACGCCGAGGAGCTg                                                                                          >  4:270920/1‑135 (MQ=255)                                            ggAGGCCGTACAGATCGTGCGCGAGCGGGTAGAAAAGCATGCTATCGCCTGCGACCTGGCCTGGGGCTACTGCGACCTGGCCAACAAACCCGCCGAGCTGCTGGGCTTTGCCCAAGACGCCGAGGAGCTGCGCgg                                                                                     >  5:66136/1‑135 (MQ=255)                                            ggAGGCCGTACAGATCGGGCGCGAGCGGGTAGAAAAGCATGCTATCGCCGGCGACCTGGCGTGGGGCTACTGCGACCTGCCCAACAAACCCGCCGAGCTGCTGGGCTTTGCCCAAGACGCCGAGGAGCTGCGCgg                                                                                     <  2:314567/135‑1 (MQ=255)                                               ggCGGTACTGATGGGGCGCGTGCGGGTAGAAAAGCATGCTATCGCCTGCGCCCTGGCCTGGGGCTACTGCGACCTGGCAAACAACCCCGCCGAGCTGCTGGGCTTTGCCCAAGACGCCGAGGAGCTGCGCGGTCt                                                                                  <  2:110988/135‑1 (MQ=255)                                               ggCCTTACCCATCGTGCTCGAGCGTGTATAAAAGATTGCTATCGCCTGCGACCGGGCGTGGGGCTACTGCCACCTGCACAACAACCCCGCCTAGCTGCTGGGCTTTGCCCAACACGCCGAGGAGCTGCGCGGTCt                                                                                  <  6:115720/135‑1 (MQ=255)                                                     aCAGATCGTGCGCGAGCGGGTAGAAAAGCATGCTATCGCCTGCGACCTGGCCTGGGGCTACTGCGACCTGGCCAACAAACCCGCCGAGCTGCTGGGCTTTGCCCAAGACGCCGAGGAGCTGCGCGGTCTGGGCTa                                                                            >  8:43425/1‑135 (MQ=255)                                                     aCAGATCGTGCGCGAGCGGGTAGAAAAGCATGCTATCGCCTGCGACCTGGCCTGGGGCTACTGCGACCTGGCCAACAAACCCGCCGAGCTGCTGGGCTTTGCCCAAGACGCCGAGGAGCTGCGCGGTCTGGGCTa                                                                            >  5:133028/1‑135 (MQ=255)                                                       aGATCGTGCGCGAGCGGGTAGAAAAGCATGCTATCGCCTGCGACCTGGCCTGGGGCTACTGCGACCTGGCCAACAAACCCGCCGAGCTGCTGGGCTTTGCCCAAGACGCCGAGGAGCTGCGCGGTCTGGGCTAt                                                                           >  3:126056/1‑134 (MQ=255)                                                            gTGCGCGAGCGGGTAGAAAAGCATGCTATCGCCTGCGACCTGGCCTGGGGCTACTGCGACCTGGCCAACAAACCCGCCGAGCTGCTGGGCTTTGCCCAAGACGCCGAGGAGCTGCGCGGTCTGGGCTATGcaca                                                                      >  4:194060/1‑134 (MQ=255)                                                             tGCGCGAGCGGGTAGAAAAGCATGCTATCGCCTGCGACCTGGCCTGGGGCTACTGCGACCTGGCCAACAAACCCGCCGAGCTGCTGGGCTTTGCCCAAGACGCCGAGGAGCTGCGCGGTCTGGGCTATGcaca                                                                      >  7:57151/1‑133 (MQ=255)                                                               cgcgAGCGGGTAGAAAAGCATGCTATCGCCTGCGACCTGGCCTGGGGCTACTGCGACCTGGCCAACAAACCCGCCGAGCTGCTGGGCTTTGCCCAAGACGCCGAGGAGCTGCGCGGTCTGGGCTATGCACATGAg                                                                  >  3:222543/1‑135 (MQ=255)                                                                 cgAGCGGGTAGAAAAGCATGCTATCGCCTGCGACCTGGCCTGGGGCTACTGCGACCTGGCCAACAAACCCGCCGAGCTGCTGGGCTTTGCCCAAGACGCCGAGGAGCTGCGCGGTCTGGGCTATGCACATGAGCt                                                                >  4:23960/1‑135 (MQ=255)                                                                 cgAGCGGGTAGAAAAGCATGCTATCGCCTGCGACCTGGACTCGGGCTACTGCGCCCTGGCCAACAAACACGCCGAGCTGCTGGGCTTTGCCCAAGCCGCCGAGGAGCTGCGCGGTCTGGGCTATGCACATGAGCt                                                                >  3:200175/1‑135 (MQ=255)                                                                  gAGCGGGTAGAAAAGCATGCTATCGCCTGCGACCTGGCCTGGGGCTACTGCGACCTGGCCAACAAACCCGCCGAGCTGCTGGGCTTTGCCCAAGACGCCGAGGAGCTGCGCGGTCTGGGCTATGCACATGAGCTg                                                               >  4:90540/1‑135 (MQ=255)                                                                  gAGCGGGTAGAAAAGCATGCTATCGCCTGCGACCTGGCCTGGGGCTACTGCGACCTGGCCAACAAACCCGCCGAGCTGCTGGGCTTTGCCCAAGACGCCGAGGAGCTGCGCGGTCTGGGCTATGCACATGAGCTg                                                               >  6:51671/1‑135 (MQ=255)                                                                             aacgCATGCTACCGCCTGCGACCGGGCCGGGGGCTACTGCCACCTGGCCAACAACCCCGCCGCGCGGCTGGGCTTTGCCCAAGACGCCGAGGAGCTGCGCGGTCTGGGCTATGCACATGAGCTGCGGCTGGTCgg                                                    <  3:291272/132‑1 (MQ=255)                                                                                 cATGATACCGCGTGCGACCTGGCGTGGGGCTACTGCGACCGGGCCAACCACCCCGCCGCGCTGCTGGGCTTGGCCCAAGACGCCGAGGAGCTGCGCGGTCTGGGCTATGCACATGAGCTGCGGCTGGTCGGCAAg                                                <  4:376628/135‑1 (MQ=255)                                                                                  aggCTAGCGCCGGCGACCTGCCGTGGGGCTACTGCGACCGGCCCAACACACCCCCCGAGCGGCTGGGCTTTGCCCAAGACGCCGAGGAGCTGCGCGGTCTGGGCTATGCACATGAGCTGCGGCTGGTCGGCAAgg                                               <  6:243825/133‑1 (MQ=255)                                                                                     cTATCGCCTGCGACCTGGCCTGGGGCTACTGCGACCTGGCCAACAAACCCGCCGAGCTGCTGGGCTTTGCCCAAGACGCCGAGGAGCTGCGCGGTCTGGGCTATGCACATGAGCTGCGGCTGGTCGGCAAGGACg                                            >  5:346188/1‑135 (MQ=255)                                                                                     cTATCGCCTGCGACCTGGCCTGCGGCTACTGCGACCTGGCCAACAAACCCGCCGAGCTGCTGGGCTTTGCCCAAGACGCCGAGGAGCTGCGCGGTCTGGGCTATGCACATGAGCTGCCGCTGGTCGGCAAGGACg                                            >  1:86817/1‑135 (MQ=255)                                                                                      tATCGCCTGCGACCTGGCCTGGGGCTACTGCGACCTGGCCAACAAACCCGCCGAGCTGCTGGGCTTTGCCCAAGACGCCGAGGAGCTGCGCGGTCTGGGCTATGCACATGAGCTGCGGCTGGTCGGCAAGGACg                                            >  5:243825/1‑134 (MQ=255)                                                                                         cGCCTGCGACCTGGCCTGGGGCTAGTGCTACCTGGCCAACAACCCCGCCGAGCGGCTGGGCTTTGCCCAAGACGCCGAGGAGCTGCGCGGTCTGGGCTATGCACATGAGCTGCGGCTGGTCGGCAAGGACGat                                          <  7:238868/133‑1 (MQ=255)                                                                                         cGCCTGCGACCTGGCCTGGGGCTACTGCGACCTGGCCAACAAACCCGCCGAGCTGCTGGGCTTTGCCCAAGACGCCGAGGAGCTGCGCGGTCTGGGCTATGCACATGAGCTGCGGCTGGTCGGCAAGGACGatat                                        >  3:376628/1‑135 (MQ=255)                                                                                         cGCCTGCGACCTGGCCTGGGGCTACTGCGACCTGGCCAACAAACCCGCCGAGCTGCTGGGCTTTGCCCAAGACGCCGAGGAGCTGCGCGGTCTGGGCTATGCACATGAGCTGCGGCTGGTCGGCAAGGACGatat                                        <  6:147416/135‑1 (MQ=255)                                                                                          gcgtGCGCCCTGGCGTGGGGGTAGTGCGCCCGGGCAAACAAACCGGCGGAGCTGCTGGGCTTTCCCCAAGACGCCGAGGAGCTGCGCGGTCTGGGCTATGCACATGAGCTGCGGCTGGTCGGCAAGGACGatat                                        <  4:222543/131‑1 (MQ=255)                                                                                          gcgtGCGACCTGGCGTGGGGCTACTGCGACCTGGCCAACAAACCCCCCTAGGTGCTGGGCTTTGCCAAAGACGCCGAGGAGCTGCGCGGTCTGGGCTATGCACATGAGCTGCGGCTGGTCGGCAAGGACGATATc                                       <  8:57151/132‑1 (MQ=255)                                                                                          gCCTGCGACGTGGCCTGGGCCTAGTGCGACGGGGCAAACAACCCCCCCGAGCTGCGGGGCTTTGCCCAAGACGCCGAGGAGCTGCGCGGTCTGGGCTATGCACATGAGCTGCGGCTGGTCGGCAAGGACGATATc                                       <  7:96965/135‑1 (MQ=255)                                                                                          gCCTGCCCCATGGCCGGGGGACACTGCGCCCGGCCCAAAACCCCCGCCGCGCTGCTGGGGTTTCCCCAAGACGCCGAGGAGCTGCGCGGTCTGGGCTATGCACATGAGCTGCGGCTGGTCGGCAAGGACGATATc                                       <  1:311510/135‑1 (MQ=255)                                                                                           ccTGCGACCTGGCCTGGGGCTACTGCGACCTGGCCAACAAACCCGCCGAGCTGCTGGGCTTTGCCCAAGACGCCGAGGAGCTGCGCGGTCTGGGCTATGCACATGAGCTGCGGCTGGTCGGCAAGGACGATATc                                       <  3:90540/134‑1 (MQ=255)                                                                                              gCGACCTGGCCTGGGGCTACTGCGACCTGGCCAACAACCCCGCCGAGCTGCTGGGCTTTGCCCAAGCCGCCGCGGAGCTGCGCGGTCTGGGCTATGCCCAGGAGCTGCGGCTGGTCGGCAAGGACGATATCcaca                                   >  5:30685/1‑135 (MQ=255)                                                                                                  cctggcctggGGCTACTGCGACCTGGCCAACAAACCCGCCGAGCTGCTGGGCTTTGCCCAAGACGCCGAGGAGCTGCGCGGTCTGGGCTATGCACATGAGCTGCGGCTGGTCGGCAAGGACGATATCCACAGCgt                               >  5:294582/1‑135 (MQ=255)                                                                                                  cctggcctggGGCTACTGCGACCTGGCCAACAAACCCGCCGAGCTGCTGGGCTTTGCCCAAGACGCCGAGGAGCTGCGCGGTCTGGGCTATGCACATGAGCTGCGGCTGGTCGGCAAGGACGATATCCACAGCgt                               >  4:162574/1‑135 (MQ=255)                                                                                                        cgggggCTACTGCGACCGGGCCAACAACCCCGCCGAGCTGCTGGGCTTTGCCCAAGACGCCGAGGAGCTGCGCGGTCTGGGCTATGCACATGAGCTGCGGCTGGTCGGCAAGGACGATAT‑CACAGCGTggtgggt                        <  3:342986/133‑1 (MQ=255)                                                                                                         tggGGAGACTGCGACCTGGCCAACACCCCCGCCGAGCGGCGGGGCTTTGCCCAGGACGCCGAGGAGCTGCGCGGTCTGGGCTATGCACATGAGCTGCGGCTGGTCGGCAAGGACGATATCCACAGCGTggtgggt                        <  6:248303/135‑1 (MQ=255)                                                                                                           gggCTACTGCGACCTGGCAAACAAACCCGCCGAGCTGCGGGGCTTTGCCCAAGACGCCGAGGAGCTGCGCGGTCTGGGCTATGCACATGAGCTGCGGCTGGTCGGCAAGGACGATATCCACAGCGTGGTGGGTGc                      <  2:214200/135‑1 (MQ=255)                                                                                                            ggCTACTGCGACCTGGCCAACAACCCCGCCGAGCGGCTGGGCTTGGCCCAAGACGCCGAGGAGCTGCGCGGTCTGGGCTATGCACATGAGCTGCGGCTGGTCGGCAAGGACGATATCCACAGCGTGGTGGGTGcc                     <  6:164066/135‑1 (MQ=255)                                                                                                              cacagccGACCTGCCCAAAAACCCCCCCTCGCGGGGGGGCTTGGCCCAGGACGCCGAGGAGCTGCGCGGTCTGGGCTATGCACATGAGCTGCGGCTGGTCGGCAAGGACGATATCCACAGCGTGGTGGGTGCCGa                   <  2:168025/129‑1 (MQ=255)                                                                                                              cTACTGCGACGTGGCAAACAACCCCGCCTATCTGATGGGGTTTGCCAAAGACGCCGAGGAGCTGCGCGGTCTGGGCTATGCACATGAGCTGCGGCTGGTCGGCAAGGACGATATCCACAGCGTGGTGGGTGCCGa                   <  6:353148/135‑1 (MQ=255)                                                                                                              cTACTGCGACCTGGCCAACAAACCCGCCGAGCTGCTGGGCTTTGCCCAAGACGCCGAGGAGCTGCGCGGTCTGGGCTATGCACATGAGCTGCGGCTGGTCGGCAAGGACGATATCCACAGCGTGGTGGGTGCCGa                   >  7:362246/1‑135 (MQ=255)                                                                                                              cTACTGCGACCTGGCCAACAAACCCGCCGAGCTGCGGGGCTTTGCCCAAGACGCCGAGGAGCTGCGCGGTCTGGGCTATGCACATGAGCTGCGGCTGGTCGGCAAGGACGATATCCACAGCGTGGTGGGTGCCGa                   <  7:43425/135‑1 (MQ=255)                                                                                                                  tGCGACCTGCACAACAAACCCGCCGAGCGGCGGGGCTTTGCCCAAGACGCCGAGGAGCTGCGCGGTCTGGGCTATGCACATGAGCTGCGGCTGGTCGGCAAGGACGATATCCACAGCGTGGTGGGTGCCGACTg                <  6:13486/134‑1 (MQ=255)                                                                                                                   gCGACCTGGCCAACAAACCCGCCGAGCTGCTGGGCTTTGCCCAAGACGCCGAGGAGCTGCGCGGTCTGGGCTATGCACATGAGCTGCGGCTGGTCGGCAAGGACGATATCCACAGCGTGGTGGGTGCCGACTGCt              >  7:346607/1‑135 (MQ=255)                                                                                                                    cGACCTGGCCAACAAACCCGCCGAGCTGCTGGGCTTTGCCCAAGACGCCGAGGAGCTGCGCGGTCTGGGCTATGCACATGAGCTGCGGCTGGTCGGCAAGGACGATATCCACAGCGTGGTGGGTGCCGACTGCTa             >  5:326755/1‑135 (MQ=255)                                                                                                                         tGGCCAACAAACCCGCCGAGCTGCTGGGCTTTGCCCAAGACGCCGAGGAGCTGCGCGGTCTGGGCTATGCACATGAGCTGCGGCTGGTCGGCAAGGACGATATCCACAGCGTGGTGGGTGCCGACTGCTATGTgg        >  8:205807/1‑135 (MQ=255)                                                                                                                         tGGCCAACAAACCCGCCGAGCTGCTGGGCTTTGCCCAAGACGCCGAGGAGCTGCGCGGTCTGGGCTATGCACATGAGCTGCGGCTGGTCGGCAAGGACGATATCCACAGCGTGGTGGGTGCCGACTGCTATGTgg        >  2:41020/1‑135 (MQ=255)                                                                                                                         tGGCCAACAAACCCGCCGAGCTGCTGGGCTTTGCCCAAGACGCCGAGGAGCTGCGCGGTCTGGGCTATGCACATGAGCTGCGGCTGGTCGGCAAGGACGATATCCACAGCGTGGTGGGTGCCGACTGCTATGTgg        >  1:284354/1‑135 (MQ=255)                                                                                                                              aacaaACCCGCCGAGCTGCTGGGCTTTGCCCAAGACGCCGAGGAGCTGCGCGGTCTGGGCTATGCACATGAGCTGCGGCTGGTCGGCAAGGACGATATCCACAGCGTGGTGGGTGCCGACTGCTATGTggggggg   >  4:351406/1‑135 (MQ=255)                                                                                                                               aaaaaCCCGCCGAGCTGCGTGTCTTTGCCCTAGACGCCGAGGAGCTGCGCGGTCTGGGCTAGGCACATGAGCTGCGGCTGGTCGGCAAGGACGATATCCACAGCGTGGTGGGTGCCGACTGCTATGTGGGGGGGc  <  8:350724/133‑1 (MQ=255)                                                                                                                                   |                                                                                                                                   TCATCGGAGAAGAAGGCGTGCGCACCATGAAACTGATGGGCCTGGAGGCCGTACAGATCGTGCGCGAGCGGGTAGAAAAGCATGCTATCGCCTGCGACCTGGCCTGGGGCTACTGCGACCTGGCCAACAAACCCGCCGAGCTGCTGGGCTTTGCCCAAGACGCCGAGGAGCTGCGCGGTCTGGGCTATGCACATGAGCTGCGGCTGGTCGGCAAGGACGATATCCACAGCGTGGTGGGTGCCGACTGCTATGTGGGGGGGC  >  ACB122/6014717‑6014977

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