Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	5,688,637	A→C	58.0%	S236R (AGC→CGC)	pyrC' →	Dihydroorotase‑like protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	5,688,637	0	A	C	58.0%	‑3.2 / 17.1	39	S236R (AGC→CGC)	pyrC'	Dihydroorotase‑like protein
Reads supporting (aligned to +/- strand):  ref base A (6/9);  new base C (8/14);  total (14/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.81e-01

TGGCCCACGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGCAGC  >  ACB122/5688503‑5688769                                                                                                                                       |                                                                                                                                     tGGCCCACGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGGGCGGGGGCCCGTCAGCCCCATAACCa                                                                                                                                      >  4:98796/1‑135 (MQ=255)           aGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCCCCCCCGCCCGTgg                                                                                                                            >  2:56449/1‑135 (MQ=255)                   tGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGGGGGCGGCCGTCCGCCACATCACCCGCCCCCCAGGCGCCCGGc                                                                                                                    >  1:88799/1‑135 (MQ=255)                             ccTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGGGCGGGGGCCCGTCAGCCCCCTCACCCCCGCCCCGGGCGCGCGGGGGCTCgcccc                                                                                                          >  7:34190/1‑134 (MQ=255)                              cgcggccggccggcCAGCCCGCCGCCCGCCTAAACGGGGGCCCGGGCGCGCAACCTGCGGGGGGTGCAACAGCGCGGCGGGCGGGCGCACTTCAGCATGACCCCCCGCGCCCGGGGCGCGCGGCTGATCgcccag                                                                                                         <  4:87631/133‑1 (MQ=255)                              cgcggccggcGGGCCAGCCGGTAGAGCGCCGAAACCGGGGCCGGGGGGCGCACCCTGGTGGGGGGGCACCAGCGCGGCGTGGGTGCGCACTCCCGCCTCAACCCCCGCGCCCGGGGCGCGCGGCTGATCgcccag                                                                                                         <  4:109173/133‑1 (MQ=255)                                  gCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCCCTTGAGCCCGATCACCCGCGCCCGGGGCGCGCGGCGGGTCgccccggccc                                                                                                     >  8:191776/1‑135 (MQ=255)                                          ggccgcccGGCGAGCGCCGAACCGGGGCCCCGGGCGCGCACCGTGGTGCGGGGGTACCCGCGGGGCGGGCGGGCGCAGTCCCGCCAGAGCCCCCGCGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTg                                                                                             <  2:103506/130‑1 (MQ=255)                                                     gagCGCCGAAACCGTGGCCCGGGCGCGCAACCGGGTGCGGGGGGAACAGAGCGGGGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGt                                                                                  <  8:205090/135‑1 (MQ=255)                                                        cgcgGAACCCGGGGCCGTGGCGCGCAACCTGCTGGGGGTGTAACTGAGCGGGGTGCGGGCGCAGTGCAGCCAGTTCACCAGGGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAAc                                                                               <  6:140821/135‑1 (MQ=255)                                                        cgcCTAAACCGTGGCCGTGGCGCGCAACCTGGTGCTGGTGGAACATAGGGGCGTGCGTGCGCACTTCAGCATGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGt                                                                                  <  8:205089/132‑1 (MQ=255)                                                        ccccGAACCCGGGGCCGGGGGGCGCAACGTGCTGCTGGGGCAACTGCGGGGGGGGCGGGGGCACTGCCCCCACAGCCCCGGCGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAAc                                                                               <  2:78232/133‑1 (MQ=255)                                                          ccGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGGGGACCAGAGCGGCGTGCGTGCGCACTTCCGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAAccg                                                                             <  4:143303/135‑1 (MQ=255)                                                           cGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGGGGTGGGACTTGAGCCCGCGACCCCACGCCCCTGGCGCGCGGCGGGTCGCCCCGGCCCCGGGCGTGGGGCGGGCGGtttacgc                                                                            >  2:231691/1‑129 (MQ=255)                                                           cGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGGGGGGGCCCATACGCCCCATCACCCCCCCCCCGGGGGCGCCGCTGCTAGCCCACGCCCCCGCCCTGGGCCGGCCCGTGAccgc                                                                            >  8:41911/1‑135 (MQ=255)                                                           cGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGGGCGCACTTGAGCCCGATAACCCGCGCCCCTGGCGCGCCGCGGATCGCCCCCGCCCCCGCCGTGGGCCTGCCGGGAAccgc                                                                            >  2:198935/1‑135 (MQ=255)                                                           cGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGGGCGCACTTCCGCCAGATCACCCGCGCCCGTGGGGCGCGGCGGCTCGCCCCGGCCCCGGCCGGGGGCCTGGCGGTAAccgc                                                                            >  6:148736/1‑135 (MQ=255)                                                                cgggggCCCGGGCGCGCAACCTGCTGCTGGTGCAACAGAGCGGCGTGCGGGCGCACTTCAGCCAGACCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATg                                                                       <  7:199312/131‑1 (MQ=255)                                                                  ggggCCGTGGCGCTCACCCTGCGGGTGGGGCAACAGGGGGGCGGGCGGGCGCACTACCCCCAGAACCCCAGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTa                                                                     <  8:117575/133‑1 (MQ=255)                                                                  ggggCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGGGGCGGGCGTGCGCACTTCAGCCAGATCACCCGCGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTa                                                                     <  7:191776/133‑1 (MQ=255)                                                                  ggggCCCGGGCGCGAACCCTGCGGGGGGGGCAACAGAGCGGCGGGCGGGCGCACTACCGCCAGAGCACCCGCGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTa                                                                     <  8:91759/133‑1 (MQ=255)                                                                  gcggCCCGGGCCGGCAACCTCTTGGTGGTGAAACCTTGCGGGGTGCTTCCGCACTACAGCCTGATCCACCGCGCCGGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTa                                                                     <  7:238635/133‑1 (MQ=255)                                                                  gTGGCCCTGGGGCGCAACCTGCTGGTGTTGCCACTGAGCGGCGTGCGGGCGAAATTCTGCAAGATCACCAGCGCGCGTGGCGCGCGGCTGATCGCCCAGGCCCAGTCCTTGGGCCTGCCGGTAACCGCCGATGTa                                                                     <  8:87698/135‑1 (MQ=255)                                                                  gTGGCCCGGGCGCGCAACCTGTTGCTGGTGGAACAGAGCGGGTGGCGTGGTCAGTTCAGCCAGATAACGAGCGCCCGTGGAGCGCGGCTGATCGCCGAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTa                                                                     <  2:85170/135‑1 (MQ=255)                                                                      cccGGGCGCTCACCGGGCTGCGGGGGCAACAGCGCGGGGGGGGGGCGCACTCCACCCATACCCCCCGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAgcgc                                                                 <  3:253510/135‑1 (MQ=255)                                                                             gcgcAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCCGCTGATCGCCCACGCCCCGGCCTTGGGCCTGCCGGGAACCGCCGATGTAGCGCGGTACCa                                                          >  7:75481/1‑135 (MQ=255)                                                                                                      cgcggggGGCGTGCGCACTGCGGCCCGACCCCCCGGGCGGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTg                                 <  7:152514/134‑1 (MQ=255)                                                                                                        ggcggggcgtgcgAAATCCAGCCCTAACCCCTGCGCCGGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTgcg                               <  6:4533/129‑1 (MQ=255)                                                                                                        cggcggGCGGGCGCACTCCAGCCAGAGCACCCGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTgcg                               <  8:172778/135‑1 (MQ=255)                                                                                                          gcgggcggGCGCACTCCCCCCAGACCACCCGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGcgcg                             <  8:17566/135‑1 (MQ=255)                                                                                                             ggcgtgcgCACTTCTGCCACAGCACCCGCGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGt                          <  2:88799/134‑1 (MQ=255)                                                                                                             ggcgtgcgCACTTCAGCCAGAACCCCCGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGt                          <  3:60583/134‑1 (MQ=255)                                                                                                              gcgtgcgCACTTCAGCCAGATCACCCGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGtt                         <  2:15784/135‑1 (MQ=255)                                                                                                                  gcgcACTGCAGCCAGATCACCCGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTcc                     <  4:73248/135‑1 (MQ=255)                                                                                                                        ttCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGCGTGCTCCAGCCt                >  6:206565/1‑134 (MQ=255)                                                                                                                           aGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTAc            >  6:244996/1‑135 (MQ=255)                                                                                                                                gATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGc                                   <  7:209728/107‑1 (MQ=255)                                                                                                                                gATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGc                                   >  8:209728/1‑107 (MQ=255)                                                                                                                                    aCCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTATCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGCAg   >  8:76203/1‑135 (MQ=255)                                                                                                                                     ccAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGCAgc  >  6:142660/1‑135 (MQ=255)                                                                                                                                       |                                                                                                                                     TGGCCCACGAAGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGCAGC  >  ACB122/5688503‑5688769

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