Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	1,605,593	A→C	23.6%	S112R (AGC→CGC)	ung →	Uracil‑DNA glycosylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	1,605,593	0	A	C	23.6%	97.0 / 5.1	56	S112R (AGC→CGC)	ung	Uracil‑DNA glycosylase
Reads supporting (aligned to +/- strand):  ref base A (17/24);  new base C (11/4);  total (28/28)
Fisher's exact test for biased strand distribution p-value = 6.83e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGGTCAGGACCCGTACCACGGGCCGGGCCAGGCACATGGCTTGTGCTTCTCGGTGCAACCGGGCGTGGCCACGCCACCGTCGCTGGTCAATATCTACAAGGAGCTGCAGCGCGACCTGAACATCCCGATCCCCAGCCACGGCTACCTCCAGAGCTGGGCGGAACAAGGCGTTCTGTTGCTCAACACGACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCGACAGGA  >  ACB122/1605459‑1605716                                                                                                                                       |                                                                                                                            tCGGTCAGGACCCGTACCACGGGCCGGGCCAGGCACATGGCTTGTGCTTCTCGGTGCAACCGGGCGTGGCCACGCCACCGTCGCTGGTCAATATCTACAAGGAGCTGCAGCGCGACCTGAACATCCCGATCCCCAg                                                                                                                            >  4:268395/1‑136 (MQ=255)  cGGTCAGGACCCGTACCCCGGGCCGGGCCAGGCACAGGGCTTGTGCTGCTCGGTGCACCCGGGCGTGGCCACGCCACCGGCGCTGGTCAATATCTACAAGGAGCTGCAGCGCGACCTGAACATCCCGATCCCCAg                                                                                                                            <  1:162031/135‑1 (MQ=255)      cAGGTCCAGTACCTCGGGCGGGGCCAGGCACATGGGTTGTGGTTGTCGTTTCATCCGGGGGTGGCCAGGCCACCGTCGCTGGTCAATATCTACTAGGAGCTGCAGCGCGACCTGAACATCCCGATCCCCAGCCACg                                                                                                                       <  8:43951/136‑1 (MQ=255)      cAGGACCCGTACCACGGGCCGGGCCAGGCACATGGCTTGTGCTTCTCGGTGCACCGGGGCGTGGCCACGCCACCGTCGCTGGTCAATATCTACAAGGAGCTGCAGCGCGACCTGAACATCCCGATCCCCAGCCACg                                                                                                                       <  7:193704/136‑1 (MQ=255)      cAGGACCCGTACCACGGGCCGGGCCAGGCACATGGCTTGTGCTTCTCGGTGCAACCGGGCGTGGCCACGCCACCGTCGCTGGTCAATATCTACAAGGAGCTGCAGCGCGACCTGAACATCCCGATCCCCAGCCACg                                                                                                                       <  4:142178/136‑1 (MQ=255)       tggccccGTACCCGGGGGCGGGCCTGGCCTTTGGGTTGTGCTGCTCGGGGCAACGGGGGGGGGCCACGCCACCGTCGCTGGTCAATATCTACAAGGAGCTGCAGCGCGACCTGAACATCCCGATCCCCAGCCACg                                                                                                                       <  5:372394/131‑1 (MQ=255)       cggACCGGTACCCGGGGCGGGGCCAGGATCTGGGCTTGTGCTTCGGGGTGCACGGGGGCGTGGCCACGCCACCGTCGCTGGTCAATATCTACAAGGAGCTGCAGCGCGACCTGAACATCCCGATCCCCAGCCACg                                                                                                                       <  5:8637/134‑1 (MQ=255)       aGGACCCGTACCACGGGCCGGGCCAGGCACATGGCTTGTGCTTCTCGGTGCAACCGGGCGTGGCCACGCCACCGTCGCTGGTCAATATCTACAAGGAGCTGCAGCGCGACCTGAACATCCCGATCCCCAGCCACg                                                                                                                       <  3:268395/135‑1 (MQ=255)       aGGACCCGTACCACGGGCCGGGCCAGGCACATGGCTTGTGCTTCTCGGTGCAACCGGGCGTGGCCACGCCACCGTCGCTGGTCAATATCTACAAGGAGCTGCAGCGCGACCTGAACATCCCGATCCCCAGCCACg                                                                                                                       <  6:80011/135‑1 (MQ=255)        ggACCCGTACCACGGGCCGGGCCAGGCACTTGGCTTGTGCTTCTCGGTGCAACCGGGCGGGGCCACGCCACCGTCGCTGGTCAATATCTACAAGGAGCTGCAGCGCGACCTGAACATCCCGATCCCCAGCCACGGc                                                                                                                     <  3:140059/136‑1 (MQ=255)        ggACCCGTACCACGGGCCGGGCCAGGCACATGGCTTGTGCTGCTCGGTGCAACCGGGCGTGGCCACGCCACCGTCGCTGGTCAATATCTACAAGGAGCTGCAGCGCGACCTGAACATCCCGATCCCCAGCCACg                                                                                                                       <  5:48227/134‑1 (MQ=255)                 ccAGGGGCCGGGCCAGGCACATGGCTTGTGCTTCTCGGTGCAACCGGGCGTGGCCACGCCACCGTCGCTGGTCAATATCTACAAGGAGCTGCAGCGCGACCTGAACATCCCGATCCCCAGCCACGGCTACCTCCag                                                                                                            <  8:121520/136‑1 (MQ=255)                  cACGGGCCGGGCCAGGCACATGGCTTGTGCTTCTCGGTGCAACCGGGCGTGGCCCCGCCCCCGGCGCGGGGGAATATATACAACGAGCGTCAGCGCGGCCTGAACATCCCGCTCCCCCCCCACCGCGACCGCCag                                                                                                            >  5:72009/1‑135 (MQ=255)                  cACGGGCCGGGCCAGGCACATGGCTTGTGCTTCTCGGTGCAACCGGGCGTGGCCCCGCCACCGTCGCTGGTCAATATCTACAAGGAGCTGCAGCGCGACCTGCACATCCCCCTCCCCAGCCCCCGCTACCTCCag                                                                                                            >  2:58559/1‑135 (MQ=255)                  cACGGGCCGGGCCAGGCACATGGCTTGTGCTTCTCGGTGCAACCGGGCGTGGCCACGCCCCCCTCGCGGGGGCATATATAAAAAGCGCTGGAGCGCGACCGGACCAACCCCCTCCCCCCCCACCGCGACCGCCag                                                                                                            >  6:8637/1‑135 (MQ=255)                  cACGGGCCGGGCCAGGCACATGGCTTGTGCTTCTCGGTGCAACCGGGCGTGGCCACGCCACCGTCGCTGGGCAATATCTACAAGGAGCTGCCGCGCGCCCTGAACATCCCGATCCCCCGCCCCGGCTACCTCCa                                                                                                             >  8:215160/1‑134 (MQ=255)                      ggccgggccTGGCACATGGCTTGGGGTTCTCGGTGCAACCGTGCGTGGCCACGCCACCGTCGCTGGTCAATATCTACAAGGAGCTGCAGCGCGACCTGAACATCCCGATCCCCAGCCACGGCTACCTCCAGAGCt                                                                                                        <  1:99087/135‑1 (MQ=255)                      ggccgggccAGGCACATGGCTTGTGCTTCTCGGTGCAACCGGGCGTGGCCCCGCCACCGTCGCGGGTGAATATATTCAACGAGCTGGAGAGCCACCGGAACATCCCGATCCCCACCCACGGCTACCACCAGAGcgg                                                                                                       >  8:57268/1‑134 (MQ=255)                         cgggccAGGCACATGGCTTGTGCTTCTCGGTGCAACCGGGCGTGGCCCCGCCCCCCTCGCGGGGGCAAATCTACAAAGAGCGGCCGCGCGCCCGGCAACCCCCCCTCCCCCCCCACCGCGACCACCACCGCTggg                                                                                                     >  6:19377/1‑135 (MQ=255)                          gggccAGGCACATGGCTTGTGCTTCTCGGTGCAACCGGGCGTGGCCACGCCACCGTCGCGGGTCAATATCTACAAGGAGCTGCAGCGCGCCCGGAACATCCCCCTCCCCAGCCACGGCTACCACCAGAGCTGGGc                                                                                                    >  5:128258/1‑135 (MQ=255)                               aGGCACATGGCTTGTGCTTCTCGGTGCAACCGGGCGTGGCCACGCCACCGTCGCTGGTCAATATCTACAAGGAGCTGCAGCGCGACCTGAACATCCCGATCCCCAGCCACGGCTACCTCCAGAGCTGGGCGGAAc                                                                                               <  2:85203/135‑1 (MQ=255)                               aGGCACATGGCTTGTGCTTCTCGGTGCAACCGGGCGTGGCCACGCCACCGTCGCTGGGCAATATCTTCAAGGCGCTGCCGCGCGCCCGGCACATCCCGCTCCCCCGCCCCGGCTACCTCCAGAGCTGGGCGGaa                                                                                                >  6:291450/1‑134 (MQ=255)                               aGGAACATGGCTTGTGCTTCTCGGTGCAACCGGGCGTGGCCACGCCACCGTCGCTGGTCAATATCTACAAGGAGCTGCAGCGCGACCTGAACATCCCGATCCCCAGCCACGGCTACCTCCAGAGCTGGGCGGAAc                                                                                               <  2:285930/135‑1 (MQ=255)                                  cacaTGGCTTGTGCTTCGCGGTGCAACCGGGCGGGGCCACGCCACCGTCGCTGGTCAATATCTACAAGGAGCTGCAGCGCGACCTGAACATCCCGATCCCCAGCCACGGCTACCTCCAGAGCTGGGCGGAACAAg                                                                                            <  1:58559/135‑1 (MQ=255)                                         cTTGTGGTTCTCGGTGCACCGGGGCGTGGCCAGGCCACCGTCGATGGTCAATATCTACAAGGAGCTGCAGCGCGACCTGAACATCCCGATCCCCAGCCACGGCTACCTCCAGAGCTGGGCGGAACAAGGCGTTCTg                                                                                    <  3:34616/136‑1 (MQ=255)                                         cTTGTGCTTCTCGGTGCAACCGGGCGTGGCCCCGCCCCCGTCGCTGGGCAATATCTACAAGGAGCTGCCGCGCGACCGGTACATCCCGCTCCCCCGCCCCGGCGACCTCCAGAGCTGGGCGGGACAAGGCGTTCt                                                                                     >  5:167857/1‑135 (MQ=255)                                         cTTGTGCTTCTCGGTGCAACCGGGCGTGGCCACGCCACCGTCGCTGGTCAATATCTACAAGGAGCTGCAGCGCGACCTGAAAATCCCGATCCCCCGCCACGGCTACCTCCCGAGCTGGGCGGAACAAGGCGTTCTg                                                                                    >  4:322403/1‑136 (MQ=255)                                         cTTGTGCTTCTCGGTGCAACCGGGCGTGGCCACGCCACCGTCGCGGGTCAATATCTACAAGGAGCTGCAGCGCGCCCTGAACATCCCGCTCCCCACCCACCGCGACCTCCAGCGCGGGGCGGAACCAGGCGTTCTg                                                                                    >  4:288465/1‑136 (MQ=255)                                          tggtgCTGCTCGGTGCAACCGGGCGTGGCCACGCCACCGTCGCTGGTCAATATCTCCAAGGAGCTGCAGCGCGACCTGAACATCCCGATCCCCAGCCACGGCTACCTCCAGAGCTGGGCGGAACAAGGCGTTCTg                                                                                    <  1:288853/133‑1 (MQ=255)                                               cTTCTCGGTGCAACCGGGCGTGGCCACGCCCCCGTCGCTGGTCAATATCTACAAGGCGCTGCAGCGCGACCTGCACCTCCCGATCCCCCGCCACGGCTACCTCCAGAGCTGGGCGGAACAAGGCGTTCTGTTGcac                                                                              >  4:128275/1‑134 (MQ=255)                                                ttCTCGGTGCAACCGGGCGTGGCCACGCCACCGTCGCTGGTCAATATCTACAAGGAGCTGCAGCGCGCCCTGAACATCCCGATACCCAGCCACGGCTACCTCCAGAGCTGGGCGGAACAAGGCGTTCTGTTGCTc                                                                              >  8:93843/1‑135 (MQ=255)                                                     ggTGCAACCGGGCGTGGCCCCGCCACCGTCGCTGGTCAATATCTAAAAGGAGCTGCAGCGCGCCCTGAACATCCCGATCCCCCGCCCCGGCGACCTCCAGAGCGGGGCGGCACAAGGCGTTCTGTTGCTCAACACg                                                                        >  7:246731/1‑136 (MQ=255)                                                     ggTGCAACCGGGCGTGGCCACGCCCCCGTCGCTGGTCAATATCTACAAGGAGCTGCAGCGCGACCTGAACATCCCGATCCCCAGCCACGGCTACCGCCAGAGCGGGGCGGAACAAGGCGTTCTGTTGCTCAacac                                                                         >  2:245139/1‑135 (MQ=255)                                                                 cGTGCCCCCGCAACCGTCGCTGTTAAATATCTACAAGGAGCTGCAGCGCGACCTGAACATCCCGATCCCCACCCACGGCTACCTCCAGAGCTGGGCGTAACACGGCGTTCTGTTGCTAAACACGACCATGACCGTc                                                            <  4:151141/136‑1 (MQ=255)                                                                       cccgccacCGTCGCGGGTGAATATGTTCAAGGAGCTGCAGCGCGGCCGGAACACCCCCCTCCCCCACCCCGGCGACCACCCGACCGGGGCGGCACAAGGCGGTCGTGTGGGCAACACGCCCCAGACCGTCGAgcg                                                       >  2:7226/3‑135 (MQ=255)                                                                       cacgccacCGTCGCGGGTCAATATATAAAATGAGCAGCAGCGCGCCCAGAAAATCCCGATCCCCCCCCACGGGGAACTCCAGAGCTGGGCGGAACAAAGCCGTCTCTTGCGCAACACGACCAATACCCTcgcgcg                                                       >  2:208786/1‑135 (MQ=255)                                                                             acCGTCGCTGGTCAATATCTACAAGGAGCTGCAGCGCGACCTGAACATCCCGATCCCCAGCCACGGCTACCTCCAGAGCTGGGCGGAACAAg                                                                                            >  8:358098/1‑92 (MQ=255)                                                                             acCGTCGCTGGTCAATATCTACAAGGAGCTGCAGCGCGACCTGAACATCCCGATCCCCAGCCACGGCTACCTCCAGAGCTGGGCGGAACAAg                                                                                            <  7:358098/92‑1 (MQ=255)                                                                             acCGTCGCTGGTCAATATCTACAAGGAGCAGCCGCGGGACCTGAACAACCCGATACCCAGCCCCGGCTACCTCCAGAGCTGGGCGGAACACTGCGTTCTGTTGGTCAACACGACCATGACCGTCGCGCGCCCCac                                                 >  5:123180/1‑134 (MQ=255)                                                                             acCGTCGCTGGTCAAGATCGACAAGGAGCTGCAGCGCGACCTGAACATCCCGATCCCCAGCCACGGCTACCTCCAGCGCTGGGCGGAACGAGGCGTTCTGGTGCTCAACACGCCCATGACCGTCGCGCGCGCCaa                                                 >  6:153304/1‑135 (MQ=255)                                                                                  cGCTGGTCAATATCTACAAGGAGCTGCAGCGCGCCCTGAACATCCCGATCCCCAGCCACGGCTACCTCCAGAGCTGGGCGGAACAAGGCGTTCTGTTGCTCAACACGACCATGACCGTCGAGCGCGCCAATTCcg                                            >  6:268336/1‑135 (MQ=255)                                                                                   gCTGGTCAATATCTACAAGCAGCTGCAGCGCGACCTGAACATCCCGAACCCCAGCCACGGCTACCTCCAGCGCTGGGCGGAACAAGGCGTTCTGTTGCTCAACACGACCATGACCGTCGAGCGCGCCAATGCcgcg                                          <  7:79357/136‑1 (MQ=255)                                                                                    cTGGTCAATATCAACAAGGAGCGCCAGCGCGACCTCAACAGCCCGATCCCCAGCCCCGCCTACCACCAGAGCTGGGCGGAACAAGGGGTTGTGTTGCTCAACACGACCATGACCGTCGAGCGCGCCAATGCcgcg                                          <  5:268336/135‑1 (MQ=255)                                                                                       gTCAATATCTACAAGGATCTGCAGCGCGCCCTGAACATCCCGCTCCCCAGCCACGGCTACCTCCAGAGCTGGGCGGCCCAAGGCGGTCTGTTGCTCACCACGTCCATGACCGTCGAGCGCGCCCATGCCGCGTCg                                       >  2:329065/1‑135 (MQ=255)                                                                                                  cAAGGAGCGGCAGCGCGACCTGAACATCCCGAGCCCCAGCCACGGCTACCTCCAGAGCTGGGCGTAACAAGGCGTTCTGTTGCTCAACACGACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCaa                               <  3:322403/132‑1 (MQ=255)                                                                                                    aGGAGCTGCAGCGCGACCTGAACATCCCGATCCCCAGCCACGGCTACCTCCAGAGCTGGGCGGAACAAGGCGTTCTGTTGCTCAACACGACCATGACCGTCGATCGCGCCAATGCCGCGTCGCATGCCAAGAAAgg                         >  8:389717/1‑136 (MQ=255)                                                                                                    aGGAGCTGCAGCGCGACCTGAACATCCCGATCCCCAGCCACGGCTACCTCCAGAGCTGGGCGGAACAAGGCGTTCTGTTGCTCAACACGACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCAGGCCAAGAAAgg                         >  4:100779/1‑136 (MQ=255)                                                                                                       cgCTGCAGCGCGACCTGAACATCCCGATCCCCCGCCACGGCTACCTCCAGAGCTGGGCGGAACAAGGCGTTCTGTTGCTCAACACGACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTg                      <  8:44113/135‑1 (MQ=255)                                                                                                         atGCAGCGCGACCTGAACATCCCGATCCCCAGCCACGGCTACCTCCAGAGCTGGGCGGAACAAGGCGTTCTGTTGCTCAACACGACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTgg                     >  6:240248/2‑135 (MQ=255)                                                                                                            cAGCGCGACCTGAACATCCCGATCCCCAGCCACGGCTACCTCCAGAGCTGGGCGGAACAAGGCGTTCTGTTGCTCAACACGACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGaa                  >  5:78371/1‑135 (MQ=255)                                                                                                                   tcgtGAAAATCCCGAGCCCCCGCCACGGGTACCGCACGAGGTGGGCGGTACTAGGGGTTGTTTTGCTCACGACGTCTATGACCGTCGAGCTCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTag          <  4:333322/133‑2 (MQ=255)                                                                                                                   aCCTGAACAGCCCGATCCCCCGCCACGGCTACCTCCAGAGCTGGGCGGAACACGGCGTGCTGTTGCTCAACACGACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTAc          <  3:91452/136‑1 (MQ=255)                                                                                                                        aaCCTCCCGACCCCCCGCCAGGGCTACCGCCCGCGGGGGGCGGAACACGGGGTGTGGTGGCCCCCCACGACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCGACa     <  7:252346/136‑1 (MQ=255)                                                                                                                         acacCCCGACCCCCAGCCCCGGCTACCCCCAGCGCGGGGCGGAACCCGGCGGGCTGTGGCCCAACACGACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCGa       <  1:376554/133‑1 (MQ=255)                                                                                                                         aCATCCCGATCCCCAGCCACGGCTACCTCCAGAGCTGGGCGGAACAAGGCGTTCTGTTGCTCAACACGACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCGACa     <  1:245139/135‑1 (MQ=255)                                                                                                                            tCCCGATCCCCAGCCACGGCTACCTCCAGAGCTGGGCGGAACAAGGCGTTCTGTTGCTCAACACGACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCGACAGGa  >  6:260445/1‑135 (MQ=255)                                                                                                                                       |                                                                                                                            TCGGTCAGGACCCGTACCACGGGCCGGGCCAGGCACATGGCTTGTGCTTCTCGGTGCAACCGGGCGTGGCCACGCCACCGTCGCTGGTCAATATCTACAAGGAGCTGCAGCGCGACCTGAACATCCCGATCCCCAGCCACGGCTACCTCCAGAGCTGGGCGGAACAAGGCGTTCTGTTGCTCAACACGACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCGACAGGA  >  ACB122/1605459‑1605716

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