Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,401,1950TG56.2% 9.3 / 16.4 16V77G (GTA→GGA) panFsodium/pantothenate symporter
Reads supporting (aligned to +/- strand):  ref base T (1/6);  new base G (9/0);  total (10/6)
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.80e-04
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GGCGATGACGCTCACCGCGACCTATATCAGTGCCAGTTCGTTTATCGGCGGGCCAGGAGCTGCTTATAAATACGGGCTGGGCTGGGTATTGCTGGCGATGATTCAGCTTCCTGCAGTCTGGCTTTCACTCGGTATTCTCGGCAAGAAGTTTGCGATTCTTGCGCGCCGCTAC  >  NZ_CP009273/3401109‑3401280
                                                                                      |                                                                                     
ggCGATGACGCTCACCGCGACCTATATCAGTGCCAGTTCGTTTATCGGCGGGCCAGGGGCTGCTTATAAATACGGGCTGGGCTGGGGAtt                                                                                    >  2:13470/1‑90 (MQ=255)
          cTCACCGCGACCTATATCAGTGCCAGTTCGTTTATCGGCGGGCCAGGAGCTGCTTATAAATACGGGCTGGGCTGGGGATTGCTGGCgatg                                                                          >  2:245106/1‑90 (MQ=255)
                  gACCTATATCAGTGCCAGTTCGTTTATCGGCGGGCCAGGAGCTGCTTATAAATACGGGCTGGGCTGGGTATTGCTGGCGATGATTCAGCt                                                                  <  2:631120/90‑1 (MQ=255)
                        tatCAGTGCCAGTTCGTTTATCGGCGGGCCAGGAGCTGCTTATAAATACGGGCTGGGCTGGGTATTGCTGGCGATGATTCAGCTTCCTGc                                                            <  1:588780/90‑1 (MQ=255)
                             gTGCCAGTTCGTTTATCGGCGGGCCCGGGGCTGCTTATAAAAACGGGCTGGGGTGGGGATTTGTGGGGGTGATTTAGCTTCCTTCCGTCt                                                       >  2:588780/1‑90 (MQ=255)
                             gTGCCAGTTCGTTTATCGGCGGGCCAGGGGGTGCTTTTAAATACGGGCTGGGGTGGGGATTTCTGGGGGTGATTTAGCTTCCTGCAGTCt                                                       >  2:601884/1‑90 (MQ=255)
                               gCCAGTTCGTTTATCGGCGGGCCAGGAGCTGCTTTTAAATACGGGCTGGGGTGGGGATTTCTGGGGGTGATTTAGCTTCCTTCAGTCTgg                                                     >  2:395463/1‑90 (MQ=255)
                               gCCAGTTCGTTTATCGGCGGGCCAGGAGCTGCTTATAAATACGGGCTGGGCTGGGGATTGCTGGCGATGATTCAGCTTCCTGCAGTCTgg                                                     >  2:106814/1‑90 (MQ=255)
                                 cAGTTCGTTTATCGGCGGGCCAGGAGCTGCTTATAAATACGGGCTGGGCTGGGGATTGCTGGCGATGATTCAGCTTCCTGCAGTCTGGCt                                                   >  2:38705/1‑90 (MQ=255)
                                   gTTCGTTTATCGGCGGGCCAGGAGCTGCTTATAAATACGGGCTGGGCTGGGTATTGCTGGCGATGATTCAGCTTCCTGCAGTCTGGCttt                                                 <  2:264977/90‑1 (MQ=255)
                                       gTTTATCGGCGGGCCAGGAGCTGCTTATAAATACGGGCTGGGGTGGGGATTGCTGGCGATGATTCAGCTTCCTGCAGTCTGGCTTTCACt                                             >  2:612150/1‑90 (MQ=255)
                                       gTTTATCGGCGGGCCAGGAGCTGCTTATAAATACGGGCTGGGCTGGGGATTGCTGGCGATGATTCAGCTTCCTGCAGTCTGGCTTTCACt                                             >  2:189965/1‑90 (MQ=255)
                                                cggGCCAGGAGCTGCTTATAAATACCTGCTTGGTTGGGTATTGCTGGCGATCATTCACCTTCCTCCAGTCTGGCTTTCCCTCGGTATtct                                    <  1:23441/90‑1 (MQ=255)
                                                                              gggctgggTATTGCTGGCGATGATTCAGCTTCCTGCAGTCTGGCTTTCACTCGGTATTCTCGGCAAGAAGTTTGCGATTCTTGCGcgccg      >  2:282881/1‑90 (MQ=255)
                                                                                  tgggTATTGCTGGCGATGATTCAGCTTCCTGCAGTCTGGCTTTCACTCGGTATTCTCGGCAAGAAGTTTGCGATTCTTGCGCGCCGCTAc  <  1:13470/90‑1 (MQ=255)
                                                                                  tgggTATTGCTGGCGATGATTCAGCTTCCTGCAGTCTGGCTTTCACTCGGTATTCTCGGCAAGAAGTTTGCGATTCTTGCGCGCCGCTAc  <  1:565938/90‑1 (MQ=255)
                                                                                      |                                                                                     
GGCGATGACGCTCACCGCGACCTATATCAGTGCCAGTTCGTTTATCGGCGGGCCAGGAGCTGCTTATAAATACGGGCTGGGCTGGGTATTGCTGGCGATGATTCAGCTTCCTGCAGTCTGGCTTTCACTCGGTATTCTCGGCAAGAAGTTTGCGATTCTTGCGCGCCGCTAC  >  NZ_CP009273/3401109‑3401280

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