Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,131,452 C→A 12.3% 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,4520CA12.3% 133.1 / 9.9 49L249L (CTG→CTT) ‡fimDtype I pili subunit FimD
Reads supporting (aligned to +/- strand):  ref base C (25/18);  new base A (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 = 4.40e-01

CTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCTTGCCCTGGCTGCCC‑‑TGGCGCCAGGACTGGGTACTGCGCAGGCGCCATTCGCCCAGGTTGATACCACTGTTCAG  >  NC_002947/2131304‑2131589
                                                                                                                                                    |                                                                                                                                           
cTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTc                                                                                                                                             <  5:126381/149‑1 (MQ=255)
cTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTc                                                                                                                                             <  3:68616/149‑1 (MQ=255)
cTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTc                                                                                                                                             <  3:1278/149‑1 (MQ=255)
      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                                                                                                                >  4:9581/1‑125 (MQ=255)
                                                     ttGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCATCAGATTGGCGAGGGTACCCGGTATATcgcg                                                                                                                <  3:9581/125‑1 (MQ=255)
                                                       gTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTCTCGCCAATCTGACGATTGGCGaaacc                                                                                                                               >  5:122145/1‑104 (MQ=38)
                                                       gTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTCTCGCCAATCTGACGATTGGCGaaacc                                                                                                                               <  6:122145/108‑5 (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                                                <  2:251612/127‑1 (MQ=255)
                                                                                                                 aaCACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCTTGCCCTGGCTGCCC‑‑TGGCGCCAGGACt                                                >  1:251612/1‑127 (MQ=255)
                                                                                                                      cTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCTTGCCCTGGCTGCCC‑‑TGGCGCCAGGACTGGGTACTGCGCAGGCGCCATTCGcc                       >  5:309478/1‑147 (MQ=255)
                                                                                                                      cTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCTTGCCCTGGCTGCCC‑‑TGGCGCCAGGACTGGGTACTGCGCAGGCGCCATTCGcc                       >  2:101671/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                 >  5:42349/1‑134 (MQ=255)
                                                                                                                                         tCGCCAATCTGACGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCTTGCCCTGGCTGCCC‑‑TGGCGCCAGGACTGGGTACTGCGCAGGCGCCATTCGCCCAGGtt                 <  6:42349/134‑1 (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)
                                                                                                                                                    |                                                                                                                                           
CTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCTTGCCCTGGCTGCCC‑‑TGGCGCCAGGACTGGGTACTGCGCAGGCGCCATTCGCCCAGGTTGATACCACTGTTCAG  >  NC_002947/2131304‑2131589

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