Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1225,688,6370AC54.8% 2.1 / 16.9 31S236R (AGC→CGC) pyrC'Dihydroorotase‑like protein
Reads supporting (aligned to +/- strand):  ref base A (10/4);  new base C (1/16);  total (11/20)
Fisher's exact test for biased strand distribution p-value = 3.51e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.86e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGCAGC  >  ACB122/5688514‑5688769
                                                                                                                           |                                                                                                                                    
ggTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCCGATCACCAGCGCCCGTGgc                                                                                                                           >  4:150431/1‑135 (MQ=255)
    ccATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCACCCAGATCACCCCCCCCCCGGgcgcgc                                                                                                                       >  2:296909/1‑135 (MQ=255)
            aGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGGGTGCGCCCGTCAGCCACCTCACCAGCGCCCGGGGCGCGCgggaggtc                                                                                                               >  5:205555/1‑129 (MQ=255)
            aGCTGCCGCGGCCTCCCGGGCAGCCCGGAGAGCGCCGAACCCGTGGCCATGGCGCGCAACCTCCTGCTGGGGCAACAGCGCGGCGTGCGTGCGCACTCCCGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATc                                                                                                               <  5:64436/135‑1 (MQ=255)
                  cTCGGCCTGCGGGTCATCCCGGAGCGCGCCGACACGGGGGCGCGGGCGCGCAACCTGGGGGTGGTCGAACAGAGCGGCGTGCGTGCTCAGTTCACCCAGATCCCCCGCGCGCGGGGCGCGCGGCTGATCgcccagg                                                                                                        <  8:15285/136‑1 (MQ=255)
                     ggCGTGCGGGGCACCCCGGAGCGCGCCGAAACCGGGGCCCGGGCGCGCAACAGGCTGCGGGTTTAACAGCGCGGCGGGGGGGCGCACTCCAGCCACACCACCCGCGCCCGGGGCGCGCGGCTGATCgcccaggc                                                                                                       <  8:6431/134‑1 (MQ=255)
                          gCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCCGGCGCACGTAACCCACAACACCACCGCCCCGGTCGCGCGGCTGCACgcccccgcccccgcc                                                                                                >  8:76942/1‑131 (MQ=255)
                                           gcgcCGAAACCGGGGCCCGGGCGCGCAACGTGCTGCGGGGGGAACAGAGCGGCGTGGGTGCGCACTGCCGCCAGACCACCAGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAAc                                                                               <  7:207512/136‑1 (MQ=255)
                                               cGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCCCTTCAGCCCCATCACCAGCCCCCGTGGCGCGCGGCGGATCGCCCCCGCCCCGGGCGTGGGCCAGCCGGGAAacgc                                                                            >  2:33443/1‑135 (MQ=255)
                                                      ggggCCCTGGCGCGCACCCGGCTGGGGGGGGACCCGCGCGGCGGGCGGGCGCACTGCCCCCAGAGCCCCCGCGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTa                                                                     <  1:63832/133‑1 (MQ=255)
                                                      ggggCCCTGGCGCGCAACCTGCTGCGGGTGGAAGAGCGCGGCGGGCGGGCGCACTTCACCCTGATCACCCGCGCCGGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTa                                                                     <  3:359651/133‑1 (MQ=255)
                                                      gTGGCCCTGGCGCGCACCCTGCTGCTGGTGGAACAGAGGGGCGTGCGTGGGCACTTCAGCCAGATCACCAGCGCCGGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTa                                                                     <  3:364028/135‑1 (MQ=255)
                                                        ggCCGTGGCGCCCAACGTGCTGCTGGTGCAACAGCGCGGGGGGCGTGCGCACTCCCCCCTGAGCCCCCGCGCCGGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAgc                                                                   <  2:299297/135‑1 (MQ=255)
                                                          cccTGGGGCGCAACCTCCTGAGGGGGCAACCGAGGGGCGTGCGGGCGCATTCCAGCCAGATCCCCCGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAgcgc                                                                 <  6:439745/135‑1 (MQ=255)
                                                          cccTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGGGTGCGCCCTTCAGCCAGATCACCAGCGCCCGGGGCGCGCGGCGGCTCGCCCCGGCCCCGGCCGGGGGCCTGCCCGGGACCGCCGCTGGAgcgc                                                                 >  3:79101/1‑135 (MQ=255)
                                                          cccTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGGGCCCTTCTGCCCGATCTCCACCGCCCGTGGCGCGCGGGGGTGCGCCCAGGCCCAGGCCCTGGGCCGGGCGCTATCCGCCGCGGTAtcgc                                                                 >  1:74217/1‑135 (MQ=255)
                                                                cgcgcAACCTGCTGCGGGGGGACCTGCGCGGGGTGCGTGCGCACTTCAGCCACATCCCCCGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTAcc                                                           <  6:444299/135‑1 (MQ=255)
                                                                   gcACCCTGCGGCGGGGGGAACAGCGCGGCGGGCGTGCGCACTTCCGCCACATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCt                                                       <  8:270365/136‑1 (MQ=255)
                                                                    cAACCTGCTGGGGGGGTACCTGAGCGGCGTGCGTGCGCACTTCCGCCATATCACCCGCGCCGGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCt                                                       <  6:182105/135‑1 (MQ=255)
                                                                             gctgGTGGAACAGAGCGGCGTGGGGGCGCCCGTCAGCCAGATCACCAGCCCCCCGGGCGCGCCGCGGGTCGCCCCCGCCCCCGGCGTGGGGCAGCCGGTAACCCCCCATGTAGCGCGGTTCCCGCTGCTCCacac                                              >  6:180583/1‑135 (MQ=255)
                                                                                             ggcgggcgggcgCACGACCGCCATACCCCCCGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTgcgc                              <  6:205555/135‑1 (MQ=255)
                                                                                                 ggcgggccCACTTCCGCCAGAGCCCCCGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGt                          <  3:71333/134‑1 (MQ=255)
                                                                                                  gcgggcgCACTTCCGCCACACCACCCGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGtt                         <  5:180583/135‑1 (MQ=255)
                                                                                                  gcgggcgCACTGCCGCCCCAGCCCCCGCGCCCGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTtc                        <  7:243983/136‑1 (MQ=255)
                                                                                                   ggggcgcACTCCCGCCAGATCCCCCGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTtc                        <  3:352322/132‑1 (MQ=255)
                                                                                                      gcgcACTGCCCCCAGACCCCCCGCGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTcc                     <  6:377416/135‑1 (MQ=255)
                                                                                                           cTTCAGCCACAGCACCCGCGCCGGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCt                <  3:273380/135‑1 (MQ=255)
                                                                                                                     aTCACCCGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGc     <  7:447139/136‑1 (MQ=255)
                                                                                                                       cACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGCa    >  1:262540/1‑135 (MQ=255)
                                                                                                                         ccAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCCGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGGTCTCCCGCCTGTACCACGTGCAgc  >  4:320968/1‑135 (MQ=255)
                                                                                                                         ccAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGCAgc  >  1:118317/1‑135 (MQ=255)
                                                                                                                           |                                                                                                                                    
GGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGCAGC  >  ACB122/5688514‑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: ‑