Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB143	5,693,183	T→G	58.2%	R238R (CGT→CGG)	pyrC' →	Dihydroorotase‑like protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB143	5,693,183	0	T	G	58.2%	‑5.5 / 14.6	43	R238R (CGT→CGG)	pyrC'	Dihydroorotase‑like protein
Reads supporting (aligned to +/- strand):  ref base T (8/10);  new base G (4/21);  total (12/31)
Fisher's exact test for biased strand distribution p-value = 8.25e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.15e-01

CCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGCAGCCGCCGCTGCG  >  ACB143/5693056‑5693317                                                                                                                                |                                                                                                                                       ccATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCCGCTCACCAGCGCCCCTGGcgcgcg                                                                                                                                >  6:291536/1‑136 (MQ=255)  cATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCACCGCCCGGGGCGCGCgg                                                                                                                               >  2:93506/1‑136 (MQ=255)        cAGCTTCCTCGGCCGGCCGGGCATCCCGGAGAGCGCCTAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATc                                                                                                                         <  3:91979/136‑1 (MQ=255)         aGCTTCCTCGGCCGGCGGGGCATCCCGGAGAGCGCCTAAACGGTGGCCCTGGCGCGCAACCTGCTGCTGGGGTAACAGAGCGGGGTGCGTGCGCACTTCAGCCAGATCACCCGCGCCCGGGGCGCGCGGCTGATCg                                                                                                                        <  1:106926/136‑1 (MQ=255)              ccTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCACCCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGGGGCCGTGCGCCACAACACCCGCCCCCGGGGCGCGCGGCTGGTCgccccg                                                                                                                   >  5:363770/1‑134 (MQ=255)              ccTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCCCTTCAGCCCGATCACCAGCGCCCGGGGCGCGCGGCTGATCgccccg                                                                                                                   >  8:224253/1‑134 (MQ=255)               cTCGGCCTGCCGGGCATCCCGGCGAGCGCCGAACCCGGGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCCGCGCCCGTGGCGCGCGGCTGATCgcccagg                                                                                                                  <  8:95572/136‑1 (MQ=255)               cTCGGCCTGCCGGGCATCCCGGAGAGCGCCTAAACCGGGGCCCTGGCGCGCAACCTGCTGCGGGTGGAACAGCGCGGGGGGCGTGCGCACTTCAGCCAGATCACCCGCGCCCGGGGCGCGCGGCTGATCgcccagg                                                                                                                  <  8:34055/136‑1 (MQ=255)                  ggCCTGCCGGGCATCCCGGAGAGCGCCGAACCCGTGGCCCTGGCGCGCAACCTGCGGCTGGTGGAACAGAGCGGCGGGCGGGCGCACTTCAGCCAGATCACCAGCGCCCGGGGCGCGCGGCTGATCgcccaggc                                                                                                                 <  5:348343/134‑1 (MQ=255)                          gggCATCCCGGAGAGCGCCTAAACCGTGGCCGGGGCGCGCAACCTGGGTCTGGGGTAACTGAGGGGGGTGCGTGCGCACTTCAGCCAGATCACCAGGGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTg                                                                                                       <  2:39050/136‑1 (MQ=255)                                        gcgcCTAACCCGGGGCCCGGGCGCGCAACCTGCTGCGGGTGCAACAGAGCGGCGTGCGGGCGCACTTCAGCCAGATCCCCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAAc                                                                                         <  3:368897/136‑1 (MQ=255)                                        ccgcCTAAACCTGGGCCCTGGCGCGCAACCTGCTGCTGGGGGAACAGAGCGGCGTGCGTGGGCACTTCTGCCAGATCACCAGCGCCCGTGGGGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAAc                                                                                         <  7:316547/135‑1 (MQ=255)                                            cGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGGGCACTTCCGCCAGCTCACCCGCGCCCCTGGCCCGCGGCGGATCGCCCCCGCCCAGGCCGTGGGCCAGCCGGTAAccgcc                                                                                     >  5:165127/1‑136 (MQ=255)                                             gAAACCGGGGCCCGGGCGCGAAACCTGCTGCGGGTGCAACCGCGCGGCGTGGGGGCGCACTACCCCAACAACACCCGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAAccgccg                                                                                    <  2:110177/136‑1 (MQ=255)                                              aaaCCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGGGCGCCCTTCAGCCCCATAACCCCCCCCCCTGGCGCCCGGCGGCTAGCCCCCGCCCCCGCCGTGGGCCTGCCGGGTAccgccc                                                                                    >  1:183385/1‑134 (MQ=255)                                               aaCCGTGGCCCTGGCGCGAAACCGGGTGGTGGTGCAACAGCGCGGCGGGGGTGCGCACTTCAGCCAGAACACCCGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAAccgccg                                                                                    <  2:91802/134‑1 (MQ=255)                                                  cggggCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGACCACCCGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTa                                                                               <  6:204257/133‑1 (MQ=255)                                                  cGTGGCCGTGCCGCGCAACCTGCTGCTGGTGGAACAGGGCGGCGGGCGGGCGCACTTCACCCAGACCACCCGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTa                                                                               <  2:224678/136‑1 (MQ=255)                                                                    ccTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGGACTTCGGCCCGATCACCACCGCCCCTGGCGCGCGGCTGCTAGCCCAGGCCCCGGCCGTGGGGCGGCCGGTAACCCCCGATGTAGCGCTGTACCAGCTGATc                                                             >  6:84757/1‑136 (MQ=255)                                                                    ccTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCCCTTCAGCCAGATCACCCGCGCCCGTGGCGCGCGGCTGATCGCCCCCGCCCACGCCGTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGAta                                                             >  8:82526/1‑135 (MQ=255)                                                                      ggctgcggGGGCAACAGAGCGGCGGGGGGGCGAAGTTCAGCCAGATCCCCCGCGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCt                                                           <  6:125962/135‑1 (MQ=255)                                                                              ggggAACAGCGCGGCGGGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACg                                                   <  7:82526/134‑1 (MQ=255)                                                                                 gAACAGAGCGGCGTGCGTGCGCACTTAAGCCAGATAACCAGCGCCCGGGGCGCGCGGGTGCTCGCCCCCGCCCCCGCCGTGGGGCAGCCGGGAACCGCCGATGTAGCGCTGTACCCGCTGATACTCCCGGACCAg                                                 >  2:278662/1‑135 (MQ=255)                                                                                      gcgcggcgTGCGTGCGCACTTCCGCCAGATCACCCGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTg                                           <  2:24961/134‑1 (MQ=255)                                                                                        ggggcgTGCGTGCGCACTACTCCCCGACCCCCCGCGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTgcg                                         <  5:331167/134‑1 (MQ=255)                                                                                        gcggcgGGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTgcg                                         <  4:256117/136‑1 (MQ=255)                                                                                         cggcgTGCGGGCGCACTTCCGCCAGATCACCAGCTCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTgcgc                                        <  5:154436/136‑1 (MQ=255)                                                                                             gtgctggcgcACTTCCCCCACATCCCCACCGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGt                                    <  3:182321/130‑1 (MQ=255)                                                                                             gtgcgtgcGCACTTCAGCCAGATCACCCGCGCCCGGGGGGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGt                                    <  5:291858/136‑1 (MQ=255)                                                                                             gggcgtgcgCACTTCAGCCCGAACACCAGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGt                                    <  2:299108/134‑1 (MQ=255)                                                                                             gggcgtgcgCACTTCACCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGt                                    <  5:167979/134‑1 (MQ=255)                                                                                             gggcgtgcgCACTCCCCCCAGAGCACCCGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGt                                    <  4:178412/134‑1 (MQ=255)                                                                                             gggcgggcgCCCTGCAGCCAGATCACCCGCGCCCGTGGCGCGCGGCTGAGCGCCCATGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGt                                    <  8:294399/134‑1 (MQ=255)                                                                                             gggcgggcgCACTGCCGCCTGACCCCCCGCGCGCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGt                                    <  2:146467/134‑1 (MQ=255)                                                                                               gcgtgcgCACTTCAGCCAGATAACCAGCGCCCGTGGCGCGCGGCTGATAGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGGtc                                  >  1:323148/1‑136 (MQ=255)                                                                                               gcgggcgAACTTCAGCCATATACCCTGCGCCCGGGGGGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTtc                                  <  8:156288/136‑1 (MQ=255)                                                                                                         ttCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTg                         >  4:154723/1‑135 (MQ=255)                                                                                                           cAGCCAGAGCCCCCGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTa                       <  2:183385/135‑1 (MQ=255)                                                                                                                    cACCCGGGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGCAg             <  8:275196/136‑1 (MQ=255)                                                                                                                     ccccgcgcCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGCAgc            <  3:209166/132‑1 (MQ=255)                                                                                                                       ccgggcCCGGGGGGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGCAgccg          <  4:285292/132‑1 (MQ=255)                                                                                                                       ccccccccGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGCAgccg          <  5:316258/131‑1 (MQ=255)                                                                                                                               gTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCCCGTGCAGCCGCCGCTgcg  >  7:209940/1‑136 (MQ=255)                                                                                                                                |                                                                                                                                       CCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGCAGCCGCCGCTGCG  >  ACB143/5693056‑5693317

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