Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1225,688,6450TG51.4% 10.8 / 24.1 35R238R (CGT→CGG) pyrC'Dihydroorotase‑like protein
Reads supporting (aligned to +/- strand):  ref base T (8/9);  new base G (5/13);  total (13/22)
Fisher's exact test for biased strand distribution p-value = 3.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.13e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGCAGCCGCCGCTGC  >  ACB122/5688514‑5688778
                                                                                                                                   |                                                                                                                                     
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)
                                                                                                                              gcCCGGGGCCCGCGGCGGCTCGCCCAGGCCCAGGCCGTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGCAGccgccg      <  1:388526/135‑1 (MQ=255)
                                                                                                                                ccGGGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGCAGCCGCCGCt    <  4:251807/135‑1 (MQ=255)
                                                                                                                                  gTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGCCGCCGCCGCTgc  >  1:165234/1‑135 (MQ=255)
                                                                                                                                  gTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGGTGTCCAGCCTGTACCCCGTGCAGCCGCCGCAgc  >  5:54693/1‑135 (MQ=255)
                                                                                                                                   |                                                                                                                                     
GGTGCCATGGCCAGCTTCCTCGGCCTGCCGGGCATCCCGGAGAGCGCCGAAACCGTGGCCCTGGCGCGCAACCTGCTGCTGGTGGAACAGAGCGGCGTGCGTGCGCACTTCAGCCAGATCACCAGCGCCCGTGGCGCGCGGCTGATCGCCCAGGCCCAGGCCTTGGGCCTGCCGGTAACCGCCGATGTAGCGCTGTACCAGCTGATCCTCACGGACGAGTCGCTGCGCGAGTTCTCCAGCCTGTACCACGTGCAGCCGCCGCTGC  >  ACB122/5688514‑5688778

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