Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1223,440,3480TG54.2% 1.4 / 12.6 28A19A (GCT→GCGCJ019_03266hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (9/2);  new base G (8/5);  total (21/7)
Fisher's exact test for biased strand distribution p-value = 3.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAA  >  ACB122/3440217‑3440473
                                                                                                                                   |                                                                                                                             
cAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTggcg                                                                                                                           >  7:435661/1‑136 (MQ=255)
                 gtAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCGCGCgccgccgc                                                                                                           >  5:400623/1‑133 (MQ=255)
                    gCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCgggggg                                                                                                       >  7:177086/1‑134 (MQ=255)
                    gCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCCCCCCGGCAcccagcgg                                                                                                        >  6:154799/1‑130 (MQ=255)
                    gCAACGGTCTGGAGGATCGCAAGGCACTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCAGGCGGCGCCCGCGCCGCCCggcg                                                                                                        >  1:279668/1‑133 (MQ=255)
                    ccAACCGACTGGAGGATCGCAAGCCCCTCTACAAGCGAGCGCTCGAGCAGCTGCAGAGAGCGCATCCGCCACCAGCCTGGTCGCCGTGGCGCTGCTGCTGGACACCACCGCAGCCCGCGCCCCGGCCCCCTgaggg                                                                                                       >  8:248292/2‑132 (MQ=255)
                         ggTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACACGCCTTGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCAGGGGGCGCCCGGGCGACCCGGGGGTgggc                                                                                                  >  8:324718/1‑134 (MQ=255)
                             tGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCTGCCGGCGCCCCCGCACCCGGGGGGGTGCAGGc                                                                                               >  6:34529/1‑135 (MQ=255)
                                       cAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGGGGCGCCCGGGCCGCCTGGGGGTGGGAGGCCCCCGccga                                                                                      >  3:117738/1‑132 (MQ=255)
                                        aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCCGGCCGCCGGGGGGTGGCGGGCCCACGCCGATaa                                                                                   >  8:136122/1‑136 (MQ=255)
                                        aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCAGGCCGCGCCCCCGCACCCTGGGTGTGGGAGGCCCCCGccga                                                                                      >  3:192430/1‑131 (MQ=255)
                                                ctacccgcgcgcgcTCGCGGTGCGGCAGTGAGCGCGGGGGCCCCCCGCCTGGGCGCCGTGGCGCTGCTGCGGGCCGCCAGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGc                                                                           <  7:67463/130‑1 (MQ=255)
                                                         agcgCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGGCCGGGCCGCCGGGGGGTGGGGGGCCCACGCCTATAAGGAGCCGCTCGcgggt                                                                   >  1:255267/1‑131 (MQ=255)
                                                         agcgCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCCCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGGGTGGCCGGCCCACCCCGATAAGGCGCCgcgcgcggct                                                                   >  8:338517/1‑131 (MQ=255)
                                                            gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGGAGGCCCCCGCCGATAAAGCGCAGCACGCGGCTCac                                                                >  1:312275/1‑134 (MQ=255)
                                                            gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCGGGGTGTGGCAGGCCAACGCCTAGAAGGCGCAGCTCGCTGATCAgg                                                               >  7:51453/1‑136 (MQ=255)
                                                              tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGGCGCCCGCGCGGGCGGCGCCCCCGCCGCCCGGGGGTGGCGGGGCCACCCCAATAACGCGCCGCTCGCGGGTCAGGc                                                              >  1:68860/1‑135 (MQ=255)
                                                              tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGGGGGGCGGGGCCACGCCGATAAGGCGCAGCGCGCGGCTCAGGc                                                              >  1:242162/1‑135 (MQ=255)
                                                              tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCGGGCGGCGCCCCGGCCGCCTGGGGGTGGGCCGCCCCCCCCTATAATGCGCCGCTCGCTGGTCAgggc                                                             >  8:178051/1‑134 (MQ=255)
                                                              tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCCCCCGCGCGCGCGGCGCCCGGGCCGCCGGGGTGGGGCCGGCCCCCGCCGAGAATGCGCAGCGGGCGGATCAggg                                                              >  5:436981/1‑134 (MQ=255)
                                                                                 gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGGCGCCGGGGTGTGGCAGGCCCACCCCGATAATGCGCCGCTCGCGGGTCAGGCCGCTGACGTAGCCAAGCAg                                           >  2:354980/1‑135 (MQ=255)
                                                                                              gCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGTGTGGGAGGCCCACGCCGATAAGGCGCCGCGCGCGGATCAGGCCGAGGGCTTCGGCGCGCCGTTGGGGGAGaa                                >  5:173801/1‑133 (MQ=255)
                                                                                                    gggcggTGGGGCTGCGGCGGGCCGCGCGGGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc                       <  7:338517/131‑1 (MQ=255)
                                                                                                     gggggTGGCGCTGCGGCGGGCGGCCGGGGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc                       <  5:154799/131‑1 (MQ=255)
                                                                                                     cGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGTCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc                       <  1:284812/135‑1 (MQ=255)
                                                                                                                     ccggCCGCCCGCGCGGCCGCCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCgg       <  1:109035/133‑1 (MQ=255)
                                                                                                                       ggCCGCGCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGga     <  1:394100/135‑1 (MQ=255)
                                                                                                                          cgccAGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGaa  <  4:417784/135‑1 (MQ=255)
                                                                                                                                   |                                                                                                                             
CAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAA  >  ACB122/3440217‑3440473

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