Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,548,401 T→C 12.1% C39C (TGT→TGC) ‡ fimG → type 1 fimbriae minor subunit FimG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,548,4010TC12.1% 109.6 / 8.3 50C39C (TGT→TGC) ‡fimGtype 1 fimbriae minor subunit FimG
Reads supporting (aligned to +/- strand):  ref base T (23/20);  new base C (4/2);  total (27/23)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.40e-01

AATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGTTGACTAATTGTCCGGTGGGAACGTCGAGGGTCACTGCCA  >  NC_000913/4548262‑4548543
                                                                                                                                           |                                                                                                                                              
aaTATCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCGCACAAGTGCGAAGATACAGGCATCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACg                                                                                                                                             <  1:784765/143‑1 (MQ=255)
   atCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGtt                                                                                                                                          >  2:938619/1‑143 (MQ=255)
     cAGTAACTGGAGATGGGCAAAAAAGGGTGCAAACGTGGGTATGTATTGGCGGCAATAGTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTACACGGtttg                                                                                                                                        <  1:825363/143‑2 (MQ=255)
     cAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTACACGGtttg                                                                                                                                        <  1:458615/143‑2 (MQ=255)
        tAACTGGAGATGCTCATGAACTGGTGTATACGTGAGAATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTccac                                                                                                                                     <  1:955157/143‑1 (MQ=255)
          aCTGGGGGTGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTccacca                                                                                                                                   <  1:316978/143‑1 (MQ=255)
            tGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGCACGGTAAGGTAGTCGCCAAACCGTGTACGGTTTccaccacc                                                                                                                                 >  1:384036/1‑143 (MQ=255)
                   gCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCa                                                                                                                          >  1:703963/1‑143 (MQ=255)
                     tCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGGCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTAGGGTCTCCACAACCAACGGCAcc                                                                                                                        >  2:109911/1‑142 (MQ=255)
                       aTGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATGGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGt                                                                                                                      <  2:319618/143‑1 (MQ=255)
                           aaTGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGAt                                                                                                                  <  1:476491/143‑1 (MQ=255)
                           aaTGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGGGAACGGTAAGGTCGTCGCCCAACCGTGTACCGGTTCCACCACCACAGCCACGGTTgct                                                                                                                  >  1:33821/1‑141 (MQ=255)
                                                tatTGGCGGCAATATTGGCGCTCGCAAGAGCGACGATACAGGCAGCCGATGTCACCATCACGGTGGACGGTGAGGTCGTCGCCAAACGGTGTACGGGTTCCACCACCAATGCCACGGTAGACATCGGCGGTATTTATCCtcta                                                                                             >  2:91881/1‑140 (MQ=255)
                                                   tGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGt                                                                                          <  2:703963/143‑1 (MQ=255)
                                                    ggcggcAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCGCCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTCTCCACCACCGATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTc                                                                                         >  1:99844/1‑143 (MQ=255)
                                                     gcggcAATATTGGCGCTCGCAAGTGCGACGATACAGGCCGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGGTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCCGTcc                                                                                        >  2:405099/1‑142 (MQ=255)
                                                       ggcAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTa                                                                                      <  1:767811/143‑1 (MQ=255)
                                                       ggcAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGCCGATCTTTATTCTTTCAGtgcga                                                                                      >  2:173590/1‑139 (MQ=255)
                                                       ggcAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTGGATCTCGGCGATCTTTATTCTTTCAGTCTTa                                                                                      >  1:576988/1‑143 (MQ=255)
                                                             atTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTAGGGTTTCCACCACCAATGCCAAGGATGATCTCGTCGATCGTTATTCTTTCAGCCTTATGTCTg                                                                                >  1:540478/1‑143 (MQ=255)
                                                                 gcgcTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTGTCCACCACCAAGGCCACGGTTGATCTCGGCGATCTTTATTCTTTTAGTCTTATGTCTGCCgg                                                                            >  1:535797/1‑143 (MQ=255)
                                                                  tcttCGCAAGTCGGACGATACAGCCAGCCGATGTCGCCAGCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCggg                                                                           <  2:756256/140‑1 (MQ=255)
                                                                    ctggCAAGTGCGGCGAGACAGGCAGTCGATGTCTCCAGCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGggc                                                                         <  2:402860/140‑1 (MQ=255)
                                                                    ctggCAAGTGCGCCGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGggc                                                                         <  1:1036342/140‑1 (MQ=255)
                                                                    cTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCGCGGCGAGGTTTGTTCTTTCAGTCTTATGTCTGCCGGggc                                                                         >  1:310911/1‑143 (MQ=255)
                                                                     tCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCCTGTACGGTTTCCACCACCAATGCCCCACTTGATCGTTGCGAACACAATTCTTTGAGTCTTATGTCTGCCGGggcg                                                                        >  1:488507/1‑143 (MQ=255)
                                                                       gCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGggcggc                                                                      <  1:204609/143‑1 (MQ=255)
                                                                               gacgaTACAGGCCGCCCATGCCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGGTGATCTCGGCGATCTTTATTCTTTCAGCCTTATGTCTGCCGGGGCGGCATCGGCCt                                                              <  1:45277/143‑1 (MQ=255)
                                                                                        tgCACCCGACGTGCCCATCCCGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTGCCACCACCAACGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCatgatg                                                     <  2:783530/142‑1 (MQ=255)
                                                                                            gCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTgc                                                 >  1:454400/1‑143 (MQ=255)
                                                                                             gtgATGTGATCATAACCGTGAACGGTAAGGGCGTCGCCAGAGCGTGGACGGTCTCCACCACCAAGGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTgcg                                                <  1:463892/141‑1 (MQ=255)
                                                                                             ccGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTCAGTCTGCCGGGGCGGCCTCGGCCTGGCATGATCTTgcg                                                >  2:133530/1‑143 (MQ=255)
                                                                                                  gtgcccATCACGGTGAGCGGTAAGGTCGGCGGCAAGCCGTGTACGGTTTCCACCGCAGGCGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGa                                           <  2:122490/139‑1 (MQ=255)
                                                                                                    cACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGGTGTTGCGCTTGAgg                                         >  1:103763/1‑142 (MQ=255)
                                                                                                     aCCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCAACGGTTGATCTCTGCGATCTTTATTCTTTCGGACTTATGTCTGCCGGGGAGACATCGGCCGGGCTTGATGTGGGGCTTggctt                                        >  1:343765/1‑139 (MQ=255)
                                                                                                        aTCACGGTGAACGGTAAGGTCGTCGCCAAACCTTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGTTGAc                                     <  2:707007/143‑1 (MQ=255)
                                                                                                                      tAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGTTGACTaa                                  <  1:924926/132‑1 (MQ=255)
                                                                                                                      tAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGTTGACTaa                                  >  2:924926/1‑132 (MQ=255)
                                                                                                                          gtcgtcGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTGTCAGTCATATGTATGCAGGGGCGGCACCGTCCTGGCATGATGCTGCGCTTGAGTTGACTCATTGTCCGGTGGGagc                   >  2:252126/1‑141 (MQ=255)
                                                                                                                              tcgCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGACTCATGTCTGCCGGGGCGGCATCGGCCTGGCATGCTGTTGCGCTTGGGTTGACTAATTGTCCGGTGGGAAGGTCg               >  2:337778/1‑143 (MQ=255)
                                                                                                                               tgCCAAACCGTGTAGGGTTGCCACCATCAATCCCGCTGTTGATCTCGGCGATCCTTATTCTCTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGAGGTTGCGCTTGAGTTGACTAATTGTCCGGTGCGAACGTCGa              <  2:602819/142‑1 (MQ=255)
                                                                                                                               gtggaaaCCGTACACGGTTTCCACCACCAATGCCACGGTTGATCCCGGCGATATTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGTTGACTAATTGTCCGGTGGGAACGTCGa              >  1:130788/5‑143 (MQ=255)
                                                                                                                                tggaaaCCGTACACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGTTGACTAATTGTCCTGTGGTAACGTCGAg             >  1:657702/4‑143 (MQ=255)
                                                                                                                                tggaaaCCGTACACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGTTGACTAATTGTCCGGTGGGCACGTCGAg             >  1:132062/4‑143 (MQ=255)
                                                                                                                                 ggaaaCCGTACACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGTTGACTAATTGTCCGGTGGGAACGTCGAgg            >  1:178869/3‑143 (MQ=255)
                                                                                                                                 ccAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCGTCGGCCTGGCATGATGTTGCGCTTGAGTTGACTAATTGGCCGGTGGGAACGACGAgg            >  1:159493/1‑143 (MQ=255)
                                                                                                                                 ccAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGTTGACTAATTGTCCGGTGGGAACGTCGAgg            <  1:938619/143‑1 (MQ=255)
                                                                                                                                 ccAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTGACAGGCATATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGTTGACTAATTGTCCGGTGGGAACGTCGAgg            <  1:243652/143‑1 (MQ=255)
                                                                                                                                 cagaggcGTGTGCGGTTGCCACCGCCGATGCCACGGGTGATCTCGGCGATCTTTATTCGTTCAGTCTTATGTCTGCCGGGGCGGCATCGGGCTGGCATGATGTTGCGCTTGAGTTGACTAATGGTCCGGTGGGAACGTCGAgg            <  2:361727/137‑1 (MQ=255)
                                                                                                                                    aaCCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGCTGTTGCGCTTGAGTTGGCTCATTGTCCGGTGGGAACGTCGGGGGTc         >  2:1001304/1‑143 (MQ=255)
                                                                                                                                      gtcggtAGGGTTGCCACGGCCAGTGCCGCGGTTGATCTCGGCGATCGTTATTCTTTCGGTCTTATGTCTGCCGGGGCGGCATCGGGCTGGCATGATGTTGCGCTTGAGTTGACTAATTGTCCGGTGGGAACGTCGAGGGTCAc       <  2:100744/139‑1 (MQ=255)
                                                                                                                                           tACGGTTTCCACGCCCGGTGCCACGGGTGAGCTCGGCGATCTGTATTCTTTCAGTCTTATGTGTGCCGGGGCGGCATCGGCCGGGCATGATGTTGCGCTTGAGGTGACTAATTGTCCGGTGGGAACGTCGAGGGTCACTGCCa  <  2:717897/143‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                              
AATATCAGTAACTGGAGATGCTCATGAAATGGTGCAAACGTGGGTATGTATTGGCGGCAATATTGGCGCTCGCAAGTGCGACGATACAGGCAGCCGATGTCACCATCACGGTGAACGGTAAGGTCGTCGCCAAACCGTGTACGGTTTCCACCACCAATGCCACGGTTGATCTCGGCGATCTTTATTCTTTCAGTCTTATGTCTGCCGGGGCGGCATCGGCCTGGCATGATGTTGCGCTTGAGTTGACTAATTGTCCGGTGGGAACGTCGAGGGTCACTGCCA  >  NC_000913/4548262‑4548543

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑