Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273497,0310AC51.8% 30.6 / 38.5 27V555G (GTG→GGG) ybaLKef family K(+) transporter
Reads supporting (aligned to +/- strand):  ref base A (13/0);  new base C (0/14);  total (13/14)
Fisher's exact test for biased strand distribution p-value = 4.99e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.40e-06
Rejected as consensus: Frequency below/above cutoff threshold.

GCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTC  >  NZ_CP009273/496960‑497119
                                                                       |                                                                                        
gcgggCGCGTTTTACCATGCCTACTTTTTTTTCATATTTTACATCCGGCACCCCCCGTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCg                                                                        <  1:1782877/87‑1 (MQ=255)
gCTGGCGCTTCTTATCATCCCTACATTTTTTTCATTTTTTACACCCGGCACCCCCCTTTTCCCCCGCCCCCCCCTCCCCCGCCGGTGGCg                                                                        <  2:304923/90‑1 (MQ=255)
     cgcgTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTcc                                                                   >  1:1011956/1‑90 (MQ=255)
        ttCTTATCATCCCTAATTTTTTTTCATATTTTCCCCCCGGCACCCCCCTTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCcagc                                                                <  2:1525506/89‑1 (MQ=255)
        gTCTTATCATGCCTCCATTTTTTTCATTTTTTACACCCGGCACCCCCCTTTTCCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCcagc                                                                <  2:375585/90‑1 (MQ=255)
        gTCTTATCATGCCTACATTTTTTTCATATTTTACATCCGGCAACCACCGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCcagc                                                                <  1:674670/90‑1 (MQ=255)
        acagtATCATGCCTACTTTTTTTTCATTTTTTCCACCGGACACCCCCCTTTCCCCCCCCCCCCCCCCCCCCCGCCGGTGGCGTTTCcagc                                                                <  2:1192526/86‑1 (MQ=255)
                  gCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCa                                                      >  1:464418/1‑90 (MQ=255)
                   ccTACATTTTTTTCATTTTTTACATCCGGCACCCCCCTTTTCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCAt                                                     <  1:988004/90‑1 (MQ=255)
                   ccTACATTTTTTTAATTTTTTACACCCGGCACCCCCCTTTTCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCAt                                                     <  1:1307759/90‑1 (MQ=255)
                   ccTACATTTTTTTAATTTTTTACACCCGGCACCCCCCTTTTCCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCAt                                                     <  1:1455019/90‑1 (MQ=255)
                   ccTACATTTTTTTAATTTTTTACACCCGGCACCCCCCGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCAt                                                     <  1:1442509/90‑1 (MQ=255)
                      aCATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCactgt                                                  >  1:1192526/1‑86 (MQ=255)
                      aCATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGt                                                  >  1:1579783/1‑90 (MQ=255)
                            ttttCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGc                                            >  2:1696743/1‑90 (MQ=255)
                                  tatTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTc                                      >  2:5294/1‑90 (MQ=255)
                                          acccGCCACCCCCCTTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGc                              <  1:1684581/88‑1 (MQ=255)
                                           cccGGCACCCCCCTTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGcc                             <  1:801346/89‑1 (MQ=255)
                                              ggCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCAt                          >  2:1248649/1‑90 (MQ=255)
                                              ggCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCAt                          >  1:463174/1‑90 (MQ=255)
                                                     cccGTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTAcc                   <  1:1350276/89‑1 (MQ=255)
                                                        gTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGa                >  2:1871866/1‑90 (MQ=255)
                                                          ttcccccGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGAtt              <  1:1267663/87‑1 (MQ=255)
                                                          ttACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGAtt              >  2:1674203/1‑90 (MQ=255)
                                                               ccGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGcac         >  1:988127/1‑90 (MQ=255)
                                                                cGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGcacc        >  2:1499153/1‑90 (MQ=255)
                                                                cGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGcacc        >  2:1750669/1‑90 (MQ=255)
                                                                      caccTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTc  <  1:1499153/90‑1 (MQ=255)
                                                                       accTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTc  >  2:1242135/1‑89 (MQ=255)
                                                                       |                                                                                        
GCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTC  >  NZ_CP009273/496960‑497119

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