Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913559,1810GA31.0% 33.7 / 10.7 29E70K (GAA→AAA) sfmCputative fimbrial chaperone SfmC
Reads supporting (aligned to +/- strand):  ref base G (12/8);  new base A (0/9);  total (12/17)
Fisher's exact test for biased strand distribution p-value = 3.29e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTATGATGACTAAAATAAAGTTATTGATGCTCATTATATTTTATTTAATCATTTCGGCCAGCGCCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACC‑AATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGC‑AGCGGTGTAAAAGAAAA‑GT‑CATTCATCA‑TTACACC‑GCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTA  >  NC_000913/558936‑559417
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                
ttATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTATGCTGACTAAAATAAAGTTATTGCTGCTCCTTATATTTTATTTAATCATTTCGTCCAGCGCCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACC‑AATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGC‑AGCGGTGTaaaagaaa                                                                                                                                                                                                                                               >  2:321258/1‑249 (MQ=255)
ttATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTATGATGACTAAAATAAAGTTATTGATGCTCATTATATTTTATTTAATCATTTCGGCCAGCGCCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACC‑AATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGC‑AGCGGTGTcacagaaa                                                                                                                                                                                                                                               >  2:88977/1‑249 (MQ=255)
                      acaAGGCGTATTAACTATGATGACAAAAATAAAGTTATTTATGCCCTTTATTTTTTATTTAATAGTTCCGGCCAGCGCCCATGCTGCCGGAGGGATCGCTTTAGTTGCCACGGGTTTTATTTACCCCGCTGATGCTAAACAGACTGCGGTGTGGTTAAGAAATAGCCATCCC‑AATGAGCGCTTTCGGGTCAATGCGTTGAGTGAAAACAGC‑AGCGTTGTAAAAGAAAG‑GT‑CTTTCATGA‑TTACACC‑GCCAc                                                                                                                                                                                                                   <  2:442859/251‑1 (MQ=255)
                       caAGGCGTATTAACTATGATGACTAAAATAAAGTTATTGATGCTCATTATATTTTATTTAATCATTTCGGCCAGCGCCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACC‑AATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGC‑AGCGGTGTAAAAGAAAA‑GT‑CATTCATCA‑TTACACC‑GCCACt                                                                                                                                                                                                                  >  1:362591/1‑251 (MQ=255)
                             gTATTAACTATGATGACTAAAATAAAGTTATTGATGCTCATTATATTTTATTTAATCATTTCGGCCAGCGCCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACC‑AATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGC‑AGCGGTGTAAAAGAAAA‑GT‑CATTCATCA‑TTACACC‑GCCACtgtttgt                                                                                                                                                                                                            >  1:143679/1‑251 (MQ=255)
                                                aaaaaaatttATTGATGTCCATTATATGTTATTTAATCTTTTCGGCCAGCGCCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGGTTATTATTTATCCCGCTGATGCAAAACAGACTGCGGTAGGGATTAGAAATACCAATACC‑AATGAGCGTTTTCTGGTCAATTCGTGGATTGAAAACAGC‑AGGGGTGTAAAAGAAAA‑GT‑CATTCATCA‑TTACACC‑GCCACTGTTTGTTAGTGAACCCAAAAGCGaa                                                                                                                                                                                         <  1:338385/243‑1 (MQ=255)
                                                                      tataTTTTATTTAATCATTTCGGCCAGCGCCCATGCTGCCGGAGGGACCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGCCTGGGGTATGTATTAGAAATAGCCATACC‑AATGAGCGTTTTGTGGTCAATTCGTGTATTGAAAACAGC‑AGCGGTGTAAAAGAAAA‑GT‑AATTCATCA‑TTCCACC‑GCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACGTTGCGTATTATTTac                                                                                                                                                                    <  1:377387/250‑1 (MQ=255)
                                                                        tattttatttAATCATTTCGGCCAGCGCCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACC‑AATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGC‑AGCGGTGTAAAAGAAAA‑GT‑CATTCATCA‑TTACACC‑GCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACAcc                                                                                                                                                                 >  2:366361/1‑251 (MQ=255)
                                                                                     cATTTCGGCCAGCGCCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACC‑AATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGC‑AGCGGTGTAAAAGAAAA‑GT‑CATTCATCA‑TTACACC‑GCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTgg                                                                                                                                                    >  2:114323/1‑251 (MQ=255)
                                                                                                                       gCATTAGGTGCCACGCGTCTTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGCTTAGAAATAGCCATACC‑AATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGC‑AGCGGTGTAAAAGAAAA‑GT‑CATTCATCA‑TTACACC‑GCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGCCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGt                                                                                                                  >  2:330828/1‑251 (MQ=255)
                                                                                                                                 ccACGCGTATTATTTATCCCGCTGATGCAAAACAGACGGCGGTATGGATTAGAAATAGCCATACC‑AATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGC‑AGCGGTGAAAAAGAAAA‑GT‑CATTCATCA‑TTCCACC‑GCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATc                                                                                                        <  1:366361/251‑1 (MQ=255)
                                                                                                                                      cgTATTATTTACCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCC‑TACCTAATGAGCGCTTTCTGGTCAATTCTTGGATT‑AAAACCGCCAGCGGTTAAAAAAAAAA‑TT‑CATTAATCA‑TACCTCC‑GCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCtt                                                                                                    <  1:104143/250‑1 (MQ=255)
                                                                                                                                                        gattaaaaccAGACTGCGTTATGGAGTAAAAAAATCAATACA‑AAGAAGCCCTTTCTGGTCAATTCGTGGATTAAAAACAGC‑ATCGTTGTAAAAAAAAA‑CT‑CATTCATTA‑TT‑CACCCGCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTTAGTCTCTGGTCTGGATGAATGTGAAGACGACCTCTTCGGTAGATAAAAATGCAtt                                                                                 <  2:47781/242‑1 (MQ=255)
                                                                                                                                                             aaaaCAGACTGGGGTATGAATAAGAAATAGCAAAACC‑AATGAGCGTTTTCTGGTCATTTCGTGGATTGAAAACAGC‑AGCGGTGTAAAAAAAAA‑GT‑CATCCTTCA‑TTACACC‑GCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACg                                                                            <  1:114323/250‑1 (MQ=255)
                                                                                                                                                                    aCTGCGGTATGGATTAGAAATAGCCATACC‑AATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGC‑AGCGGTGTAAAAGAAAA‑GT‑CATTCATCA‑TTACACC‑GCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGa                                                                      >  1:448230/1‑250 (MQ=255)
                                                                                                                                                                      tGCGGTATGGATTAGAAATAGCCATACC‑AATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGC‑AGCGGTGTAAAAGAAAA‑GT‑CATTCATCA‑TTACACC‑GCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATg                                                                   >  1:417968/1‑251 (MQ=255)
                                                                                                                                                                      ggggttgtGGATTAAAAATAGCCATACA‑AATGGGCGCTTGCTGTTAATTTCGTGGATTGAAAACAGC‑AGGGGGGTAAAAAAAAA‑GT‑CATTCATAA‑TTACACC‑GCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATg                                                                   <  2:393721/244‑1 (MQ=255)
                                                                                                                                                                           tATGGATTAGAAATAGCCATACC‑AATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGC‑AGCGGTGTAAAAGAAAA‑GT‑CATTCATCA‑TTACACC‑GCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGaa                                                                     >  1:150248/1‑244 (MQ=255)
                                                                                                                                                                           aaTGGATTAGAAATAGCCATACC‑AATGAGCGCTTTTTGGTCAATTGGTGGATGAAAAACAGC‑AGCGGTGAAAAAGAAAA‑GT‑CATTCATCA‑TTACACC‑GCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGaa                                                                     <  2:150248/243‑1 (MQ=255)
                                                                                                                                                                             gggATTAGAAATACATATCCC‑AAGGAGCGCTTGCTGGTAAATTCGTTGATTAAAAACACC‑TGCGGTGTAAAAGAAAA‑GT‑CAG‑CAACAATTACACC‑GCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGc                                                             <  2:171738/249‑1 (MQ=255)
                                                                                                                                                                               gATTAGAAATAGCCATACC‑AATGCGCGCTTTCTGGTCAATTCGTGGATTTAAACCAGC‑AGCGGTGTAAAAGCAAA‑GT‑CCTTCATCA‑TTACACC‑GCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAACAGATCGTGAGTCTCTGTTCTGGATGAGTGTTAAGACTATCCCTTCGGTAAATAAAAATACATTGAACGGCCGGAATGTGTTTCAAc                                                          >  2:447218/1‑251 (MQ=255)
                                                                                                                                                                                  tAGAAATAGCCTTCCC‑AATGAGCGTTTTCTGGTCAATTCGTGGATGAAAAACAGC‑AGCGGTTAAAAAAAAAA‑GT‑CATTCATCA‑TTACACC‑GCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTg                                                        <  1:107314/250‑1 (MQ=255)
                                                                                                                                                                                  tAGAAAAAGACATACA‑AATGAGCTTTTTCTGGTAATTTGGTGGATTGAAAACAGC‑AGGGGT‑TAAAAAAAAA‑ATACATTAATCA‑TTACCCC‑GCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTg                                                        <  2:162097/250‑1 (MQ=255)
                                                                                                                                                                                        aaGCCATACC‑AATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGC‑AGCGGTGAAAAAGAAAA‑GT‑CATTCATCA‑TTACACC‑GCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGAttt                                                 <  1:439941/250‑1 (MQ=255)
                                                                                                                                                                                                     ggAGCGCTGTCGGGTAAATTCGTGAATGGAAAAC‑GG‑AGGGGTGTAAAAAAAAAATT‑CATTCATCA‑TTCCTCC‑GCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGa                                     <  2:77649/250‑1 (MQ=255)
                                                                                                                                                                                                      gAGCGTTTTCGGGTCAATTCGGGGATTGAAAACAGC‑AGGGGTGTAAAAGAAAA‑AT‑AATTCATCA‑TTACACC‑GCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTGACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGGTCTGGATTAATGTTAAGGCGATCCCGTCTGTAGATAAAAATGCAGTGAACGGCAGGAATGTTTTTCAACTGGCGATTTTagggctcatgaa                                    <  2:444487/251‑12 (MQ=255)
                                                                                                                                                                                                               cTGGTCAATTCGTGGATTGAAAACAGC‑AGCGGTGTAAAAGAAAA‑GT‑CATTCATCA‑TTACACC‑GCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGCTAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTtct                           >  2:213311/1‑251 (MQ=255)
                                                                                                                                                                                                                tGGTCAATTCGTGGATTGAAAACAGC‑AGCGGTGTAAAAGAAAA‑GT‑CATTCATCA‑TTACACC‑GCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTtctc                          >  1:432725/1‑251 (MQ=255)
                                                                                                                                                                                                                       ttCGTGGATTGAAAACAGC‑AGCGGTGTAAACGAAAA‑GT‑CATTCATCA‑TTACACC‑GCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCca                    >  1:204965/1‑250 (MQ=255)
                                                                                                                                                                                                                                      cagc‑agGGGTGTAAAAAAAAA‑GT‑CATCCATCA‑TTACACC‑GCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAAt    <  1:84285/251‑1 (MQ=255)
                                                                                                                                                                                                                                        gc‑agcGGTTTAAAAAAAAA‑GT‑CTTTCATCA‑TTACACC‑GCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTa  <  1:294027/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                
TTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTATGATGACTAAAATAAAGTTATTGATGCTCATTATATTTTATTTAATCATTTCGGCCAGCGCCCATGCTGCCGGAGGGATCGCATTAGGTGCCACGCGTATTATTTATCCCGCTGATGCTAAACAGACTGCGGTATGGATTAGAAATAGCCATACC‑AATGAGCGCTTTCTGGTCAATTCGTGGATTGAAAACAGC‑AGCGGTGTAAAAGAAAA‑GT‑CATTCATCA‑TTACACC‑GCCACTGTTTGTTAGTGAACCCAAAAGCGAAAATACTTTGCGTATTATTTACACCGGTCCACCGCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTA  >  NC_000913/558936‑559417

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