Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP00927369,9000AC56.2% 10.4 / 17.4 16V370G (GTG→GGG) thiPthiamine/thiamine pyrophosphate ABC transporter permease ThiP
Reads supporting (aligned to +/- strand):  ref base A (6/1);  new base C (0/9);  total (6/10)
Fisher's exact test for biased strand distribution p-value = 8.74e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.07e-03
Rejected as consensus: Frequency below/above cutoff threshold.

GATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAA  >  NZ_CP009273/69815‑69979
                                                                                     |                                                                               
gATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTg                                                                             >  2:3963/1‑90 (MQ=255)
      gttgttGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGc                                                                       >  2:387586/1‑90 (MQ=255)
                aGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTc                                                             >  2:591669/1‑90 (MQ=255)
                            ccGGTACCCAGCCCACCCCCGGGCAGGGGGAGGATCAACATGCCCCTCATCCCCAGCCCCTGACCCGCCAGCATTTTCTGCCGCGCCCGc                                                 <  1:497054/90‑1 (MQ=255)
                                         caATCCCCGGCATGGCGGGGACAAACATGCCCCTCATCCCCAGCCCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGc                                    <  1:581405/90‑1 (MQ=255)
                                           aTCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTc                                  >  2:205933/1‑90 (MQ=255)
                                                 ggCAGGGGGAGGATAAACATGCCGCCCATCCCCAGCCCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATagc                            <  2:519978/90‑1 (MQ=255)
                                                     tGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCCCCCGCAGTTCGCGACTGCTCCATAgcagca                        >  1:222329/1‑90 (MQ=255)
                                                        ggAGGATCAACATGCCGCTCATCCCCAGCCCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGcatca                     <  1:209100/89‑1 (MQ=255)
                                                             aTCAACATGCCGCTCATCCCCAGCCCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTc                <  2:189400/90‑1 (MQ=255)
                                                               cAACATCCCCCTCATCCCCACCCCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAg              <  1:3963/90‑1 (MQ=255)
                                                                    tcccGCTCACCCCCAGCCCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTa         <  1:222889/88‑1 (MQ=255)
                                                                    tGCCGCTCACCTCCAGCCCCTGCCCTGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTa         <  1:222888/90‑1 (MQ=255)
                                                                      ccGCTCATCCCCAGCCCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACg       <  1:247676/90‑1 (MQ=255)
                                                                       cGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGc      <  2:167472/90‑1 (MQ=255)
                                                                           cATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGGCAGCACTACGCACaa  >  2:310293/1‑90 (MQ=255)
                                                                                     |                                                                               
GATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAA  >  NZ_CP009273/69815‑69979

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: