Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,132,962 G→A 9.1% T17M (ACG→ATG)  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,9620GA 9.1% 118.8 / 3.3 44T17M (ACG→ATG) PP_1890type I pili chaperone protein FimC
Reads supporting (aligned to +/- strand):  ref base G (22/18);  new base A (2/2);  total (24/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.35e-01

GGTTTCGAACAACGACGGATGTACTGCGTTTGTCGCCGTTGAAAACGAGGCGCGTGTTGTTTAATGTCAAGGCGGCATGACTGTTTACCGAGAGGGTAAGCATCAATAAGGTCAGGCAGACGCAAACAACGACGCGGGCGGGTATCGTGAATGTTGACGGATTGGTCATCGGGTAACTCGCTGCATCCAGACGCACAACTACGCCTGGATGGTTTAACCGTGCAAGGACTGCACAGAATTTATGCGTAATTATTGATAGTCGACGGTGTAACGAACCTGGGCGGATGCAAGG  >  NC_002947/2132816‑2133107
                                                                                                                                                  |                                                                                                                                                 
ggTTTCGAACAACGACGGATGTACTGCGTTTGTCGCCGTTGAAAACGAGGCGCGTGTTGTTTAATGTCAAGGCGGCATGACTGTTTACCGAGAGGGTAAGCATCAATAAGGTCAGGCAGACGCAAACAACGACGCGGGCGGGTATCGt                                                                                                                                                  <  3:194845/148‑1 (MQ=255)
 gTTTCGAACAACGAGGGATGTACTGCGTTTGTCGCCGTTGAAAACGAGGCGCGTGTTGTTTAATGTCAAGGCGGCATGACTGTTTACCGAGAGGGTAAGCATCAATAAGGTCAGGCAGACGCAAACAACGACGCGGGCGGGTATCGTg                                                                                                                                                 >  7:347839/1‑148 (MQ=255)
  tttCGAACAACGACGGATGTACTGCGTTTGTCGCCGTTGAAAACGAGGCGCGTGTTGTTTAATGTCAAGGCGGCATGACTGTTTACCGAGAGGGTAAGCAACAATAAGGTCAGGCAGACGCAAACAACGACGCGGGCGGGTATCATGaa                                                                                                                                               >  8:379586/1‑149 (MQ=255)
             gacgGATGTACTGCGTTTGTCGCCGTTGAAAACGAGGCGCGTGTTGTTTAATGTCAAGGCGGCATGACTGTTTACCGAGAGGGTAAGCATCAATAAGGTCAGGCAGACGCAAACAACGACGCGGGCGGGTATCGTGAATGTTGACGGAt                                                                                                                                    <  2:155701/149‑1 (MQ=255)
                       cTGCGTTTGTCGCCGTTGAAAACGAGGCGCGTGTTGTTTAATGTCAAGGCGGCATGACTGTTTACCGAGAGGGTAAGCATCAATAAGGTCAGGCAGACGCAAACAACGACGCGGGCGGGTATCGTGAATGTTGACGGATTGGTCATCgg                                                                                                                          >  5:248439/1‑149 (MQ=255)
                       cTGCGTTTGTAGCCGTTGAAAACGAGGCGCGTGTTGTTTAATGTCAAGGCGGCATGACTGTTTACCGAGAGGGTAAGCATCAATAAGGTCAGGCAGACGCAAACAACGACGCCGGACGGTATCGTGAATGTTGACGGATAGGTCATCgg                                                                                                                          >  8:229404/1‑149 (MQ=255)
                        tGCGTTTGTCGCCGTTGAAAACGAGGCGCGTGTTGTTTAATGTCAAGGCGGCATGACTGTTTACCGAGAGGGTAAGCATCAATAAGGTCAGGCAGACGCAAACAACGACGCGGGCGGGTATAGTGAATGTTGACGGATTgg                                                                                                                                 >  5:340379/1‑141 (MQ=255)
                        tGCGTTTGTCGCCGTTGAAAACGAGGCGCGTGTTGTTTAATGTCAAGGCGGAATGACTGTTTACCGAGAGGGTAAGCATCAATAAGGTCAGGCAGACGCAAACAACGACGCGGGCGGGTATCGTGAATGTTGACGGATTgg                                                                                                                                 <  6:340379/141‑1 (MQ=255)
                           gTTTGTCGCCGTTGAAAACGAGGCGCGTGTTGTTTAATGTCAAGGCGGCATGACTGTTTACCGAGAGGGTAAGCATCAATAAGGTCAGGCAGACGCAAACAACGACGCGGGCGGGTATCGTGAATGTTGACGGATTGGTCATCGGGTaa                                                                                                                      >  1:27255/1‑149 (MQ=255)
                                 cGCCGTTGAAAACGAGGCGCGTGTTGTTTAATGTCAAGGCGGCATGACTGTTTACCGAGAGGGTAAGCATCAATAAGGTCAGGCAGACGCAAACAACGACGCGGGCGGGTATCGTGAATGTTGACGGATTGGTCATCGGGTAACTCGCt                                                                                                                <  8:347839/149‑1 (MQ=255)
                                     gTTGAAAACGAGGCGCGTGTTGTTTAATGTCAAGGCGGCATGACTGTTTACCGAGAGGGTAAGCATCAATAAGGTCAGGCAGACGCAAACAACGACGCGGGCGGGTATCGTGAATGTTGACGGATTGGTCATCGGGTAACTCGCTGCAt                                                                                                            <  5:169145/149‑1 (MQ=255)
                                      ttGAAAACGAGGCGCGTGTTGTTTAATGTCAAGGCGGCATGACTGTTTACCGAGAGGGTAAGCATCAATAAGGTCAGGCAGACGCAAACAACGACGCGGGCGGGTATCGTGAATGTTGACGGATTGGTCATCGGGTAACTCGCTGCATc                                                                                                           >  8:337809/1‑149 (MQ=255)
                                                 gcgcGTGTTGTTTAATGTCAAGGCGGCATGACTGTTTACCGAGAGGGTAAGCATCAATAAGGTCAGGCAGACGCAAACAACGACGCGGGCGGGTATCGTGAATGTTGACGGATTGGTCATCGGGTAACTCGCTGCATCCAGACGCACaa                                                                                                <  2:27255/149‑1 (MQ=255)
                                                 gcgcGTGTTGTTTAATGTCAAGGCGGCATGACTGTTTACCGAGAGGGTAAGCATCAATAAGGTCAGGCAGACGCAAACAACGACGCGGGCGGGTATCGTGAATGTTGACGGATTGGTCATCGGGTAACTCGCTGCATCCAGACGCACaa                                                                                                >  5:329983/1‑149 (MQ=255)
                                                        ttgttTAATGTCAAGGCGGCATGACTGTTTACCGAGAGGGTAAGCATCAATAAGGTCAGGCAGACGCAAACAACGACGCGGGCGGGTATCGTGAATGTTGACGGATTGGTCATCGGGTAACTCGCTGCATCCAGACGCACAACTACGcc                                                                                         >  5:353897/1‑149 (MQ=255)
                                                          gttTAATGTCAAGGCGGCATGACTGTTTACCGAGAGGGTAAGCATCAATAAGGTCAGGCAGACGCAAACAACGACGCGGGCGGGTATCATGAATGTTGACGGATTGGTCATCGGGTAACTCGCTGCATCCAGACGCACAACTACGCCTg                                                                                       >  8:227593/1‑149 (MQ=255)
                                                           tttAATGTCAAGGCGGCATGACTGTTTACCGAGAGGGTAAGCATCAATAAGGTCAGGCAGACGCAAACAACGACGCGGGCGGGTATCGTGAATGTTGACGGATTGGTCATCGGGTAACTCGCTGCATCCAGACGCACAACTACGCCTg                                                                                       >  3:90178/1‑148 (MQ=255)
                                                             tAATGTCAAGGCGGCATGACTGTTTACCGAGAGGGTAAGCATCAATAAGGTCAGGCAGACGCAAACAACGACGCGGGCGGGTATCGTGAATGTTGACGGATTGGTCATCGGGTAACTCGCTGCATCCAGACGCACAACTACGCCTGGAt                                                                                    >  4:48024/1‑149 (MQ=255)
                                                                    aaGGCGGCATGACTGTTTACCGAGAGGGTAAGCATCAATAAGGTCAGGCAGACGCAAACAACGACGCGGGCGGGTATCGTGAATGTTGACGGATTGGTCATCGGGTAACTCGCTGCATCCAGACGCACAACTACGCCTGGATGGTTTaa                                                                             >  6:13328/1‑149 (MQ=255)
                                                                     aGGCGGCATGACTGTTTACCGAGAGGGTAAGCATCAATAAGGTCAGGCAGACGCAAACAACGACGCGGGCGGGTATCATGAATGTTGACGGATTGGTCATCGGGTAACTCGCTGCATCCAGACGCACAACTACGCCTGGATGGTTTAAc                                                                            <  7:227593/149‑1 (MQ=255)
                                                                       gcggcATGACTGTTTACCGAGAGGGTAAGCATCAATAAGGTCAGGCAGACGCAAACAACGACGCGGGCGGGTATCGTGAATGTTGACGGATTGGTCATCGGGTAACTCGCTGCATCCAGACGCACAACTACGCCTGGATGGTTTAACCg                                                                          <  7:337809/149‑1 (MQ=255)
                                                                         ggcATGACTGTTTACCGAGAGGGTAAGCATCAATAAGGTCAGGCAGACGCAAACAACGACGCGGGCGGGTATCGTGAATGTTGACGGATTGGTCATCGGGTAACTCGCTGCATCCAGACGCACAACTACGCCTGGATGGTTTAACCGTg                                                                        <  3:48024/149‑1 (MQ=255)
                                                                             tGACTGTTTACCGAGAGGGTAAGCATCAATAAGGTCAGGCAGACGCAAACAACGACGCGGGCGGGTATCGTGAATGTTGACGGATTGGTCATCGGGTAACTCGCTGCATCCAGACGCACAACTACGCCTGGATGGTTTAACCGTGCAAg                                                                    <  2:339635/149‑1 (MQ=255)
                                                                             tGACTGTTTACCGAGAGGGTAAGCATCAATAAGGTCAGGCAGACGCAAACAACGACGCGGGCGGGTATCGTGAATGTTGACGGATTGGTCATCGGGTAACTCGCTGCATCCAGACGCACAACTACGCCTGGATGGTTTAACCGTGCAAg                                                                    <  6:329983/149‑1 (MQ=255)
                                                                               aCTGTTTACCGAGAGTGTAAGCATCAATAAGGTCAGGCAGACGCAAACAACGACGCGGGCGGGTATCGTGAATGTTGACGGATTGGTCATCGGGTAACTCGCTGCATCCAGACGCACAACTACGCCTGGATGGTTTAACCGTGCAAgg                                                                   <  5:200554/148‑1 (MQ=255)
                                                                                      aCCGAGAGGGTAAGCATCAATAAGGTCAGGCAGACGCAAACAACGACGCGGGCGGGTATCGTGAATGTTGACGGATTGGTCATCGGGTAACTCGCTGCATCCAGACGCACAACTACGCCTGGATGGTTTAACCGTGCAAGGACTGcac                                                            <  6:353897/148‑1 (MQ=255)
                                                                                       ccGAGAGGGTAAGCATCAATAAGGTCAGGCAGACGCAAACAACGACGCGGGCGGGTATCGTGAATGTTGACGGATTGGTCATCGGGTAACTCGCTGCATCCAGACGCACAACTACGCCTGGATGGTTTAACCGTGCAAGGACTGca                                                             >  1:25478/1‑146 (MQ=255)
                                                                                        cGAGAGGGTAAGCATCAATAAGGTCAGGCAGACGCAAACAACGACGCGGGCGGGTATCGTGAATGTTGACGGATTGGTCATCGGGTAACTCGCTGCATCCAGACGCACAACTACGCCTGGATGGTTTAACCGTGCAAGGACTGCACAg                                                          >  5:390670/1‑148 (MQ=255)
                                                                                             gggTAAGCATCAATAAGGTCAGGCAGACGCAAACAACGACGCGGGCGGGTATCGTGAATGTTGACGGATTGGTCATCGGGTAACTCGCTGCATCCAGACGCACAACTACGCCTGGATGGTTTAACCGTGCAAGGACTGCACAGAAttt                                                     <  6:390670/148‑1 (MQ=255)
                                                                                              ggTAAGCATCAATAAGGTCAGGCAGACGCAAACAACGACGCGGGCGGGTATCGTGAATGTTGACGGATTGGTCATCGGGTAACTCGCTGCATCCAGACGCACAACTACGCCTGGATGGTTTAACCGTGCAAGGACTGCACAGAATTTAt                                                   <  6:248439/149‑1 (MQ=255)
                                                                                                tAAGCATCAATAAGGTCAGGCAGACGCAAACAACGACGCGGGCGGGTATCGTGAATGTTGACGGATTGGTCATCGGGTAACTCGCTGCATCCAGACGCACAACTACGCCTGGATGGTTTAACCGTGCAAGGACTGCACAGAATTTATGc                                                 >  6:351127/1‑149 (MQ=255)
                                                                                                        aaTATGGTCATGCAGACGCAAACAACGACGCTGGCGGGTATCGTGAATGTTGACGGATTGGTCATCGGGTAACTCGCTGCATCCAGACGCACAACTACGCCTGGATGGTTTAACCGTGCAAGGACTGCACAGAATTTATGCGTAattat                                         <  8:189129/149‑1 (MQ=255)
                                                                                                            aGGTCAGGCAGACGCAAACAACGACGCGGGCGGGTATCATGAATGTTGACGGATTGGTCATCGGGTAACTCGCTGCATCCAGACGCACAACTACGCCTGGATGGTTTAACCGTGCAAGGACTGCACAGAATTTATGCGTAATTATTGa                                      <  7:379586/148‑1 (MQ=255)
                                                                                                             ggTCAGGCAGACGCAAACAACGACGCGGGCGGGTATCGTGAATGTTGACGGATTGGTCATCGGGTAACTCGCTGCATCCAGACGCACAACTACGCCTGGATGGTTTAACCGTGCAAGGACTGCACAGAATTTATGCGTAATTATTGa                                      >  1:29051/1‑147 (MQ=255)
                                                                                                                 aGGCAGACGCAAACAACGACGCGGGCGGGTATCGTGAATGTTGACGGATTGGTCATCGGGTAACTCGCTGCATCCAGACGCACAACTACGCCTGGATGGTTTAACCGTGCAAGGACTGCACAGAATTTATGCGTAATTATTGATAGTc                                 <  5:13328/148‑1 (MQ=255)
                                                                                                                   gCAGAAGCAAACAACGACGCGGGCGGGTATCGTGAATGTTGACGGATTGGTCATCGGGTAACTCGCTGCATCCAGACGCACAACTACGCCTGGATGGTTTAACCGTGCAAGGACTGCACAGATTTTATGCGTAATTATTGATAGTCGAc                              >  2:152220/1‑149 (MQ=255)
                                                                                                                         gCAAACAACGACGCGGGCGGGTATCGTGAATGTTGACGGATTGGTCATCGGGTAACTCGCTGCATCCAGACGCACAACTACGCCTGGATGGTTTAACCGTGCAAGGACTGCACAGAATTTATGCGTAATTATTGATAGTCGACGgtgt                         >  3:376792/1‑148 (MQ=255)
                                                                                                                         gCAAACAACGACGCGGGCGGGTATCGTGAATGTTGACGGATTGGTCATCGGATAACTCGCTGCATACAGACGCACAACTACGCCTGGATGGTTTAACCGTGCAAGGACAGCACAGAATTTATGCGTAATTATTGATAGTCGACGgtgt                         >  1:99039/1‑148 (MQ=255)
                                                                                                                           aaaCAACGACGCGGGCGGGTATCGTGAATGTTGACGGATTGGTCATCGGGTAACTCGCTGCATCCAGACGCACAACTACGCCTGGATGGTTTAACCGt                                                                         >  6:191475/1‑98 (MQ=255)
                                                                                                                           aaaCAACGACGCGGGCGGGTATCGTGAATGTTGACGGATTGGTCATCGGGTAACTCGCTGCATCCAGACGCACAACTACGCCTGGATGGTTTAACCGt                                                                         <  5:191475/98‑1 (MQ=255)
                                                                                                                               aacGACGCGGGCGGGTATCGTGAATGTTGACGGATTGGTCATCGGGTAACTCGCTGCATCCAGACGCACAACTACGCCTGGATGGTTTAACCGTGCAAGGACTGCACAGAATTTATGCGTAATTATTGATAGTCGACGGTGTAACGAAc                  >  8:6130/1‑149 (MQ=255)
                                                                                                                                    cgcgGGCGGGTATCGTGAATGTTGACGGATTGGTCATCGGGTAACTCGCTGCATCCAGACGCACAACTACGCCTGGATGGTTTAACCGTGCAAGGACTGCACAGAATTTATGCGTAATTATTGATAGTCGACGGTGTAACGAACCTggg             <  2:25478/149‑1 (MQ=255)
                                                                                                                                           gggTATCGTGAATGTTGACGGATTGGTCATCGGGTAACTCGCTGCATCCAGACGCACAACTACGCCTGGATGGTTTAACCGTGCAAGGACTGCACAGAATTTATGCGTAATTATTGATAGTCGACGGTGTAACGAACCTGGGCGGATg       >  8:61439/1‑148 (MQ=255)
                                                                                                                                               aTCGTGAATGTTGACGGATTGGTCATCGGGTAACTCGCTGCATCCAGACGCACAACTACGCCTGGATGGTTTAACCGTGCAAGGACTGCACAGAATTTATGCGTAATTATTGATAGTCGACGGTGTAACGAACCTGGGCGGATGCAAgg  >  6:379567/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GGTTTCGAACAACGACGGATGTACTGCGTTTGTCGCCGTTGAAAACGAGGCGCGTGTTGTTTAATGTCAAGGCGGCATGACTGTTTACCGAGAGGGTAAGCATCAATAAGGTCAGGCAGACGCAAACAACGACGCGGGCGGGTATCGTGAATGTTGACGGATTGGTCATCGGGTAACTCGCTGCATCCAGACGCACAACTACGCCTGGATGGTTTAACCGTGCAAGGACTGCACAGAATTTATGCGTAATTATTGATAGTCGACGGTGTAACGAACCTGGGCGGATGCAAGG  >  NC_002947/2132816‑2133107

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