Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	5,688,645	T→G	55.5%	R238R (CGT→CGG)	pyrC' →	Dihydroorotase‑like protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	5,688,645	0	T	G	55.5%	‑3.3 / 16.4	41	R238R (CGT→CGG)	pyrC'	Dihydroorotase‑like protein
Reads supporting (aligned to +/- strand):  ref base T (8/10);  new base G (5/18);  total (13/28)
Fisher's exact test for biased strand distribution p-value = 1.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

AGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGCAGCC  >  ACB122/5688513‑5688770                                                                                                                                     |                                                                                                                              aGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGGGCGCCCCTCCGCCCGATAACCAGCGCCCGTgg                                                                                                                             >  2:126442/1‑135 (MQ=255)         tGGCCAGCTTCCTCGGCCGGCCGGGCATCCCGGAGCGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGGGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGc                                                                                                                     <  3:382744/135‑1 (MQ=255)              aGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGGGCGCCCTTCAGCCAGATCACCCGCGCCCGGGGCGCGCGGCGGATCg                                                                                                               >  7:212660/1‑136 (MQ=255)                   ccTCGGCCTGCCGGGCATCCAGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGT‑GAACAGAGCGGCGTGGGGGGGCACGTGCGCCATAACCCCCCAGCCCCTGGCGCCCCGCGGATCgccccg                                                                                                          >  6:70234/1‑133 (MQ=255)                     tCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCGGCGGGTGGAACAGAGCGGCGGGCGTGCGCACTTCCGCCAGATCACCCGCGCCCGGGGCGCGCGGCTGATCgcccagg                                                                                                         <  3:186788/135‑1 (MQ=255)                       ggCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGGGCGCCCGTCCCCCACATCACCCCCGCCCGGGGCGCCCGGGGGGTCgccccggcc                                                                                                       >  2:246654/1‑135 (MQ=255)                             ccGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGGGGGGGCCCTGGAGCCACACCCCCACCGCCCCGGGCGCGCGGGGGCTCgcccccgcccccgcc                                                                                                 >  3:253619/1‑130 (MQ=255)                                        gcagagCGCCTAAACGGGGGCCCGGGCGCGAAACCTGGTGCTGGTGGAACCGAGCGGGGTGCGTGCGCACTTCCGCCAGATCCCCCGCGCCGGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGcc                                                                                      <  3:123163/133‑1 (MQ=255)                                             gcgcCGAACCGGGGGCCGGGGGGCGCACCGTGGGGGTGGTGGAACCGCGGGGGGGGGGGGCGCACTTCCGCCACAGCCCCCGCGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAAc                                                                                <  8:364602/136‑1 (MQ=255)                                             gcgcCGAACCGGGGGCCCTGGCGCGCAACCGGCTGCTGGTGTAACCGAGCGGCGTGGGTGCGCACTTCAGCCCTAACCCCCGCGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAAc                                                                                <  7:315396/136‑1 (MQ=255)                                               gcCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCCCTTCAGCCAGATCACCCGCCCCCCTGGCGCGCGGCTGCTCGCCCCCGCCCCGGCCGTGGGGCTGCCgggtacc                                                                               >  1:32848/1‑130 (MQ=255)                                                 cGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGGGGGGGGCCTGGAGCCACCTAACCCGCGCCCGGGGCGCGCGGCGGCTCGCCCCGGCCCCGGCCGTGGGGCGGCCGGGGAccgg                                                                             >  5:187535/1‑134 (MQ=255)                                                 cGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGGGGGGGGCCGTCAGCCAGATCAACCCCCCCCCGGGCGCGCCGCGGGACGCCCCCGCCCCGGGCGTGGGCCTGCCCGGAAcccc                                                                             >  6:359857/1‑133 (MQ=255)                                                 cGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGGGCACTTCAGCCCGAACACCCCCCCCCCTGGGGCGCGGCGGCTAGCCCCCGCCCACGCCGTGGGCCAGCCCGGAAccgc                                                                             >  6:112585/1‑135 (MQ=255)                                                   aaaCCGGGGCCCTGGCGCGCAACCTGCTGCTGGGGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGAGCCCCCGCGCCGGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAAccgccg                                                                           <  3:53954/135‑1 (MQ=255)                                                    aaCGGTGGCCGTGGCGGGCAACCGGGTGCGGGTGCACCAGAGCGGCGGGCGGGCGCAGTTCAGCCAGACCACCCGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGa                                                                          <  4:215227/135‑1 (MQ=255)                                                       gggggCCCGGGCGCGAAACCTGCGGGGGGGGCAACAGAGCGGCGGGCGGGCGCACTTCAGCCAGATCCCCAGCGCCGGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTa                                                                      <  7:257385/133‑1 (MQ=255)                                                       cggggCCCTGGCGCGCAACCTGCTGCGGGGGGAACAGAGCGGCGGGCGTGCGCACTTCCGCCACAGCCCCCGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTa                                                                      <  8:370936/133‑1 (MQ=255)                                                        gTGGCCCTGGCGCGCAACCTGCTGCGGGGGGAACAGCGCGGCGTGCGTGCGCACTTCAGCCAGATCCCCCGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTa                                                                      <  2:32848/135‑1 (MQ=255)                                                        gTGGCCCTGGCGCGCAACCGGGTGGGGGGGGAACTGAGGGGGGTGCGTGCGCACTGCTGCCAGAGCCCCGGCGCCCGTGGGGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTa                                                                      <  4:314060/135‑1 (MQ=255)                                                                     gcAACCTGGGGCTGGGGAACCAGCGCGGCGGGCGGGCGCACTTCCGCCAGATCACCCGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCt                                                        <  7:322633/136‑1 (MQ=255)                                                                     gaaaCCTGCGGCTGGGGGAACCGGGCGGGGGGCGGGGGCACTTCTGCCAGAGCACCCGGGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGc                                                         <  1:246654/133‑1 (MQ=255)                                                                         ccTGCTGCTGGTGGAACAGAGCGGCGTGCGAGCGCCCTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGCTCGCCCCGGCCCCGGCCGTGGGCCTGCCCGTAACCGCCGATGTAGCGCTGTACCAGCTGAt                                                     >  1:57948/1‑135 (MQ=255)                                                                                         cAGAGCGGCGTGCGTGCGCACATCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCCGGCCGTGGGCCTGCCCGTAACCCCCCATGTAGCGCTGTACCAGCTGCTACTCACGGCCGAGTCg                                     >  8:319947/1‑135 (MQ=255)                                                                                            cgcggcgGGCGTGCGCACTTCCGCCAGAGCACCCGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTg                                  <  1:3372/134‑1 (MQ=255)                                                                                             gcggcgTGCGGGCGCACTCCCCCCAGACCACCCGCGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTgcg                                <  7:201757/136‑1 (MQ=255)                                                                                             gcggcgGGCGGGCGCACCCCCGCCACACCCCCGCCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTgcg                                <  7:346479/136‑1 (MQ=255)                                                                                              cggcgTGCGTGCGCACTTCAGCCCCAGCCCCCGCGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTgcg                                <  2:57948/135‑1 (MQ=255)                                                                                                  gggcgggcgCACTTCCGCCATATCACCCGCGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGt                           <  7:319947/134‑1 (MQ=255)                                                                                                  gggcgggcgCACTCCAGCCAGATCCCCAGCGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGt                           <  8:79315/134‑1 (MQ=255)                                                                                                   tgggggcgcACTCCAGCCAGACCCCCCGCGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGt                           <  5:112585/130‑1 (MQ=255)                                                                                                   tgcgtgcgAACTGCAGCCAGAGCACCCGGGCCGGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGt                           <  6:125761/135‑1 (MQ=255)                                                                                                   ggcgggcgCACTCCCCCCAGAACACCCGCGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGt                           <  6:305021/134‑1 (MQ=255)                                                                                                    gcgggcgCACTTCCGCCAGAGCACCCGCGCCCGTGCGGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGtt                          <  4:127569/135‑1 (MQ=255)                                                                                                    gcgggcgCACTCCCGCCAGACCACCCGCGCCGGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTtc                         <  8:366354/136‑1 (MQ=255)                                                                                                     cgtgccAACTGCAGCCAGAGCCCCCGCGCCCGGGCCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTATCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTtc                         <  6:190352/135‑1 (MQ=255)                                                                                                     cgggcgcACTTCAGCCACTTCACCCGGGCCGGTGGGGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATGCTCACGGACGAGTCGCTGCGCGAGTtc                         <  6:304293/132‑1 (MQ=255)                                                                                                       tgcgCAGTTCAGCCAGATCACCAGCGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTcc                      <  8:84320/136‑1 (MQ=255)                                                                                                                           ccAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTtc                         <  3:312524/113‑1 (MQ=255)                                                                                                                           ccAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTtc                         >  4:312524/1‑113 (MQ=255)                                                                                                                           ccAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGCAgcc  >  7:126775/1‑136 (MQ=255)                                                                                                                                     |                                                                                                                              AGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGCAGCC  >  ACB122/5688513‑5688770

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