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

CAGGCACCTTTGCCGTTGTTGAGGAAAGAGTACAGCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACGCAGGCTTCCTCGATCGAACGCCCAGCGATCAGCAGGC  >  ACB143/2308168‑2308438
                                                                                                                                        |                                                                                                                                      
cAGGCACCTTTGCCGTTGTTGAGGAAAGAGTACAGCTTCGATCAGCCCAGGCAGTCGACATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTCGGGGTGGAAATCCAGTCGTGGGc                                                                                                                                        >  1:240730/1‑137 (MQ=255)
 cggCACCTTTGCCGGTGTTGAGGAAAGAGTACAGCTGCGCTCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGcc                                                                                                                                       <  4:236050/136‑1 (MQ=255)
      ccTTTGCCGTTGTTGAGGAAAGAGTACAGCTTCGATCAGCGCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCttt                                                                                                                                  <  7:252904/137‑1 (MQ=255)
      ccTTTGCCGTTGTTGAGGAAAGAGTACAGCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCttt                                                                                                                                  >  6:261801/1‑137 (MQ=255)
      ccTTTGCCGTTGTTGAGGAAAGAGTACAGCTTCGATCAGCCCACGCAGTCGCCATGGGCACGCAAGGCACGGCGCCCGTAGTAGCTGAAGGCAGCACCGTGAAGTTAGGGGGTGGAAATCCAGTCGAGGGCCTCttt                                                                                                                                  >  7:185870/1‑137 (MQ=255)
                    gAGGAAAGAGTACAGCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCTGCACCGTGACGTTTGGGGGTGTAAATCCAGTCGTGGGGCTCTTTTGCCAGTGGTTAgg                                                                                                                    >  8:216093/1‑137 (MQ=255)
                    gAGGAAAGAGTACAGCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCCAGGCTTCgg                                                                                                                    >  8:221269/1‑137 (MQ=255)
                    gAGGAAAGAGTACAGCTCCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGGCGTGGGCCTCTTTGCCCAGTGCTTCgg                                                                                                                    <  6:201442/137‑1 (MQ=255)
                               aCAGCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCGCTTTGGCCAGTGCTTCGGGGATGGGGCgg                                                                                                         >  8:2741/1‑137 (MQ=255)
                                      cGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCTCGTTGCCCAGGGCTTCGGGGATGGGGCGGCTGTCAc                                                                                                  >  6:109749/1‑137 (MQ=255)
                                         tCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCGCGTTGCCCAGTGCGGCGGGGATGGGGGGGGTGTGACCgg                                                                                               >  1:41199/1‑137 (MQ=255)
                                           aGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCACTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGgcg                                                                                             >  5:334522/1‑137 (MQ=255)
                                                aGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCAt                                                                                        <  8:155121/137‑1 (MQ=255)
                                                aGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCGGTTTGCCCCGTGCGTGGGGGTTGGGGGGGGTGTCTCCGGCGCCCAt                                                                                        >  4:235538/1‑137 (MQ=255)
                                                       cGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCcagc                                                                                 <  7:168740/137‑1 (MQ=255)
                                                         ccATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGTGTTGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCcagcag                                                                               <  8:202735/137‑1 (MQ=255)
                                                          cATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGGTCGGTTGGCAGTGCTTCGGGGGTGGGGCGGGTGTTAACGGCCCCCCAGGCcagcaga                                                                              >  1:165213/1‑136 (MQ=255)
                                                          cATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCTCTTTGGCCAGGGGTTCGGGGATTGGGCGGATGTCACCGGCGCCCCGGGCCAGCAGt                                                                              >  8:168740/1‑137 (MQ=255)
                                                          cATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCGCTTTGGCCAGGGCGTGGGGGCTGGGGCGGGTGTGACCGGCGCCCCAGGCccgccct                                                                              >  5:117149/1‑134 (MQ=255)
                                                                     aaGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCGCGTTGGCCAGGGCTTAGGGGCTGGGGCGGCTGTCACCGGCGCCCCAGGCCAGCACTTGCAGCTTtgg                                                                   >  3:270749/1‑136 (MQ=255)
                                                                       ggCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCGGGTTGCCCAGTGGGTGGGGGATGGGGGGGCTGTCAACGGCGCCCATGGCCCGCAGTTGGCTCTTtgctg                                                                 >  2:346898/1‑137 (MQ=255)
                                                                         cACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTgcg                                                               <  2:227774/137‑1 (MQ=255)
                                                                            ggCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTgcgcgc                                                            <  5:109749/137‑1 (MQ=255)
                                                                                cgcgTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGGTCGTTGTCCACTGCTTGGGGGGTGGGGCGGATGTCACCGGCCCCCAAGGGCAGCAGTTGCATCTTTGCTTCGCGCTCgc                                                        >  4:226803/1‑136 (MQ=255)
                                                                                 gcgTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGaa                                                       <  7:221269/137‑1 (MQ=255)
                                                                                 gcgTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGaa                                                       <  7:216093/137‑1 (MQ=255)
                                                                                     agtagCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCGCTTTGCCCAGGGGTTGGGGGGTGGGGCGGGTGTCACCGGCGCCCCTGGGCCGCAGGTGGATCTTTTCTTCGCGGGCGCAcagc                                                   >  6:4002/1‑137 (MQ=255)
                                                                                         gCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCagag                                               <  4:142540/137‑1 (MQ=255)
                                                                                              aGGCAGCTCCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCTCTTTGCCGAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGGAGTTGCATCTTTGCTGCGCGCTCGAACAGCAgagcgagc                                          <  5:193345/137‑1 (MQ=255)
                                                                                                 cagcaCCGTGACGTTTGGGGGTGGAAATCCAGTCGTGCGCCGCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACg                                       <  5:258964/137‑1 (MQ=255)
                                                                                                       cGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGGCCAGTGCTTCCGGGTTGGGGCGGCTGTCCCCGGCGCCCCAGGCCCCCACTTTGATCATTGCTGCGCGCTCGAAAACAACAGCCAGCCCGCAGGCt                                 >  8:381019/1‑137 (MQ=255)
                                                                                                            cGTTTGGGGGTGGAAATCCAGTCGTGGGCCGCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCCGCGCCCAAGGCCCGCAGTTGGATCTTTTCTGCGCGCTCGAACACCAGAGCGAGCACGCAGGCTTCCtc                            >  2:100595/1‑137 (MQ=255)
                                                                                                            cGTTTGGGGGTGGAAATCCAGACGTGGGCCTCTTTGCCCAGAGCTTCGGGGATGGGGCGGATGTAACCGGCGCCCATGGCCAGCAGTTGCATCTTTTCTGCGCGCGCGAACACCACAGCGAGCACCCAGGCTTCCtc                            >  6:61051/1‑137 (MQ=255)
                                                                                                               ttGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACGCAGGCTTCCtcgat                         <  4:75284/137‑1 (MQ=255)
                                                                                                                      tGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTGGGGGGTGGGGCGGGTGTCACCGGCGCCCCTGGCCAGCAGTTGCATCTTTTCTGCGCGCTCGCACCGCAGAGCGAGCACGCCGGCTTCCGCGATCGAACGc                  >  8:123564/1‑137 (MQ=255)
                                                                                                                      gggAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACGCAGGCTTCCTCGATCGAACGc                  <  2:165213/136‑1 (MQ=255)
                                                                                                                              cAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACGCAGGCTTCCTCGATCGAACGCCCAGCGAt          <  7:72008/137‑1 (MQ=255)
                                                                                                                               aGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACGCAGGCTTCCTCGATCGAACGCCCAGCGATc         <  1:100595/137‑1 (MQ=255)
                                                                                                                                 tCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGGGGATGTCACCGGCGCCCAAGGCCACCACTTGGATCTTTGCTGCGCCCACGCAACGCAGAGCGAGCACGCAGGCTTCCTCGATCGAACGCCCAGCGATcag       >  6:183915/1‑137 (MQ=255)
                                                                                                                                  cGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACGCAGGCTTCCTCGATCGAACGCCCAGCGATcagc      <  3:241346/137‑1 (MQ=255)
                                                                                                                                    tGGGGCTGTTTGGCCAGTGGGTCGGGGGTGGGGCGGATGTCACCGGCGCCCATGGCCCGCAGTTGCATCTTTGCTGCGCGCGAGACCAGCAGAGCGCGCACGCAGGCTTGCTCGATCGAACGCCCCGCGCTcagcag    >  2:88285/1‑137 (MQ=255)
                                                                                                                                      ggCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACGCAGGCTTCCTCGATCGAACGCCCAGCGATCAGCAGGc  <  4:184297/137‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                      
CAGGCACCTTTGCCGTTGTTGAGGAAAGAGTACAGCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACGCAGGCTTCCTCGATCGAACGCCCAGCGATCAGCAGGC  >  ACB143/2308168‑2308438

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: