Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1223,316,6550TG52.2% 9.0 / 19.8 23T99P (ACG→CCG) CJ019_03117hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (6/5);  new base G (0/12);  total (6/17)
Fisher's exact test for biased strand distribution p-value = 4.58e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GCTGCAAGTGCTCAGCCAGTTCCTCTGGCAACGGCACGAAGCGTACAAACGCGTGCAAGTGATGCGCCTCGCGGGTTACCTGCTTGATGCGCCGTTGCAGCTCACTGCCCAGCCGGTCACCGGCCAGCATTGCCGTGCGGTCGCCATGGGCGACGCGCCACAGCACTTCATACAGCAGGTTCCAGCGCTGCTCGCCACGGTAGCGAGAGGCCTGCTCCAACTGGCTCAGCAGTGCGGCCGGG  >  ACB122/3316520‑3316761
                                                                                                                                        |                                                                                                          
gCTGCAAGTGCTCAGCCAGTTCCTCTGGCAACGGCACGAAGCGTACAAACGCGTGCAAGTGATGCGCCTCGCGGGTTACCTGCTTGATGCGCCGTTGCAGCTCACTGCCCAGCCGGTCACCGGCCAGCATTGCCGTg                                                                                                           >  3:127099/1‑137 (MQ=255)
 cTGCAAGTGCTCAGCCAGTTCCTCTGGCAACGGCACGAAGCGTACAAACGCGTGCAAGTGATGCGCCTCGCGGGTTACCTGCTTGATGCGCCGTTGCAGCTCACTGCCCAGCCGGTCACCGGCCAGCATTGCCGTGc                                                                                                          <  6:41212/137‑1 (MQ=255)
                    tCCTCTGGCAACGGCACGAAGCGTACAAACGCGTGCAAGTGATGCGCCTCGCGGGTTACCTGCTTGATGCGCCGTTGCAGCTCACTGCCCAGCCGGTCACCGGCCAGCATTGCCGGGCAGTCGCCATGGGCGAcgcg                                                                                       <  6:223292/137‑1 (MQ=255)
                      ctctGGCAACGGCACGAAGCGTACAAACGCGTGCAAGTGATGCGCCTCGCGGGTTACCTGCTTGATGCGCCGTTGCAGCTCACTGCCCAGCCGGTCACCGGCCAGCATTGCCGTGCGGGGGCCCTGGGCGCCCCGcc                                                                                     >  8:246719/1‑137 (MQ=255)
                      ctctGGCAACGGCACGAAGCGTACAAACGCGTGCAAGTGATGCGCCTCGCGGGTTACCTGCTTGATGCGCCGTTGCAGCTCACTGCCCAGCCGGTCACCGGCCAGCATTGCCGGGCGGGCGCCATGGGCGACGCGcc                                                                                     <  6:284133/137‑1 (MQ=255)
                              aCGGCACTAAGGGTACAACGGCGTGAATGTGATGCGCCTCGGGGGTTACCTGCTTGATGCGCCTTTGCAGGGCACTCCCCCGCGGGTCACCGGCCTGCATTGCCGGGGGGGCGCCTTGGGCGACGCGCCACAGCACt                                                                             <  6:310411/137‑1 (MQ=255)
                                                      gCAAGTGATGCGCCTCGCGGGTTACCTGCTTGATGCGCCGTTGCAGCTCACTGCCCAGCCGGTCACCGGCCAGCATTGCCGTGGGGTAGCCCGCGGGGCCCCCCCCCCCCACTTCATAAACCACGTTACAgcgcata                                                     >  7:83725/1‑134 (MQ=255)
                                                         aGTGATGCGCCCCGCGGGTTACCTGCTTGATGCGCCGTTGCAGCTCACTGCCCAGCCGGTCACCGGCCAGCATTGCCGTGCGGGCGCCATGGGCGACGCGCCACAGCACTTCATACAGCAGGTTCCAGCGCTGCTCg                                                  <  5:218989/137‑1 (MQ=255)
                                                         aGTGAGGCGCCTCGCGGGTTACCTGCTTGATGCGCCGTGGCAGCTCACTCCCCCGCCGGGCCCCGGCCAGCATTGCCGGGCGGGCGCCATGGGCGACGCGCCACAGCACTTCATACAGCAGGTTCCAGCGCTGCTCg                                                  <  2:53215/137‑1 (MQ=255)
                                                           tGATGCGCCTCGGGGGTTACCTGCTTGATGCGCCGTTGAAGCTCACTGCCCCGCCGGGCACCGGCCAGCTTTGCCGGGGGGTCGCCATGGGCTACGCGCCACAGCACTTCATACAGCAGGTTCCAGCGCTGCTCGcc                                                <  5:122943/137‑1 (MQ=255)
                                                             aTGCGCCTCGCGGGTTACCTGCTTGATGCGCCGTTGCAGCTCACTGCCCAGCCGGTCACCGGCCAGCATTGCCGTGCGGGCGCCATGGGCGACGCCCCCCACCCCTTCAGACAGCACGTTCCAGCGCTGCTCGCCAc                                              >  8:271353/1‑137 (MQ=255)
                                                                  ccgcgggggTTACCTGCTTGATGCGCCGTTGCAGCTCACTGCCCAGCCGGTCCCCGGCCAGCATTGCCGGGGGGGCGCCATGGGCGACGCGCCACAGCACTTCATACAGCAGGTTCCAGCGCTGCTCGCCACGGTAg                                         <  8:193170/131‑1 (MQ=255)
                                                                     cgcgGGTTACCTGCTTGATGCGCCGTTGCAGCTCACTGCCCAGCCGGTCACCGGCCAGCATTGCCGTGGGGGGGCCCGGGGGGGCGCGCCCCAGCCCGGGATACAGCCGGTTGCAGCGCGGCTCGCCACGGTGGCgg                                      >  7:299535/1‑136 (MQ=255)
                                                                                 gCTTGTGGCGCGGTTGCAGCTCAGGCCCCCGCGGGTCCCCGCCCAGCTTGCCGGGGGGGGGGCCATGGGCGACGCGCCACAGCACTTCATACAGCAGGTTCCAGCGCTGCTCGCCACGGTAGCGAGAGGCCTGCTcc                          <  3:54407/137‑1 (MQ=255)
                                                                                 gCTTGAGGCGCCGTGGCAGCTCACTCCCCAGCCGGTCACGGGCCAGCATTGCCGTGCGGTCGCCATGGGCGACGCGCCACAGCACTTCATACAGCAGGTTCCAGCGCTGCTCGCCACGGTAGCGAGAGGCCTGCTcc                          <  7:271353/137‑1 (MQ=255)
                                                                                             gTTGCAGCTCACTGCCCCGCCGGTCACCGCCCAGCATGGCCGGGCGGTCGCCATGGGCGACGCGCCACAGCACTTCATACAGCAGGTTCCAGCGCTGCTCGCCACGGTAGCGAGAGGCCTGCTCCAACTGGCTcagc              <  7:246719/137‑1 (MQ=255)
                                                                                              ttGCAGCTCACTGCCCAGCCGGTCACCGGCCAGCATTGCCGTGCGGTCGCCATGGGCGACGCGCCACAGCACTTCATACAGCAGGTTCCAGCGCTGCTCGCCACGGTAGCGAGAGGCCTGCTCCAACTGGCTcagca             <  4:222951/137‑1 (MQ=255)
                                                                                              ttGCAGCTCACTGCCCAGCCGGCCACCGGCCAGCATTGCCGGGCGGTCGCCATGGGCGACGCGCCACAGCACTTCATACAGCAGGTTCCAGCGCTGCTCGCCACGGTAGCGAGAGGCCTGCTCCAACTGGCTcagca             <  1:20198/137‑1 (MQ=255)
                                                                                                    cccACTCCCCCGCCGTGCCCCGCCAAGCATGGCCGTGCGGTCCCCAGGGCCTACGCGCCACAGCACTTCATACAGCAGGTTCCAGCGCTGCTCGCCACGGTAGCGAGAGGCCTGCTCCAACTGGCTCAGCAGTGCgg       <  2:56468/135‑1 (MQ=255)
                                                                                                    cTCACCGCCCCCCCGTCCACCGGCCAGCATTGCCGGGCGGGCCCCAGGGGCGACGCGCCACAGCACTTCATACAGCAGGTTCCAGCGCTGCTCGCCACGGTAGCGAGAGGCCTGCTCCAACTGGCTCAGCAGTGCgg       <  2:60995/137‑1 (MQ=255)
                                                                                                       aCTGCCCAGCCGGTCACCGGCCAGCATTGCCGTGGGGGGGCCCTGGGGGCCGCGCCCCCCCCCATCATACACCACGTTGCCGCGCGGCGCGCCCCCGGGGCGCGAGGGCGGCGCCCCCTGGCGCAGCCGTGCgggcc    >  2:56597/1‑134 (MQ=255)
                                                                                                         tGCCCCGCCGGTCACCGGCCTGCATGGCCGGGGGGGCGCCTGGGGCGACGCGCCACAGCACTTCATACAGCAGGTTCCAGCGCTGCTCGCCACGGTAGCGAGAGGCCTGCTCCAACTGGCTCAGCAGTGCGGCCggg  <  4:22176/137‑1 (MQ=255)
                                                                                                         ggCCCCGCCGGTCACCGGCCAGCATTGCCGGGCGGTCGCCATGGGCGACGCGCCACAGCACTTCATACAGCAGGTTCCAGCGCTGCTCGCCACGGTAGCGAGAGGCCTGCTCCAACTGGCTCAGCAGTGCGGCCggg  <  8:82288/136‑1 (MQ=255)
                                                                                                                                        |                                                                                                          
GCTGCAAGTGCTCAGCCAGTTCCTCTGGCAACGGCACGAAGCGTACAAACGCGTGCAAGTGATGCGCCTCGCGGGTTACCTGCTTGATGCGCCGTTGCAGCTCACTGCCCAGCCGGTCACCGGCCAGCATTGCCGTGCGGTCGCCATGGGCGACGCGCCACAGCACTTCATACAGCAGGTTCCAGCGCTGCTCGCCACGGTAGCGAGAGGCCTGCTCCAACTGGCTCAGCAGTGCGGCCGGG  >  ACB122/3316520‑3316761

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: