Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1223,440,3670TG53.3% 5.2 / 21.5 31W26G (TGG→GGG) CJ019_03266hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (10/4);  new base G (16/0);  total (27/4)
Fisher's exact test for biased strand distribution p-value = 3.65e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.44e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

AGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAGCGCTCGATCAGCTGGC  >  ACB122/3440236‑3440503
                                                                                                                                   |                                                                                                                                        
aGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCCCCAGCGCTGGCGGCGCCCCGGCAGCCTGGgtg                                                                                                                                     >  6:244107/1‑137 (MQ=255)
 gCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGgtgt                                                                                                                                    >  1:8938/1‑137 (MQ=255)
 gCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCCCAGGCGGCGCCCGGGCCGCCGGGgtgt                                                                                                                                    >  7:91918/1‑137 (MQ=255)
                    cAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGCCAACCCCTaaaa                                                                                                                 >  3:121686/1‑137 (MQ=255)
                     aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCCGGCGGCCCCCGCCCAGCCTGGGTGTGGCAGGCCAACCCCAAAAAg                                                                                                                >  3:84488/1‑137 (MQ=255)
                                 aaGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGCTCGCCGTGGCGCTGCTGCTGGCCCCCCCCGCTGCCGCCGCCCGGGCAGCCTCGGTGTGGCAGGCCCACGCCTATAAGGCCCAGCTcccg                                                                                                    >  7:56501/1‑134 (MQ=255)
                                      agcgCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCTGGGTGTTGCAGGCCAACGCCGATAAGGCGCAGCGAGCTGATCa                                                                                               >  4:249985/1‑137 (MQ=255)
                                         gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCGGGGTGTGGCAGGCCAACGCCGATAATGCGCAGCTCGCTGATCAGGc                                                                                            >  3:160233/1‑137 (MQ=255)
                                         gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCGGGGGGTGGCAGGCCCACGCCTATAATGCGCAGCGCGCTGATCAGGc                                                                                            >  6:318005/1‑137 (MQ=255)
                                         gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCTTGCGGCGCCCCGGCAGCCTGGGTTGGGCCGCCCAACCGCTAAAAGGCGCAGCTCCCAGAtaacgc                                                                                            >  6:73827/1‑132 (MQ=255)
                                         gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCGGGGGGCGCCCGGGCCGCCGGGGTGTGGGAGGCCCACGCCGATAAAGCGCAGCGCGCGGGTCAGGc                                                                                            >  3:35847/1‑137 (MQ=255)
                                         gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCACCGCTGGCGGCGCCCGGGCAGCCGGGGTGTGGGCGGCCAACGCCGTTAAGGCGCAGCTCGCTGATAAggg                                                                                            >  2:250170/1‑136 (MQ=255)
                                         gcTCGAGGTGCTGCAGTGAGCGCTTGGGCAACCCGCCTGGTCGCCGTGGCGCTGCTGATGGCCGCCACCGCTCGAGGCGCCCGGCCACCCGGGGTGTGGCACGGCACACCCTAGAAGGCGCCGCGCGCTGATGAGGc                                                                                            >  5:298829/1‑137 (MQ=255)
                                           tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCTGGGGGTGGCCGGCCCACGCCGATAAGGCGCCGCGCGCTGCTCAGGCCg                                                                                          >  7:258674/1‑137 (MQ=255)
                                           tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGGGGCGCCCGGGCCGCGGGGGGGGGGCCGGCCAACGCCGATAAGGCGCCGCGCGGCGCTGAGGCCg                                                                                          >  8:303070/1‑137 (MQ=255)
                                           tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCCGCCGGGGGGTGGCAGGCCACCGCCGAGAAGGCGCAGCTCGCGGATCAGGCCg                                                                                          >  3:323987/1‑137 (MQ=255)
                                           tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCCCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCg                                                                                          >  6:138734/1‑137 (MQ=255)
                                                      cAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCCCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGTTTGGCAGGCCAACGCCTATAAGGCGCAGCGCGCGGGTAAGGCCGAGGCCTTgggg                                                                               >  4:230623/1‑136 (MQ=255)
                                                              gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGGCGCCCGCCCTGGCGGCGCCCGGGCAGCCAGGGGGGGGGAGGCCAACCCCAATAAGGCGCAGCTCGCTGATCACGCCGAAGAAGTCGGCCAACAgat                                                                       >  4:91854/1‑135 (MQ=255)
                                                              gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGGGCCCCGGCCGCCGGGGTGTGGCAGGCCCACCCCGATTAGGCGCAGCGCGCGGATCAGGCCGATGACGTGGGGAAGCAggt                                                                       >  4:261693/1‑135 (MQ=255)
                                                              gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCCGCCGGGGGGTGGCAGGGCCACGCCTATAATGCGCAGCTCGCGGATCAGGCCGCTGACTTCGGCCAGCCGtt                                                                       >  2:321241/1‑137 (MQ=255)
                                                                   ggCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGCCAACCCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTAGGCAACCAGTTGgggg                                                                  >  2:49053/1‑137 (MQ=255)
                                                                            ccTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGGCCGGGCCGGCGGGGGGGGGCGGGCCCCCGCCGTTAAGTGGCCGCGGGCGGCTCAGGCGGCTGAGATCGGCAAGCAGTGGGCGGAGAAGGAtt                                                         >  8:25503/1‑136 (MQ=255)
                                                                            ccTGGTCGCCGTGGCGCTGCTGCTGGCCCCCCCCGCGGGCGGCGCCCCGGCACCGGGGGTGTTGGAGGGCAACCCCGTAAAGTCGCAGCTGGCTGCTCAGGCCGATTACGTAGGCCACATGTTGGCGGGGCAgggtg                                                         >  6:207146/1‑134 (MQ=255)
                                                                                gTCGCCGTGGGGGTGCGGCTCGCCGCGAGGGGGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc                                                     <  5:118097/137‑1 (MQ=255)
                                                                                                     cgcgccAGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGaa                                <  5:351756/135‑1 (MQ=255)
                                                                                                            gcgcGGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGCTCAGGCCGCTGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCCCGGTCTGGAGGGGGAGAAGGCTgcc                         >  1:191422/1‑136 (MQ=255)
                                                                                                             ggcgggcggcGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTgcag                        <  2:8938/133‑1 (MQ=255)
                                                                                                                    ggcgCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCTGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGAAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCagcagcagc                 >  3:134393/1‑137 (MQ=255)
                                                                                                                     gcgcCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAgcg                >  8:262342/1‑137 (MQ=255)
                                                                                                                           gggCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAGCGCTCGAt          <  4:35847/137‑1 (MQ=255)
                                                                                                                                   tGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAGCGCTCGATCAgctggc  <  8:259222/137‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                        
AGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAGCGCTCGATCAGCTGGC  >  ACB122/3440236‑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: ‑