Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1222,308,3210CG35.4% 68.5 / 15.3 48A230P (GCC→CCC) fucAL‑fuculose phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base C (5/26);  new base G (17/0);  total (22/26)
Fisher's exact test for biased strand distribution p-value = 6.20e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AGGCACCTTTGCCGTTGTTGAGGAAAGAGTACAGCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACGCAGGCTTCCTCGATCGAACGCCCAGCGATCA  >  ACB122/2308186‑2308449
                                                                                                                                       |                                                                                                                                 
cggCACCTGTGCCGTTGTTGAGGAAAGAGTACAGCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGcc                                                                                                                                  <  2:286429/136‑1 (MQ=255)
aGGCACCTTTGCCGTTGTTGAGGAAAGAGTACAGCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGcc                                                                                                                                  <  5:374746/137‑1 (MQ=255)
     ccTTTGCCGTTGTTGAGGAAAGAGTACAGCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCttt                                                                                                                            <  6:181027/137‑1 (MQ=255)
     ccTTTGCCGTTGTTGAGGAAAGAGTACAGCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCttt                                                                                                                            <  6:1862/137‑1 (MQ=255)
       tttGCCGTTGTTGAGGAAAGAGTACAGCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGc                                                                                                                          <  6:205680/137‑1 (MQ=255)
       tttGCCGTTGTTGAGGAAAGAGTACAGCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGc                                                                                                                          <  5:346641/137‑1 (MQ=255)
           ccGTTGTTGAGGAAAGAGTACAGCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCGCTTTGCCCac                                                                                                                      >  6:58267/1‑136 (MQ=255)
             gttgttGAGGAAAGAGTACAGCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTg                                                                                                                    <  4:48737/137‑1 (MQ=255)
                   gAGGAAAGAGTACAGCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCgg                                                                                                              <  6:244301/137‑1 (MQ=255)
                              aCAGCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCGGGTTGACCCGGTCTTCTGGGATGgggggg                                                                                                   >  4:27198/1‑137 (MQ=255)
                              aCAGCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCGCGTTGTCCAGTGCTTGGGGGATTgggggg                                                                                                   >  5:231939/1‑137 (MQ=255)
                              aCAGCTTCGATCAGCCCAGGCAGTCGCCATCGCCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTTGCCACTGCTACGGGGATTGGGCgg                                                                                                   >  8:155929/1‑137 (MQ=255)
                                 gCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCGCGTTTCCCAGTGCTTGGGGGGTGggggggttt                                                                                                >  3:81453/1‑134 (MQ=255)
                                   ttCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTc                                                                                              <  7:305337/137‑1 (MQ=255)
                                     cGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCTCTTTGGCCCGTGCGTCGGGGGTGGGGGGGGTGTCAc                                                                                            >  5:62138/1‑137 (MQ=255)
                                       aTCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCGCGTTGCGCAGTGGGTCGGGGATGGGGGGGGTGTGAccc                                                                                          >  5:56321/1‑136 (MQ=255)
                                        tCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCTCTTTGGCCAGTGCTTGGGGGTTGGGGGGGGTGTCACCgg                                                                                         >  3:285837/1‑137 (MQ=255)
                                        tCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCTCTTTGCCCCGGGGTTGGGGGTTGGGGGGGGTGTGACCgg                                                                                         >  5:67251/1‑137 (MQ=255)
                                        tCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCTCTTTGCCCAGTGCGTCGGGGTTGGGGGGGGTGTCACCgg                                                                                         >  6:189413/1‑137 (MQ=255)
                                        tCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCGCTTTGGGCAGTGGGTGGGGGTTGGGGCGGGTGGGACCgg                                                                                         >  7:250206/1‑137 (MQ=255)
                                               aGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCTCTTTGGCCAGTGGTTAGGGGTTTGGGCGGCTGTGACCGGCGCCCaa                                                                                  >  7:162624/1‑136 (MQ=255)
                                                  cAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGc                                                                               <  1:112151/137‑1 (MQ=255)
                                                      cGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGTGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCcagc                                                                           <  2:173500/137‑1 (MQ=255)
                                                        ccATGTGCACGCAAGGCACTGCGCTCGTAGTATCTGAAGGCAGCACCGTGACGTTTGGGTGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCcagcag                                                                         <  5:210/137‑1 (MQ=255)
                                                        ccATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGGGCTTTGGCCAGTGGGTGGGGGGTGGGGGGGATGTGACCGGCGCCCATGGCccgcag                                                                         >  7:104330/1‑137 (MQ=255)
                                                        ccATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCGCTTTGCCCAGTGGGTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCcaccag                                                                         >  2:229569/1‑137 (MQ=255)
                                                        ccATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCcagcag                                                                         <  2:159471/137‑1 (MQ=255)
                                                         cATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGGGCTGTGGCCAGTGGGTTGGGGGTTGGGGGGCTGTGGACGGCGCCCATGGCcacccgt                                                                        >  3:331259/1‑132 (MQ=255)
                                                         cATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCGCGTTGGCCCGTGGGTCGGGGCTGGGGGGGGTGTGACCGGCGCCCATGGGCCGCAGt                                                                        >  5:219077/1‑137 (MQ=255)
                                                         cATGGGCACCCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCGCCATTACCACGGCGACCGGGAGGCGGCGGGTATCACCCGCGCCCCTGGCccgcaga                                                                        >  6:13048/1‑136 (MQ=255)
                                                          aTGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCTGTTTGGCCCGTTCTTAGGGGATAGGGGGGTTGTAACCCGCGCCCATGGCcagaaatt                                                                       >  6:210/1‑132 (MQ=255)
                                                                   cacggCACCGCGCGCGTACTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGTGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTTCATCTTtg                                                              <  8:157377/134‑1 (MQ=255)
                                                                        cACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTgcg                                                         <  2:265425/137‑1 (MQ=255)
                                                                              gcgcgTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCGCGTTTGCCCGTGGGGGGGGGGTGGGGGGGCTGTGACCGGCGCCCCAGGCCCGCCGGTGGCTCTTTTGTgcgcgcggg                                                   >  1:107090/1‑134 (MQ=255)
                                                                                gcgTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGaa                                                 <  3:27198/137‑1 (MQ=255)
                                                                                       agCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGGCGTGTGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCaga                                          <  3:114000/137‑1 (MQ=255)
                                                                                             aGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAgagcgagc                                    <  5:189413/137‑1 (MQ=255)
                                                                                                cagcaCCGTGACGTTTGGGGGTGGAAATCCAGTCGGGGTCCTCTTGGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACg                                 <  8:104771/137‑1 (MQ=255)
                                                                                                           cGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTTACCAGTGGTTCGGGGATGGGGCGGATGTCACCGGCGCCCAAGGCCAGCAGTTGCATCTTTTCTTCGCGCTCGAACAGCAGAGCGCGCACCCAGGCGTACtc                      >  6:166854/1‑137 (MQ=255)
                                                                                                           cGTTTGGGGGTGGAAATCCAGTCGTGGGCCGCTTTTGCCAGTGCGTCGGGGGTTGGGGGGATTTGACCCGCGCCCAAGGCCAGCAGTTGCATATTTTCTTCGCGCGCGAACACCAGACCGAGCCCGCCCGCTTccac                      >  4:191691/1‑135 (MQ=255)
                                                                                                              ttGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACGCAGGCTTCCtcgat                   <  6:275619/137‑1 (MQ=255)
                                                                                                                   ggTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACGCAGGCTTCCTCGATCGAAc              <  4:199816/137‑1 (MQ=255)
                                                                                                                   ggTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACGCAGGCTTCCTCGATCGAAc              <  5:162778/137‑1 (MQ=255)
                                                                                                                   ggTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACGCAGGCTTCCTCGATCGAAc              <  1:159477/137‑1 (MQ=255)
                                                                                                                   ggTGGAAATCCAGTCGTGGGCCGCTTTTGCCAGTGGTTAGGGGATTGGGGGGATGTCACCCGCGCCCCAGGCCAGCAGTTGGAACTTTTATGCGCGCGCGCAAACCAGAGCGAGCACCCAGGCTTACTCGATAGaaa              >  4:170965/1‑136 (MQ=255)
                                                                                                                       taaaTCCAGGCGTGGGCCTCTTTGGCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACGCAGGCTTGCTCGATCGAACGccc          <  7:81073/136‑1 (MQ=255)
                                                                                                                        aaaTGCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACGCAGGCTTCCTCGATCGAACGCCCa         <  1:106156/137‑1 (MQ=255)
                                                                                                                             cAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACGCAGGCTTCCTCGATCGAACGCCCAGCGAt    <  3:191691/137‑1 (MQ=255)
                                                                                                                              aGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCATCAGTTGGATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACGCAGGCTTCCTCGATCGAACGCAGAGCGATc   <  5:352530/137‑1 (MQ=255)
                                                                                                                               gTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACGCAGGCTTCCTCGATCGAACGCCCAGCGATca  <  5:58267/137‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                 
AGGCACCTTTGCCGTTGTTGAGGAAAGAGTACAGCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACGCAGGCTTCCTCGATCGAACGCCCAGCGATCA  >  ACB122/2308186‑2308449

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: