Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB143	5,693,175	A→C	48.0%	S236R (AGC→CGC)	pyrC' →	Dihydroorotase‑like protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB143	5,693,175	0	A	C	48.0%	3.6 / 11.2	41	S236R (AGC→CGC)	pyrC'	Dihydroorotase‑like protein
Reads supporting (aligned to +/- strand):  ref base A (8/12);  new base C (4/16);  total (12/29)
Fisher's exact test for biased strand distribution p-value = 3.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGCCCACGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGCAG  >  ACB143/5693041‑5693306                                                                                                                                       |                                                                                                                                    tGGCCCACGAAGGTGCCATGCCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTCCTGCTCGTGGACCAGAGCGGCGTGGCTGCGCACATCCGCCAGATCTCCAg                                                                                                                                    >  1:360960/1‑136 (MQ=255)    cccACGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCTAAACCGGGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGTGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAgcgc                                                                                                                                 <  2:291165/136‑1 (MQ=255)                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)                                                                                                                                       |                                                                                                                                    TGGCCCACGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGCAG  >  ACB143/5693041‑5693306

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