Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,262,076 C→T 7.1% L888L (CTC→CTT)  fimV → motility protein FimV

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,262,0760CT 7.1% 165.7 / 2.0 56L888L (CTC→CTT) fimVmotility protein FimV
Reads supporting (aligned to +/- strand):  ref base C (20/32);  new base T (2/2);  total (22/34)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.73e-01

GCGCAACCAGAGCCTGAGCCTTTGGACGACCTGGACTTCGACTTCTTCTCCGGCAGTGACGAAGTCGCCACCAAGCTCGACCTGGCCCGTGCCTATATCGACATGGGTGACAGCCAGGGCGCGCGCGACATCCTCGACGAAGTGGTCAAGGATGGCGACGACAGCCAGCGTCAGGAAGCCGAGGACATGCTCTCCCGGCTGGTCTGATGCAAGCAGGTCAGGCCTCTTGGCGGGCTGACCGCGCCCGCGAAGAGGCCCGCGCATGCAGTGCGCAATCA  >  NC_002947/2261942‑2262219
                                                                                                                                      |                                                                                                                                               
gcgcAACCAGAGCCTGAGCCTTTGGACGACCTGGACTTCGACTTCTTCTCCGGCAGTGACGAAGTCGCCACCAAGCTCGACCTGGCCCGTGCCTATATCGACATGGGTGACAGCCAGGGCGCGCGCGACATCCTCGACGAAGTGGTc                                                                                                                                     >  6:104018/1‑147 (MQ=255)
 cgcAACCAGAGCCTGAGCCTTTGGACGACCTGGACTTCGACTTCTTCTCCGGCAGTGACGAAGTCGCCACCAAGCTCGACCTGGCCCGTGCCTATATCGACATGGGTGACAGCCAGGGCGCGCGCGACATCCTCGACGAAGTGGTCAAg                                                                                                                                  <  2:350054/149‑1 (MQ=255)
  gcAACCAGAGCCTGAGCCTTTGGACGACCTGGACTTCGACTTCTTCTCCGGCAGTGACGAAGTCGCCACCAAGCTCGACCTGGCCCGTGCCTATATCGACATGGGTGACAGCCAGGGCGCGCGCGACATCCTCGACGAAGTGGTCAAg                                                                                                                                  >  2:52250/1‑148 (MQ=255)
     aCCAGAGCCTGAGCCTTTGGACGACCTGGACTTCGACTTCTTCTCCGGCAGTGACGAAGTCGCCACCAAGCTCGACCTGGCCCGTGCCTATATCGACATGGGTGACAGCCAGGGCGCGCGCGACATCCTTGACCACTTCGTCGAGGATg                                                                                                                              <  4:176362/149‑1 (MQ=255)
     aCCAGAGCCTGAGCCTTTGGACGACCTGGACTTCGACTTCTTCTCCGGCAGTGACGAAGTCGCCACCAAGCTCGACCTGGCCCGTGCCTATATCGACATGGGTGACAGCCAGGGCGCGCGCGACATCCTCGACGAAGTGGTCAAGGATg                                                                                                                              <  2:55554/149‑1 (MQ=255)
          agcctgagcctTTGGACGACCTGGACTTCGACTTCTTCTCCGGCAGTGACGAAGTCGCCACCAAGCTCGACCTGGCCCGTGCCTATATCGACATGGGTGACAGCCAGGGCGCGCGCGACATCCTCGACGAAGTGGTCAAGGATGGcgac                                                                                                                         <  2:224535/149‑1 (MQ=255)
            cctgagcctTTGGACGACCTGGACTTCGACTTCTTCTCCGGCAGTGACGAAGTCGCCACCAAGCTCGACCTGGCCCGTGCCTATATCGACATGGGTGACAGCCAGGGCGCGCGCGACATCCTCGACGAAGTGGTCAAGGATGGcgacga                                                                                                                       >  4:354268/1‑149 (MQ=255)
            cctgagcctTTGGACGACCTGGACTTCGACTTCTTCTCCGGCAGTGACGAAGTCGCCACCAAGCTCGACCTGGCCCGTGCCTATATCGACATGGGTGACAGCCAGGGCGCGCGCGACATCCTCGACGAAGTGGTCAAGGATGGcgacga                                                                                                                       <  5:208693/149‑1 (MQ=255)
                agcctTTGGACGACCTGGACTTCGACTTCTTCTCCGGCAGTGACGAAGTCGCCACCAAGCTCGACCTGGCCCGTGCCTATATCGACATGGGTGACAGCCAGGGCGCGCGCGACATCCTCGACGAAGTGGTCAAGGATGGCGACGAcagc                                                                                                                   <  1:205037/149‑1 (MQ=255)
                agcctTTGGACGACCTGGACTTCGACTTCTTCTCCGGCAGTGACGAAGTCGCCACCAAGCTCGACCTGGCCCGTGCCTATATCGACATGGGTGACAGCCAGGGCGCGCGCGACATCCTCGACGAAGTGGTCAAGGATGGCGACGAcag                                                                                                                    <  2:308465/148‑1 (MQ=255)
                       ggACGACCTGGACTTCTACTTCATCTCCGGCAGTTACGAAGTCGCCACCAATCTCGACCTGTCCCGTTCATATATCGACATGGGTGACAGCCAGGGCGCGCGCGACATCCTCGACGAAGTGGTCAAGGATGGCGACTACAGCCAGCGt                                                                                                             <  2:271009/148‑1 (MQ=255)
                        gacgccCGGGACTTCGACTTCTTCTCCGGCAGTGACGAAGTCGCCACCAAGCTCGACCTGGCCCGTGCCTATATCGACATGGGTGACAGCCAGGGCGCGCGCGACATCCTCGACGAAGTGGTCAAGGATGGCGACGACAGCCAGCGTCa                                                                                                           <  2:135876/149‑1 (MQ=255)
                        gacgacCTGGACTTCGACTTCTTCTCCGGCAGTGACGAAGTCGCCACCAAGCTCGACCTGGCCCGTGCCTATATCGACATGGGTGACAGCCAGGGCGCGCGCGACATCCTCGACGAAGTGGTCAAGGATGTCGACGACAGCCAGCGTCa                                                                                                           <  2:131421/149‑1 (MQ=255)
                        gacgacCTGGACTTCGACTTCTTCTCCGGCAGTGACGAAGTCGCCACCAAGCTCGACCTGGCCCGTGCCTATATCGACATGGGTGACAGCCAGGGCGCGCGCGACATCCTCGACGAAGTGGTCAAGGATGGCGACGACAGCCAGCGTc                                                                                                            <  5:104018/148‑1 (MQ=255)
                          cgacCTGGACTTCGACTTCTTCTCCGGCAGTGACGAAGTCGCCACCAAGCTCGACCTGGCCCGTGCCTATATCGACATGGGTGACAGCCAGGGCGCGCGCGACATCCTCGACGAAGTGGTCAAGGATGGCGACGACAGCCAGCGTCAgg                                                                                                         <  4:143146/149‑1 (MQ=255)
                          cgacCTGGACTTCGACTTCTTCTCCGGCAGTGACGAAGTCGCCACCAAGCTCGACCTGGCCCGTGCCTATATCGACATGGGTGACAGCCAGGGCGCGCGCGACATCCTCGACGAAGTGGTCAAGGATGGCGACGACAGCCAGCGTCAgg                                                                                                         <  2:165262/149‑1 (MQ=255)
                          cgacCTGGACTTCGACTTCTTCTCCGGCAGTGACGAAGTCGCCACCAAGCTCGACCTGGCCCGTGCCTATATCGACATGGGTGACAGCCAGGGCGCGCGCGACATCCTCGACGAAGTGGTCAAGGATGGCGACGACAGCCAGCGTCAg                                                                                                          >  3:340913/1‑148 (MQ=255)
                                ggACTTCGACTTCTTCTCCGGCAGTGACGAAGTCGCCACCAAGCTCGACCTGGCCCGTGCCTATATCGACATGGGTGACAGCCAGGGCGCGCGCGACATCCTCGACGAAGTGGTCAAGGATGGCGACGACAGCCAGCGTCAGGAAGcc                                                                                                    >  2:195424/1‑148 (MQ=255)
                                 gacttcgacttcTTCTCCGGCAGTGACGAAGTCGCCACCAAGCTCGACCTGGCCCGTGCCTATATCGACATGGGTGACAGCCAGGGCGCGCGCGACATCCTCGACGAAGTGGTCAAGGATGGCGACGACAGCCAGCGTCAGGAAGCCGa                                                                                                  <  2:329124/149‑1 (MQ=255)
                                  acttcgacttcTTCTCCGGCAGTGACGAAGTCGCCACCAAGCTCGACCTGGCCCGTGCCTATATCGACATGGGTGACAGCCAGGGCGCGCGCGACATCCTCGACGAAGTGGTCAAGGATGGCGACGACAGCCAGCGTCAGGAAGCCGAg                                                                                                 <  1:52250/149‑1 (MQ=255)
                                   cttcgacttcTTCTCCGGCAGTGACGAAGTCGCCACCAAGCTCGACCTGGCCCGTGCCTATATCGACATGGGTGACAGCCAGGGCGCGCGCGACATCCTCGACGAAGTGGTCAAGGATGGCGACGACAGCCAGCGTCAGGAAGCCGAgg                                                                                                <  6:197574/149‑1 (MQ=255)
                                   cttcgacttcTTCTCCGGCAGTGACGAAGTCGCCACCAAGCTCGACCTGGCCCGTGCCTATATCGACATGGGTGACAGCCAGGGCGCGCGCGACATCCTCGACGAAGTGGTCAAGGATGGCGACGACAGCCAGCGTCAGGAAGCCGAgg                                                                                                <  4:344982/149‑1 (MQ=255)
                                    ttcgacttcTTCTCCGGCAGTGACGAAGTCGCCACCAAGCTCGACCTGGCCCGTGCCTATATCGACATGGGTGACAGCCAGGGCGCGCGCGACATCCTCGACGAAGTGGTCAAGGATGGCGACGACAGCCAGCGTCAGGAAGCCGAGGa                                                                                               <  2:378731/149‑1 (MQ=255)
                                    ggcgacttcGTCACTGCCGGAGAAGAAGTCGCCACCAAGCTCGACCTGGCCCGTGCCTATATCGACATGGGTGACAGCCAGGGCGCGCGCGACATCCTCGACGAAGTGGTCAAGGATGGCGACGACAGCCAGCGTCAGGAAGCCGAGGa                                                                                               >  1:192411/3‑149 (MQ=255)
                                      cgacttcTTCTCCGGCAGTGACGAAGTCGCCACCAAGCTCGACCTGGCCCGTGCCTATATCGACATGGGTGACAGCCAGGGCGCGCGCGACATCCTCGACGAAGTGGTCAAGGATGGCGACGACAGCCAGCGTCAGGAAGCCGAGGACa                                                                                             >  3:235888/1‑149 (MQ=255)
                                         cttcttCTCCGGCAGTGACGAAGTCGCCACCAAGCTCGACCTGGCCCGTGCCTATATCGACATGGGTGACAGCCAGGGCGCGCGCGACATCCTCGACGAAGTGGTCAAGGATGGCGACGACAGCCAGCGTCAGGAAGCCGAGGACATGc                                                                                          <  1:195424/149‑1 (MQ=255)
                                               ctcCGGCAGTGACGAAGTCGCCACCAAGCTCGACCTGGCCCGTGCCTATATCGACATGGGTGACAGCCAGGGCGCGCGCGACATCCTCGACGAAGAGGTCAAGGAAGGCGACGACAGCCAGCGTCAGGAAGCCGAGGACATGCTCTccc                                                                                    >  8:169558/1‑149 (MQ=255)
                                                tcCGGCAGTGACGAAGTCGCCACCAAGCTCGACCTGGCCCGTGCCTATATCGACATGGGTGACAGCCAGGGCGCGCGCGACATCCTCGACGAAGTGGTCAAGGATGGCGACGACAGCCAGCGTCAGGAAGCCGAGGACATGCTCTCCCg                                                                                   >  7:176410/1‑149 (MQ=255)
                                                 ccGGCAGTGACGAAGTCGCCACCAAGCTCGACCTGGCCCGTGCCTATATCGACATGGGTGACAGCCAGGGCGCGCGCGACATCCTCGACGAAGTGGTCAAGGATGGCGACGACAGCCAGCGTCAGGAAGCCGAGGACATGCTCTCCCgg                                                                                  >  5:340629/1‑149 (MQ=255)
                                                  cGGCAGTGACGAAGTCGCCACCAAGCTCGACCTGGCCCGTGCCTATATCGACATGGGTGACAGCCAGGGCGCGCGCGACATCCTCGACGAAGTGTTCAAGGATGGCGACGACAGCCAGCGTCAGGAAGCCGAGGACATGCTCTCCCGGc                                                                                 <  6:264714/149‑1 (MQ=255)
                                                  cGGCAGTGACGAAGTCGCCACCAAGCTCGACCTGGCCCGTGCCTATATCGACATGGGTGACAGCCAGGGCGCGCGCGACATCCTCGACGAAGTGGTCAAGGATGGCGACGACAGCCAGCGTCAGGAAGCCGAGGACATGCTCTCCCGGc                                                                                 >  7:27971/1‑149 (MQ=255)
                                                   ggCAGTGACGAAGTCGCCACCAAGCTCGACCTGGCCCGTGCCTATATCGACATGGGTGACAGCCAGGGCGCGCGCGACATCCTCGACGAAGTGGTCAAGGATGGCGACGACAGCCAGCGTCAGGAAGCCGAGGACATGCTCTCCCGGCt                                                                                <  3:65977/149‑1 (MQ=255)
                                                      aGTGACGAAGTCGCCACCAAGCTCGACCTGGCCCGTGCCTATATCGACATGGGTGACAGCCAGGGCGCGCGCGACATCCTCGACGAAGTGGTCAAGGATGGCGACGACAGCCAGCGTCAGGAAGcc                                                                                                    >  1:19097/1‑126 (MQ=255)
                                                      aGTGACGAAGTCGCCACCAAGCTCGACCTGGCCCGTGCCTATATCGACATGGGTGACAGCCAGGGCGCGCGCGACATCCTCGACGAAGTGGTCAAGGATGGCGACGACAGCCAGCGTCAGGAAGcc                                                                                                    <  2:19097/126‑1 (MQ=255)
                                                           cGAAGTCGCCACCAAGCTCGACCTGGCCCGTGCCTATATCGACATGGGTGACAGCCAGGGCGCGCGCGACATCCTCGACGAAGTGGTCAAGGATGGCGACGACAGCCAGCGTCAGGAAGCCGAGGACATGCTCTCCCGGCTGGTCTGAt                                                                        >  6:167410/1‑149 (MQ=255)
                                                           agAAGTCGCCACCAAGCTCGACCTGGCCCGTGCCTATATCGACATGGGTGACAGCCAGGGCGCGCGCGACATCCTCGACGAAGTGGTCAAGGATGGCGACGACAGCCAGCGTCAGGAAGCCGAGGACATGCTCTCCCGGCTGGTCTGAt                                                                        <  2:192411/148‑1 (MQ=255)
                                                                   ccaccaAGCTCGACCTGGCCCGTGCCTATATCGACATGGGTGACAGCCAGGGCGCGCGCGACATCCTCGACGAAGTGGTCAAGGATGGCGACGACAGCCAGCGTCAGGAAGCCGAGGACATGCTCTCCCGGCTGGTCTGATGCAAGCAg                                                                >  6:209720/1‑149 (MQ=255)
                                                                            tCGACCTGGCCCGTGCCTATATCGACATGGGTGACAGCCAGGGCGCGCGCGACATCCTCGACGAAGTGGTCAAGGATGGCGACGACAGCCAGCGTCAGGAAGCCGAGGACATGCTCTCCCGGCTGGTCTGATGCAAGCAGGTCAGGCct                                                       >  7:155251/1‑149 (MQ=255)
                                                                              gACCTGGCCCGTGCCTATATCGACATGGGTGACAGCCAGGGCGCGCGCGACATCCTCGACGAAGTGGTCAAGGATGGCGACGACAGCCAGCGTCAGGAAGCCGAGGACATGCTCTCCCGGCTGGTCTGATGCAAGCAGGTCAGGCctc                                                      <  8:176410/148‑1 (MQ=255)
                                                                                ccTGGCCCGTGCCTATATCGACATGGGTGACAGCCAGGGCGCGCGCGACATCCTCGACGAAGTGGTCAAGGATGGCGACGACAGCCAGCGTCAGGAAGCCGAGGACATGCTCTCCCGGGTGGTCTGATGCAAGCAGGTCAGGCCTCTTg                                                   <  1:34990/149‑1 (MQ=255)
                                                                                 cTGGCCCGTGCCTATATCGACATGGGTGACAGCCAGGGCGCGCGCGACATCCTCGACGAAGTGGTCAAGGATGGCGACGACAGCCAGCGTCAGGAAGCCGAGGACATGCTCTCCCGGCTGGTCTGATGCAAGCAGGTCAGGCCTCTTgg                                                  >  1:285801/1‑149 (MQ=255)
                                                                                      ccGTGCCTATATCGACATGGGTGACAGCCAGGGCGCGCGCGACATCCTCGACGAAGTGGTCAAGGATGGCGACGACAGCCAGCGTCAGGAAGCCGAGGACATGCTCTCCCGGCTGGTCTGATGCAAGCAGGTCAGGCCTCTTGGCGGGc                                             >  4:102714/1‑149 (MQ=255)
                                                                                         tGCCTATATCGACATGGGTGACAGCCAGGGCGCGCGCGACATCCTCGACGAAGTGGTCAAGGATGGCGACGACAGCCAGCGTCAGGAAGCCGAGGACATGCTCTCCCGGCTGGTCTGATGCAAGCAGGTCAGGCCTCTTGGCGGGCTGa                                          <  5:163641/149‑1 (MQ=255)
                                                                                            cTATATCGACATGGGTGACAGCCAGGGCGCGCGCGACATCCTCGACGAAGTGGTCAAGGATGGCGACGACAGCCAGCGTCAGGAAGCCGAGGACATGCTCTCCCGGCTGGTCTGATGCAAGCAGGTCAGGCCTCTTGGCGGGCTGACcg                                       >  3:97293/1‑149 (MQ=255)
                                                                                                 tCGACATGGGTGACAGCCAGGGCGCGCGCGACATCCTCGACGAAGTGGTCAAGGATGGCGACGACAGCCAGCGTCAGGAAGCCGAGGACATGCTCTCCCGGCTGGTCTGATGCAAGCAGGTCAGGCCTCTTGGCGGGCTGACCGCGcc                                   <  8:27971/148‑1 (MQ=255)
                                                                                                   gACATGGGTGACAGCCAGGGCGCGCGCGACATCCTCGACGAAGTGGTCAAGGATGGCGACGACAGCCAGCGTCAGGAAGCCGAGGACATGCTCTCCCGGCTGGTCTGATGCaa                                                                    >  3:203419/1‑113 (MQ=255)
                                                                                                   gACATGGGTGACAGCCAGGGCGCGCGCGACATCCTCGACGAAGTGGTCAAGGATGGCGACGACAGCCAGCGTCAGGAAGCCGAGGACATGCTCTCCCGGCTGGTCTGATGCaa                                                                    <  4:203419/113‑1 (MQ=255)
                                                                                                         ggTGACAGCCAGGGCGCGCGCGACATCCTCGACGAAGTGGTCAAGGATGGCGACGACAGCCAGCGTCAGGAAGCCGAGGACATGCTCTCCCGGCTGGTCTGATGCAAGCAGGTCAGGCCTCTTGGCGGGCTGACCGCGCCCGCGAagag                          <  6:340629/149‑1 (MQ=255)
                                                                                                                  cAGGGCGCGCGCGACATCCTCGACGAAGTGGTCAAGGATGGCGACGACAGCCAGCGTCAGGAAGCCGAGGACATGCTCTCCCGGCTGGTCTGATGCAAGCAGGTCAGGCCTCTTGGCGGGCTGACCGCGCCCGCGAAGAGGCCCGCGCa                 >  1:106108/1‑149 (MQ=255)
                                                                                                                    gggCGCGCGCGACATCCTCGACGAAGTGGTCAAGGATGGCGACGACAGCCAGCGTCAGGAAGCCGAGGACATGCTCTCCCGGCTGGTCTGATGCAAGCAGGTCAGGCCTCTTGGCGGGCTGACCGCGCCCGCGAAGAGGCCCGCGCATg               <  4:340913/149‑1 (MQ=255)
                                                                                                                               aCATCCTCGACGAAGTGGTCAAGGATGGCGACGACAGCCAGCGTCAGGAAGCCGAGGACATGCTCTCCCGGCTGGTCTGATGCAAGCAGGTCAGGCCTCTTGGCGGGCTGACCGCGCCCGCGAAGAGGCCCGCGCATGCAGTGCGCAAt    <  3:102714/149‑1 (MQ=255)
                                                                                                                                cATCCTTGACCACTTCGTCGAGGATGGCGACGACAGCCAGCGTCAGGAAGCCGAGGACATGCTCTCCCGTCTGGTATGATGCAAGCAGGTCCGGCCTCTTGGCGGGATGACCGCGACCTCGAAGAGGCCCGCCCATGCAGTGCGCAAt    >  7:176148/1‑148 (MQ=255)
                                                                                                                                cATCCTTGACCACTTCGTCGAGGATGGCGACGACAGCCAGCGTCAGGAAGCCGAGGACATGCTCTCCCGGCTGGTCTGATGCAAGCAGGTCAGGCCTCTTGGCGGGCTGACCGCGCCCGCGAAGAGGCCCGCGCATGCAGTgcgc       >  3:222813/1‑145 (MQ=255)
                                                                                                                                cATCCTTGACCACTTCGTCGAGGATGGCGACGACAGCCAGCGTCAGGAAGCCGAGGACATGCTCTCCCGGCTGGTCTGATGCAAGCAGGTCAGGCCTCTTGGCGGGCTGACCGCGCCCGCGAAGAGGCCCGCGCATGCAGTgcgc       <  4:222813/145‑1 (MQ=255)
                                                                                                                                 atacTCGACGAAGTGGTCAAGGATGGCGACGACATCCAGCGTCAGGAAGCCTAGGACATGCGCTCCCGGCTGGTCTGATGCAAGCAGGTCAGGCCTCTTGGCGGGCTGACCGCGCCCGCGAAGAGGCCCGCGCATGCAGTGCGCAATCa  <  4:11138/146‑1 (MQ=255)
                                                                                                                                 aTCCTCGACGAAGTGGTCAAGGATGGCGACGACAGCCAGCGTCAGGAAGCCGAGGACATGCTCTCCCGGCTGGTCTGATGCAAGCAGGTCAGGCCTCTTGGCGGGCTGACCGCGCCCGCGAAGAGGCCCGCGCATGCAGTGCGCAATCa  <  3:354268/149‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                               
GCGCAACCAGAGCCTGAGCCTTTGGACGACCTGGACTTCGACTTCTTCTCCGGCAGTGACGAAGTCGCCACCAAGCTCGACCTGGCCCGTGCCTATATCGACATGGGTGACAGCCAGGGCGCGCGCGACATCCTCGACGAAGTGGTCAAGGATGGCGACGACAGCCAGCGTCAGGAAGCCGAGGACATGCTCTCCCGGCTGGTCTGATGCAAGCAGGTCAGGCCTCTTGGCGGGCTGACCGCGCCCGCGAAGAGGCCCGCGCATGCAGTGCGCAATCA  >  NC_002947/2261942‑2262219

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: ‑