Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1432,308,3040CG31.4% 57.6 / 10.6 35A230P (GCC→CCC) fucAL‑fuculose phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base C (6/18);  new base G (11/0);  total (17/18)
Fisher's exact test for biased strand distribution p-value = 2.97e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCACCTTTGCCGTTGTTGAGGAAAGAGTACAGCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACGCAGGCTTCCTCGATCGAACGCCCAGCGATCAGC  >  ACB143/2308171‑2308434
                                                                                                                                     |                                                                                                                                  
gCACCTTTGCCGTTGTTGAGGAAAGAGTACAGCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGcc                                                                                                                                   <  8:272622/135‑1 (MQ=255)
  aCCTTTGCCGTTGTTGAGGAAAGAGTACAGCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCctc                                                                                                                                 >  7:185526/1‑135 (MQ=255)
     tttGCAGTTGTTGAGGAAAGAGTACAGCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCttt                                                                                                                              <  4:78165/135‑1 (MQ=255)
        gCCGTTGTTGAGGAAAGAGTACAGCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTg                                                                                                                             <  3:62452/133‑1 (MQ=255)
        gCCGTTGTTGAGGAAAGAGTACAGCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGc                                                                                                                            <  4:143007/134‑1 (MQ=255)
        gCCGTTGTTGAGGAAAGAGTACAGCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGc                                                                                                                            <  5:143381/134‑1 (MQ=255)
         ccGTTGTTGAGGAAAGAGTACAGCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGcc                                                                                                                           <  3:81243/134‑1 (MQ=255)
                 gAGGAAAGAGTACAGCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCTCTTTTGCCAGTGCGTc                                                                                                                  >  3:329466/1‑135 (MQ=255)
                        gagTACAGCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGTCCTCTTTGCCCAGTGCTTCGGGGATg                                                                                                           <  4:251601/135‑1 (MQ=255)
                            aCAGCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCTCTTTGGCCAGTGCTTCGGGGATTGGGc                                                                                                       >  3:77773/1‑135 (MQ=255)
                                      tCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCTCTTTGGCCAGTGCTTGGGGGATGGGGCGGCTGTGAcc                                                                                             >  4:172499/1‑135 (MQ=255)
                                        aGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGGGGGTTTCCCAGGGGTTGGGGGAAGGGGCGGGTGTGAcccc                                                                                           >  3:182620/1‑133 (MQ=255)
                                        aGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCGAGTTGGGCAGTGCGGGGGGGTTGGGGCGGCTGGGACCgg                                                                                           >  3:87372/1‑135 (MQ=255)
                                        aGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCGTTGGCCACGGCGTCGGGGATGGGGCGGCTGTCACCgg                                                                                           >  6:78065/1‑135 (MQ=255)
                                          cccAGGCAGTCGCCATGGGCACGCAAGGCACGGCTCGCGTATTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGACCAGTGCTTCGGGGATGGGGCGGATGTCACCGgcg                                                                                         <  4:277076/135‑1 (MQ=255)
                                             aGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCGGGGTGCCCCGTGGTGGGGGGAGGGGGGGGGTGTGACCGGGGcc                                                                                       >  3:202847/1‑134 (MQ=255)
                                                 aGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATgg                                                                                  <  7:292949/135‑1 (MQ=255)
                                                  gTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCGGTTTGCCCAGTGCTTGGGGGCTGGGGGGGGGGTGACCCGCGccccaggc                                                                                 >  5:291247/1‑130 (MQ=255)
                                                      ccATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCGCGTTTCCCAGTGCGTCGGGGATGGGGGGGATGTGACCGGCGCCCATGGCccgc                                                                             >  4:142250/1‑135 (MQ=255)
                                                      ccATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCGCTTTGCCCCGTGGGTAGGGGGTGGGGGGGCTGTCACCGGCGCCCCAGGCcagc                                                                             >  8:302916/1‑135 (MQ=255)
                                                       cATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCTCTTTGGCCAGTGCTTGGGGGGTGGGGGGGTTGTCAACGGCGCCCAAGGCcagca                                                                            >  2:215465/1‑135 (MQ=255)
                                                                   aGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTtg                                                                <  3:2356/135‑1 (MQ=255)
                                                                           cgcgcgTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTgcgcgc                                                        <  7:201125/135‑1 (MQ=255)
                                                                           cgcgcgTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTgcgcgc                                                        <  8:303371/135‑1 (MQ=255)
                                                                             cgcgTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCGCTTTGGCCAGTGCGTGGGGGGTGGGGCGGATGTCACCCGCGCCCCTGGCCAGCAGGTGGCTCGTTTATgcgcgcgc                                                      >  2:86269/1‑133 (MQ=255)
                                                                              gcgTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCg                                                     <  5:114749/135‑1 (MQ=255)
                                                                                             gcagcaCCGTGACGTTTGTGGGTGGACATCCAGTGGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCAGGGCCAGCAGTTGGATCGTTGCTGCGCGCTCGAACAGCAgagcgagc                                      <  5:33981/135‑1 (MQ=255)
                                                                                             gcagcaCCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAgagcgagc                                      <  8:328990/135‑1 (MQ=255)
                                                                                               agcaCCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCCGGGCTTGGGGGGTGGGGCGGGTGTCACCGGCGCCCCTGGGCCGCCGTTGGCTCATTGCTGCGCGCGCGACCCGCAgcccgagccc                                    >  2:338511/1‑133 (MQ=255)
                                                                                                    cGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCGCTTTGCCCACTGCGTCGGGGATGGGGGGGCTGTCAACGGCGCCCATGGCCAGCCGTTGCATCTTTGCTGCGCGCGCCAAAACCACAGCGAGCACCCAgg                               >  3:227910/1‑135 (MQ=255)
                                                                                                         cGTTTGGGGGTGGAAATCCAGTCGTGGGCCACGTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGCCGCCCAAGGCCCGCAGTTGCATCTTTTCTTCGCCCGCCAAAACCAGATCGAGCACCCACGCGTcc                          >  8:236803/1‑135 (MQ=255)
                                                                                                              gggggTGGAAATCCAGTCGTGGGGCTCTTTGCCCAGGGCGTGGGGGGTGGGGGGGCTGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTTCTGCGCGCTCGCACAGCAGAGCGAGCACGCCGGCGTCCtcgat                     >  1:228696/1‑135 (MQ=255)
                                                                                                                   tGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACGCAGGCTTCCTCGATCGAAc                <  4:144322/135‑1 (MQ=255)
                                                                                                                             gTCGTGGGCCTCTTTGCCCTGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACGCAGGCTTCCTCGATCGAACGCCCAGCGAt      <  6:179049/135‑1 (MQ=255)
                                                                                                                             gTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACGCAGGCTTCCTCGATCGAACGCCCAGCGAt      <  8:65310/135‑1 (MQ=255)
                                                                                                                                 tGGGCCTCTTTGGCCAGTGCTTCTGGGATGGGGCTGAGGTCACCGGCGCCCATGGCCAGCAGTTGGATATTTGCTGGGCGCTCGAACAGCAGAGCGAGCAGGCAGGATTCCTCGATCGAACGCCCAGCGATctgc  <  2:162394/135‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                  
GCACCTTTGCCGTTGTTGAGGAAAGAGTACAGCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACGCAGGCTTCCTCGATCGAACGCCCAGCGATCAGC  >  ACB143/2308171‑2308434

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