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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1435,693,1750AC40.7% 13.0 / 5.7 36S236R (AGC→CGC) pyrC'Dihydroorotase‑like protein
Reads supporting (aligned to +/- strand):  ref base A (10/11);  new base C (4/11);  total (14/22)
Fisher's exact test for biased strand distribution p-value = 3.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGG‑CCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGCAGCCG  >  ACB143/5693051‑5693309
                                                                                                                            |                                                                                                                                       
aGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTCGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGGGCGGGGGGCCGTGAGCCACGACACCCGCGCCCGGGgc                                                                                                                              >  5:281654/1‑136 (MQ=255)
  gTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCTAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGgcg                                                                                                                             <  7:271037/135‑1 (MQ=255)
      cATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGGGCCCTTCAGCCCGATAACCAGCGCCCCTGGCGCGCgg                                                                                                                        >  8:271312/1‑136 (MQ=255)
        tGGCCCGCTCCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGGGGCCCTGGCGCGCAACCTGCTGCGGGGGGACCAGCGCGGCGTGCGTGCGCACTTCAGCCACATCACCAGCGCCCGTGGCGCGCGGCt                                                                                                                      <  5:72742/136‑1 (MQ=255)
        gggCCTGCTTCCTCGGCCTGCCGGGCAACCCGGCGCGCCCCGAAACCGGGCCCCTGGCGCGCAACCTGCTGCGGGTGGAACTCTGCGGCGTGCGTGCGCACTTCCGCCAGATCACCCGCGCCCGGGGCGCGCGGCt                                                                                                                      <  2:178455/135‑1 (MQ=255)
            cAGCTTCCTCGGCCTGCCGGGCAGCCCGGAGAGCGCGGAAACCGTGGCCCTGCCGCGCAACCTGCTGCGGGTGGAACAGTGCGGGGTGCGGGCGCACTTCAGCCAGATCCCCAGCGCCCGGGGCGCGCGGCTGATc                                                                                                                  <  1:246852/136‑1 (MQ=255)
             cgCTGCCTGGCCCGGCCGGGCACCCCGCCGAGCGCCTAAACCGGGGCCCTGGGGCGAAACCGGCTGCGGGGGCAACATAGCGGGGGGCGGGCGCACTACCGCCCTACCACCCGCGCCCGTGGCGCGCGGCTGATCg                                                                                                                 <  4:131217/135‑1 (MQ=255)
             aGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCCCGTAAGGCAGATCACCAGCGCCCGTGGCGCGCGGCTGCTCg                                                                                                                 >  5:250253/1‑136 (MQ=255)
                  ccTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCGGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCgcccag                                                                                                            <  3:199148/136‑1 (MQ=255)
                   cTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCCGCGCCCGTGGCGCGCGGCTGATCgcccagg                                                                                                           <  5:249729/136‑1 (MQ=255)
                   cTCGGCCGGCCGGGCATCCGGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCATATCACCAGCGCCCGTGGCGCGCGGCTGATCgcccagg                                                                                                           <  1:12616/136‑1 (MQ=255)
                       gCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGGGCGGACTTCAGCCCGATCACCAGCGCCCGGGGCGCGCGGCTGATCgccccgg‑cccc                                                                                                      >  8:303671/1‑135 (MQ=255)
                                     ccGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGGGGGCCCTTCAGCCAGATAACCCCCGCCCGTGGCGCGCGGCTGATCGCCCCGG‑CCCCGGCCGTGGGCCGGc                                                                                        >  1:60833/1‑136 (MQ=255)
                                      cGGCGAGCGCCGAAACCGTGGCCCTGGCGCGAAACCTGCTGCTGGTGGAACAGAGCGGCGTGGGTGCGCACTTCAGCCAGATCACCAGCGCCGGTGGCGCGCGGCTGATCGCCCAGG‑CCCAGGCCTTGGGCCTGcc                                                                                       <  2:198054/136‑1 (MQ=255)
                                            gcgcCTAAACGGTGGCCGTGGCGCGCAACGTGGTGCTGGTGTAACAGAGCGGGGTGCGTGCGCACTTCCGCCAGATCACCCGCGCCCGTGGCGCGCGGCTGATCGCCCAGG‑CCCAGGCCTTGGGCCTGCCGGTAAc                                                                                 <  3:132783/136‑1 (MQ=255)
                                            gcgcCGAAACCGTGGCCCGGGCGCGCAACCTGCTGCTGGTGCAACAGCGCGGCGTGCGGGCGCACTTCCGCCAGACCCCACGCGCCCGGGGCGCGCGGCTGATCGCCCAGG‑CCCAGGCCTTGGGCCTGCCGGTAAc                                                                                 <  1:40490/136‑1 (MQ=255)
                                                cGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCCATTGAGCCAGATACCCCGCGCCCCGGGCGCGCGGCTGATCGCCCCGGCCCCAGGCC‑TGGGGCTGCCGGGAAccgcc                                                                             >  8:209830/1‑136 (MQ=255)
                                                 gAAACCGTGGCCCGGGCGCGCAACCTGCGGCTGGTGCAACAGAGCGGCGTGCGGGCGCACTTCAGCCAGAACACCCGCGCCGGTGGCGCGCGGCTGATCGCCCAGG‑CCCAGGCCTTGGGCCTGCCGGTAAccgccg                                                                            <  5:168995/136‑1 (MQ=255)
                                                      ggTGGCCCTGGGGGGCAACCTGCGGGTGGTGTAACAGAGGGGGGTGCTTGGGCACTTCAGCCAGATCACCAGCGCCCGTGGGGCGCGGCTGATCGCCCAGG‑CTCAGTAATTGGGCCTGCCGGTAACCGCCGATGTa                                                                       <  8:23965/135‑1 (MQ=255)
                                                                  gcgcAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCCCTTCAGCCAGCTACCCCGCGCCCGTGGCGCGCGGCTGATCGCCCCCG‑CCCCGGCCGTGGGGCTGCCGGTAACCGCCGATGTAGCGCTGTACCAg                                                           >  1:26804/1‑136 (MQ=255)
                                                                   cgcAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCCCTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCCGG‑CCCCGGCCGTGGGCCAGCCGGTTACCGCCGCTGTAGCGCTGTACCCGc                                                          >  6:59934/1‑136 (MQ=255)
                                                                         ctgcggctgGTGGAACAGAGGGGGGGGCGTGCGCACTTCAGCCAGACCACCCGCGCCCGGGGCGCGCGGCTGATCGCCCAGG‑CCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATcc                                                    <  3:289432/136‑1 (MQ=255)
                                                                                            gcggcgTGCGGGCGCACTTCCCCCAGATCCCCCGCGCCCGGGGCGCGCGGCTGATCGCCCAGG‑CCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTgcg                                 <  2:222178/136‑1 (MQ=255)
                                                                                             cggcgTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCgcccagg                                                                                                           >  8:84244/1‑62 (MQ=255)
                                                                                             cggcgTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGGGGCGCGCGGCTGATCgcccagg                                                                                                           <  7:84244/62‑1 (MQ=255)
                                                                                              ggcgTGGGGGCGCACATCAGCCCGATCACCACCGCCCCTGGGGCGCGGCGGCTAGCCCCCG‑CCCACGCCGTGGGGCGGCCGGGTACCGCCGATGTAGCGCTGTACCACCTGATACGCACCGCCCAGTCGCGGcgcg                               >  8:21910/1‑136 (MQ=255)
                                                                                               gcgtgcgtGCGCACTTCAGCCACATAACCAGCGCCCGTGGCGCGCGGCTGCTACCCCAGG‑CCCAGGCCATTGGCCTGCCGGTAACCGCCGATGTTGCGCAGTACCAGCAGATACTCACCGACCAGTAGCAgcgc                                >  6:133403/1‑134 (MQ=255)
                                                                                                 gggcgtgcgCACTTCCGCCACAACACCCGCGCCCGTGGCGCGCGGCTGATCGCCCAGG‑CCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGt                            <  1:289330/134‑1 (MQ=255)
                                                                                                 gggcgggcgCACTCCCCCCAGAGCACCCGCGCCCGGGGCGCGCGGCTGATCGCCCAGG‑CCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGt                            <  3:307708/134‑1 (MQ=255)
                                                                                                  ggggggcgcACTTCCCCCAGATCCCCCGCGCCCGTGGCGCGCGGCTGATCGCCCAGG‑CCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGtt                           <  2:60833/131‑1 (MQ=255)
                                                                                                     gtgcgCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCCGG‑CCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCAGTACCAGCTGATCCTCACGGCCGAGTCGCTGCGCGCGGtctc                        >  7:274257/1‑136 (MQ=255)
                                                                                                            cTTCAGCCACAGCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGG‑CCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTg                 <  2:242128/136‑1 (MQ=255)
                                                                                                                          ccAGCGCCCGTGGCGCGCGGCTGATCGCCCAGG‑CCCAGGCCTTGGGCCTGCCGGTAAccgccg                                                                            <  4:178433/63‑1 (MQ=255)
                                                                                                                          ccAGCGCCCGTGGCGCGCGGCTGATCGCCCAGG‑CCCAGGCCTTGGGCCTGCCGGTAAccgccg                                                                            >  3:178433/1‑63 (MQ=255)
                                                                                                                          ccAGCGCCCGTGGCGCGCGGCTGATCGCCCAGG‑CCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCCCGTGCAgcc   >  4:86109/1‑136 (MQ=255)
                                                                                                                           cAGCGCCCGTGGCGCGCGGCTGATCGCCCAGG‑CCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGCAgccg  <  7:228242/136‑1 (MQ=255)
                                                                                                                            |                                                                                                                                       
AGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGG‑CCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGCAGCCG  >  ACB143/5693051‑5693309

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