Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273744,1030TC57.1% 1.8 / 13.6 14G112G (GGA→GGGybgOfimbrial protein
Reads supporting (aligned to +/- strand):  ref base T (3/3);  new base C (4/4);  total (7/7)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.67e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CATATGAATATTAGCGGGATCGCCCGAGCAAATTTTGTTCTGCCATCCCATCGCGATGATCTGCTCAGGGGTATAAATATCAATGTTCTTACTGTCGATACACTGGTTAGTGTCGATGCGACTTTTTCCCAGACTGGCATCGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTACCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGT  >  NZ_CP009273/743977‑744230
                                                                                                                              |                                                                                                                               
cATATGAATATTAGCGGGATCGCCCGAGCAAATTTTGTTCTGCCATCCCATCGCGATGATCTGCTCAGGGGTATAAATATCAATGTTCTTACTGTCGATACACTGGTTAGTGTCGATGCGACTTTTCCCCAGACTGGCa                                                                                                                     <  2:245440/139‑1 (MQ=255)
            aGCGGGATCGCCCGAGCAAATTTTGTTCTGCCATCCCATCGCGATGATCTGCTCAGGGGTATAAATATCAATGTTCTTACTGTCGATACACTGGTTAGTGTCGATGCGACTTTTCCCCAGACTGGCATCGTAATCGaca                                                                                                         >  2:458669/1‑139 (MQ=255)
             gCGGGATCGCCCGAGCAAATTTTGTTCTGCCATCCCATCGCGATGATCTGCTCAGGGGTATAAATATCAATGTTCTTACTGTCGATACACTGGTTAGTGTCGATGCGACTTTTTCCCAGACTGGCATCGTAATCGaca                                                                                                         <  2:424527/138‑1 (MQ=255)
                                                  tcGCGATGATCTGCTCAGGGGTATAAATATCAATGTTCTTACTGTCGATACACTGGTTAGTGTCGATGCGACTTTTCCCCAGACTGGCATCGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTg                                                                   <  1:39659/139‑1 (MQ=255)
                                                       atgatCTGCTCAGGGGTATAAATATCAATGTTCTTACTGTCGATACACTGGTTAGTGTCGATGCGACTTTTTCCCAGACTGGCATCGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCAccc                                                              >  1:160120/1‑139 (MQ=255)
                                                       atgatCTGCTCAGGGGTATAAATATCAATGTTCTTACTGTCGATACACTGGTTAGTGTCGATGCGACTTTTCCCCAGACTGGCATCGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCAccc                                                              <  1:284254/139‑1 (MQ=255)
                                                               cTCAGGGGTATAAATATCAATGTTCTTACTGTCGATACACTGGTTAGTGTCGATGCGACTTTTTCCCAGACTGGCATCGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATa                                                       <  2:160120/138‑1 (MQ=255)
                                                                         taAATATCAATGTTCTTACTGTCGATACACTGGTTAGTGTCGATGCGACTTTTTCCCAGACTGGCATCGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTAcc                                             <  2:416797/138‑1 (MQ=255)
                                                                                                        ggTTAGTGTCGATGCGACTTTTTCCCAGACTGGCATCGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTACCCAGGCGTAGACATGCTCGCTTTCAAAGTTGcc             >  1:52500/1‑139 (MQ=255)
                                                                                                                tCGATGCGACTTTTTCCCAGACTGGAATCGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTACCCAGGAGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTAt     >  1:228695/1‑139 (MQ=255)
                                                                                                                 cGATGCGACTTTTCCCCAGACTGGCATCGTAATCGaca                                                                                                         <  1:79534/38‑1 (MQ=255)
                                                                                                                 cGATGCGACTTTTCCCCAGACTGGCATCGTAATCGaca                                                                                                         >  2:79534/1‑38 (MQ=255)
                                                                                                                   aTGCGACTTTTCCCCAGACTGGCATCGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTACCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTAttgt  >  1:373786/1‑139 (MQ=255)
                                                                                                                   aTGCGACTTTTCCCCAGACTGGCATCGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTACCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTAttgt  >  1:397521/1‑139 (MQ=255)
                                                                                                                              |                                                                                                                               
CATATGAATATTAGCGGGATCGCCCGAGCAAATTTTGTTCTGCCATCCCATCGCGATGATCTGCTCAGGGGTATAAATATCAATGTTCTTACTGTCGATACACTGGTTAGTGTCGATGCGACTTTTTCCCAGACTGGCATCGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTACCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGT  >  NZ_CP009273/743977‑744230

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