Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,131,453 A→C 12.0% L249R (CTG→CGG) ‡ 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,4530AC12.0% 126.5 / 9.2 49L249R (CTG→CGG) ‡fimDtype I pili subunit FimD
Reads supporting (aligned to +/- strand):  ref base A (25/18);  new base C (3/3);  total (28/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 = 8.75e-01

CCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCTTGCCCTGGCTGCCC‑‑TGGCGCCAGGACTGGGTACTGCGCAGGCGCCATTCGCCCAGGTTGATACCACTGTTCAGGTACAGGTCGT  >  NC_002947/2131310‑2131600
                                                                                                                                               |                                                                                                                                                     
cccGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTg                                                                                                                                                  >  8:103027/1‑149 (MQ=255)
         gggCCACGCCGCGAATGATCGGAGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTa                                                                                                                                         >  3:172936/1‑149 (MQ=255)
          ggCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGa                                                                                                                                              <  6:134427/143‑1 (MQ=255)
          ggCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGa                                                                                                                                              >  5:134427/1‑143 (MQ=255)
          ggCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTCTCGCCAATCTGACGATTGGCGaaacc                                                                                                                                          <  4:150055/147‑5 (MQ=255)
           gCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTAcc                                                                                                                                       >  1:176425/1‑149 (MQ=255)
               cgccgcGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGt                                                                                                                                   <  5:152300/149‑1 (MQ=255)
                gccgcGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGta                                                                                                                                  >  6:137429/1‑149 (MQ=255)
                  cgcgAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAg                                                                                                                                             <  1:223668/136‑1 (MQ=255)
                  cgcgAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAg                                                                                                                                             >  2:223668/1‑136 (MQ=255)
                             ggggCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGc                                                                                                                     <  7:103027/149‑1 (MQ=255)
                               agCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGt                                                                                                                   <  4:172936/148‑1 (MQ=255)
                                              gTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTCGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCgg                                                                                                     <  1:229445/148‑1 (MQ=255)
                                               ttGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCATCAGATTGGCGAGGGTACCCGGTATATcgcg                                                                                                                           <  3:9581/125‑1 (MQ=255)
                                               ttGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCATCAGATTGGCGAGGGTACCCGGTATATcgcg                                                                                                                           >  4:9581/1‑125 (MQ=255)
                                                 gTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTCTCGCCAATCTGACGATTGGCGaaacc                                                                                                                                          <  6:122145/108‑5 (MQ=38)
                                                 gTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTCTCGCCAATCTGACGATTGGCGaaacc                                                                                                                                          >  5:122145/1‑104 (MQ=38)
                                                   cGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCAt                                                                                               <  2:176425/149‑1 (MQ=255)
                                                    ggCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCtgt                                                                                                                       <  3:71964/124‑1 (MQ=255)
                                                    ggCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCtgt                                                                                                                       >  4:71964/1‑124 (MQ=255)
                                                    ggCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCAtt                                                                                              <  5:137429/149‑1 (MQ=255)
                                                          cATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCt                                                                                        >  6:168343/1‑149 (MQ=255)
                                                               ccATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCtt                                                                                       <  3:5012/145‑1 (MQ=255)
                                                               ccATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCTtgccc                                                                                   >  2:216632/1‑149 (MQ=255)
                                                               ccATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGAGCTGTGCGTAGGTCTGTGCCCTGGTCCATTAAGTCtt                                                                                       >  4:5012/1‑145 (MQ=255)
                                                                    tCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCTtgccctgg                                                                                <  4:1075/147‑1 (MQ=255)
                                                                    tCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCTtgccctgg                                                                                >  3:1075/1‑147 (MQ=255)
                                                                      cGAGGCCCCGCGCACGCCTGAAATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCAGTGCGTAGGTCTGTGCACGGGTCCATTCAGTCAtgccctggctga                                                                            >  4:112670/1‑148 (MQ=255)
                                                                      cGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGt                                                                                                               <  6:212340/114‑1 (MQ=255)
                                                                      cGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGt                                                                                                               >  5:212340/1‑114 (MQ=255)
                                                                            cACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAAGCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCtgt                                                                                                                       >  8:55738/1‑100 (MQ=255)
                                                                                     gCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCTTGCCCTGGCTGCCC‑‑TGGCGCCAGGACt                                                           >  3:76858/1‑149 (MQ=255)
                                                                                       cTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCTTGCCCTGGCTGCCC‑‑TGGCGCCAGGACTgg                                                         <  5:168343/149‑1 (MQ=255)
                                                                                          atatCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCTTGCCCTGGCTGCCC‑‑TGGCGCCAGGACTggg                                                        <  3:112670/147‑1 (MQ=255)
                                                                                            atCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCTTGCCCTGGCTGCCC‑‑TGGCGCCAGGACTGGGTACt                                                    >  8:251094/1‑149 (MQ=255)
                                                                                                      aaCGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCTTGCCCTGGCTGCCC‑‑TGGCGCCAGGACTGGGTACTGCGCAGgcgc                                          >  3:138327/1‑149 (MQ=255)
                                                                                                         ggAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCTTGCCCTGGCTGCCC‑‑TGGCGCCAGGACTGGGTACTGCGCAGGCGCCAt                                       >  2:136819/1‑149 (MQ=255)
                                                                                                           aaCACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCTTGCCCTGGCTGCCC‑‑TGGCGCCAGGACt                                                           >  1:251612/1‑127 (MQ=255)
                                                                                                           aaCACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCTTGCCCTGGCTGCCC‑‑TGGCGCCAGGACt                                                           <  2:251612/127‑1 (MQ=255)
                                                                                                                cTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCTTGCCCTGGCTGCCC‑‑TGGCGCCAGGACTGGGTACTGCGCAGGCGCCATTCGcc                                  >  2:101671/1‑147 (MQ=255)
                                                                                                                cTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCTTGCCCTGGCTGCCC‑‑TGGCGCCAGGACTGGGTACTGCGCAGGCGCCATTCGcc                                  >  5:309478/1‑147 (MQ=255)
                                                                                                                         gatctaGGTTTCGCCACTCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCTTGCCCTGGCTGCCC‑‑TGGCGCCAGGACTGGGTACTGCGCAGGCGCCATTCGCCCAGGTTGATAc                       <  1:15304/144‑1 (MQ=255)
                                                                                                                               ggTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCTTGCCCTGGCTGCCC‑‑TGGCGCCAGGACTGGGTACTGCGCAGGCGCCATTCGCCCAGGTTGATAc                       >  2:15304/1‑143 (MQ=255)
                                                                                                                                  ttCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCTTGTAATAATAGCCCTATCGGGCAAGGAGTGGTAAGTGCGCAGGCGCCTTTCGCCCAGGTTGATCCCACTGt                 >  1:269244/1‑148 (MQ=255)
                                                                                                                                   tCGCCAATCTGACGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCTTGCCCTGGCTGCCC‑‑TGGCGCCAGGACTGGGTACTGCGCAGGCGCCATTCGCCCAGGtt                            <  6:42349/134‑1 (MQ=255)
                                                                                                                                   tCGCCAATCTGACGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCTTGCCCTGGCTGCCC‑‑TGGCGCCAGGACTGGGTACTGCGCAGGCGCCATTCGCCCAGGtt                            >  5:42349/1‑134 (MQ=255)
                                                                                                                                   tCGCCAATCTGACGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCTTGCCCTGGCTGCCC‑‑TGGCGCCAGGACTGGGTACTGCGCAGGCGCCATTCGCCCAGGTTGATACCACTGTTCAg             >  6:148807/1‑149 (MQ=255)
                                                                                                                                             tCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCTTGCCCTGGCTGCCC‑‑TGGCGCCAGGACTGGGTACt                                                    <  7:117308/100‑1 (MQ=255)
                                                                                                                                             tCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCTTGCCCTGGCTGCCC‑‑TGGCGCCAGGACTGGGTACt                                                    >  8:117308/1‑100 (MQ=255)
                                                                                                                                               aGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCTTGCCCTGGCTGCCC‑‑TGGCGCCAGGACTGGGTACTGCTCAGGCGCCATTCGCCCAGGTTGATACCACTGTTCAGGTACAGGTCGt  <  6:309478/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
CCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCTTGCCCTGGCTGCCC‑‑TGGCGCCAGGACTGGGTACTGCGCAGGCGCCATTCGCCCAGGTTGATACCACTGTTCAGGTACAGGTCGT  >  NC_002947/2131310‑2131600

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