Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,131,452 C→A 10.0% L249L (CTG→CTT) ‡ 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,4520CA10.0% 122.1 / 4.1 40L249L (CTG→CTT) ‡fimDtype I pili subunit FimD
Reads supporting (aligned to +/- strand):  ref base C (19/17);  new base A (2/2);  total (21/19)
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.15e-01

CGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCTTGCCCTGGCTGCCCTGGCGCCAGGACTGGGTACTGCGCAGGCGCCATTCGCCCAGGTTGATACCACTGTTCAGGTACAGGTCG  >  NC_002947/2131312‑2131599
                                                                                                                                            |                                                                                                                                                   
cGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGc                                                                                                                                             >  3:260960/1‑149 (MQ=255)
      tGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAggg                                                                                                                                        >  1:58512/1‑148 (MQ=255)
      tGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAggg                                                                                                                                        >  4:25135/1‑148 (MQ=255)
              gccgcGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGta                                                                                                                               <  2:113695/149‑1 (MQ=255)
                 gcgAATGATCGGGGCATAACCGCTCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTatat                                                                                                                            <  6:247347/149‑1 (MQ=255)
                           ggggCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGc                                                                                                                  <  4:255875/149‑1 (MQ=255)
                                             ttGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGccc                                                                                                    >  7:120339/1‑145 (MQ=255)
                                             ttGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGccc                                                                                                    <  8:120339/145‑1 (MQ=255)
                                                tCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCa                                                                                             <  4:260960/149‑1 (MQ=255)
                                                 cGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCAt                                                                                            >  6:62487/1‑149 (MQ=255)
                                                  ggCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCTGACGATTGGCGaaac                                                                                                                                        >  5:80144/1‑101 (MQ=255)
                                                  ggCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCTGACGATTGGCGaaac                                                                                                                                        <  6:80144/104‑4 (MQ=255)
                                                     cAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACGGAACGGAACACCTCCGTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATc                                                                                                                           >  1:192128/1‑114 (MQ=255)
                                                     cAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCTTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATc                                                                                                                           <  2:192128/114‑1 (MQ=255)
                                                                  tCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCTtgccctggct                                                                           >  8:33671/1‑149 (MQ=255)
                                                                  tCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCTtgccctggct                                                                           <  5:204079/149‑1 (MQ=255)
                                                                          cACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCTtgccctggctgccctggc                                                                   >  7:1466/1‑149 (MQ=255)
                                                                              cgcACGCCTGATATCGCCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCTTGCCCTGGCTGCCCTGGCGCCa                                                               <  3:35029/149‑1 (MQ=255)
                                                                                 aCGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCTTGCCCTGGCTGCCCTGGCGCCAGGa                                                            <  8:1466/149‑1 (MQ=255)
                                                                                 aCGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCTTGCCCTGGCTGCCCTGGCGCCAGGa                                                            <  2:58512/149‑1 (MQ=255)
                                                                                  cGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCgg                                                                                                  <  2:58621/110‑1 (MQ=255)
                                                                                   gCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCTTGCCCTGGCTGCCCTGGCGCCAGGACt                                                          >  5:112749/1‑149 (MQ=255)
                                                                                   gCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCTTGCCCTGGCTGCCCTGGCGCCAGGACt                                                          >  8:163127/1‑149 (MQ=255)
                                                                                    ccTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCTTGCCCTGGCTGCCCTGGCGCCAGGACTg                                                         <  5:62487/149‑1 (MQ=255)
                                                                                               cACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCTTGCCCTGGCTGCCCTGGCGCCAGGACTGGGTACTGCGCa                                              <  3:25135/149‑1 (MQ=255)
                                                                                                   gAACGGAACAACTACCTGTCCGTCAATGTTTCGCCAATCGTCACATTGGCTAGGGTACCCGGTATATCGCTCTGTGAGTAGGTCTTTGACCGGGTCCAATCAGTCTTGCACTTGCTGCCCTGTCGCCAGGACTGGGTACTGCGCAGtcg                                          >  4:167579/1‑149 (MQ=255)
                                                                                                          acacCTCCCTGTCCGAGAAGGTTTCGCCAATAGACAGATTGGCGAGGGAACCCGGTATATAGCGCTCTGCGTAGGTCTGTGCCCGGGTCCATACACACATGCACTGGCTGCCCTGGCGCCAGGACTGCGTACTGCGCAGGCGCCAATCg                                   >  1:117728/1‑149 (MQ=255)
                                                                                                                    gTCCGGGATGGTTTCGCCAATCGTCAGATTGGCGTGGGTACCGGGTATATCGCGCTGTGCGTAGGTCTGTGCCCgg                                                                                                  <  1:15885/76‑1 (MQ=255)
                                                                                                                     tCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCTTGCCCTGGCTGCCCTGGCGCCAGGACTGGGTACTGCGCAGGCGCCATTCGCCCAGGTTGa                         >  7:204021/1‑148 (MQ=255)
                                                                                                                        gTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCTTGCCCTGGCTGCCCTGGCGCCAGGACTGGGTACTGCGCAGGCGCCATTCGCCCAGGTTGATAcc                     >  7:176478/1‑149 (MQ=255)
                                                                                                                           aaGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCCTAGGTCTGTGCCCGGGTCCATTCAGTCTTGCCCTGGCTGCCCTGGCGCCAGGACTGGGTACTGCGCAGGCGCCATTCGCCCAGGTTGATACCACt                  >  8:220425/1‑149 (MQ=255)
                                                                                                                                ctCGCCAATCTGACGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCTTGCCCTGGCTGCCCTGGCGCCAGGACTGGGTACTGCGCAGGCGCCATTCGCCCAgg                             <  6:47794/132‑1 (MQ=255)
                                                                                                                                ctCGCCAATCTGACGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCTTGCCCTGGCTGCCCTGGCGCCAGGACTGGGTACTGCGCAGGCGCCATTCGCCCAgg                             >  5:47794/2‑133 (MQ=255)
                                                                                                                                 tCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCTTGCCCTGGCTGCCCTGGCGCCAGGACTGGGTACTGCGCAGGCGCCATTCGCCCAGGTTg                          >  7:35978/1‑135 (MQ=255)
                                                                                                                                 tCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCTTGCCCTGGCTGCCCTGGCGCCAGGACTGGGTACTGCGCAGGCGCCATTCGCCCAGGTTg                          <  8:35978/135‑1 (MQ=255)
                                                                                                                                 tCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCTTGCCCTGGCTGCCCTGGCGCCAGGACTGGGTACTGCGCAGGCGCCATTCGCCCAGGTTGATACCACTGTTCa             <  8:176478/148‑1 (MQ=255)
                                                                                                                                    ccAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCTTGCCCTGGCTGCCCTGGCGCCAGGACTggg                                                       <  6:6684/103‑1 (MQ=255)
                                                                                                                                    ccAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCTTGCCCTGGCTGCCCTGGCGCCAGGACTggg                                                       >  5:6684/1‑103 (MQ=255)
                                                                                                                                     cAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCTTGCCCTGGCTGCCCTGGCGCCAGGACTGGGTACTGCGCAGGCGCCATTCGCCCAGGTTGATACCACTGTTCAGGTa         >  7:156393/1‑148 (MQ=255)
                                                                                                                                           tCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCTTGCCCTGGCTGCCCTGGCGCCAGGACTGGGTACTGCGCAGGCGCCATTCGCCCAGGTTGATACCACTGTTCAGGTACAGGTCg  >  8:222709/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                                   
CGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCTTGCCCTGGCTGCCCTGGCGCCAGGACTGGGTACTGCGCAGGCGCCATTCGCCCAGGTTGATACCACTGTTCAGGTACAGGTCG  >  NC_002947/2131312‑2131599

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: