Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1223,440,3670TG68.0% ‑4.7 / 14.5 26W26G (TGG→GGG) CJ019_03266hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (4/4);  new base G (17/0);  total (22/4)
Fisher's exact test for biased strand distribution p-value = 5.53e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.20e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAGCGCTCGATCAGCTGGC  >  ACB122/3440237‑3440503
                                                                                                                                  |                                                                                                                                        
gCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGgtgt                                                                                                                                    >  5:248898/1‑137 (MQ=255)
gCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCAGGgtgt                                                                                                                                    >  6:449134/1‑137 (MQ=255)
          ggaggaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGCCaa                                                                                                                          >  3:378758/1‑137 (MQ=255)
                   cAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCCGGCACCCGGGGTGTGGCAGGCCACCGCCTATaa                                                                                                                 >  8:8440/1‑137 (MQ=255)
                   cAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCCCGGGCGGCGCCCCGGCAGCCGGGGGGTGGCAGGCCAACGCCTATaa                                                                                                                 >  4:125807/1‑137 (MQ=255)
                    aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGGCGCCCGCCCAGGCGGCGCCCCGGCCGCCGGGGGGTGGCCGGCCCCAGCCGATAAg                                                                                                                >  2:123152/1‑137 (MQ=255)
                    aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCTGGGGGTGGCAGGCCCACCCCGATaaa                                                                                                                >  2:123148/1‑136 (MQ=255)
                    aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCCCGGGCGGCGCCCGGGCCGCCGGGGTGTGGGAGGCCAACCCCGATaaa                                                                                                                >  3:242415/1‑136 (MQ=255)
                    aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCACCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGGCACCGCCGATAAg                                                                                                                >  5:232086/1‑137 (MQ=255)
                             tACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCGGGGTGTGGCAGGCCAACGCCTATAATGCGCAGCTc                                                                                                       >  6:72931/1‑137 (MQ=255)
                              aCAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGCTGGCCAACGCCGAAAAGGCGCAGCTCg                                                                                                      >  5:30583/1‑137 (MQ=255)
                                aaGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCTGGCGGCGCCCGCGCCGCCGGGGTGTGGCACGGCAACGCCGATAATGCGCAgcgcgcg                                                                                                    >  1:329454/1‑136 (MQ=255)
                                 cccgcgcgcTCGTGGTGCGGCAGTTCCGGCTTGGGCCCCCCGGGGGGGCGCGGGGGCGCTGCTGCGGGCCCCCCGGGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTg                                                                                                   <  8:194619/132‑1 (MQ=255)
                                        gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCGGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGc                                                                                            >  8:136680/1‑137 (MQ=255)
                                        gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGACGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCACCGCCTATAAGGCGCAGCTCGCTGATCAGGc                                                                                            >  1:127766/1‑137 (MQ=255)
                                             aGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCCGCCGGGGGGTGGGAGGCCCACGCCGATAAGGCGCCGCTCGCTGATCAGGCCGCTg                                                                                       >  6:398540/1‑137 (MQ=255)
                                                tgctgcAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGTGTGGGCGGCCAACGCCGATAAGGCGCAGCGCGGGGATCAGGCCGCTGAcg                                                                                    >  4:191546/1‑136 (MQ=255)
                                                             gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGCCGGCCAACCCCGATAAGGCGCCGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGtt                                                                       >  8:108020/1‑137 (MQ=255)
                                                             gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGGGTGGCCGGCCCACCCCGATAAGGCGCCGCTCGCTGCTCAGGCCGCTGACGTCGGCAAGCAGtt                                                                       >  6:225777/1‑137 (MQ=255)
                                                                       cccGCCTGGTCGCCGTGGCGCTGCTGCTGGGCGCCCGCCCTGGCGGGGCCCCGGCCCCCGGGGGGTGGGGCGCCCACGCCTATAACGCGCAGCACGCCGCTCAGGCCGCGGTCATCGGCCAGCACTTGGGGGGGaaa                                                             >  4:349608/1‑136 (MQ=255)
                                                                          gCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGGGTGGGCGGCCCACGCCGATAAGGCGCCGCTCGCTGCTCAGGCCGCTGACTTCGGCAAGCAGTTGGGGGGGAAGGAt                                                          >  6:182620/1‑137 (MQ=255)
                                                                               ggcgcgGTGGCGCTGCTGCGGGCCCCCCGGGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc                                                     <  8:232619/135‑1 (MQ=255)
                                                                                     gTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCTGGGGGTGGGCGGCCCACGCCGCTAAGGCGCCGCGCGGGGGTCAGGCCGAGGGCTTGGGCGAGCGGTTGGGGGGGAAGGGTcgcgcgggcgg                                               >  2:50988/1‑132 (MQ=255)
                                                                                          gccgcggctgGCCGGGGGGGCCGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTgg                                          <  1:92086/131‑1 (MQ=255)
                                                                                                 ggggCCGCCCGCGGTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGga                                   <  6:373163/135‑1 (MQ=255)
                                                                                                                          gggCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGAAGCAGCGCTCGAt          >  5:91857/1‑137 (MQ=255)
                                                                                                                                  tGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAGCGCTCGATCAgctggc  <  7:136680/137‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                        
GCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAGCGCTCGATCAGCTGGC  >  ACB122/3440237‑3440503

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: