Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273497,0360AC55.6% 5.6 / 16.4 18G553G (GGT→GGGybaLKef family K(+) transporter
Reads supporting (aligned to +/- strand):  ref base A (4/4);  new base C (0/10);  total (4/14)
Fisher's exact test for biased strand distribution p-value = 2.29e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.08e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAA  >  NZ_CP009273/496919‑497163
                                                                                                                     |                                                                                                                               
cACGCCGCATCGGGCATTGGGCTCTCGTTGCTTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATTTTTTACACCCGGCAACCACCTTTTCCCCCGTCCCCACCCCCCCCGCCGGTGGCGTTTCcag                                                                                                             <  2:287842/138‑1 (MQ=37)
        aTCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACTCCCGGCACCCCCCTTTTCCCCCGCCCCCACCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAg                                                                                                    <  2:381046/139‑1 (MQ=37)
             gCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCCCCGTTTCCCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCa                                                                                                  <  1:254164/136‑1 (MQ=255)
                            tGCCTGAGGCGCCGCTGGCGCGTCTTATCATTCCTACATATTTTTCATATTTTACATCCGGCACCCACCGTTTCCCCCGTCACCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCAc                                                                                <  1:316568/139‑1 (MQ=255)
                            tGCCTGAGGCGCCCCTGGCGCGTCTTACCATCCCTACATATTTTTCATATTTTACACCCGGCACCCCCCTTTTCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCAc                                                                                <  1:254081/139‑1 (MQ=255)
                               cTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATTTTTTACACCCGGCACCCACCGTTTCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCt                                                                             <  1:195361/139‑1 (MQ=38)
                                  aTGCGCCGCTGGCGCGTCTTATCATCCCTACATATTTTTCATATTTTCCATCCGGCAACCCCCTTTTCCCCCGCCCCCACCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCg                                                                           <  1:34991/138‑1 (MQ=255)
                                    gCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGcc                                                                         >  2:428049/1‑138 (MQ=255)
                                      gACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCAt                                                                      >  2:347906/1‑139 (MQ=255)
                                                   cTTATCATGCCTACATATTTTTCATATTTTACATCCGCCAACCCCCGTTTCCCCCGTCACCACCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGAtt                                                          <  2:369781/138‑1 (MQ=255)
                                                         aTGCCTACATATTTTTCATATTTTACATCCGCCAACCACCGTTTCCCCCGTCACCACCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGcacc                                                    <  1:92956/138‑1 (MQ=255)
                                                            ccTACATTTTTTTCATATTTTACATCCGGCACCCCCCGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACg                                                 <  2:206622/138‑1 (MQ=255)
                                                                cATATTTTTCATATTTTACACCCGGCACCCACCGTTTCCCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCg                                             <  2:436694/138‑1 (MQ=255)
                                                                                      cGGCAACCCCCTTTTCCCCCGCCACCACCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGt                       <  1:178615/138‑1 (MQ=255)
                                                                                       ggCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCa                     >  1:354948/1‑139 (MQ=255)
                                                                                                cGTTTCCCCCGTCACCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGgcg            <  2:210662/139‑1 (MQ=255)
                                                                                                   ttcccccGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGgcgc           <  2:354948/134‑1 (MQ=255)
                                                                                                          gTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaataa  >  1:91971/1‑139 (MQ=255)
                                                                                                                     |                                                                                                                               
CACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAA  >  NZ_CP009273/496919‑497163

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