Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB143994,2350TG52.0% 6.3 / 19.2 25R54R (CGT→CGGCJ019_00938hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (7/5);  new base G (0/13);  total (7/18)
Fisher's exact test for biased strand distribution p-value = 1.65e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.07e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TGGATGTGGCCGGCTTGCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTCGGGCGTTGCCTGGAGCATGGCTGGGGTGGCGAGGTGAATTTGGTGC  >  ACB143/994105‑994365
                                                                                                                                  |                                                                                                                                   
tggATGTGGCCGGCTTGCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTgggct                                                                                                                                >  4:327240/1‑133 (MQ=255)
       ggCCGGCTTGCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCCTCcagcag                                                                                                                         >  1:121973/1‑136 (MQ=255)
             cTTGCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCCTCCCGCAGGCTGcc                                                                                                                  >  6:76937/1‑136 (MQ=255)
                   ggCCAGGCCCGAGGACAGCCCCGGAAAGGCCGCCCATGCGATTCTGGCGGCAGCGGGGCTGGGTGTGTTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGTTGGGCGTGGCATCCAGCTGGATGCCACGGtt                                                                                                            <  8:241649/136‑1 (MQ=255)
                          cgcGAGCAAACCCCCGGCACGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGGCGCGGCGCAGTTGCTGCGGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCgg                                                                                                      <  2:51531/133‑1 (MQ=255)
                                 aaccccccGGCAAGGCGGCCCCTGCTTTTGTGGGGGCCGGGGGGCAGGGGGTGGCGGCGGCGCTGTGGGTGGGGGGGAATACCCGTGGGGTGGGGCGGGGACCCCAGCTGGATGCCACGGTTGCCCGGCGTTGGtt                                                                                              <  8:318772/131‑1 (MQ=255)
                                     tCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATTCCACGGTTGCCCGGCGTTGGGTCggg                                                                                          >  2:128624/1‑135 (MQ=255)
                                      cccGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGGCTCCCGCAGGGTGACCCCGGTGGCCCGgggtgggtggggt                                                                                          >  7:182301/1‑131 (MQ=255)
                                         ggCACGGCCGCCCAGGCGATGGTGGCGGCAGCGGGGCAGGGGGTGGTCGTGGCGCAGTTGCTGGTGGGGCTGATCCTGCTGGATGGGCGGGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCg                                                                                      <  1:23808/136‑1 (MQ=255)
                                           cccggCCGCCCCGGCTATTCTGGCGGCTGGGGGGCAGGGTGGGGTCGTGGCGCTGTGGCTGCGGGGGCATAGCGTGGGGGTGGGGCGGGGCATCCTGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCt                                                                                    <  1:72311/133‑1 (MQ=255)
                                           cacggCCGCCCCGCGGTTGCTGGGGGCAGCGGGGCTGGGTGGGGTCGCGCCGCAGTTGCTGCGGGGGCATATCCTGGTGGATGGGCGGGGCATCCTGCCGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCt                                                                                    <  6:53937/133‑1 (MQ=255)
                                                ccgccCAGGCGATTCGGGCGGCAGGGGGGCTGGGGGTGGGCGAGGCGCAGTTGCTGCTGGGGCATAGCCTGCTGGATGGGCGGGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGcaca                                                                               <  3:354151/136‑1 (MQ=255)
                                                      aGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGGCTCCAGCACGCTGCCCCCGGTGGCCGGGGGTGGGTGGGGGTCGGGGCCCAGGGTgg                                                                         >  8:263736/1‑136 (MQ=255)
                                                                ggcgggAGGGGGGCAGGGTGGGGGGGGGGGGCCGGGCGTGGTGGGTCAGATCCTGCGGGTTGGGCGTGGCATCCCGAAGCATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATg                                                               <  3:48216/136‑1 (MQ=255)
                                                                  gggcAGGGGGGCAGGGTGGGGTCGAGGCGCAGTGGCTGCTGGGGCAGAGCCTGCTGGATGGGCGGGGCACCCAGCCGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGgc                                                             <  4:22093/135‑1 (MQ=255)
                                                                    gcAGCGGGTCAGGGTTGGGTCGTGGCGCTGTTGCTGGTGGGGAATATCCTGGTGTCTGGGCGGGGCAGCCAGAAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGgcgc                                                           <  4:60285/136‑1 (MQ=255)
                                                                      aggggggCCGGGTGGGGCCGCGCCGCCGTGGCGGCTGGGGCAGATCGTGGTGGTTGGGCGGGCCACCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGcac                                                         <  3:173487/133‑1 (MQ=255)
                                                                                 gggtgGGCGCGGCGCAGTTGGGGCTGGGTCAGATCCGGTTGTAGGGGGGTGCCATACCGCTGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTCgg                                              <  7:331934/134‑1 (MQ=255)
                                                                                   gtgGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATTCCACGGTTGCCCGGCGTTGGTTGGGTATCGCTTCACAAGGTGGCAGTGGCATGGCGCACAACATGCTCGGGc                                            >  8:261143/1‑136 (MQ=255)
                                                                                      ggcGAGGCGCAGTTGCGGCTGGGGCAGATCCTGCGGGAGGGGCGGGGCACCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTCGGGCGtt                                         <  2:121973/134‑1 (MQ=255)
                                                                                        cgcgccgcATTGGCGGCGGGGGCAGTTCCTGCTGGTTGGGGGGGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTCGGGCGTTGc                                       <  6:267033/131‑1 (MQ=255)
                                                                                                            ggCAGATCCTGCTGGTTGGGGGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTCGGGCGTTGCCTGGAGCATGGCTGGGGTgg                   <  7:261143/136‑1 (MQ=255)
                                                                                                             gcttaTCCTGCTGTTTGGGCGGGGCAGCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTCGGGCGTTGCCTGGAGCATGGCTGGGGTGGc                  <  6:283623/132‑1 (MQ=255)
                                                                                                                      gctgGAGGGGCGGGGCATCCAGAGGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTCGGGCGTTGCCTGGAGCATGGCTGGGGTGGCGAGGTGAAt         <  4:32575/136‑1 (MQ=255)
                                                                                                                             gggCGGGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTCGGGCGTTGCCTGGAGCATGGCTGGGGTGGCGAGGTGAATTTGGTGc  <  3:113653/136‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                   
TGGATGTGGCCGGCTTGCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATGGCGCACAACATGCTCGGGCGTTGCCTGGAGCATGGCTGGGGTGGCGAGGTGAATTTGGTGC  >  ACB143/994105‑994365

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: