Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273822,9550TG55.3% 35.4 / 49.1 38T495P (ACC→CCC) ybhFATP‑binding cassette domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (4/13);  new base G (21/0);  total (25/13)
Fisher's exact test for biased strand distribution p-value = 4.40e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.92e-06
Rejected as consensus: Frequency below/above cutoff threshold.

CTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCT  >  NZ_CP009273/822871‑823038
                                                                                    |                                                                                    
cTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGgggggggg                                                                                 >  2:580172/1‑90 (MQ=255)
cTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAggg                                                                                 >  2:543319/1‑90 (MQ=255)
               gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGGCAACGCCGGGAg                                                                 >  2:299759/1‑90 (MQ=255)
               gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGTGGGGGCTACGCCGGGAg                                                                 >  1:451737/1‑90 (MQ=255)
                  tggtgACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAAAACCGGAAGGCg                                                              >  2:696750/1‑90 (MQ=255)
                      gACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTtc                                                          <  1:125850/90‑1 (MQ=255)
                         cATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGGAGGCGGgtcgtc                                                       >  1:508803/1‑90 (MQ=255)
                           tCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCga                                                     <  1:53615/90‑1 (MQ=255)
                              ccGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGGACGGCCGGAAGCGGGTTGTCgagga                                                  >  1:1042420/1‑88 (MQ=255)
                                gTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGTCGGGTCGTCGAGaaaa                                                >  2:53615/1‑89 (MQ=255)
                                   aCGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGa                                             <  2:1150490/90‑1 (MQ=255)
                                        tttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGGCGGGTCGTCGGGAAACAGAAAGTc                                        >  2:144076/1‑90 (MQ=255)
                                        tttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCCGAAGGCGGGTCGTCGAGAAAAAGaaaggc                                        >  1:990402/1‑86 (MQ=255)
                                               aCCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGGCGGGTCGTCGAGAAAAAGAATGTCCGGGTCa                                 >  2:63797/1‑90 (MQ=255)
                                                  aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGc                              <  2:269101/90‑1 (MQ=255)
                                                  aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAAAGCCGGAAGGCGGGTTGTCGGGAAAAAAAAAGTCCGGGTCATGc                              >  2:1109342/1‑90 (MQ=255)
                                                  aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGTCGGGTTGTCGGGAAAAAGAAAGTCCGGGTCATGc                              >  2:781360/1‑90 (MQ=255)
                                                  aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGGCGGGTCGTCGAGAAAAAAAAAGGCCGGGTCATGc                              >  1:1163/1‑90 (MQ=255)
                                                  aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGTCGGTTCGTCGGGAAAAAGAATGTCCGGGTCCTgg                              >  1:750442/1‑89 (MQ=255)
                                                  aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGGCGGGTCGTCGAGAAAAAGAATGTCCGGGTCATGc                              >  2:83312/1‑90 (MQ=255)
                                                  aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCCGGAGTCGGGTCGTCGAGAAAAAGAAAGTCCGGGTCATGc                              >  1:205245/1‑90 (MQ=255)
                                                    gCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGTCGGGTCGTCGAGAAACAGAATGTCCGGGTCATGCAt                            >  2:228081/1‑90 (MQ=255)
                                                     cTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATc                           <  2:897853/90‑1 (MQ=255)
                                                       gTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAg                         <  1:543319/90‑1 (MQ=255)
                                                            tgtgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc                    <  1:781360/90‑1 (MQ=255)
                                                            tgtgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc                    <  2:750442/90‑1 (MQ=255)
                                                            tgtgCAGCCAAAATTCACGGCGGGGGAGGGGGGCAAAGCCGGAAGTCGGGTCGTCGCGAAAAAGAAAGTCCGGGTTCTGCATCAGCGaaa                    >  2:822987/1‑89 (MQ=255)
                                                               gCAGCCAAAATTCACGGCGGGTGAGGGGGGCAACGCCGGAAGTCGGTTCGTCGAGAAAAAGAATGTCCGGTTCATGCATCAGCGAACAgg                 >  1:284489/1‑90 (MQ=255)
                                                               gCAGCCAAAATTCACGGCGGGTGAGGGGGGCAACGCCGGAAGTCGGGTCGTCGAGAAAAAGAAAGTCCGGGTCATGCATCAGCGAACAgg                 >  1:320062/1‑90 (MQ=255)
                                                               gCAGCCAAAATTCACGGCGGGGGGGGGGGGGGACCCCCCAAGGCGGTTCGCTGAGAAACAAAATATCCGGGTCATGGATCAGCGCACAgg                 >  1:1140847/1‑90 (MQ=255)
                                                               gCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGTAGGCGGGTCGGCGTAAAACAGAAAAGGCGGGTTATGCAGCAGCGGACaag                 >  1:897853/1‑88 (MQ=255)
                                                               gCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCCGAAGTCGGTTTGTCGAGAAACAGAATGTCCGGGTTCTTCATCAGCGCACAgg                 >  2:911656/1‑90 (MQ=255)
                                                                 aGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGgcc               <  1:918416/90‑1 (MQ=255)
                                                                 aGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGGAGTCGGGTCGTCGGGAAAAAGAAAGTGCGGGGCATGCCTCAGCGGAAAGgcc               >  1:931361/1‑90 (MQ=255)
                                                                 aGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGTCGGGTCGTCGAGAAAAAGAATGTCCGGGTTATGCATCAGCGAACAGgcc               >  2:1042041/1‑90 (MQ=255)
                                                                             cggcggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGc   <  1:940332/90‑1 (MQ=255)
                                                                             cggcggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGc   <  1:83312/90‑1 (MQ=255)
                                                                             cggcggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGc   <  1:1042987/90‑1 (MQ=255)
                                                                              ggcggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCt  <  2:725932/90‑1 (MQ=255)
                                                                                    |                                                                                    
CTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCT  >  NZ_CP009273/822871‑823038

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