Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB122790,0070TG53.8% 1.5 / 15.4 26intergenic (+159/‑31)CJ019_00751/CJ019_00752hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (7/5);  new base G (14/0);  total (21/5)
Fisher's exact test for biased strand distribution p-value = 1.20e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.38e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCTTTTTTAGTCACCTGATGAGTCCTGAACTTCATGGAACGTATCCTCGAAAACGCGATGTATGCCTCGCGCTGGCTGCTCGCACCGATCTACTTCGGCTTGTCCCTCGGCCTGCTGGCCCTGGCGCTGAAATTTT  >  ACB122/789895‑790140
                                                                                                                |                                                                                                                                      
ccACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCTTTTTTAGTGACCTGATGCGTcc                                                                                                                  >  2:71439/1‑135 (MQ=255)
     ggTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCTTTTTTAGTCACCTGATGAGTCCTGAAc                                                                                                            >  2:235130/1‑135 (MQ=255)
             aTCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCGTTTTTAGTCACCTGATGAGTACTGAACGTCaatg                                                                                                     >  3:129307/1‑131 (MQ=255)
               cGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCACCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCGTTTTTAGTAAACTGATGAGTCCTGCACTTGATTGAAc                                                                                                  >  7:318974/1‑135 (MQ=255)
                 gCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCTTTTTTAGTCACCTGATGAGTCCTGATCTTCATGGGACgg                                                                                                >  6:107879/1‑134 (MQ=255)
                  cTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCTTTTTTAGTGCCCTGATGAGGCCTGCACTTGATGGacccga                                                                                               >  2:46334/1‑130 (MQ=255)
                                    ccTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCTTTTTTTGTGACCAGATGATTGCTGAACATCATTGAAAGTATACTCCAAAACCCGATGt                                                                             >  8:346294/1‑135 (MQ=255)
                                    ccTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCGTTTTTTTTCACCCGATGAGTCCTGCACGTCATGGAACGGATACTAGAAAACGCGAttt                                                                             >  7:265512/1‑133 (MQ=255)
                                      ttGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCTTTTTTATTCACCTGATGTGTCCGGAACTTGATTGAACGTTTACTCGAAAACGCGATGttt                                                                           >  8:208251/1‑133 (MQ=255)
                                      ttGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCGTTTTTAGTCACCGGATGAGTCCGGAACATCCAGGCACCTATCCGCCAAAACCCCCTGTAt                                                                           >  8:213767/1‑135 (MQ=255)
                                         tgtCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCGTTTTTTGTCAACGGATGAGTCGTGAAAGTCATGGAACGTATCCTAGAAAACCCGATGTATGcc                                                                        >  4:138323/1‑135 (MQ=255)
                                           tCGCGAACAGCCTGCGGAACAGCCACCGGGAGATTGCTGTAACAGGTGGCTCTCTATACTCGCCGTCGCGTTTTTATACACCAAAACAATCATGAACAACATGGAAACTATCCACCAAAACACGATGTaaaccac                                                                      >  4:91757/1‑129 (MQ=255)
                                            cgcgATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCGTTTTTTTTGACCAGATTAGTCCAGGAAGTCAATGGAAGTTTACACGAAAAAGCCATGTATggccgg                                                                     >  8:346371/1‑130 (MQ=255)
                                              cgttctggCTGGGGAGCCGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCTTTTTTAGTCACCTGATGAGTCCTGAACTTCATGGAACGTATCCTCGAAAACGCGATGTATGCCTcgcg                                                                   <  6:54637/129‑1 (MQ=255)
                                                 tccggCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCGGTTTTAGTGCCCGGTTGAGTCCTGGACTTCATGGAACGTATACTCGAAAACCCGATGTATTCCTACCgctg                                                                >  8:102111/4‑135 (MQ=255)
                                                 tCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCGTTTTTATTCACCGGATGGGTCCGGAACTTCATGGAACGTATCCTCGAAAACGCGATGTTTGCCGCGCgctg                                                                >  2:332203/1‑135 (MQ=255)
                                                             gcagcCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCTTTTTTAGTCACCTGATGAGTCCTGAACTTCATGGAACGTATCCTCGAAAACGCGATGTATGCCTCGCGCTGGCTGCTCGCAcc                                                    <  1:235130/135‑1 (MQ=255)
                                                             gcagcCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCTTTTTTAGTCACCTGATGAGTCCTGAACTTCATGGAACGTATCCTCGAAAACGCGATGTATGCCTCGCGCTGGCTGCTCGCAcc                                                    <  7:223318/135‑1 (MQ=255)
                                                                        gCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCTTTTTTAGTCACCTGATGAGTCCTGAACTTCATGGAACGTATCCTCGAAAACGCGATGTATGCCTCGCGCTGGCTGCTCGCACCGATCTACTTCg                                         <  2:239989/135‑1 (MQ=255)
                                                                              gCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCTTTTTTAGTCACCAGATGAGTCCGGAAATTCATGGAACGTATCCTAGAAAACGCGATGTATGCCTCCAGCTGGCTGCTCGCACCGATCTACTTCGGCgtgt                                   >  8:122883/1‑135 (MQ=255)
                                                                              gCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCGTTTTTTGTAAACAGATGAGTCCTGAACGTAATGGAAAGTATCCTAGAAAACGCGAGTTATGCCTCCCGCTGGCTGCGCGCCCCCCTCTACTTGGGCgtgt                                   >  6:205566/1‑135 (MQ=255)
                                                                              gCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCGTTTTTAGGCCCCTGATGAGTCCTGAACTTCATGGCACGGATACTCGAAAACGCGATGTATGCCTCGCGCTGGCGGCTCGCCCCCATCGACTTAGGCTTGt                                   >  1:261022/1‑135 (MQ=255)
                                                                                       ggTGGCTGGCTATACTCGCCGTCGCGTTTTTATTCACCTGATGAGTCCTGAACTTCATGGAACGTATCCTCGAAAACGCGATGTATGCCTCGCGCTGGCTGCTCGCACCGATCTACTTCGGCTTGTCCCTCGGcc                          >  6:352899/1‑135 (MQ=255)
                                                                                            ctggctATACTCGCCGTCGCGTTTTTTGTCCACCGGTGAGTCCTGCACATGATGGTACCTATCCGGGCAAACCCGGTTGATGCCGCGCGCGGGCTGCTGGCACCGGTATACATCGGCGTGTCCCACGGGcagctg                     >  2:87574/1‑135 (MQ=255)
                                                                                               gctATACTCGCCGTCGCGTTTTTAGTCACCAGATGAGTCCTGGACTTCATGGGACGGTTGCACGAAAACGCGATGTATGCCTCGCGCGGGCTGCGCGCACCGCTATACTTCGGCGTGTCCCAcggccggcgtgcc                  >  4:125699/1‑130 (MQ=255)
                                                                                                               cTTTTTTAGTCACCTGATGAGGCCTGAACTTCATGGAACGTCTCCTCGAAAACGCGATGTATGCCTCGCGCTGGCTGCTCGCACCGATCTACTTCGGCTTGTCCCTCGGCCTGCTGGCCCTGGCGCTGAAAtttt  <  7:334051/135‑1 (MQ=255)
                                                                                                                |                                                                                                                                      
CCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCTTTTTTAGTCACCTGATGAGTCCTGAACTTCATGGAACGTATCCTCGAAAACGCGATGTATGCCTCGCGCTGGCTGCTCGCACCGATCTACTTCGGCTTGTCCCTCGGCCTGCTGGCCCTGGCGCTGAAATTTT  >  ACB122/789895‑790140

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: