Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,548,9320AT32.1% 49.6 / 10.1 29Y42F (TAT→TTT) fimHminor component of type 1 fimbriae
Reads supporting (aligned to +/- strand):  ref base A (2/17);  new base T (8/1);  total (11/18)
Fisher's exact test for biased strand distribution p-value = 1.19e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCC‑TGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCG‑CCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAAT‑CTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATG‑TCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAA  >  NC_000913/4548713‑4549165
                                                                                                                                                                                                                             |                                                                                                                                                                                                                                           
agCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCC‑TGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCG‑CCAACGTTTTTTTAAAACTTGCcccccg                                                                                                                                                                                                                           >  1:1033708/1‑236 (MQ=255)
                     aGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCC‑TGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGGGGGAGCGCCCAACGTTT‑TGTAAACCATGCGCCCCTCGCGATTATGGGGGACAAccag                                                                                                                                                                                                     >  1:932582/1‑237 (MQ=255)
                                     aTTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCC‑TGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGAGGCAGCG‑CCAAAGTTTATGTAAACCTTTGGCCCGCCGTGGAAGTGGGGGAAAACCATGTCGGGGAGCttttg                                                                                                                                                                                      >  1:890435/1‑236 (MQ=255)
                                                taTTAACATCACCTATACCTACAGCTGAACCAGAAGAGATTATTATAAGAAAACGTGTTGTCCCCCTTGTTGCTGTACGGCTGTTGGGCTGGTGGAAAATTGCCTGGT‑ATTCCCCTTGTAAACCCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCG‑CCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCa                                                                                                                                                                                 <  2:1074995/231‑1 (MQ=255)
                                                      cATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCC‑TGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCG‑CCCAAGGTTTTGTAAAACTTTCGCCCGTCGTGGATGGGGGGGAAAACCTGGTCGGGGGTCTTTTGTCGCCAATACTTTTcccc                                                                                                                                                                    >  1:1074995/1‑238 (MQ=255)
                                                                     cAGCTGAACCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTAATGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCC‑TGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGGGGCGGGG‑CCCAAGTTTTTGtc                                                                                                                                                                                                                                         >  1:1430475/1‑154 (MQ=255)
                                                                          gAACCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCC‑TGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCG‑CCAATGTTTATTTAAAc                                                                                                                                                                                                                                      >  1:204856/1‑153 (MQ=255)
                                                                           aaCCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCC‑TGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGGAACG‑CCCATGTTTTTg                                                                                                                                                                                                                                           >  1:1293680/1‑147 (MQ=255)
                                                                                        tgtgttAAAGAAACGTTTTCTTACCCTGTTTGGTGTACTGCTGAGGGGTTGTGCGAAAAATGCTGGGCAATCCCCC‑TGAAAACCCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCG‑CCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAAT‑CTTTTGCCATAACGATTATCCGGAAACCATTACa                                                                                                                                           <  2:1212884/224‑1 (MQ=255)
                                                                                                        tttattcCCTTTTTGGTTGTCCTGTTGGGGGGCTGGTCGGAAAATGCCTGGTCTTTCGCC‑TGAAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCG‑CCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAAT‑CTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATG‑TCACACTGCAACGAGGCt                                                                                                                 <  2:285351/237‑1 (MQ=255)
                                                                                                             taCCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCC‑TGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGGGGGAGGG‑GCCAAGTTTTTGTaaa                                                                                                                                                                                                                                       >  1:1356473/1‑117 (MQ=255)
                                                                                                                     ttGCTGTACTGCTGATGGCCTGGTCGGAAAATGCCTGTCCTTTCGCG‑TGAAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCG‑CCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAAT‑CTTTTGCCATAACGATTATCCGGAAACCATTACAGCCTATG‑TCACACTGCAACGAGGCTCGGCTTATggcggc                                                                                                   <  1:83489/239‑1 (MQ=255)
                                                                                                                      tGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCC‑TGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGGAGAG‑CCAATGTTTTTGTAAAACTTgcgc                                                                                                                                                                                                                               >  1:908181/1‑116 (MQ=255)
                                                                                                                      ggCTGTACGGCTGAGGGGGGGGTCGGAAAGTCCCGGTTCTTTCGCT‑TGAAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCG‑CCAATGTTTATGTAAACCTTGCGCCCGTCGTGAAtgt                                                                                                                                                                                                                  <  2:192189/128‑1 (MQ=255)
                                                                                                                             ctgtggaggggCTGGTCGGAAATTGCCTGTTCTTTCGCC‑TGAAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCG‑CCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAAT‑CTTTTGCCATAACGATTATCCGGaaa                                                                                                                                                   <  1:572520/176‑1 (MQ=255)
                                                                                                                                ctgATGGGCTGGTCGGTAAATGCCTGGTCATTCGCC‑TGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCC‑CCCATGTTTTTGTAAACCTTTCGCCCGGCGTGAATGTGGGGCAAAAc                                                                                                                                                                                                        >  1:192189/1‑129 (MQ=255)
                                                                                                                                ctgAGGGGCGGGTCGTTAAATGCCTGTTCATCCCCC‑TATAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCG‑CCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAAT‑CTTTTGCCATAACGATTATCCGGAAACCATTACa                                                                                                                                           <  1:1082436/189‑1 (MQ=255)
                                                                                                                                            tgggaaaaTGCCGGGTTATTCCCC‑GGAAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCG‑CCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAAT‑CTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATG‑TCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                             <  1:450837/233‑1 (MQ=255)
                                                                                                                                            tCGGTAAATGCCTGGTCATTCGCC‑TTTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCG‑CCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAAT‑CTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATG‑TCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                             <  1:86402/238‑1 (MQ=255)
                                                                                                                                            gcgtaaactGCGTGGTCATCGGCT‑TAAAAAGCCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCG‑CCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAAT‑CTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATG‑TCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                             <  2:1033708/230‑1 (MQ=255)
                                                                                                                                                tAAATGCCTGGTCATTCGCC‑TGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCGGCG‑CCCATATTTTTTTTAACCCTTCCCCCCGCGTGGATTGGGGGGGAAAACAGGGCGGGGGTTTTTTCTCGCCAAT‑CTTTTTTCCTAACAAATTTTCCGAAAACATTACAGACTTTGTTCACACT‑CCACAAGGGGCGGCGTCTTGCGGCGTGTGTTATCAttttttttgggc                                                                            >  2:450837/1‑228 (MQ=255)
                                                                                                                                                                      taaaagCCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCG‑CCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAAT‑CTTTTGCCATAACGATTATCCGGAAACCTTTACAGACTATG‑TCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTAt                                                   <  2:1036302/237‑1 (MQ=255)
                                                                                                                                                                                    ggTACCGCTATCCCTATTGGCGGTGGCAGCG‑CCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAAT‑CTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATG‑TCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCa                                                           <  1:36259/217‑1 (MQ=255)
                                                                                                                                                                                           cTATCCCTATTGGCGGTGGCAGCG‑CCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAAT‑CTTTTGCCATAACGATTATCCGGAAACCATTACAGa                                                                                                                                         <  1:645326/133‑1 (MQ=255)
                                                                                                                                                                                                    ttGGCGGTGGCAGCG‑CCAATGTTTTTGTAAGCCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAAT‑CTTTTGCCATAACGATTTTCCGGAAACCATTACAGACTATG‑TCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCcgcgc                     <  1:464060/239‑1 (MQ=255)
                                                                                                                                                                                                     tGGCGGGCGGAGCG‑GAAATGTGTATGTAAGCCTTGCGGCGGTCGTGAAGGTGGGGCAAAGCCTGGTCGTAGATGTGTCGACGGCAAT‑CTTTTGCCAGAAAGATTATCGGGAAGGCGTTGGAGGCGATG‑TCAGGGGGGGGCGGGGGGGGCCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTTGCTGTAGCTATCCATTTCCTACCACCAGCGAAACGCcgcggt                    <  2:246675/239‑3 (MQ=255)
                                                                                                                                                                                                        cGGTGGCAGCG‑CCAATGTTTATGTAAACCTTTCGCCCTTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAAT‑CTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATG‑TCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCgtt                  <  2:542754/238‑1 (MQ=255)
                                                                                                                                                                                                                      aaTGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAAT‑CTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATG‑TCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCgag    <  1:453529/239‑1 (MQ=255)
                                                                                                                                                                                                                        tGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAAT‑CTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATG‑TCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGaa  <  1:232110/239‑1 (MQ=255)
                                                                                                                                                                                                                             |                                                                                                                                                                                                                                           
AGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCC‑TGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCG‑CCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAAT‑CTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATG‑TCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAA  >  NC_000913/4548713‑4549165

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑