Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	4,863,130	C→T	9.0%	N305N (AAC→AAT)	pucD →	putative xanthine dehydrogenase subunit D

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	4,863,130	0	C	T	9.0%	129.2 / 3.4	44	N305N (AAC→AAT)	pucD	putative xanthine dehydrogenase subunit D
Reads supporting (aligned to +/- strand):  ref base C (14/26);  new base T (2/2);  total (16/28)
Fisher's exact test for biased strand distribution p-value = 6.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCGCGGTAATTGCGCGCCTGACCGGCCAGCCAACCAAGATGCGCCTGCCGCGGGTCGAAGACATGACCATGACCGGCAAGCGCCACCCGTTCTACGTCGAGTATGACGTGGGCTTCGACGACAACGGCCGGCTGCACGGTATCAATTTCGACCTGGCTGGCAACTGTGGCTATTCGCCCGACCTGTCCGGTTCGATCGTCGACCGTGCCATGTTCCACTCCGACAACGCTTACTACCTGGGCGATGCCACCGTCCACGG  >  ACB122/4863005‑4863264                                                                                                                              |                                                                                                                                       tgCGCGGTAATTGCGCGCCTGACCGGCCAGCCAACCAAGATGCGCCTGCCGCGGGTCGAAGACATGACCATGACCGGCAAGCGCCACCCGTTCTACGTCGAGTATGACGTGGGCTTCGACGACAACGGCCGGCt                                                                                                                                <  6:106598/134‑1 (MQ=255)  gcgcGGTAATTGCGCGCCTGACCGGCCAGCCAACCAAGATGCGCCTGCCGCGGGTCGAAGACATGACCATCACCGGCAAGCGCCACCCGGTCTACGTCGAGTATGACGTGGGCTTCGACGACAACGGCCGGCTGc                                                                                                                              <  8:32447/135‑1 (MQ=255)      ggTAATTGCGCGCCTGACCGGCCAGCCAACCAAGCTGCGCCTGCCGCGGGTCGAAGACATGCCCATGACCGGCAAGCGCCACCCGTTCTCCGTCGAGTATGACGTGGGCTTCGACGACCACGGCCGGCAGCACg                                                                                                                           >  4:244940/1‑134 (MQ=255)      ggTAATTGCGCGCCTGACCGGCCAGCCAACCAAGATGCGCCTGCCGCGGGTCGAAGACATGACCATGACCGGCAAGCGCCACCCGTTCTACGTCGAGTATGACGTGGGCTTCGACGACAACGGCCGGCt                                                                                                                                >  4:210414/1‑129 (MQ=255)      ggTAATTGCGCGCCTGACCGGCCAGCCAACCAAGATGCGCCTGCCGCGGGTCGAAGACATGACCATGACCGGCAAGCGCCACCCGTTCTACGTCGAGTATGACGTGGGCTTCGACGACAACGGCCGGCt                                                                                                                                <  3:210414/129‑1 (MQ=255)               gcgcCTGACCGGCCAGCCAACCAAGATGCGCCTGCCGCGGGTCGAAGACATGACCATGACCGGCAAGCGCCACCCGTTCTACGTCGAGTATGACGTGGGCTTCGACGACAACGGCCGGCTGCACGGTATCAAttt                                                                                                                 <  3:121676/135‑1 (MQ=255)                   cTGACCGGCCAGCCAACCAAGATGCGCCTGCCGCGGGTCGAAGACATGACCATGACCGGCAAGCGCCACCCGTTCTACGTCGAGTATGACGTGGGCTTCGACGACAATGGCCGGCTGCACGGTATCAATTTCGAc                                                                                                             >  3:190857/1‑135 (MQ=255)                   cTGACCGGCCAGCCAACCAAGATGCGCCTGCCGCGGGTCGAAGACATGACCATGACCGGCAAGCGCCACCCGTTCTACGTCGAGGATGACGTGGGCTTCGACGACAACGGCCGGCTGCTCGGTATCAATTTCGAc                                                                                                             >  6:41351/1‑135 (MQ=255)                     gACCGGCCAGCCAACCAAGATGCGCCTGCCGCGGGTCGAAGACATGACCATGACCGGCAAGCGCCACCCGTTCTACGTCGAGTATGACGTGGGCTTCGACGACAATGGCCGGCTGCACGGTATCAATTTCGACct                                                                                                           >  3:191597/1‑135 (MQ=255)                     gACCGGCCAGCCAACCAAGATGCGCCTGCCGCGGGTCGAAGACATGACCATGACCGGCAAGCGCCACCCGTTCTACGTCGAGTATGACGTGGGCTTCGACGACAACGGCCGGCTGCACGGTATCAATTTCGACct                                                                                                           >  3:130084/1‑135 (MQ=255)                           ccagccaCCCACGAGGCGCCTGCGGCGTGTCTCCGACATTCCCATGACCTGCAAGCGCCACCCTTGATACGTCGGTTATGCCGTGGGCTGCGACGCCACCGGCCGGCTGCCCGGTATCAATTTCGACctggctgg                                                                                                     <  4:3674/135‑1 (MQ=255)                                    caaGATGCGCCTGCCGCGGGTCGAAGACATGACCATGACCGGCAAGCGCCACCCGTTCTACGTCGAGTATGACGTGGGCTTCGACGACAATGGCCGGCTGCACGGTATCAATTTCGACCTGGCTGGCAACTGTgg                                                                                            <  4:190857/135‑1 (MQ=255)                                    caaGATGCGCCTGCCGCGGGTCGAAGACATGACCATGACCGGCAAGCGCCACCCGTTCTACGTCGAGTATGACGTGGGCTTCGACGACAACGGCCGGCTGCACGGTATCAATTTCGACCTGGCTGGCAACTGTgg                                                                                            <  8:272518/135‑1 (MQ=255)                                      aGATGCGCCTGCCGCGGGTCGAAGACATGACCATGACCGGCAAGCGCCACCCGTTCTACGTCGAGTATGACGTGGGCTTCGACGACAACGGCCGGCTGCACGGTATCAATTTCGACCTGGCTGGCAACTGTGGCt                                                                                          <  3:198593/135‑1 (MQ=255)                                             ccTGCCGCGGGTCGAAGACATGACCATGACCGGCAAGCGCCACCCGTTCTACGTCGAGTATGACGTGGGCTTCGACGACAACGGCCGGCTGCACGGTATCAATTTCGACCTGGCTGGCAACTGTGGCTATTCGcc                                                                                   <  6:181221/135‑1 (MQ=255)                                                gCCGCGGGTCGAAGACATGACCATGACCGGCAAGCGCCACCCGTTCTACGTCGAGTATGACGTGGGCTTCGACGACAACGGCCGGCTGCACGGTATCAATTTCGACCTGGCTGGCAACTGTGGCTATTCGCCCGa                                                                                >  1:282939/1‑135 (MQ=255)                                                        tCGAACACATGACCATGACCGCCAAGCGCCACCCGTTCTACGGCGAGTATGACGTGGGCTTCGACGACAACGGCCGGCTGCACGGTATCAATTTCGACCTGGCTGGCAACTGTGGCTATTCGCCCGACCTGTCCg                                                                        <  1:154333/135‑1 (MQ=255)                                                                    ccatgaccGGCAAGCGCCACCCGTTCTACGTCGAGTATGACGTGGGCTTCGACGACAACGGCCGGCTGCACGGTATCAATTTCGACCTGGCTGGCAACTGTGGCTATTCGCCCGACCTGTCCGGTTCGAtcgtcg                                                            >  5:108888/1‑135 (MQ=255)                                                                    ccatgaccGGCAAGCGCCACCCGTTCTACGTCGAGTATGACGTGGGCTTCGACGACAACGGCCGGCTGCACGGTATCAATTTCGACCTGGCTGGCAACTGTGGCTATTCGCCCGACCTGTCCGGTTCGAtcgtcg                                                            >  3:100653/1‑135 (MQ=255)                                                                    ccatgaccGGCAAGCGCCACCCGTTCTACGTCGAGTATGACGTGGGCTTCGACGACAACGGCCGGCTGCACGGTATCAATTTCGACCTGGCTGGCAACTGTGGCTATTCGCCCGACCTGTCCGGTTCGAtcgtcg                                                            <  2:170035/135‑1 (MQ=255)                                                                       tgaccGGCAAGCGCCACCCGTTCTACGTCGAGTATGACGTGGGCTTCGACGACAACGGCCGGCTGCACGGTATCAATTTCGACCTGGCTGGCAACTGTGGCTATTCGCCCGACCTGTCCGGTTCGATCGTCGAcc                                                         <  8:20617/135‑1 (MQ=255)                                                                          ccGGCAAGCGCCACCCGTTCTACTTCGAGTATGACGTGGGCTTCGACGACAACGGCCGGCTGCACGGTATCAATTTCGACCTGGCTGGCAACTGTGGCTATTCGCCCGACCTGTCCGGTTCGATCGTCGACCGt                                                       <  5:254479/134‑1 (MQ=255)                                                                          ccGGCAAGCGCCACCCGTTCTACGTCGAGTATGACGTGGGCTTCGACGACAACGGCCGGCTGCACGGTATCAATTTCGACCTGGCTGGCAACTGTGGCTATTCGCCCGACCTGTCCGGTTCGATCGTCGACCGt                                                       <  7:253896/134‑1 (MQ=255)                                                                          ccGGCAAGCGCCACCCGTTCTACGTCGAGTATGACGTGGGCTTCGACGACAACGGCCGGCTGCACGGTATCAATTTCGACCTGGCTGGCAACTGTGGCTATTCGCCCGACCTGTCCGGTTCGATCGTCGACCGTg                                                      >  2:168900/1‑135 (MQ=255)                                                                                aGCGCCACCCGTTCTACGTCGAGTATGACGTGGGCTTCGACGACAACGGCCGGCTGCACGGTATCAATTTCGACCTGGCTGGCAACTGTGGCTATTCGCCCGACCTGTCCGGTTCGATCGTCGACCGTGCCATGt                                                >  5:56302/1‑135 (MQ=255)                                                                                  cgcCACCCGTTCTACGTCGAGTATGACGTGGGCTTCGACGACAACGGCCGGCTGCACGGTATCAATTTCGACCTGGCTGGCAACTGTGGCTATTCGCCCGACCTGTCCGGTTCGATCGTCGACCGTGCCATGTTc                                              <  3:110074/135‑1 (MQ=255)                                                                                        ccGTTCTACGTCGAGTATGACGTGGGCTTCGACGACAACGTCCGGCTGCACGGTATCAATTTCGACCTGGCTGGCAACTGTGGCTATTCGCCCGACCTGTCCGGTTCGATCGTCGACCGTGCCATGTTCCACTcc                                        <  2:268496/135‑1 (MQ=255)                                                                                         cGTTCTACGTCGAGTATGACGTGGGCTTCGACGAGAACGGCCGGCTGCACGGTATCAATTTCGACCTGGCTGGCAACTGTGGCTATTCGCCCGACCTGTCCGGTTCGATCGTCGACCGTGCCATGTTCCACTCCg                                       <  7:172779/135‑1 (MQ=255)                                                                                         cGTTCTACGTCGAGTATGACGTGGGCTTCGACGACAACGGCCGGCTGCAAGGTATCAATTTCGACCTGGCTGGCAACTGTGGCTATTCGCCCGACCTGTCCGGTTCGATCGTCGACCGTGCCATGTTCCACTCCg                                       <  8:244576/135‑1 (MQ=255)                                                                                          gTTCTCCGTCGAGTATGCCGGGGGCTGCGACGACAACGGCCGGCTGCACGGTATCAATTTCGACCTGGCTGGCAACTGTGGCTATTCGCCCGACCTGTCCGGTTCGATCGGCGACCGTGCCATGTTCCACTCCGa                                      <  3:244940/135‑1 (MQ=255)                                                                                          gTTCTACGTCGAGTATGACGTGGGCTTCGACGACAACGGCCGGCTGCACGGTATCAATTTCGACCTGGCTGGCAACTGTGGCTATTCGCCCGACCTGTCCGGTTCGATCGTCGACCGTGCCATGTTCCACTCCGa                                      <  6:56302/135‑1 (MQ=255)                                                                                          gTTCTACGTCGAGTATGACGTGGGCTTCGACGACAACGGCCGGCTGCACGGTATCAATTTCGACCTGGCTGGCAACTGTGGCTATTCGCCCGACCTGTCCGGTTCGATCGTCGACCGTGCCATGTTCCACTCCGa                                      <  5:100685/135‑1 (MQ=255)                                                                                                 gTCGAGTATGACGTGGGCTTCGACGACAACGGCCGGCTGCACGGTATCAATTTCGACCTGGCTGGCAACTGTGGCTATTCGCCCGACCTGTCCGGTTCGATCGTCGACCGTGCCATGTTCCACTCCGACAAc                                  <  4:80945/132‑1 (MQ=255)                                                                                                        aTGACGTGGGCTTCGACGACAACGGCCGGCTGCACGGTATCAATTTCGACCTGGCTGGCAACTGTGGCTATTCGCCCGACCTGTCCGGTTCGATCGTCGACCGTGCCATGTTCCACTCCGACAACGCTTACTAcc                        >  5:124857/1‑135 (MQ=255)                                                                                                               gggCTTCGACGACAACGGCCGGCTGCACGGTATCAATTTCGACCTGGCTGGCAACTGTGGCTATTCGCCCGACCTGTCCGGTTCGATCGTCGACCGTGCCATGTTCCACTCCGACAACGCTTACTACCTGGGCGa                 >  1:224744/1‑135 (MQ=255)                                                                                                                 gCTTCGACGACAACGGCCGGCTGCACGGTATCAATTTCGACCTGGCTGGCAACTGTGGCTATTCGCCCGACCTGTCCGGTTCGATCGTCGACCGTGCCATGTTCCACTCCGACAACGCTTACTACCTGGGCGATg               <  5:140089/135‑1 (MQ=255)                                                                                                                 gCTTCGACGACAACGGCCGGCTGCACGGTATCAATTTCGACCTGGCTGGCAACTGTGGCTATTCGCCCGACCTGTCCGGTTCGATCGTCGACCGTGCCATGTTCCACTCCGACAACGCTTACTACCTGGGCGATg               <  3:259718/135‑1 (MQ=255)                                                                                                                 gCTGCGACGACAATGGCCGGCTGAACGTTATAAATTTCGACCTGGATGGCAACTGTGGCTATTCGCCGGACCTGTCCGGTTCGAGCGTCGACCGTGCCATGTTCCACTCCTACAACGCTTACTACGTGGGCTATg               <  4:191597/135‑1 (MQ=255)                                                                                                                 gCTGCGACGACAACGGCCGGCTGCACGGTATCAATTTCGACCTGGCTGGCAACTGTGGCTATTCGCCCGACCTGTCCGGTTCGATCGTCGACCGTGCCATGATCCACTCCGACAACGCTTACTACCTGGGCGATg               <  6:254939/135‑1 (MQ=255)                                                                                                                  cTTCGACGACAACGGCCGGCTGCACGGTATCAATTTCGACCTGGCTGGCAACAGAGGCTATTCTCCCGACCTGTCCGGATCGATCGTCGACAGTTCAAGGTTCCACTCCGACAACGATTCCTACCTGGTCTATGc              >  7:224128/1‑135 (MQ=255)                                                                                                                        cgacAACGGCCGGCTGCACGGTATCAATTTCGACCTGGCTGGCAACTGTgg                                                                                            <  6:155357/51‑1 (MQ=255)                                                                                                                        cgacAACGGCCGGCTGCACGGTATCAATTTCGACCTGGCTGGCAACTGTgg                                                                                            >  5:155357/1‑51 (MQ=255)                                                                                                                        cgacAACGGCCGGCTGCACGGTATCAATTTCGACCTGGCTGGCAACTGTGGCTATTCGCCCGACCTGTCCGGTTCGATCGTCGACCGTGCCATGTTCCACTCCGACAACGc                                >  7:273079/1‑111 (MQ=255)                                                                                                                        cgacAACGGCCGGCTGCACGGTATCAATTTCGACCTGGCTGGCAACTGTGGCTATTCGCCCGACCTGTCCGGTTCGATCGTCGACCGTGCCATGTTCCACTCCGACAACGc                                <  8:273079/111‑1 (MQ=255)                                                                                                                              cggccggcTGCACGGTATCAATTTCGACCTGGCTGGCAACTGTGGCTATTCGCCCGACCTGTCCGGTTCGATCGTCGCCCGTGCCATGTTCCACTCCGACAACGCTTACTACCTGGGCGGTGCCACCGTCCACgg  >  8:269436/1‑135 (MQ=255)                                                                                                                              cggccggcTGCACGGTATCAATTTCGACCTGGCTGGCAACTGTGGCTATTCGCCCGACCTGTCCGGTTCGATCGTCGACCGTGCCATGTTc                                              >  4:360/1‑91 (MQ=255)                                                                                                                              cggccggcTGCACGGTATCAATTTCGACCTGGCTGGCAACTGTGGCTATTCGCCCGACCTGTCCGGTTCGATCGTCGACCGTGCCATGTTc                                              <  3:360/91‑1 (MQ=255)                                                                                                                              cggccggcTGCACGGTATCAATTTCGACCTGGCTGGCAACTGTGGCTATTCGCCCGACCTGTCCGGTTCGATCGTCGACCGTGCCATGTTCCACTCCGACAACGCTTACTACCTGGGCGATGCCACCGTCCACgg  >  8:116077/1‑135 (MQ=255)                                                                                                                              |                                                                                                                                       TGCGCGGTAATTGCGCGCCTGACCGGCCAGCCAACCAAGATGCGCCTGCCGCGGGTCGAAGACATGACCATGACCGGCAAGCGCCACCCGTTCTACGTCGAGTATGACGTGGGCTTCGACGACAACGGCCGGCTGCACGGTATCAATTTCGACCTGGCTGGCAACTGTGGCTATTCGCCCGACCTGTCCGGTTCGATCGTCGACCGTGCCATGTTCCACTCCGACAACGCTTACTACCTGGGCGATGCCACCGTCCACGG  >  ACB122/4863005‑4863264

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