Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,373,1590TC71.4% 26.9 / 14.4 21L340S (TTA→TCA) ycjTkojibiose phosphorylase
Reads supporting (aligned to +/- strand):  ref base T (2/4);  new base C (8/7);  total (10/11)
Fisher's exact test for biased strand distribution p-value = 6.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CGTTAGATTATGCGCTTTATCATCTGCGCATCATGACGCCTGCCCACGACGAGCGCAGCAGCATTGCGGCAAAAGGCTTAACCGGCGAAGGCTACAAAGGCCACGTTTTCTGGGATACAGAAGTATTTTTGTTACCGTTTCATCTGTTTAGCGATCCGACGGTTGCCCGAAGTTTACTGCGTTATCGCTGGCACAACTTGCCAGGCGCGCAGGAGAAAGCGCGACGCAACGGCTGGCAGGGCGCGCTATTTCCGTGGGAAAGCGCG  >  NZ_CP009273/1373027‑1373292
                                                                                                                                    |                                                                                                                                     
ggTTAGATTATGCGCTTTATCATCTGCGCATCATGACCCTTGCCCACGACGAGGGCAGCAGCATGGGGGCGAACGCCTTAACCGGCGATGGCGACGAAGGCCAGGCTTGCAGGGACAGAGTATTATTTTTGTCACCGtt                                                                                                                                 >  2:143925/2‑139 (MQ=255)
  ttAGATTATGCGCTTTATCATCTGCGCATCATGACGCCTGCCCACGACGAGCGCAGCAGCATTGCGGCAAAAGGCTTAACCGGCGAAGGCTACAAAGGCCACGTTTTCTGGGATACAGAAGTATTTTTGtt                                                                                                                                       >  2:64838/1‑131 (MQ=255)
  ttAGATTATGCGCTTTATCATCTGCGCATCATGACGCCTGCCCACGACGAGCGCAGCAGCATTGCGGCAAAAGGCTTAACCGGCGAAGGCTACAAAGGCCACGTTTTCTGGGATACAGAAGTATTTTTGtt                                                                                                                                       <  1:64838/131‑1 (MQ=255)
     gATTATGCGCTTTATCATCTGCGCATCATGACGCCTGCCCACGACGAGCGCAGCAGCATTGCGGCAAAAGGCTTAACCGGCGAAGGCTACAAAGGCCACGTTTTCTGGGATACAGAAGTATTTTTGTTACCGTTTCAt                                                                                                                             <  2:473577/138‑1 (MQ=255)
     gATTATGCGCTTTATCATCTGCGCATCATGACGCCTGCCCACGACGAGCGCAGCAGCATTGCGGCAAAAGGCTTAACCGGCGAAGGCTACAAAGGCCACGTTTTCTGGGATACAGAAGTATTTTTGTTACCGTTTCAt                                                                                                                             <  1:51974/138‑1 (MQ=255)
     gATTATGCGCTTTATCATCTGCGCATCATGACGCCTGCCCACGACGAGCGCAGCAGCATTGCGGCAAAAGGCTTAACCGGCGAAGGCTACAAAGGCCACGTTTTCTGGGATACAGAAGTATTTTTGTCACCGTTTCAt                                                                                                                             <  2:800341/138‑1 (MQ=255)
                         gcgcATCATGACGCCTGCCCACGACGAGCGCAGCAGCATTGCGGCAAAAGGCTTAACCGGCGAAGGCTACAAAGGCCACGTTTTCTGGGATACAGAAGTATTTTTGTCACCGTTTCATCTGTTTAGCGATCCGACGGtt                                                                                                        >  1:110269/1‑139 (MQ=255)
                         gcgcATCATGACGCCTGCCCACGACGAGCGCAGCAGCATTGCGGCAAAAGGCTTAACCGGCGAAGGCTACAAAGGCCACGTTTTCTGGGATACAGAAGTATTTTTGTCACCGTTTCATCTGTTTAGAGATCCGACggat                                                                                                        >  2:301245/1‑137 (MQ=255)
                              tcatGACGCCTGCCCACGACGAGCGCAGCAGCATTGCGGCAAAAGGCTTAACCGGCGAAGGCTACAAAGGCCACGTTTTCTGGGATACAGAAGTATTTTTGTTACCGTTTCATCTGTTTAGCGATCCGACGGTTGCCCg                                                                                                   <  1:81670/139‑1 (MQ=255)
                                    cGCCTGCCCACGACGAGCGCAGCAGCATTGCGGCAAAAGGCTTAACCGGCGAAGGCTACAAAGGCCACGTTTTCTGGGATACAGAAGTATTTTTGTTACCGTTTCATCTGTTTAGCGATCCGACGGTTGCCCGAAGtt                                                                                              >  2:772105/1‑138 (MQ=255)
                                                   aGCGCAGCAGCATTGCGGCAAAAGGCTTAACCGGCGAAGGCTACAAAGGCCACGTTTTCTGGGATACAGAAGTATTTTTGTCACCGTTTCATCTGTTTAGCGATCCGACGGTTGCCCGAAGTTTACTGCGTTATCGCTg                                                                              >  1:640647/1‑139 (MQ=255)
                                                                  cGGCAAAAGGCTTAACCGGCGAAGGCTACAAAGGCCACGTTTTCTGGGATACAGAAGTATTTTTGTCACCGTTTCATCTGTTTAGCGATCCGACGGTTGCCCGAAGTTTACTGCGTTATCGCTGGCACAACTTGCCAgg                                                               >  1:309887/1‑139 (MQ=255)
                                                                       aaaGGCTTAACCGGCGAAGGCTACAAAGGCCACGTTTTCTGGGATACAGAAGTATTTTTGTCACCGTTTCAt                                                                                                                             <  1:852673/72‑1 (MQ=255)
                                                                       aaaGGCTTAACCGGCGAAGGCTACAAAGGCCACGTTTTCTGGGATACAGAAGTATTTTTGTCACCGTTTCAt                                                                                                                             >  2:852673/1‑72 (MQ=255)
                                                                            cTTAACCGGCGAAGGCTACAAAGGCCACGTTTTCTGGGATACAGAAGTATTTTTGTCACCGTTTCATCTGTTTAGCGATCCGACGGTTGCCCGAAGTTTACTGCGtt                                                                                     >  2:187851/1‑107 (MQ=255)
                                                                            cTTAACCGGCGAAGGCTACAAAGGCCACGTTTTCTGGGATACAGAAGTATTTTTGTCACCGTTTCATCTGTTTAGCGATCCGACGGTTGCCCGAAGTTTACTGCGtt                                                                                     <  1:187851/107‑1 (MQ=255)
                                                                               aaCCGGCGAAGGCTACAAAGGCCACGTTTTCTGGGATACAGAAGTATTTTTGTCACCGTTTCATCTGTTTAGCGATCCGACGGTTGCCCGAAGTTTACTGCGTTATCGCTGGCACAACTTGCCAGGCGCGCAGgaga                                                    <  2:110269/137‑1 (MQ=255)
                                                                                                 aGGCCACGTTTTCTGGGATACAGAAGTATTTTTGTCACCGTTTCATCTGTTTAGCGATCCGACGGTTGCCCGAAGTTTACTGCGTTATCGCTGGCACAACTTGCCAGGCGCGCAGGAGAAAGCGCGACGCAACGGCTgg                                <  2:309887/139‑1 (MQ=255)
                                                                                                    ccACGTTTTCTGGGATACAGAAGTATTTTTGTCACCGTTTCATCTGTTTAGCGATCCGACGGTTGCCCGAAGTTTACTGCGTTATCGCTGGCACAACTTGCCAGGCGCGCAGGAGAAAGCGCGACGCAACGGCTGGCAg                             >  1:639695/1‑139 (MQ=255)
                                                                                                      aCGTTTTCTGGGATACAGAAGTATTTTTGTTACCGTTTCATCTGTTTAGCGATCCGACGGTTGCCCGAAGTTTACTGCGTTATCGCTGGCACAACTTGCCAGGCGCGCAGGAGAAAGCGCGACGCAACGGCTGGCAggg                           <  1:139420/139‑1 (MQ=255)
                                                                                                      aCGTTTTCTGGGATACAGAAGTATTTTTGTCACCGTTTCATCTGTTTAGCGATCCGACGGTTGCCCGAAGTTTACTGCGTTATCGCTGGCACAACTTGCCAGGCGCGCAGGAGAAAGCGCGACGCAACGGCTGGCAggg                           <  1:730910/139‑1 (MQ=255)
                                                                                                                ggATACAGAAGTATTTTTGTCACCGTTTCATCTGTTTAGCGATCCGACGGTTGCCCGAAGTTTACTGCGTTATCGCTGGCACAACTTGCCAGGCGCGCAGGAGAAAGCGCGACGCAACGGCTGGCAGGGCGCGCTAt                   <  2:181707/137‑1 (MQ=255)
                                                                                                                               tttGTTACCGTTTCATCTGTTTAGCGATCCGACGGTTGCCCGAAGTTTACTGCGTTATCGCTGGCACAACTTGCCAGGCGCGCAGGAGAAAGCGCGACGCAACGGCTGGCAGGGCGCGCTATTTCCGTGGGAAAgcgcg  >  2:492514/1‑139 (MQ=255)
                                                                                                                                    |                                                                                                                                     
CGTTAGATTATGCGCTTTATCATCTGCGCATCATGACGCCTGCCCACGACGAGCGCAGCAGCATTGCGGCAAAAGGCTTAACCGGCGAAGGCTACAAAGGCCACGTTTTCTGGGATACAGAAGTATTTTTGTTACCGTTTCATCTGTTTAGCGATCCGACGGTTGCCCGAAGTTTACTGCGTTATCGCTGGCACAACTTGCCAGGCGCGCAGGAGAAAGCGCGACGCAACGGCTGGCAGGGCGCGCTATTTCCGTGGGAAAGCGCG  >  NZ_CP009273/1373027‑1373292

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