Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,669,4650CA55.0% 10.8 / 27.2 20L240F (TTG→TTTyphDABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (7/2);  new base A (0/11);  total (7/13)
Fisher's exact test for biased strand distribution p-value = 4.64e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.51e-03
Rejected as consensus: Frequency below/above cutoff threshold.

GCCGTACCGCCGACCACGACCGCGGCGATGACGTCAAACTCCAGACCGTTTGCGGCACCTGCGTTACCAGAACCGAGGCGCGCCGCCAACAAAATGCCGGTCACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGAT  >  NZ_CP009273/2669376‑2669536
                                                                                         |                                                                       
gCCGTCCCGCCGCCCACGCCCGCGGCGATGACGTCAAACTCCAGACCGTTTGCGGCACCTGCGTTACCAGAACCGAGGCGCGCCGCcaac                                                                         <  1:526739/90‑1 (MQ=255)
                cGACCGCGGCGATGACGTCAAACTCCAGACCGTTTGCGGCACCTGCGTTACCAGAACCGAGGCGCGCCGCCAACAAAATGCCGGTCAccg                                                         >  1:395392/1‑90 (MQ=255)
                   ccGCGGCGATGACGTCAAACTCCAGACCGTTTGCGGCACCTGCGTTACCAGAACCGAGGCGCGCCGCCAACAAAATGCCGGTCAccgccg                                                      >  1:307287/1‑90 (MQ=255)
                     gcggcgATGACGTCAAACTCCAGACCGTTTGCGGCACCTGCGTTACCAGAACCGAGGCGCGCCGCCAACAAAATGCCGGTCACCGCCGCt                                                    >  1:816321/1‑90 (MQ=255)
                                       cccAGCCCTTTTGCGGCCCCTGCGTTCCCAAACCCGAGGCGCGCCGCCAAAAAAATGCCGGTCACCGCCGCTAATAATCCCGAAAGGGTa                                  <  2:307287/89‑1 (MQ=255)
                                       cccAGCCCTTTTGCGCCCCCTGCTTTACCAGAACCGGGGCGCGCCGCCAAAAAAATGCCGGTCACCGCCGCTAATAATCCCGAAAGGGTa                                  <  1:194532/89‑1 (MQ=255)
                                       cccAGCCCTTTGGCGGCCCCTGCGTTACCAGACCCGGGGCGCGCCGCAAAAAAAATGCCGGTCACCGCCGCTAATAATCCCGAAAGGGTa                                  <  2:472140/89‑1 (MQ=255)
                                       cccAGCCCGTTTGCGGCACCTGCGTTCCCAGAACCGAGGCGCGCCGCAAACAAAATGCCGGTCACCGCCGCTAATAATCCCGAAAGGGTa                                  <  2:28546/89‑1 (MQ=255)
                                       cccAGCCCGTTTGCGGCACCTGCGTTCCCAAAACCGAGGCGCGCCCCCAAAAAAATGCCGGTCACCGCCGCTAATAATCCCGAAAGGGTa                                  <  2:311806/89‑1 (MQ=255)
                                       cccAGCCCGTTTGCGCCCCCTGCTTTCCCAAACCCGGGGCGCGCCGCCAACAAAATGCCGGTCACCGCCGCTAATAATCCCGAAAGGGTa                                  <  2:529901/89‑1 (MQ=255)
                                       cccAGCCCGTTGGCGCCCCCGGCTTTCCCAGACCCGGGGCGCGCCGCCAAAAAAATGCCGGTCACCGCCGCTAATAATCCCGAAAGGGTa                                  <  2:887602/89‑1 (MQ=255)
                                          agcccGTTGGCGGCCCCTGCTTTCCAAGAACCGAGGCGCGCCGCCAAAAAAATGCCGGTCACCGCCGCTAATAATCCCGAAAGGGTAAAg                               <  2:658372/87‑1 (MQ=255)
                                           gACCGTTTGCGGCACCTGCGTTACCAGAACCGAGGCGCGCCGCCAACAAAATGCCGGTCACCGCCGCTAATAATCCCGAAAGGGTAAAGa                              >  2:156025/1‑90 (MQ=255)
                                              cGTTTGCGGCCCCTGCTTTCCCAAAACCGGGGCGCGCCGCAAAAAAAAGGCCGGTCCCCGCCGCTAATAATCCCGAAAGGGTAAAGATaa                           <  1:859486/90‑1 (MQ=255)
                                                tttGCGGCACCTGCGTTACCAGAACCGAGGCGCGCCGCCAACAAAATGCCGGTCACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGa                         >  2:569740/1‑90 (MQ=255)
                                                    gggCCCCTGCTTTCCCAAAACGGGGGCGCCCCCCAAAAAAAATGCCGGTCCCCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATgc                     <  1:320067/89‑1 (MQ=255)
                                                         ccTGCGTTACCAGAACCGAGGCGCGCCGCCAACAAAATGCCGGTCACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCAcc                >  1:168589/1‑90 (MQ=255)
                                                                    aGAACCGAGGCGCGCCGCCAACAAAATGCCGGTCACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGtt     >  2:213432/1‑90 (MQ=255)
                                                                     aacccGAGGCGCGCCGCCAAAAAAATGCCGGTCACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTg    <  2:248852/87‑1 (MQ=255)
                                                                       aCCGGGGCGCGCCGCAAAAAAAACCCCGGCCCCCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGAt  <  1:152789/90‑1 (MQ=255)
                                                                       aCCGAGGCGCCCCCCAAAAAAAATGCCGGTCCCCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGAt  <  2:788170/90‑1 (MQ=255)
                                                                                         |                                                                       
GCCGTACCGCCGACCACGACCGCGGCGATGACGTCAAACTCCAGACCGTTTGCGGCACCTGCGTTACCAGAACCGAGGCGCGCCGCCAACAAAATGCCGGTCACCGCCGCTAATAATCCCGAAAGGGTAAAGATAAGAATGCGCACCCGACGAACGTTGAT  >  NZ_CP009273/2669376‑2669536

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