Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1435,693,1750AC51.6% 4.6 / 19.1 31S236R (AGC→CGC) pyrC'Dihydroorotase‑like protein
Reads supporting (aligned to +/- strand):  ref base A (10/5);  new base C (3/13);  total (13/18)
Fisher's exact test for biased strand distribution p-value = 1.13e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.79e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

AGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGCAGCC  >  ACB143/5693051‑5693308
                                                                                                                            |                                                                                                                                     
aGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCCGCCAGATCACCCGCGCCCGTGgc                                                                                                                            <  2:205965/136‑1 (MQ=255)
aGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCAGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCACCCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCCCGTCCGCCAGATCCCCCGCGCCCGGGgc                                                                                                                            >  7:122020/1‑136 (MQ=255)
           ccAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCCGCGCCCGTGGCGCGCGGCAGaa                                                                                                                 >  1:279739/1‑135 (MQ=255)
           ccAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACATCAGCCAGATCACCCGCGCCCGTGGCGCGCGGCGGa                                                                                                                  >  1:274908/1‑135 (MQ=255)
             aGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGGGCGCACTTCAGCCAGATCCCCCGCGCCCGTGGCGCGCGGCTGATCg                                                                                                               <  1:289857/136‑1 (MQ=255)
                ttCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCGGATCgccc                                                                                                            >  2:289857/1‑136 (MQ=255)
                  ccTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCCCGTGAGCCAGATCACCAGCGCCCGTGGCGCGCGGCGGATCgccccg                                                                                                          >  8:46333/1‑134 (MQ=255)
                   cTCGGCCTGCGGGGCATCCCGGAGAGCGCCGAAACCGGGGCCCTGGCGCGCACCCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCCGCCAGATCCCCCGCGCCCGTGGCGCGCGGCTGATCgcccagg                                                                                                         <  4:30539/136‑1 (MQ=255)
                                            gcgcCTAAACGGTGGCCCTGGCGCGCAACCTGGTGGTGGTGTAACAGAGCGGGGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAAc                                                                                <  4:170997/136‑1 (MQ=255)
                                            gcgcCCAAACGGGGGCCCTGCGGCGCAACGGGCGGCGGGGGCAACCGTGCGGCGTGCGGGCGCACTGCCGCCACAACCCCCGCGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAAc                                                                                <  2:333915/136‑1 (MQ=255)
                                                 gAAACCGGGGCCCTGGCGCGCAACCTGCTGCTGGTGTAACAGAGCGGCGGGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAAccgccg                                                                           <  7:211677/136‑1 (MQ=255)
                                                          gCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCCGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCCGGCCGTGGGCCTGCCGGGAACCGCCGATGTAgcgc                                                                  >  2:56982/1‑136 (MQ=255)
                                                           cccTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGGGTGCGCACTTGAGCCAGATCACCAGCGCCCGTGGCGCGCGGCGGGTCGCCCACGCCCAGGCCTTGGGGCTGCCGGGTACCGCCGAGGTTGCGCt                                                                 >  8:96401/1‑136 (MQ=255)
                                                                  gcgcAACCTGCTGCGGGGGGACCAGAGGGGCGTGCGGGCGCACTTCCGCCAGATCACCAGCGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAg                                                          <  3:91818/136‑1 (MQ=255)
                                                                         ctgctgctgGTGGAACTGAGCGGCGGGCGTGCGCACTTCAGCCAGATCACCCGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATcc                                                   <  4:341459/136‑1 (MQ=255)
                                                                                            ggggcgTGGGGGCGCCCTGCAGCCAGAGCCCCCGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTgcg                                <  4:30817/134‑1 (MQ=255)
                                                                                            gcggcgGGCGTGCGCACTACAGCCAGATCCCCCGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTgcg                                <  5:148582/136‑1 (MQ=255)
                                                                                            gcgcggtgcgggcGCACTGCCCCAAGACCACCCGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTgcg                                <  1:361759/131‑1 (MQ=255)
                                                                                                 gtgcgtgcGCACTTCAGCCATATCACCCGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGt                           <  4:306507/136‑1 (MQ=255)
                                                                                                 gtgcgggcGCACTTCAGCCAGATCACCCGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGt                           <  5:115881/136‑1 (MQ=255)
                                                                                                 gggcgtgccCACTTCACCCATAACACCCGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGt                           <  6:17314/134‑1 (MQ=255)
                                                                                                 gggcgggcgCACTTCAGCCCGATCCCCCGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGt                           <  8:113780/134‑1 (MQ=255)
                                                                                                   gcgtgcgCACTTCAGCCTGAGAACCAGCGCCCTTGGCGCGCGGCTGATCGCCCACGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTtc                         <  8:181755/136‑1 (MQ=255)
                                                                                                      tgcgCACTTCAGCCATATCACCCGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTcc                      <  1:272321/136‑1 (MQ=255)
                                                                                                                 gacAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCa           <  1:319463/134‑1 (MQ=255)
                                                                                                                     gATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGCTGTAGCGCAGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTg       >  7:145979/1‑136 (MQ=255)
                                                                                                                     gATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTg       >  1:136035/1‑136 (MQ=255)
                                                                                                                        cACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGCAg    >  4:186730/1‑136 (MQ=255)
                                                                                                                        cACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGCAg    >  3:39155/1‑136 (MQ=255)
                                                                                                                        cACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGCAg    >  1:112613/1‑136 (MQ=255)
                                                                                                                          ccAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGCAgcc  >  2:198289/1‑136 (MQ=255)
                                                                                                                            |                                                                                                                                     
AGGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGCAGCC  >  ACB143/5693051‑5693308

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