Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1433,440,3460AC68.2% ‑5.5 / 14.1 22A24A (GCA→GCC) CJ019_03266hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (1/6);  new base C (15/0);  total (16/6)
Fisher's exact test for biased strand distribution p-value = 9.38e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.57e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAGCGCTCGATCA  >  ACB143/3440217‑3440480
                                                                                                                                 |                                                                                                                                      
gtAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCCGgccgccgg                                                                                                                                   >  4:154503/1‑133 (MQ=255)
 tAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCTgg                                                                                                                                  >  7:299604/1‑135 (MQ=255)
   gCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCTgggg                                                                                                                                >  3:307060/1‑134 (MQ=255)
          tCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCCCCAGCGCTGGCGGCGCCCGGGCCGCCGGGGGGTGGccg                                                                                                                         >  8:98475/1‑133 (MQ=255)
             ggaggaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGcc                                                                                                                      >  3:98531/1‑135 (MQ=255)
               aggaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCTGGCGGCGCCCGGGCCGCCTCGGTGTGGCAGGCCaa                                                                                                                    >  5:297186/1‑135 (MQ=255)
               aggaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCCCCACCGCAGGCCGCGCCCGGGCCGCCGCGGGGTGGCCGGCCac                                                                                                                    >  3:241857/1‑134 (MQ=255)
                      cAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCCGCCTGGGTGTGGCAGGCCAACGCCta                                                                                                              >  6:90850/1‑134 (MQ=255)
                      cAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCACCGCTGGCGGCCCCCCCGCCGCCCGGGTGTGGGCTGCCACCGCCtat                                                                                                             >  6:106191/1‑135 (MQ=255)
                                 aCAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCTGGGTGTGGCAGGCCCACGCCGATAAGGCGCCGCt                                                                                                  >  3:297431/1‑135 (MQ=255)
                                 aCAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGCGCCGCCGCGGTGTGGCCGGCCCCCCCCGAGAAGgcgccccg                                                                                                  >  5:107761/1‑131 (MQ=255)
                                        agcgCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGCCAACCCCTATAAGGCGCAGCTCGCTGa                                                                                            >  3:275791/1‑134 (MQ=255)
                                           gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCCGGCCGCCGGGGGGTGGCCGGCCGACGCCAATAATGCGCACCGCGCGGATCAg                                                                                        >  3:258084/1‑135 (MQ=255)
                                           gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCCGCCTGGGTGTGGCAGGCCACCGCCGATAAGGCGCCGCTCGCTGATCAg                                                                                        >  7:114393/1‑135 (MQ=255)
                                           gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAg                                                                                        >  7:12636/1‑135 (MQ=255)
                                                                             gCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCTGGCCGCGCCCGGGCCGCCTGGGTGTGGCAGGCCCACACCTAAAACGCGCAGCGAGCGGATCACGCCGATGACGACGGCTAGCACTTGGGGGAGAAgg                                                      >  8:258617/1‑135 (MQ=255)
                                                                                 ggTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGGCGCCGGGGGGTGGGCGGCCCCCCCCGATAAGGCGCCGCGCGCGGGTAAGGGCGCTGACGTCGGCAAGCAGTTGGCGGAGCAGGATcg                                                  >  4:253780/1‑135 (MQ=255)
                                                                                   gggcggTGGGGCTGCCGCGGCCCGCCCGGGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcg                                                <  3:262925/130‑1 (MQ=255)
                                                                                    cgggggggCGCGGCGGCGGCCCGCCCGGGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGTCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc                                               <  4:307060/129‑1 (MQ=255)
                                                                                    cgcgGTGGGGCTGCGGGTGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc                                               <  6:240520/135‑1 (MQ=255)
                                                                                    cGCCGTGGGGGTGGTGCTGGCCGCCAGCGGTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc                                               <  1:26760/135‑1 (MQ=255)
                                                                                              ctgcggcggGCCGCCCGGGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTg                                     <  4:138307/135‑1 (MQ=255)
                                                                                                            ccgcgcTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGc                       <  5:116872/133‑1 (MQ=255)
                                                                                                                                 aGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAGCGCTCGATCa  >  8:248208/1‑135 (MQ=255)
                                                                                                                                 |                                                                                                                                      
GTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAGCGCTCGATCA  >  ACB143/3440217‑3440480

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: ‑