Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273303,5320GT53.8% 6.5 / 16.5 13Y504* (TAC→TAAecpCfimbrial usher EcpC
Reads supporting (aligned to +/- strand):  ref base G (1/5);  new base T (7/0);  total (8/5)
Fisher's exact test for biased strand distribution p-value = 4.66e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.56e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TATTGGCGTTGCTGTCGCCGTTGTTATAGCGCTGAATACCGGCCCGCAGGCCAAGCGAACCAAAGGTACCGCTGTAGACATTTTGATAGTAATCTGCCGTGTAATAATGGCTGTTGTAACGGCGGTCATCATTGTAGCTGATGCTGAATGTGCCCAGCTTCGACCACAGTGAGTTCAGGTTGAGTGTGCCGCCGATTGCACGGTTGTCGGCATCGCTACGTCGCAATTGATTGCCAATGCGGGTTTTTTCCTGATTAACCCA  >  NZ_CP009273/303399‑303660
                                                                                                                                     |                                                                                                                                
tATTGGCGTTGCTGTCGCCGTTGTTATAGCGCTGAATACCGGCCCGCAGGCCAAGCGAACCAAAGGTACCGCTGTAGACATTTTGATAGTAATCTGCCGTGTAATAATGGCTGTTGTAACGGCGGTCATCATTGTAGCt                                                                                                                             <  1:254310/139‑1 (MQ=255)
          gCTGTCGCCGTTGTTATAGCGCTGAATACCGGCCCGCAGGCCAAGCGAACCAAAGGTACCGCTGTAGACATTTTGATAGTAATCTGCCGTGTAATAATGGCTGTTGTAACGGCGGGCATCATTTTAACTGGTGCTGa                                                                                                                     >  2:285447/1‑137 (MQ=255)
                          taGCGCTGAATACCGGCCCGCAGGCCAAGCGAACCAAAGGTACCGCTGTAGACATTTTGATAGTAATCTGCCGTGTAATAATGGCTGTTGTAACGGCGGTCATCATTGTAGCTGATGCTGAATGTGCCCAGCTTCGAcc                                                                                                   <  1:368257/139‑1 (MQ=255)
                               cTGAATACCGGCCCGCAGGCCAAGCGAACCAAAGGTACCGCTGTAGACATTTTGATAGTAATCTGCCGTGTAATAATGGCTGTTGTAACGGCGGTTATCATTTTAGCTTGTGGTGAAATTGCCCCGGTTCcaccaccag                                                                                              >  2:140934/1‑136 (MQ=255)
                                  aaTACCGGCCCGCAGGCCAAGCGAACCAAAGGTACCGCTGTAGACATTTTGATAGTAATCTGCCGTGTAATAATGGCTGTTGTAACGGCGGTCATCATTTTAGCTGATGCTGAATGTGCCCAGCTTCGGCCCCAgaggg                                                                                           >  1:85096/1‑135 (MQ=255)
                                                      gCGAACCAAAGGTACCGCTGTAGACATTTTGATAGTAATCTGCCGTGTAATAATGGCTGTTGTAACGGCGGTCATCATTGTAGCTGATGCTGAATGTGCCCAGCTTCGACCACAGTGAGTTCAGGTTGAGTGTgccgcc                                                                       <  1:168286/139‑1 (MQ=255)
                                                                        tGTAGACATTTTGATAGTAATCTGCCGTGTAATAATGGCTGTTGTAACGGCGGTCATCATTTTAGCTGATGCTGGATTTGCCCAGCTTCGCCCCCCATGGGTTTAGGTTGTGTGTGCCCCCGATTTCCCGGGTGTCGGc                                                     >  1:59672/1‑139 (MQ=255)
                                                                         gTAGACATTTTGATAGTAATCTGCCGTGTAATAATGGCTGTTGTAACGGCGGTCATCATTTTAGCTGGTGCTGAATGTGGCCCGCTTCGGCCACCGTGAGGTCAGGGTGAGGGGGCCGCCGGTTGCACGGGGGTCGGCa                                                    >  1:36820/1‑139 (MQ=255)
                                                                                 tttGATAGTAATCTGCCGTGTAATAATGGCTGTTGTAACGGCGGTCCACAATTTGGCTGATGCTGAATGTGGCCAGCGTCGAACACAGGGAGGTGTGGGTGTGTGTGGCGCCGGTTGCGCGGGGGTGGGCGTCGCTACg                                            >  1:65312/1‑139 (MQ=255)
                                                                                                             gCTGTTGTAACGGCGGTCATCATTGTAGCTGATGCTGAATGTGCCCAGCTTCGACCACAGTGAGTTCAGGTTGAGTGTGCCGCCGATTGCACGGGTGTCGGCATCGCTACGTCGCAATTGATTGCCAATGCGGGttttt                <  1:402156/139‑1 (MQ=255)
                                                                                                                           ggTCATCATTTTAGATGATGCTGAATGGGCCCAGGTTCGACCACAATGGGTTTAGGGTGGGGGGGCCGCCGGGTGTACGGGGGTCGGCCTCGCTACGGCGCAATTGATTGCCAATGCGGGGTTTTTCCTGGTTAACCca  >  2:145541/1‑139 (MQ=255)
                                                                                                                           ggTCATCATTGTAGCTGATGCTGAATGTGCCCAGCTTCGACCACAGTGAGtt                                                                                         >  1:432705/1‑52 (MQ=255)
                                                                                                                           ggTCATCATTGTAGCTGATGCTGAATGTGCCCAGCTTCGACCACAGTGAGtt                                                                                         <  2:432705/52‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                
TATTGGCGTTGCTGTCGCCGTTGTTATAGCGCTGAATACCGGCCCGCAGGCCAAGCGAACCAAAGGTACCGCTGTAGACATTTTGATAGTAATCTGCCGTGTAATAATGGCTGTTGTAACGGCGGTCATCATTGTAGCTGATGCTGAATGTGCCCAGCTTCGACCACAGTGAGTTCAGGTTGAGTGTGCCGCCGATTGCACGGTTGTCGGCATCGCTACGTCGCAATTGATTGCCAATGCGGGTTTTTTCCTGATTAACCCA  >  NZ_CP009273/303399‑303660

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: