Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273832,0420CG49.0% 42.1 / 23.0 49A79G (GCG→GGG) hcxBhydroxycarboxylate dehydrogenase HcXB
Reads supporting (aligned to +/- strand):  ref base C (9/16);  new base G (24/0);  total (33/16)
Fisher's exact test for biased strand distribution p-value = 8.30e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.77e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GGCATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCAC  >  NZ_CP009273/831957‑832130
                                                                                     |                                                                                         
ggCATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGggggcg                                                                                       >  2:1446016/1‑90 (MQ=255)
     gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGggggggggggggcc                                                                                  >  1:1573334/1‑86 (MQ=255)
     gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGggggggggggggac                                                                                  >  1:154990/1‑86 (MQ=255)
     gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGgggggggggggccc                                                                                  >  1:782067/1‑86 (MQ=255)
     gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGgggcgggtac                                                                                  >  1:1308596/1‑86 (MQ=255)
     gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGgggcggggac                                                                                  >  2:1361958/1‑86 (MQ=255)
     gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGgggcggggac                                                                                  >  1:657452/1‑86 (MQ=255)
     gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGgggcgggccc                                                                                  >  2:840698/1‑86 (MQ=255)
     gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGcggcgggccc                                                                                  >  2:596049/1‑86 (MQ=255)
     gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCAc                                                                                  >  2:991008/1‑90 (MQ=255)
     gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCAc                                                                                  >  2:1354316/1‑90 (MQ=255)
     gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCAc                                                                                  >  2:557538/1‑90 (MQ=255)
     gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGGCAc                                                                                  >  2:1612956/1‑90 (MQ=255)
                gtGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCACGCTCGATGGCg                                                                       >  2:693419/2‑90 (MQ=255)
                aTGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCgggggggg                                                                                         >  1:589924/1‑72 (MQ=255)
                aTGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGggcgg                                                                                         <  2:589924/72‑1 (MQ=255)
                aTGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCACGCCCCATGGCg                                                                       >  1:1542621/1‑90 (MQ=255)
                aTGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCACGCCCCCTGGCg                                                                       >  1:1870049/1‑90 (MQ=255)
                       cTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGGACACGCGATGGCGGTCgcgc                                                               >  1:1880066/1‑90 (MQ=255)
                       cTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCACGCCCCCTGGCGATCgcgc                                                               >  1:1847223/1‑90 (MQ=255)
                       cTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGGCACGCTCGATGGCGATCgcgc                                                               >  2:1278055/1‑90 (MQ=255)
                        tCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGGGGGCACGCTCGATGGCGATCgcgcc                                                              >  2:701830/1‑89 (MQ=255)
                          cTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCgggggggg                                                                                         >  1:1702689/1‑62 (MQ=255)
                          cTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGggcgg                                                                                         <  1:275479/62‑1 (MQ=255)
                          cTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGggcgg                                                                                         >  2:275479/1‑62 (MQ=255)
                          cTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGggcgg                                                                                         <  2:1702689/62‑1 (MQ=255)
                          cTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGGACGCCCCATGGCGATCGCGCCtt                                                            >  2:1377636/1‑90 (MQ=255)
                          cTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCACGCTCGATGGCGATCGCGCAtt                                                            >  2:728637/1‑90 (MQ=255)
                               gTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCAATCGCGCATTTGGTc                                                       <  1:1446016/90‑1 (MQ=255)
                                 cAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCCCCCGCGATAGGGGTCGCGCCTTTGGttgg                                                     >  1:1155685/1‑87 (MQ=255)
                                 cAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGGACGCCCCCTTGCGGGCGCGCACTTGGTGAg                                                     >  1:1796334/1‑90 (MQ=255)
                                 cAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGGCACGCTCGATGGCGATCGCGCATTTGGTCAg                                                     >  2:763465/1‑90 (MQ=255)
                                   ggggCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCACGCCCGATGGGGGTCGCGCATTTTGTCAggg                                                   >  2:998183/1‑89 (MQ=255)
                                   ggggCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCACGCCCGATGGCGGTCGCGCATTTGGTCAggg                                                   >  1:628333/1‑89 (MQ=255)
                                   ggggCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCACCCCCGGTGGGGAGCGCGCATTTGGTTAggg                                                   >  1:951676/1‑89 (MQ=255)
                                         ccTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCACGCCCCGTGGGGAGCGCGCCATTGGTCAGGGGGCggg                                             >  1:1112796/1‑89 (MQ=255)
                                         ccTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCACGCCCGATGGCGATCGCGCATTTGGGCAGGGCGCGGc                                             >  1:495965/1‑90 (MQ=255)
                                             cAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACAt                                         <  1:314134/90‑1 (MQ=255)
                                                   aaCCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCACGCTCGATGGCGGTCCCGCCTTTTGTCAGGGCGCGGGACAAGAAGCg                                   >  1:339350/1‑90 (MQ=255)
                                                   aaCCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGGGCCGCTCGGAGGTGGTCGCGGGATTTGGCAGGGCGCGGCACATGAAACg                                   >  2:1673134/1‑90 (MQ=255)
                                                     ccATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGAt                                 <  1:236598/90‑1 (MQ=255)
                                                     ccATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGAt                                 <  1:701830/90‑1 (MQ=255)
                                                          atGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCACCCCCCATGGCGATCCCGCCTTTTGGTAGGGCGCGGGCCATgaagggaaggggg                            >  1:912862/1‑87 (MQ=255)
                                                           tGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCt                           <  1:65765/90‑1 (MQ=255)
                                                            gCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTg                          <  2:339350/90‑1 (MQ=255)
                                                            gCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTg                          <  2:495965/90‑1 (MQ=255)
                                                            gCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTg                          <  1:991008/90‑1 (MQ=255)
                                                               aaaaCCGTTAAAGAGGCGGGGGCGGCGGGCACGCTCGATGGCGATCGCGCATTTTGTCAGGGCGCGGCACATGAAGCGATGGGGCTgggg                       >  2:184661/1‑90 (MQ=255)
                                                                     gTTAAAGAGGCGGGGGGGGGGGGCACGCTCGATGGCGATCGCGCATTTTGTCAGGGCGCGGCCCATGAAGCGATGGGGCTGGGGGTTggg                 >  2:2079800/1‑88 (MQ=255)
                                                                     gTTAAAGAGGCGGGGGGGGGGGGCACCCCCGATGGCGGTCGCGCATTTTGGCAGGGGGCGGCACATGAAACGAGGGGGGTGGGGGtttgg                 >  2:1564149/1‑87 (MQ=255)
                                                                     gTTAAAGAGGCGGGGGGGGCGGGCACGCTCCATGGGGGTCGCGCATTTTGTCAGGGCGGGGCACATGAAGCGAGGGGGCTGGGGGTTggg                 >  2:1594731/1‑88 (MQ=255)
                                                                     gTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTTGTCAGGGCGCGGCACATGAAGCGATGGCGCTGGGGGTTgag                 >  2:442037/1‑90 (MQ=255)
                                                                     gTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACAt                                         <  1:9752/66‑1 (MQ=255)
                                                                     gTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACAt                                         >  2:9752/1‑66 (MQ=255)
                                                                     gTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGAt                                 <  1:632521/74‑1 (MQ=255)
                                                                     gTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGAt                                 >  2:632521/1‑74 (MQ=255)
                                                                          agagGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAgc            <  1:184661/90‑1 (MQ=255)
                                                                             ggCGGGGGGGGCGGGCACCCCCGATGGCGATCGCGCGTTTTGTCAGGGGGCGGCACAAGAAGCGATGGGGCCGGGGGGTTAGAAAACGCa         >  1:582762/1‑90 (MQ=255)
                                                                             ggCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCa         <  2:657452/90‑1 (MQ=255)
                                                                              gCGGGGGGGGCGGGCACGCTCGATGGCGATCGCGCATTTTGTCAGGGCGCGGCACATGAAGCGGTGGGGGTGGGGGTTGAGAAAGCGCAt        >  1:1050743/1‑90 (MQ=255)
                                                                               cGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATc       <  2:428210/90‑1 (MQ=255)
                                                                               cGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATc       <  1:1354316/90‑1 (MQ=255)
                                                                                  gggCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGc    <  1:327764/90‑1 (MQ=255)
                                                                                  gggCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGc    <  2:1112796/90‑1 (MQ=255)
                                                                                    gcggcgGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCAc  <  1:998183/90‑1 (MQ=255)
                                                                                     |                                                                                         
GGCATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCAC  >  NZ_CP009273/831957‑832130

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: