Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 151,519 A→C 12.4% I27M (ATT→ATG yadK ← putative fimbrial protein YadK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913151,5190AC12.4% 124.5 / 13.6 58I27M (ATT→ATGyadKputative fimbrial protein YadK
Reads supporting (aligned to +/- strand):  ref base A (28/22);  new base C (3/4);  total (31/27)
Fisher's exact test for biased strand distribution p-value = 6.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.13e-02

GTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAG  >  NC_000913/151380‑151660
                                                                                                                                           |                                                                                                                                             
gTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATg                                                                                                                                            >  1:753975/1‑143 (MQ=255)
  ttaTAACTGTCAGCAGTAACATTACCTGAAAAATATCCCACCCCATCGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCGTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGc                                                                                                                                          <  2:318111/141‑1 (MQ=255)
   gATACCTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCt                                                                                                                                         <  2:202326/143‑1 (MQ=255)
        ccttCAGGAGTAACACTACCTTAAAAATATCCCCCCCTAACGACGCCAACATCAATATTACCCTGGTTACTGATTCCGCTTTTTCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCCATTGCTGGACa                                                                                                                                    <  1:789277/141‑1 (MQ=255)
        cctgcAGGAGTAACATTACCTGAAATATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCCATTGCTGGACa                                                                                                                                    <  2:152826/139‑1 (MQ=255)
        cTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCCATTGCTGGACa                                                                                                                                    <  1:591174/143‑1 (MQ=255)
        cTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCCATTGCTGGACa                                                                                                                                    <  2:69320/143‑1 (MQ=255)
             aGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGa                                                                                                                               <  2:549128/143‑1 (MQ=255)
              ggAGTTACAGCACCTGAAAAATATCCCACCCCAGCGACGCCAAGAGCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGAt                                                                                                                              <  2:251258/143‑1 (MQ=255)
                 gTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAACGGCTGGACAGGCGATGCa                                                                                                                           >  1:39991/1‑143 (MQ=255)
                  tcaCATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCaa                                                                                                                          <  2:236247/141‑1 (MQ=255)
                     cATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCTGCAATGGCTGGACAGGCGATGCAAAAc                                                                                                                       >  2:422153/1‑143 (MQ=255)
                      attttcTGAAAAATAACCAACCCCAACGAGGCCAAGATCAAGATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCCATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACa                                                                                                                      <  1:350387/138‑1 (MQ=255)
                            tGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTATGCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTGTTTTACATTGGCCACAAGGTCGATATCCTGTCCAGCAACGGCTGGACAGGCGATGCAAAACAGTAATg                                                                                                                >  2:637433/1‑143 (MQ=255)
                              aaaaaTATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGac                                                                                                              <  2:889843/143‑1 (MQ=255)
                                aaaTATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGCTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGCTGCAAAACAGTAATGacag                                                                                                            >  2:806962/1‑143 (MQ=255)
                                          ccccAACTACGCCAAGATCAATATTACCCTGTTTACTGATTCCGCTTTTGCTGGTGTTGTTTTTTATTTTGTCAACAAGGTCAATAGCCTGGCGAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAatt                                                                                                  <  1:621386/143‑1 (MQ=255)
                                           cccAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAatta                                                                                                 >  1:745132/1‑143 (MQ=255)
                                              aaCGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCCATGCAAAACAGTAATGACAGCAACAGAATTATtc                                                                                              >  1:692691/1‑143 (MQ=255)
                                                 gacgCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTctct                                                                                           >  2:415864/1‑143 (MQ=255)
                                                      cAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTtcatc                                                                                      >  2:161498/1‑143 (MQ=255)
                                                          aTCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGCCAGAAACAGAATTATTCTCTTCATCAGCt                                                                                  >  1:21494/1‑143 (MQ=255)
                                                             aaTATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACCGACACCGCATTATTCTCTTCATCAGCttgc                                                                               >  1:281005/1‑141 (MQ=255)
                                                                    cccTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTATTCATTAGCTTCCGTTGAGt                                                                        >  2:880831/1‑143 (MQ=255)
                                                                     ccTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTg                                                                       >  1:309361/1‑143 (MQ=255)
                                                                          ttACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGCCAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATg                                                                  <  1:702435/143‑1 (MQ=255)
                                                                            aCTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCGAGCAATGGCTGGACAGGCGATGCATAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCa                                                                <  2:444121/143‑1 (MQ=255)
                                                                            aCTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGATCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAACCAGTCATGCCAGAAACAGCATTATTCTCTTCATCAGCTTACGTTGAGTCGGCTGCa                                                                >  1:781471/1‑143 (MQ=255)
                                                                            aCTGATTCCGCTTTTGCAGGTGCTGCTGTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCa                                                                <  2:644761/143‑1 (MQ=255)
                                                                              tGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCACTATCCTGTCCAGCAATGGCTGGACAGGCGGTGCAAAACAGTAATGACAGAAACAAAATTATTCTCTTCATCAGCTTACGTTGGGTGGGATGCata                                                              >  1:819385/1‑143 (MQ=255)
                                                                                aTTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGGGGGATGCagata                                                            >  2:178790/1‑143 (MQ=255)
                                                                                     gCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATCTATAGtt                                                       >  1:623924/1‑143 (MQ=255)
                                                                                      cTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTc                                                      <  1:149401/143‑1 (MQ=255)
                                                                                       ttttGCAGGTGCGGTTGTCTACGTTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTcc                                                     <  1:26344/143‑1 (MQ=255)
                                                                                       ctttGCAGGTGCTGCGTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTcc                                                     <  1:361316/142‑1 (MQ=255)
                                                                                        tttGCAGGTGCTGTTTTTTACATTGGCAACACGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGACGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCACCCGCTACCGCTGAGTGGGAGCCATATCTAGTTCCt                                                    >  2:132552/1‑143 (MQ=255)
                                                                                            cAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGCAt                                                >  2:815750/1‑143 (MQ=255)
                                                                                            cAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAAt                                                <  2:694658/143‑1 (MQ=255)
                                                                                            cAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAAt                                                >  1:532980/1‑143 (MQ=255)
                                                                                                 gcggTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTa                                           <  2:667581/140‑1 (MQ=255)
                                                                                                  cTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTaa                                          <  1:247472/143‑1 (MQ=255)
                                                                                                     ttttttACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAACACAGTAATGACAGAAACAGAATTATTCTCTTCATGAGCTTCCGCTGAGTACGCTGCAAATATAGTTCCTGAATCCGTAAGTc                                       >  2:630735/1‑143 (MQ=255)
                                                                                                     ttttttACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTc                                       >  1:822668/1‑143 (MQ=255)
                                                                                                            cATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAACTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTgttagt                                <  2:447944/143‑1 (MQ=255)
                                                                                                              ttGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTgttcgtta                              >  1:154906/1‑143 (MQ=255)
                                                                                                              ttGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATGATTCTCTTCATCAGCTTACGTAGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGtgtgacagtga                              >  2:735100/1‑136 (MQ=255)
                                                                                                                ggCAACAAGGTCAATTTCCTGCCCAGCAATGGCGTGACAGGCGATCCAAAACAGTAATGACAGAAACAGAATTATTCTCTGCATCAGCTTACGTTGAGTTGGATGCATATATAGTTCCTGAATCTGTAAGTCTCTTAGTtata                            <  2:714815/143‑1 (MQ=255)
                                                                                                                 gCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATTTCGTTCCTGAATCTGTAAGTCTGTTAGTTATAg                           <  1:89890/143‑1 (MQ=255)
                                                                                                                     caaGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTa                       <  1:454732/143‑1 (MQ=255)
                                                                                                                      aaGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTAc                      >  1:205428/1‑143 (MQ=255)
                                                                                                                         gTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAAGCTGTAAGTCTGTTAGGTATAGGTTACTGt                   >  2:174220/1‑143 (MQ=255)
                                                                                                                            aaTATCCTGTCCAGCAAAGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTaaa                <  1:457762/143‑1 (MQ=255)
                                                                                                                              tcgccTGTCCAGCCATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAAtg              >  2:205060/4‑143 (MQ=255)
                                                                                                                                 ccTGTCCTGCGATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGACTCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTgg           <  2:109230/143‑1 (MQ=255)
                                                                                                                                 ccTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTgg           >  1:770590/1‑143 (MQ=255)
                                                                                                                                  cTGTCCAGCCATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGCTTTAGGTTACTGTCAATTTGGc          >  1:142238/1‑143 (MQ=255)
                                                                                                                                   tGTCCAGCCATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCt         >  1:477345/1‑143 (MQ=255)
                                                                                                                                          cAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGGGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAg  >  1:69440/1‑143 (MQ=255)
                                                                                                                                           |                                                                                                                                             
GTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAG  >  NC_000913/151380‑151660

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: