Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273497,0200AC46.3% 88.7 / 33.6 67*559E (TAA→GAA) ybaLKef family K(+) transporter
Reads supporting (aligned to +/- strand):  ref base A (32/4);  new base C (0/31);  total (32/35)
Fisher's exact test for biased strand distribution p-value = 5.85e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGA  >  NZ_CP009273/496931‑497105
                                                                                         |                                                                                      
ggCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTa                                                                                        >  1:911151/1‑90 (MQ=255)
   aTTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTAccc                                                                                     <  1:1636898/90‑1 (MQ=255)
    tttGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTAcccc                                                                                    >  1:47830/1‑90 (MQ=255)
       ggCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTc                                                                                 <  1:1664726/90‑1 (MQ=255)
             cGATGCCTGATGCGCCGCTGGCGCGTCTTATCATGCCTACATTTTTTTCATATTTTACATCCGGCACCCCCCGTTTCCCCCGTCaccccc                                                                          <  1:1125852/90‑1 (MQ=25)
                 gCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCAc                                                                      >  1:525491/1‑90 (MQ=255)
                 gCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCAc                                                                      >  2:2860992/1‑90 (MQ=255)
                 gCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCAc                                                                      >  2:2769197/1‑90 (MQ=255)
                     gATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTAccc                                                                                     <  1:2817389/72‑1 (MQ=255)
                     gATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTAccc                                                                                     >  2:2817389/1‑72 (MQ=255)
                          gACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTg                                                             >  1:1518830/1‑90 (MQ=255)
                           aCGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTgg                                                            >  2:2767590/1‑90 (MQ=255)
                             gCTGGCGCGTCTTATCATGCCTACATTTTTTTCATATTTTACACCCGGCAACCACCGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCg                                                          <  2:1443327/90‑1 (MQ=38)
                                ggCGCGTCTTATCATGCCTACATTTTTTTCATATTTTACATCCGGCAACCCCCGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGttt                                                       <  2:904309/90‑1 (MQ=255)
                                ggCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGttt                                                       >  2:2595432/1‑90 (MQ=255)
                                ggCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGttt                                                       >  1:83790/1‑90 (MQ=255)
                                    cgTCTTATCATGCCTACATATTTTTCATATTTTACACCCGGCACCCCCCGTTTCCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCcag                                                   <  2:905061/90‑1 (MQ=255)
                                     gTCTTATCATGCCTACATTTTTTTCATATTTTACATCCGGCACCCCCCGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCcagc                                                  <  2:757663/90‑1 (MQ=255)
                                     gTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCcagc                                                  >  2:780066/1‑90 (MQ=255)
                                     gTCTTATCATGCCTAAATTTTTTTCATTTTTTACATCCGGCAACCCCCGTTTCCCCCGTCCCCACCTCACCCGCCGGTGGCGTTTCcagc                                                  <  2:1261114/90‑1 (MQ=255)
                                     gTCTTATCATGCCTAAATTTTTTTCATATTTTACACCCGGAACCCCCCGTTTCCCCCGCCACCCCCTCCCCCGCCGGTGGCGTTTCcagc                                                  <  2:1434840/90‑1 (MQ=255)
                                      tCTTACCATCCCTACATTTTTTTCATTTTTTCCACCCGGCACCCCCCGTTCCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCcagca                                                 <  1:1440888/90‑1 (MQ=255)
                                        ttATCATCCCTACTTTTTTTTCATTTTTTACACCCGGCACCCCCCTTTTCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGt                                               <  1:1195285/90‑1 (MQ=255)
                                           tCATGCCTACATTTTTTTAATTTTTTACATCCGGCAACCACCTTTTCCCCCGTCACCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTcc                                            <  2:525491/90‑1 (MQ=255)
                                           tCATCCCTAAATTTTTTTAATTTTTTACACCCGGCACCCCCCTTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTcc                                            <  2:1868870/90‑1 (MQ=255)
                                           tCACGCCTACATTTTTTTCATTTTTTACACCCGGCACCCCCCTTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTcc                                            <  2:3020474/90‑1 (MQ=255)
                                              tcccTATATTTTTTTCATTTTTTACCTCCGACAACCCCCTTTTCCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGc                                         <  1:1517407/88‑1 (MQ=255)
                                               gCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTca                                                                                >  1:2547533/1‑51 (MQ=255)
                                               gCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTca                                                                                <  2:2547533/51‑1 (MQ=255)
                                               gCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCa                                        >  2:2786774/1‑90 (MQ=255)
                                               gCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCa                                        >  1:2251057/1‑90 (MQ=255)
                                               gCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCa                                        >  2:1150623/1‑90 (MQ=255)
                                               gCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCa                                        >  2:1156060/1‑90 (MQ=255)
                                               gCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCa                                        >  1:1870452/1‑90 (MQ=255)
                                                ccaaCATTTTTTTTATTTTTCCCCCCCGCAACCCCCCTTTCCCCCCCCCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCAt                                       <  1:1285980/87‑1 (MQ=255)
                                                ccTACATTTTTTTCATTTTTTCCACCCGGCACCCCCCTTTTCCCCCCCCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCAt                                       <  2:2118307/90‑1 (MQ=255)
                                                ccTACATTTTTTTCATTTTTTCCACCCGGCACCCCCCGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCAt                                       <  2:2251057/90‑1 (MQ=255)
                                                ccTACATTTTTTTCATTTTTTACCCCCGCCACCCCCCGTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCAt                                       <  1:110093/90‑1 (MQ=255)
                                                ccTACATTTTTTTCATATTTTACACCCGGCAACCCCCGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCAt                                       <  2:346170/90‑1 (MQ=255)
                                                ccTACATTTTTTTATTTTTTTACACCCGCCACCCCCCTTTCCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCAt                                       <  1:972522/90‑1 (MQ=255)
                                                ccTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCAt                                       >  2:1285980/1‑90 (MQ=255)
                                                  tACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATgg                                     >  2:1773673/1‑90 (MQ=255)
                                                   aCATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGt                                    >  2:1073531/1‑90 (MQ=255)
                                                   aCATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGt                                    >  2:2253265/1‑90 (MQ=255)
                                                     atttttttCATTTTTTACACCCGGCACCCCCCTTTTCCCCCCCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTAc                                  <  1:780066/87‑1 (MQ=255)
                                                       tttttccatcttttccatcCGCCCCCCCCCTTTTCCCCCCTCCCCACCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACgg                                <  2:1659209/89‑1 (MQ=255)
                                                        tttttCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACggg                               >  2:1736688/1‑90 (MQ=255)
                                                        tttttCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACggg                               >  2:367193/1‑90 (MQ=255)
                                                        tttttCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACggg                               >  1:2339073/1‑90 (MQ=255)
                                                         ttttCATTTTTTACACCCGGCACCCCCCTTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGc                              <  1:1073531/90‑1 (MQ=255)
                                                         ttttCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGc                              >  2:351763/1‑90 (MQ=255)
                                                         ttttCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGc                              >  2:2944765/1‑90 (MQ=255)
                                                             cATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGAtt                          >  2:421337/1‑90 (MQ=255)
                                                             cATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGAtt                          >  1:2607428/1‑90 (MQ=255)
                                                                    tccATCCGGCACCCACCTTTTCCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCt                   <  2:492241/88‑1 (MQ=255)
                                                                    tACATCCGGCACCCCCCGTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCt                   <  1:2247774/90‑1 (MQ=255)
                                                                     aCATCCGGCAACCCCCGTTTCCCCCGTCCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTc                  <  2:252063/90‑1 (MQ=255)
                                                                       cgccGCCAACCGCTGTTTCCCCCTCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGc                <  1:2412653/88‑1 (MQ=255)
                                                                        cccGGCACCCCCCTTTTCCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGcc               <  2:2882586/89‑1 (MQ=255)
                                                                        cccGGCACCCCCCGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGcc               <  2:927768/89‑1 (MQ=255)
                                                                         ccGGCACCCCCCTTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGccc              <  1:367193/90‑1 (MQ=255)
                                                                         ccGGCACCCCCCGTTTCCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGccc              <  2:424924/90‑1 (MQ=255)
                                                                          cGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCa             >  1:1397676/1‑90 (MQ=255)
                                                                           ggCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCAt            >  2:1145692/1‑90 (MQ=255)
                                                                           ggCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCAt            >  1:1760150/1‑90 (MQ=255)
                                                                                  accGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTAcc     >  1:2883696/1‑90 (MQ=255)
                                                                                   ccttttCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCt    <  2:1303868/87‑1 (MQ=255)
                                                                                     gTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGa  >  2:2344802/1‑90 (MQ=255)
                                                                                         |                                                                                      
GGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGA  >  NZ_CP009273/496931‑497105

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