Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,131,363 G→A 7.4% A279V (GCC→GTC)  fimD ← type I pili subunit FimD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,131,3630GA 7.4% 156.8 / 2.3 54A279V (GCC→GTC) fimDtype I pili subunit FimD
Reads supporting (aligned to +/- strand):  ref base G (21/29);  new base A (2/2);  total (23/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.38e-01

TACTCAGGTCGTCGATGGCGTAGGGGCCGGCGCTGACGTACGTCGAGTAGATCGGATAACCATGCTGCCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTT‑TCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTC  >  NC_002947/2131226‑2131506
                                                                                                                                         |                                                                                                                                                
tACTCAGGTCGTCGATGGCGTTGGGGCCGGCGCTGACGTACGTCGAGTAGATCGGATAACCATGCTGCCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGccc                                                                                                                                       <  3:186053/149‑1 (MQ=255)
     aGGTCGTCGATGGCGTAGGGGCCGGCGCTGACGTACGTCGAGTAGATCGGATAACCATGCTGCCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTc                                                                                                                                  >  2:218742/1‑149 (MQ=255)
      ggTCGTCGATGGCGTAGGGGCCGGCGCTGACGTACGTCGAGTAGATCGGATAACCATGCTGCCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTcc                                                                                                                                 >  4:86652/1‑149 (MQ=255)
       gtcgtcGATGGCGTAGGGGCCGGCGCTGACGTACGTCGAGTAGATCGGATAACCATGCTGCCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTcc                                                                                                                                 >  1:350944/1‑148 (MQ=255)
        tcgtcgATGGCGTAGGGGCCGGCGCTGACGTACGTCGAGTAGATCGGATAACCATGCTGCCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGa                                                                                                                               <  3:86652/149‑1 (MQ=255)
        tcgtcgATGGCGTAGGGGCCGGCGCTGACGTACGTCGAGTAGATCGGATAACCATGCTGCCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGa                                                                                                                               <  6:157077/149‑1 (MQ=255)
           tcgATGGCGTAGGGGCCGGCGCTGACGTACGTCGAGTAGATCGGATAACCATGCTGCCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGc                                                                                                                            <  5:276492/149‑1 (MQ=255)
             gATGGCGTAGGGGCCGGCGCTGACGTACGTCGAGTAGATCGGATAACCATGCTGCCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGACCAGCATG‑CCATGTCCGAGGCCAc                                                                                                                         >  5:172964/1‑149 (MQ=255)
             gATGGCGTAGGGGCCGGCGCTGACGTACGTCGAGTAGATCGGATAACCATGCTGCCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGACCAGCATG‑CCATGTCCGAGGCCAc                                                                                                                         >  7:228358/1‑149 (MQ=255)
                           cGGCGCTGACGTACGTCGAGTAGATCGGATAACCATGCTGCCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCa                                                                                                                          <  3:191700/135‑1 (MQ=255)
                           cGGCGCTGACGTACGTCGAGTAGATCGGATAACCATGCTGCCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCa                                                                                                                          >  4:191700/1‑135 (MQ=255)
                              cgcTGACGTACGTCGAGTAGATCGGATAACCATGCTGCCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATc                                                                                                         >  5:97330/1‑149 (MQ=255)
                                       acgtCGAGTAGATCGGATAACCATGCTGCCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGaa                                                                                                >  1:325773/1‑149 (MQ=255)
                                        cgtCGAGTAGATCGGATAACCATGCTGCCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAAc                                                                                               <  6:175331/149‑1 (MQ=255)
                                                  aTCGGATAACCATGCTGCCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACgg                                                                                             <  6:258165/141‑1 (MQ=255)
                                                  aTCGGATAACCATGCTGCCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACgg                                                                                             >  5:258165/1‑141 (MQ=255)
                                                  aTCGGATAACCATGCTGCCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTc                                                                                     <  6:147625/149‑1 (MQ=255)
                                                   tCGGATTACCATGCTGCCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCGACGCGCACGCCTGATATCGGCACCGAACGGAACACCTcc                                                                                    <  3:60111/149‑1 (MQ=255)
                                                   tCGGATAACCATGCTGCCAGATCTCGACTTTCGCCCGGCTCTGGTCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTcc                                                                                    <  6:235814/149‑1 (MQ=255)
                                                   tCGGATAACCATGCTGCCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTcc                                                                                    <  5:275278/149‑1 (MQ=255)
                                                   tCGGATAACCATGCTGCCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTcc                                                                                    <  5:355350/149‑1 (MQ=255)
                                                   tCGGATAACCATGCTGCCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTcc                                                                                    >  5:253453/1‑149 (MQ=255)
                                                   tCGGATAACCATGCTGCCAGATCTCGACTCTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTcc                                                                                    <  5:334866/149‑1 (MQ=255)
                                                        tAACCATGCTGCCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTc                                                                               <  3:83955/149‑1 (MQ=255)
                                                            cATGCTGCCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGt                                                                            >  5:108257/1‑148 (MQ=255)
                                                            cATGCTGCCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGt                                                                            >  5:183810/1‑148 (MQ=255)
                                                              tgctgcCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGaa                                                                         <  2:193917/149‑1 (MQ=255)
                                                                ctgcCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAg                                                                        <  6:108257/148‑1 (MQ=255)
                                                                  gcgaGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGtt                                                                     >  4:237333/4‑149 (MQ=255)
                                                                     aGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTT‑TCg                                                                 >  2:298599/1‑149 (MQ=255)
                                                                           cGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTT‑TCGCCAATc                                                           >  5:363878/1‑149 (MQ=255)
                                                                            gACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTT‑TCGCCaa                                                             <  2:51476/146‑1 (MQ=255)
                                                                            gACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTT‑TCGCCAATc                                                           <  6:147562/148‑1 (MQ=255)
                                                                            gACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGACCAGCATG‑CCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTT‑TCGCCAATc                                                           <  8:228358/147‑1 (MQ=255)
                                                                            gACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGACCAGCATG‑CCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTT‑TCGCCAATc                                                           <  6:172964/147‑1 (MQ=255)
                                                                                        gCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTT‑TCGCCAATCGTCAGATTGGCGa                                              <  6:253453/149‑1 (MQ=255)
                                                                                                ccACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTT‑TCGCCAATCGTCAGATTGGCGAg                                             <  4:73246/142‑1 (MQ=255)
                                                                                                ccACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTT‑TCGCCAATCGTCAGATTGGCGAg                                             >  3:73246/1‑142 (MQ=255)
                                                                                                 cACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTT‑TCGCCAATCGTCAGATTGGCGAGGGTACCCg                                     >  3:302502/1‑149 (MQ=255)
                                                                                                     ccgcGAATGATCGGGGCATAACCGCGCGGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTT‑TCGCCAATCGTCAGATTGGCGAGGGTACCCGGtat                                 <  1:218742/149‑1 (MQ=255)
                                                                                                      cgcgAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTT‑TCGCCAATCGTCAGATTGGCGAGGGTACCCGGtata                                <  3:237333/149‑1 (MQ=255)
                                                                                                         gAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTT‑TCGCCAATCGt                                                         <  7:341998/121‑1 (MQ=255)
                                                                                                         gAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTT‑TCGCCAATCGt                                                         >  8:341998/1‑121 (MQ=255)
                                                                                                         gAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTT‑TCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATcg                             >  5:315874/1‑149 (MQ=255)
                                                                                                                       tAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTT‑TCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGt               <  6:97330/149‑1 (MQ=255)
                                                                                                                         aCCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTT‑TCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCt             <  6:183810/149‑1 (MQ=255)
                                                                                                                          ccGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTT‑TCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCtg            >  4:34276/1‑149 (MQ=255)
                                                                                                                           cgcgCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTT‑TCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCtgt           >  7:158163/1‑149 (MQ=255)
                                                                                                                               cTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTT‑TCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGccc       <  4:302502/149‑1 (MQ=255)
                                                                                                                                  gTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTCTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCgg     <  1:298599/149‑1 (MQ=255)
                                                                                                                                  gTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTT‑TCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCggg    <  6:363878/149‑1 (MQ=255)
                                                                                                                                  gTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTT‑TCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCggg    >  5:244477/1‑149 (MQ=255)
                                                                                                                                     gTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTT‑TCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTc  >  4:230967/1‑148 (MQ=255)
                                                                                                                                     gTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTT‑TCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTc  <  2:350944/148‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                                
TACTCAGGTCGTCGATGGCGTAGGGGCCGGCGCTGACGTACGTCGAGTAGATCGGATAACCATGCTGCCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTT‑TCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTC  >  NC_002947/2131226‑2131506

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