Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,178,2500AC56.2% 4.4 / 12.8 16G12G (GGT→GGGpotDspermidine/putrescine ABC transporter substrate‑binding protein PotD
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 = 5.70e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAGCTGGCTGGCA  >  NZ_CP009273/1178163‑1178339
                                                                                       |                                                                                         
tCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCAcc                                                                                         <  1:112077/90‑1 (MQ=255)
         cACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAg                                                                                >  2:427213/1‑90 (MQ=255)
          ccGTCCTCGGCCCAGTTGTAAAAATACGGGGTGTTGTTGTCATCGGCGGGGGCGGCGCTCATGCCCAGTGCAAGAGCCCCCGCCGCGAGc                                                                               <  2:478857/89‑1 (MQ=255)
                  ggTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCg                                                                       >  1:149598/1‑90 (MQ=255)
                   gTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGt                                                                      >  2:215889/1‑90 (MQ=255)
                                       cGTGTTGTTGTCAGCGGGGTGGGGGGCGCTCATCCCAAGTACCAGAGCCCCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACg                                                  <  2:530524/90‑1 (MQ=255)
                                           ttgttgtCAGCGGGGTGAGCGGCGCTCATGCCCAGTGCCAGAGCCCCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTccc                                              <  2:197803/90‑1 (MQ=255)
                                            ttttgtCATCGGCGGGGGGGGCGCTCATGCCCAGTGCCAGACCCCCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTcccc                                             <  2:75385/88‑1 (MQ=255)
                                            tgttgtCATCGGGGTGGGCGGCCCTCCTGCCCAGTGCAAGAGCCCCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTcccc                                             <  2:413096/90‑1 (MQ=255)
                                                  cATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGtt                                       >  1:377191/1‑90 (MQ=255)
                                                  cATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGGGGCGTGGCCCTTTTTTCATTTAACGTCCCCTGGGtt                                       >  2:101047/1‑90 (MQ=255)
                                                        cGTGGGGGGCGCTCATGCCCAGTGCCAGAGCCCCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCttt                                 <  1:513105/90‑1 (MQ=255)
                                                             gcggcgCTCTTCCCCAGTCCCAGCCCCCCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGttt                            <  2:86653/90‑1 (MQ=255)
                                                                     cATGCCCAGTCCCAGACCCCCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGa                    <  2:149598/90‑1 (MQ=255)
                                                                                   gagCCCCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAgctggct      <  1:336704/90‑1 (MQ=255)
                                                                                   gagCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAgctggct      >  1:310528/1‑90 (MQ=255)
                                                                                       aCCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGa                    >  2:137563/1‑72 (MQ=255)
                                                                                       aCCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGa                    <  1:137563/72‑1 (MQ=255)
                                                                                       aCCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAGCTGGCTGGCa  >  1:386745/1‑90 (MQ=255)
                                                                                       aCCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAGCTGGCTGGCa  >  1:375535/1‑90 (MQ=255)
                                                                                       |                                                                                         
TCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAGCTGGCTGGCA  >  NZ_CP009273/1178163‑1178339

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