Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273892,8110TG52.9% 10.5 / 17.4 17V91G (GTG→GGG) potIputrescine ABC transporter permease PotI
Reads supporting (aligned to +/- strand):  ref base T (1/7);  new base G (9/0);  total (10/7)
Fisher's exact test for biased strand distribution p-value = 4.11e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-03
Rejected as consensus: Frequency below/above cutoff threshold.

AGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCACCGCGCCGCTGGTAATGCCAGATGTCATCAC  >  NZ_CP009273/892726‑892889
                                                                                     |                                                                              
aGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTgcgc                                                                            <  1:375815/90‑1 (MQ=255)
           tttAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGGGGGCTTTTGGAgggt                                                                 >  1:231745/1‑88 (MQ=255)
                gTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTcgcg                                                            <  1:309893/90‑1 (MQ=255)
                gTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGGGCTGGGGCGCTTTTGCAGGGTTcgcg                                                            >  1:170741/1‑90 (MQ=255)
                  cTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGGGGGGGTGGGGGGCGTTGGGAGGGTTCGCggg                                                          >  2:332655/1‑89 (MQ=255)
                    gACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGGGCTGGGGGGCTTTGGCAGGGTTTGCGGATc                                                        >  1:52832/1‑90 (MQ=255)
                    gACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGGGGGGGTGGGGGGGTTTGGCGGGGTTTGCGGGTc                                                        >  2:309893/1‑90 (MQ=255)
                              gCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAATGGCttt                                              <  2:86092/90‑1 (MQ=255)
                              gCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAATGGCttt                                              <  1:539550/90‑1 (MQ=255)
                              gCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGGGCTGGGGCGCTTTTGCAGGGTTCGCGGGTCAAAAAGCttt                                              >  1:687625/1‑90 (MQ=255)
                                  gtgCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGGGCTGGGGCGCTTTTGGAGGGTTTTCGGGTCAAAAGGGttttgcc                                          >  1:451492/1‑86 (MQ=255)
                                                 cggcgATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCACCgcgc                           <  1:332655/90‑1 (MQ=255)
                                                    cgATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCACCGcgccgc                        <  2:30917/90‑1 (MQ=255)
                                                                  aTTGCGGCGGTGGTGCTGGGGCGCTTTGGCAGGTTTCGCGGGTCAAATTGCTTTGCCTTTTTGATCACCGCGCCGCTGGGAATGCCAGAt          >  1:567807/1‑90 (MQ=255)
                                                                     gcggcgGTGGTGCTGGTGCGTTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCGCCGCGCCGCTGGTAATGCCAGATGtc       <  1:298225/90‑1 (MQ=255)
                                                                       ggcggTGGGGGTGGTGCGCTTTTGCAGGTTTCGCGGGTCAAAAGGGTTTTCCTTTATGATCACCCCGCCGCTGGGAATTCCAGGTGgcat     >  2:55062/1‑90 (MQ=255)
                                                                          ggtggtGCTGGGGCGCTTTTGCAGGGTTCCCGGATCAAAAGGCTTTGCCTTTTTGATCACCCCGCCCCTGGTAATGCCAGATGTCATCAc  >  2:658828/1‑90 (MQ=255)
                                                                                     |                                                                              
AGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCACCGCGCCGCTGGTAATGCCAGATGTCATCAC  >  NZ_CP009273/892726‑892889

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