Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1435,689,1720AC42.5% 44.9 / 10.7 40P118P (CCA→CCCCJ019_05336hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (7/16);  new base C (17/0);  total (24/16)
Fisher's exact test for biased strand distribution p-value = 6.69e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGCCCAGGACAACCAGCAAGGCAGCGGCACCCTGGACAAGAAAGCGGTGCCCAAGACCACTGAGCTGGCGCCGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCACCCCCTCCCCCGAAAAAGTCGGCAGTGGTGACGACGGCGCCCAAGCCGCACAAGGTCGAGCCCAAACCCAAGGAAAGCAAGGCGCAGCCCAAGCCAGCCGCGCCGACGCCGGACTT  >  ACB143/5689037‑5689288
                                                                                                                                       |                                                                                                                     
ggCCCAGGACAACCAGCAAGGCAGCGGCACCCTGGACAAGAAAGCGGTGCCCAAGACCACTGAGCTGGCGCCGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCAc                                                                                                                      <  5:277462/137‑1 (MQ=255)
  cccAGGACAACCAGCAAGGCAGCGGCACCCTGGAGAAGAAAGCGGTGCCCAAGCCCACTGAGCTGGCGCCGTTCCAGGACAGCAAGATCAACAAGATCACCCCCACGCCTGAGGCCAGGCCCGAGGTGGTCCCAccc                                                                                                                    >  2:212567/1‑137 (MQ=255)
   ccAGGACCACCAGCAAGGCAGCGGCACCCTGGACAAGAAAGCGGTGCCCAAGACCACTGAGCTGGCGACGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCAcccc                                                                                                                   <  3:26954/137‑1 (MQ=255)
   ccAGGACAACCAGCAAGGCAGCGGCACCCTGGACAAGAAAGCGGTGCCCAAGACCACTGAGCTGGCGCCGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTcccccccc                                                                                                                   >  7:248836/1‑137 (MQ=255)
   ccAGGACAACCAGCAAGGCAGCGGCACCCTGGACAAGAAAGCGGTGCCCAAGACCACTGAGCTGGCGCCGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCAcccc                                                                                                                   <  8:225507/137‑1 (MQ=255)
   ccAGGACAACCAGCAAGGCAGCGGCACCCTGGACAAGAAAGCGGTGCCCAAGACCACTGAGCTGGCGCCGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCAcccc                                                                                                                   <  8:220423/137‑1 (MQ=255)
        aCAACCAGCAAGGCAGCGGCACCCTGGACAAGAAAGCGGTGCCCAAGACCACTGAGCTGGCGCCGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCccacccccaccc                                                                                                              >  1:233985/1‑137 (MQ=255)
        aCAACCAGCAAGGCAGCGGCACCCTGGACAAGAAAGCGGTGCCCAAGACCACTGAGCTGGCGCCGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCCCCCCCTccc                                                                                                              >  3:102128/1‑137 (MQ=255)
        aCAACCAGCAAGGCAGCGGCACCCTGGACAAGAAAGCGGTGCCCAAGACCACTGAGCTGGCGCCGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCACCCCCGccc                                                                                                              >  8:355664/1‑137 (MQ=255)
          aaCCAGCAAGGCAGCGGCACCCTGGACAAGAAAGCGGTGCCCAAGACCACTGAGCTGGCGCCGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCacccccaccccc                                                                                                            >  7:261917/1‑137 (MQ=255)
          aaCCAGCAAGGCAGCGGCACCCTGGACAAGAAAGCGGTGCCCAAGACCACTGAGCTGGCGCCGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCACCCCCTccccc                                                                                                            <  6:58333/137‑1 (MQ=255)
          aaCCAGCAAGGCAGCGGCACCCTGGACAAGAAAGCGGTGCCCAAGACCACTGAGCTGGCGCCGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCACCCCCTccccc                                                                                                            >  3:580/1‑137 (MQ=255)
              aGCAAGGCAGCGGCACCCTGGACAAGAAAGCGGTGCCCAAGACCACTGAGCTGGCGCCGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCCCCCCCGCCCCCGcaa                                                                                                        >  4:21064/1‑137 (MQ=255)
              aGCAAGGCAGCGGCACCCTGGACAAGAAAGCGGTGCCCAAGACCACTGAGCTGGCGCCGTTCCAGGACAGAAAGATCAACAAGATAACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCCCACCCTCCCCCGaaa                                                                                                        >  6:166869/1‑137 (MQ=255)
               gCAAGGCAGCGGCACCCTGGACAAGAAAGCGGTGCCCAAGACCACTGAGCTGGCGCCGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGCGGTGGTCCCCCCCCCACCCCCGaaaa                                                                                                       >  8:289918/1‑137 (MQ=255)
                    gCAGCGGCACCCTGGACAAGAAAGCGGTGCCCAAGACCACTGAGCTGGCGCCGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCACCCCCTCCCCCGAAAAAGTCg                                                                                                  <  7:212568/137‑1 (MQ=255)
                         ggCACCCTGGACAAGAAAGCGGTGCCCAAGACCACTGAGCTGGCGCCGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCCCCCCCGCCCCCCAAAAACTCGGCAgg                                                                                             >  3:286365/1‑136 (MQ=255)
                         ggCACCCTGGACAAGAAAGCGGTGCCCAAGACCACTGAGCTGGCGCCGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCCCCCCCACCCCCCCAAAAGTCGGCAgg                                                                                             >  3:281030/1‑136 (MQ=255)
                               cTGGACAAGAAAGCGGGGCCCAAGACCACTGAGCTGGCGCCGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCACCCCCTCCCCCGAAAAAGTCGGCAGTGGTgac                                                                                       <  2:307728/137‑1 (MQ=255)
                                        aaaGCGGTGCCCAAGACCACTGAGCTGGCGCCGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCCCCCCCGCCCCCCAAAAAGTAGGCCGTGGGggcgccggcgcc                                                                             >  4:160428/1‑137 (MQ=255)
                                         aaGCGGTGCCCAAGACCACTGAGCTGGCGCCGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCACCCCCTCCCCCGAAAAAGTCGGCAGTGGTGACGACGGCGccc                                                                            <  8:292499/137‑1 (MQ=255)
                                               tGCCCAAGACCACTGAGCTGGCGCCGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCCCCCCCTCCCCCGAAAAAGTAGGCAGTGGTGAAGACGGCGGCCCACGcc                                                                       >  8:153714/1‑137 (MQ=255)
                                               tGCCCAAGACCACTGAGCTGGCGCCGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCCCCCCCGCCCCCCAAAAAGTGGGCAGTGGTGCCGACGGCGCCCAAgacc                                                                      >  1:241728/1‑134 (MQ=255)
                                                  ccAAGACCACTGAGCTGGCGCCGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCCCCCCCCCCCCCGAAAAAGACGGCAGTGGTGAACACGGCCCCCAACCCGcac                                                                   >  5:191787/1‑137 (MQ=255)
                                                     aGACCACTGAGCTGGCGCCGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCACCCCCTCCCCCGAAAAAGTCGGCAGTGGTGACGACGGCGCCCAAGCCGCACAAg                                                                <  1:216867/137‑1 (MQ=255)
                                                          aCTGAGCTGGCGCCGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCACCCCCTCTCCCGAAAAAGTCGGCAGTGGTGACGACGGCGCCCAAGCCGCACAAGGTCGa                                                           <  8:46544/137‑1 (MQ=255)
                                                              aGCTGGCGCCGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCCCCCCCACCCCCCAAAAAATCGGCAGTGGTGACCACCGCGCCCCAGCCCCCCAACGGAGAGccc                                                       >  4:150962/1‑137 (MQ=255)
                                                                      ccGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCACCCCCTCCCCCGAAAAAGGCGGCCGTGGTGACGACGGCGCCCAAGCCGCACAAGGTAGAGCcccaacccaa                                               >  6:181496/1‑137 (MQ=255)
                                                                      ccGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCACCCCCAACCCCCCAAAAGTCGGCAATGGTGTCGACGGCGCCCACGCCGCACCACGTCGAGCCCAAACCCaa                                               >  7:332838/1‑137 (MQ=255)
                                                                           ccAGGACAGCATTATCAACAAGATCACGCCACCGACTTGGGCCTGGCCCTTGGTGGTCCCACCCCCTCCCCCGAAATTGTCGTCTGTGGTGACGACGGCGACCAAGCCTCACTAGGTCTAGACCAAACCCAAGGaaa                                          <  2:152408/137‑1 (MQ=255)
                                                                           ccAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCCCCCCCGCCCCCGCAAAAGTCGGCCGTGGTGACGACGGCGCCCCACCCGCCCAACGTCGCGCCCCAACCCAAGGcaa                                          >  8:113732/1‑137 (MQ=255)
                                                                            cAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCACCCCCTCCCCCGAAAAAATCGGCACTGGTGTCCACGGCGCCCCAGCAGCACACGCACGAGCCCACACCCAAGGAAAg                                         >  4:330193/1‑137 (MQ=255)
                                                                                aCAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCACCCCCTCCCCCGAAAAAGTCGGCAGTGGTGACGACGGCGCCCAAGCCGCACAAGGTCGAGCCCAAACCCAAGGAAAGCAAg                                     <  2:233985/137‑1 (MQ=255)
                                                                                aCAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCACCCCCTCCCCCGAAAAAGTCGGCAGTGGTGACGACGGCGCCCAAGCCGCACAAGGTCGAGCCCAAACCCAAGGAAAGCAAg                                     <  8:305484/137‑1 (MQ=255)
                                                                                aCAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCACCCCCTCCCCCGAAAAAGTCGGCAGTGGTGACGACGGCGCCCAAGCCGCACAAGGTCGAGCCCAAACCCAAGGAAAGCAAg                                     <  7:276663/137‑1 (MQ=255)
                                                                                    cAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTGCCCCCCCCACCCCCCAAAAACTCGGCAGTTGTGTAGAACGCGCCCCAGCCCCCCAACGGAGAGCCCAAaaccaacgcaaccaacgccc                                 >  4:270123/1‑135 (MQ=255)
                                                                                       gATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTGCCCCCCCCCCCCCCCCAAAAATCGGCCGGGGGGTCGCCCGCCCCCCAGCCCCCCAAAGGAGCGCCCCAACCCCAGGAAAACCACgcgcccc                              >  2:19704/1‑134 (MQ=255)
                                                                                       gATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCCCCCCCTCCCCCGCAAAAGTCGGCAGTGGTGCCGACGGCGCCCAAGCCGCACAAGGTCGCGCCCAAACCCAAGGTAAGCAAGGcgccgc                              >  3:63052/1‑137 (MQ=255)
                                                                                       gATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCCCCCCCGCCCCCGAAAAAGTCGGCAGTTGTGGCCCCGGCGCCCCACCCGCCCAAAGTCGCGCCCCAACCCAAGGAAAGCCAGGCGCAGc                              >  6:74096/1‑137 (MQ=255)
                                                                                                        cACCGCCTGCGGCCAGGCCCGAGGTGGTCCCACCCCCTCCCCCGAAAAACTCGGCAGTGGTGACGCCGGCGCCCCAGCCGCACAAGGTCGCGCCCCAACCCAAGGAAAGCCAGGCGCCGCCCCAGCCCGCCGCCCCg             >  3:119967/1‑137 (MQ=255)
                                                                                                         aCCGCCGGCGGCCAGGCCCGAGGTGGTCCCACCCCCTCCCCCGAAAAAGTCGGCAGTGGTGACGACGGCGCCCAAGCCGCACAAGGTCGAGCCCAAACCCAAGGAAAGCAAGGCGCAGCCCAAGCCAGCCGCGCCGa            <  6:284936/137‑1 (MQ=255)
                                                                                                                gCGGCCAGGCCCGAGGTGGTCCACCCCCCACACCCGAAAAAGTCGGGAGTTGTGTCCACGGCCCCCACGCCCCACAATGTGGCGCCCCACCCCAAGGAAAGCAAGGCCCAGCCCAAGCCCCGCGCCCCGCCGCCggc     >  6:353875/1‑136 (MQ=255)
                                                                                                                   gCCAGGCCCGAGGTGGTCCCACCCCCTCCCCCTAAAAAGTCGGCAGTGGTGACGACGGCGCCCAAGCCGCACAAGGTCGAGCCCAAACCCAAGGAAAGCAAGGCGCAGCCCAAGCCAGCCGCGCCGACGCCGGACtt  <  6:191787/137‑1 (MQ=255)
                                                                                                                                       |                                                                                                                     
GGCCCAGGACAACCAGCAAGGCAGCGGCACCCTGGACAAGAAAGCGGTGCCCAAGACCACTGAGCTGGCGCCGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCACCCCCTCCCCCGAAAAAGTCGGCAGTGGTGACGACGGCGCCCAAGCCGCACAAGGTCGAGCCCAAACCCAAGGAAAGCAAGGCGCAGCCCAAGCCAGCCGCGCCGACGCCGGACTT  >  ACB143/5689037‑5689288

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: