Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,132,762 T→C 8.9% N84D (AAT→GAT)  PP_1890 ← type I pili chaperone protein FimC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,132,7620TC 8.9% 124.0 / 2.7 45N84D (AAT→GAT) PP_1890type I pili chaperone protein FimC
Reads supporting (aligned to +/- strand):  ref base T (22/19);  new base C (2/2);  total (24/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GCGGAATCTCCTGCACGTTGAAGAAGAACAGCGACTCGCGGTCGTGTGGCAGTGTATTGGGCAGTCCAATTATTTTTAGAAGTTGTTCACTGTTCGGGTTGATCTTGAACAGCGGTGGGGACGTTGCAAAGGGTACGCTGGTCGTATTGTCGTCGGCGTCCGTATTTATCCACGCCTGCACGGCATAGATGCTTTTCGATGGGTTTCGAACAACGACGGATGTACTGCGTTTGTCGCCGTTGAAAACGAGGCGCGTGTTGTTTAATGTCAAGGCGGCATGACTGTTTAC  >  NC_002947/2132615‑2132903
                                                                                                                                                   |                                                                                                                                             
gCGGAATCTCCTGCACGTTGAAGAAGAACAGCGACTCGCGGTCGTGTGGCAGTGTATTGGGCAGTCCAATTATTTTTAGAAGTTGTTCACTGTTCGGGTTGATCTTGAACAGCGGTGGGGACGTTGCAAAGGGTACGCTAGTCGTAtt                                                                                                                                               <  8:31252/148‑1 (MQ=255)
          cTGCACGTTGAAGAAGAACAGCGACTCGCGGTCGTGTGGCAGTGTATTGGGCAGTCCAATTATTTTTAGAAGTTGTTCACTGTTCGGGTTGATCTTGAACAGCGGTGGGGACGTTGCAAAGGGTACGCTGGTCGTATTGTCGTCGGCg                                                                                                                                     <  2:191442/148‑1 (MQ=255)
          cTGCACGTTGAAGAAGAACAGCGACTCGCGGTCGTGTGGCAGTGAATTGGGCAGTCCAATTATTTTTAGAAGTTGTTCACTGTTCGGGTTGATCTTGAACAGCGGTGGGGCCGTTGCAAAGGGTACCCAGGTCGTATTGACGTCGGCGt                                                                                                                                    >  8:5891/1‑149 (MQ=255)
            gCACGTTGAAGAAGAACAGCGACTCGCGGTCGTGTGGCAGTGTATTGGGCAGTCCAATTATTTTTAGAAGTTGTTCACTGTTCGGGTTGATCTTGAACAGCGGTGGGGACGTTGCAAAGGGTACGCTGGTCGTATTGTCGTCGGCGTcc                                                                                                                                  >  3:97880/1‑149 (MQ=255)
             cACGTTGAAGAAGAACAGCGACTCGCGGTCGTGTGGCAGTGTATTGGGCAGTCCAATTATTTTTAGAAGTTGTTCACTGTTCGGGTTGATCTTGAACAGCGGTGGGGACGTTGCAAAGGGTACGCTGGTCGTATTGTCGTCGGCGTCCg                                                                                                                                 >  3:34300/1‑149 (MQ=255)
                   gaagaagaaCAGCGACTCGCGGTCGTGTGGCAGTGTATTGGGCAGTCCAATTATTTTTAGAAGTTGTTCACTGTTCGGGTTGATCTTGAACAGCGGTGGGGACGTTGCAAAGGGTACGCTGGTCGTATTGTCGTCGGCGTCCGTATTTa                                                                                                                           >  8:226720/1‑149 (MQ=255)
                      gaagaaCAGCGACTCGCGGTCGTGTGGCAGTGTATTGGGCAGTCCAATTATTTTTAGAAGTTGTTCACTGTTCGGGTTGATCTTGAACAGCGGTGGGGACGTTGCAAAGGGTACGCTGTTCGTATCGTCGTCGGCGTCCGTATTTATc                                                                                                                         >  3:55576/1‑148 (MQ=255)
                      gaagaaCAGCGACTCGCGGTCGTGTGGCAGTGTATTGGGCAGTCCAATTATTTTTAGAAGTTGTTCACTGTTCGGGTTGATCTTGAACAGCGGTGGGGACGTTGCAAAGGGTACGCTGGTCGTATCGTCGTCGGCGTCCGTATTTATc                                                                                                                         >  3:50654/1‑148 (MQ=255)
                        agaaCAGCGACTCGCGGTCGTGTGGCAGTGTATTGGGCAGTCCAATTATTTTTAGAAGTTGTTCACTGTTCGGGTTGATCTTGAACAGCGGTGGGGACGTTGCAAAGGGTACGCTGGTCGTATTGTCGTCGGCGTCCGTATTTATCCAc                                                                                                                      <  4:97880/149‑1 (MQ=255)
                           aCAGCGACTCGCGGTCGTGTGGCAGTGTATTGGGCAGTCCAATTATTTTTAGAAGTTGTTCACTGTTCGGGTTGATCTTGAACAGCGGTGGGGACGTTGCAAAGGGTACGCTGGTCGTATTGTCGTCGGCGTCCGTATTTATCCACGcc                                                                                                                   <  2:13344/149‑1 (MQ=255)
                                     gcgGTCGTGTGGCAGTGTATTGGGCAGTCCAATTATTTTTAGAAGTTGTTCACTGTTCGGGTTGATCTTGAACAGCGGTGGGGACGTTGCAAAGGGTACGCTGGTCGTATTGTCGTCGGCGTCCGTATTTATCCACGCCTGCACGGCa                                                                                                          >  4:13361/1‑148 (MQ=255)
                                     gcgGTCGTGTGGCAGTGTATTGGGCAGTCCAATTATTTTTAGAAGTTGTTCACTGTTCGGGTTGATCTTGAACAGCGGTGGGGACGTTGCAAAGGGTACGCTGGTCGTATTGTCGTCGGCGTCCGTATTTATCCACGCCTGCACGGCa                                                                                                          <  7:30236/148‑1 (MQ=255)
                                        gTCGTGTGGCAGTGTATTGGGCAGTCCAATTATTTTTAGAAGTTGTTCACTGTTCGGGTTGATCTTGAACAGCGGTGGGGACGTTGCAAAGGGTACGCTGGTCGTATTGTCGTCGGCGTCCGTATTTATCCACGCCTGCACGGCa                                                                                                          >  4:253028/1‑145 (MQ=255)
                                           gtttgGCAGTGTAGTGGGCAGTCCATTTATTTATAGAATTTGTTCACTGTTCGGGTTGATCTTGAACAGGGGTGGGGACGTTGCAAAGGGTACGCTGGTCGTATTTTCGTCGGCGTCCGTATTTATCCACGCCTGCACGGCATAGATGc                                                                                                   <  6:29337/146‑1 (MQ=255)
                                                  aGTGTATTGGGCAGTCCAATTATTTTTAGAAGTTGTTCACTGTTCGGGTTGATCTTGAACAGCGGTGGGGACGTTGCAAAGGGTACGCTGGTCGTATTGTCGTCGGCGTCCGTATTTATCCACGCCTGCACGGCATAGATGCTTTTCGa                                                                                            <  7:331046/149‑1 (MQ=255)
                                                      tATTGGGCAGTCCAATTATTTTTAGAAGTTGTTCACTGTTCGGGTTGATCTTGAACAGCGGTGGGGACGTTGCAAAGGGTACGCTGGTCGTATTGTCGTCGGCGTCCGTATTTATCCACGCCTGCACGGCATAGATGCTTTTCGATggg                                                                                        >  6:317539/1‑149 (MQ=255)
                                                                tCCAATTATTTTTAGAAGTTGTTCACTGTTCGGGTTGATCTTGAACAGCGGTGGGGACGTTGCAAAGGGTACGCTGGTCGTATTGTCGTCGGCGTCCGTAt                                                                                                                              >  6:11696/1‑101 (MQ=255)
                                                                tCCAATTATTTTTAGAAGTTGTTCACTGTTCGGGTTGATCTTGAACAGCGGTGGGGACGTTGCAAAGGGTACGCTGGTCGTATTGTCGTCGGCGTCCGTAt                                                                                                                              <  5:11696/101‑1 (MQ=255)
                                                                    attattTTTAGAAGTTGTTCACTGTTCGGGTTGATCTTGAACAGCGGTGGGGACGTTGCAAAGGGTACGCTGGTCGTATTGTCGTCGGCGTCCGTATTTATCCACGCCTGCACGGCATAGATGCTTTTCGATGGGTTTCGAACAacgac                                                                          >  1:1433/1‑149 (MQ=255)
                                                                        tttttAGAAGTTGTTCACTGTTCGGGTTGATCTTGAACAGCGGTGGGGACGTTGCAAAGGGTACGCTGGTCGTATTGTCGTCGGCGTCCGTATTTATCCACGCCTGCACGGCATAGATGCTTTTCGATGGGTTTCGAACAACGACGGa                                                                       <  3:13361/148‑1 (MQ=255)
                                                                           ttAGAAGTTGTTCACTGTTCGGGTTGATCTTGAACAGCGGTGGGGACGTTGCAAAGGGTACGCTGGTCGTATTGTCGTCGGCGTCCGTATTTATCCACGCCTGCACGGCATAGATGCTTTTCGATGGGTTTCGAACAACGACGGATga                                                                    >  1:135440/1‑147 (MQ=255)
                                                                             aGAAGTTGTTCACTGTTCGGGTTGATCTTGAACAGCGGTGGGGACGTTGCAAAGGGTACGCTGGTCGTATTGTCGTCGGCGTCCGTATTTATCCACGCCTGCACGGCATAGATGCTTTTCGATGGGTTTCGaa                                                                                 <  8:3469/133‑1 (MQ=255)
                                                                             aGAAGTTGTTCACTGTTCGGGTTGATCTTGAACAGCGGTGGGGACGTTGCAAAGGGTACGCTGGTCGTATTGTCGTCGGCGTCCGTATTTATCCACGCCTGCACGGCATAGATGCTTTTCGATGGGTTTCGaa                                                                                 >  7:3469/1‑133 (MQ=255)
                                                                                 gttgttCACTGTTCGGGTTGATCTTGAACAGCGGTGGGGACGTTGCAAAGGGTACGCTGGTCGTATTGTCGTCGGCGTCCGTATTTATCCACGCCTGCACGGCATAGATGCTTTTCGATGGGTTTCGAACAACGACGGATGTACTGCGt                                                             <  7:226720/149‑1 (MQ=255)
                                                                                            ttCGGGTTGATCTTGAACAGCGGTGGGGACGTTGCAAAGGGTACGCTGGTCGTATTGTCGTCGGCGTCCGTATTTATCCACGCCTGCACGGCATAGATGCTTTTCGATGGGTTTCGAACAACGACGGATGTACTGCGTTTGTCGCCGtt                                                  <  2:265031/149‑1 (MQ=255)
                                                                                            ttCGGGTTGATCTTGAACAGCGGTGGGGACGTTGCAAAGGGTACGCTGGTCGTATTGTCGTCGGCGTCCGTATTTATCAACGCCTGCACGGCATAGATGCTTTTAGATGGGTTTCGAACAACGACCGATGTACTGCGTTTGTCGCCGtt                                                  >  1:376388/1‑149 (MQ=255)
                                                                                                 gTTGATCTTGAACAGCGGTGGGGACGTTGCAAAGGGTACGCTGGTCGTATTGTCGTCGGCGTCCGTATTTATCCACGCCTGCCCGGCATAGATGCTTTTCGATGGGTTTCGAAAAACGAGGGATGTCCTGCGTTTGTCGCCGTTGaaac                                             >  6:292984/1‑148 (MQ=255)
                                                                                                 gTTGATCTTGAACAGCGGTGGGGACGTTGCAAAGGGTACGCTGGTCGTATTGTCGTCGGCGTCCGTATTTATCCACGCCTGCACGGCATAGATGCTTTTCGATGGGTTTCGAACAACGACGGATGTACTGCGTTTGTCGCCGTTGaaaa                                             >  6:32578/1‑149 (MQ=255)
                                                                                                  ttGATCTTGAACAGCGGTGGGGACGTTGCAAAGGGTACGCTGGTCGTATTGTCGTCGGCGTCCGTATTTATCCACGCCTGCACGGCATAGATGCTTTTCGATGGGTTTCGaa                                                                                 >  8:286432/1‑112 (MQ=255)
                                                                                                  ttGATCTTGAACAGCGGTGGGGACGTTGCAAAGGGTACGCTGGTCGTATTGTCGTCGGCGTCCGTATTTATCCACGCCTGCACGGCATAGATGCTTTTCGATGGGTTTCGaa                                                                                 <  7:286432/112‑1 (MQ=255)
                                                                                                       cTTGAACAGCGGTGGGGACGTTGCAAAGGGTACGCTGGTCGTATTGTCGTCGGCGTCCGTATTTATCCACGCCTGCACGGCATAGATGCTTTTCGATGGGTTTCGaacaa                                                                              >  5:193474/1‑110 (MQ=255)
                                                                                                       cTTGAACAGCGGTGGGGACGTTGCAAAGGGTACGCTGGTCGTATTGTCGTCGGCGTCCGTATTTATCCACGCCTGCACGGCATAGATGCTTTTCGATGGGTTTCGaacaa                                                                              <  6:193474/110‑1 (MQ=255)
                                                                                                        ttGAACAGCGGTGGGGACGTTGCAAAGGGTACGCTGGTCGTATTGTCGTCGGCGTCCGTATTTATCCACGCCTGCACGGCATAGATGCTTTTCGATGGGTTTCGAACAACGACGGATGTACTGCGTTTGTCGCCGTTGAAAACGAggag                                      >  3:246791/1‑147 (MQ=255)
                                                                                                        ttGAACAGCGGTGGGGACGTTGCAAAGGGTACGCTGGTCGTATTGTCGTCGGCGTCCGTATTTATCCACGCCTGCACGGCATAGATGCTTTTCGATGGGTTTCGAACAACGACGGATGTACTGCGTTTGTCGCCGTTGAAAACGAGgcg                                      >  3:6770/1‑149 (MQ=255)
                                                                                                                     gggACGTTGCAAAGGGTACGCTGGTCGTATTGTCGTCGGCGTCCGTATTTATCCACGCCTGCACGGCATAGATGCTTTTCGATGGGTTTCGAACAACGACGGATGTACTGCGTTTGTCGCCGTTGAAAACGAGGCGCGTGTTGTTTAAt                         <  2:1433/149‑1 (MQ=255)
                                                                                                                      ggACGTTGCAAAGGGTACGCTGGTCGTATTGTCGTCGGCGTCCGTATTTATCCACGCCTGCACGGCATAGATGCTTTTCGATGGGTTTCGAACAACGACGGATGTACTGCGTTTGTCGCCGTTGAAAACGAGGCGCGTGTTGTTTAATg                        <  7:5891/149‑1 (MQ=255)
                                                                                                                       gACGTTGCAAAGGGTACGCTGGTCGTATCGTCGTCGGCGTCCGTATTTATCCACGCCTGCACGGCATAGATGCTTTTCGATGGGTTTCGAACAACGACGGATGTACTGCGTTTGTCGCCGTTGAAAACGAGGCGCGTGTTGTTTAATGt                       <  4:55576/149‑1 (MQ=255)
                                                                                                                       gACGTTGCAAAGGGTACGCTGGGCGTATCGTCGTCGGCGTCCGTATTTATCCACGCCTGCACGGCATAGATGCTTTTCGATGGGCTTCGAACAACGACGGATGTACTGCGTTTGTCGCCGTTGAAAACGAGGCGCGTGTTGTTTAATGt                       <  4:50654/149‑1 (MQ=255)
                                                                                                                           ttGCAAAGGGTACGCTGGTCGTATTGTCGTCGGCGTCCGTATTTATCCACGCCTGCACGGCATAGATGCTTTTCGATGGGTTTCGAACAACGACGGATGTACTGCGTTTGTCGCCGTTGAAAACGAGGCGCGTGTTGTTTAATGTCAAg                   >  8:115790/1‑149 (MQ=255)
                                                                                                                             gCAAAGGGTACGCTGGTCGTATTGTCGTCGGCGTCCGTATTTATCCACGCCTGCACGGCATAGATGCTTTTCGATGGGTTTCGAACAACGACGGATGTACTGCGTTTGTCGCCGTTGAAAACGAGGCGCGTGTTGTTTAATGTCAAggc                 >  8:159745/1‑149 (MQ=255)
                                                                                                                             gCAAAGGGTACGCTGGTCGTATTGTCGTCGGCGTCCGTATTTATCCACGCCTGCACGGCATAGATGCTTTTCGATGGGTTTCGAACAACGACGGATGTACTGCGTTTGTCGCCGTTGAAAACGAGGCGCGTGTTGTATAATGTCAAggc                 >  4:179675/1‑149 (MQ=255)
                                                                                                                               aaaGGGTACGCTGGTCGTATTGTCGTCGGCGTCCGTATTTATCCACGCCTGCACGGCATAGATGCTTTTCGATGGGTTTCGAACAACGACGGATGTACTGCGTTTGTCGCCGTTGAAAACGAGGCGCGTGTTGTTTAATGTCAAggcgg               <  5:317539/149‑1 (MQ=255)
                                                                                                                                aaGGGTACGCTGGTCGTATTGTCGTCGGCGTCCGTATTTATCCACGCCTGCACGGCATAGATGCTTTTCGATGGGTTTCGAACAACGACGGATGTACTGCGTTTGTCGCCGTTGAAAACGAGGCGCGTGTTGTTTAATGTCAAggcggc              <  3:324204/149‑1 (MQ=255)
                                                                                                                                  gggTACGCTGGTCGTATTGTCGTCGGCGTCCGTATTTATCCACGCCTGCACGGCATAGATGCTTTTCGATGGGTTTCGAACAACGACGGATGTACTGCGTTTGTCGCCGTTGAAAACGAGGCGCGTGTTGTTTAATGTCAAGGCGGCAt            <  7:159745/149‑1 (MQ=255)
                                                                                                                                   ggTACGCTGGTCGTATTGTCGTCGGCGTCCGTATTTATCCACGCCTGCACGGCATAGATGCTTTTCGATGGGTTTCGAACAACGACGGATGTACTGCGTTTGTCGCCGTTGAAAACGAGGCGCGTGTTGTTTAATGTCAAGGCGGCAt            >  2:65176/1‑148 (MQ=255)
                                                                                                                                            gTCGTATTGTCGTCGGCGTCCGTATTTATCCACGCCTGCACGGCATAGATGCTTTTCGATGGGTTTCGAACAACGACGGATGTACTGCGTTTGTCGCCGTTGAAAACGAGGCGCGTGTTGTTTAATGTCAAGGCGGCATGACTGTTTAc  <  1:65176/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
GCGGAATCTCCTGCACGTTGAAGAAGAACAGCGACTCGCGGTCGTGTGGCAGTGTATTGGGCAGTCCAATTATTTTTAGAAGTTGTTCACTGTTCGGGTTGATCTTGAACAGCGGTGGGGACGTTGCAAAGGGTACGCTGGTCGTATTGTCGTCGGCGTCCGTATTTATCCACGCCTGCACGGCATAGATGCTTTTCGATGGGTTTCGAACAACGACGGATGTACTGCGTTTGTCGCCGTTGAAAACGAGGCGCGTGTTGTTTAATGTCAAGGCGGCATGACTGTTTAC  >  NC_002947/2132615‑2132903

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