Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 151,519 A→C 7.2% 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,5190AC 7.2% 150.2 / 3.8 56I27M (ATT→ATG) yadKputative fimbrial protein YadK
Reads supporting (aligned to +/- strand):  ref base A (25/27);  new base C (2/2);  total (27/29)
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.11e-01

GTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGA  >  NC_000913/151380‑151661
                                                                                                                                           |                                                                                                                                              
gTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATg                                                                                                                                             >  1:1191940/1‑143 (MQ=255)
  tGATAACTTTCAGGAGTAACATTACCGGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGc                                                                                                                                           <  2:882575/143‑1 (MQ=255)
    aTAACTTTCAGTAGTAAAATTACCTGAAAAATATCCCACCCGAACGACGCCAAGATCAATATTACCCCGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCCACAAGGTCAATATCCTGTCCAGCAATGGCTg                                                                                                                                         <  1:945203/143‑1 (MQ=255)
         cttCAGGAGGAACGTGACCTGAAAAATAGCCCACCCGAACGACGCCAGGACCAATAGTACCCTGGTTACTGACTCCGCTTTGGCATGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCCATTGCTGGACAg                                                                                                                                    <  1:278317/142‑1 (MQ=255)
           tCAGGAGTAACATTACCTGAAAACTATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCCATTGCTGGACAgga                                                                                                                                  <  1:210431/143‑2 (MQ=255)
           tCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGc                                                                                                                                  >  1:330506/1‑143 (MQ=255)
             aGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCCCTTGCAGGTACTGTTTTCTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGa                                                                                                                                >  2:332140/1‑143 (MQ=255)
              ggATTAACATTACCTTAAAAAGATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTGTTTTACAGTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGAt                                                                                                                               <  2:551531/143‑1 (MQ=255)
                aGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGc                                                                                                                             >  2:820139/1‑143 (MQ=255)
                      atcaCCGGAAAAACCTCCCACCCCCACGACGCCAACCTCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACa                                                                                                                       <  1:1049141/140‑1 (MQ=255)
                          ccTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTaa                                                                                                                   <  2:1074433/143‑1 (MQ=255)
                           cTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAAt                                                                                                                  >  1:1106530/1‑143 (MQ=255)
                            tGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATg                                                                                                                 >  2:183643/1‑143 (MQ=255)
                              aaaaaTATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTaa                                                                                                                   >  1:346930/1‑139 (MQ=255)
                              aaaaaTATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTaa                                                                                                                   <  2:346930/139‑1 (MQ=255)
                                 aaTATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAATAATGaaaga                                                                                                            >  2:581873/1‑143 (MQ=255)
                                     tCCCACCGCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacagaaaca                                                                                                        <  2:267316/143‑1 (MQ=255)
                                        cACCCCAACGACGCCAAGATCAATATTACCCTTGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGTTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacagaaacagaa                                                                                                     <  2:859080/143‑1 (MQ=255)
                                        cACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacagaaacagaa                                                                                                     <  1:982682/143‑1 (MQ=255)
                                             cAACGACGCCAAGATCAAGATTACCCTGGTTAGTGATTCCGCTTTTGCAGGTGCTGTGTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAattatt                                                                                                <  2:1113705/143‑1 (MQ=255)
                                              aaCGACG‑CAAGATCAATATTACCCTGCTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATtct                                                                                              <  1:417931/143‑1 (MQ=255)
                                               acgacgCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATtct                                                                                              <  2:1132795/143‑1 (MQ=255)
                                                 gacgCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTctct                                                                                            >  1:22689/1‑143 (MQ=255)
                                                  acgCCAAGATCAATATGACCCTGGTTACTGATTCCGCTGTTGCAGGTGCTGTTTTTTACATGGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCtt                                                                                           <  1:575911/143‑1 (MQ=255)
                                                    gCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTtca                                                                                         <  1:451803/143‑1 (MQ=255)
                                                    gCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTGTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTtca                                                                                         >  1:115737/1‑143 (MQ=255)
                                                      cAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTtcatc                                                                                       >  1:1242820/1‑143 (MQ=255)
                                                      cAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTtcatc                                                                                       >  1:789194/1‑143 (MQ=255)
                                                         gATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGAAAAACAGTAATGACAGAAACAGAATTAGTCTCTtcatcacc                                                                                    >  2:893064/1‑141 (MQ=255)
                                                                      cTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTgg                                                                       <  2:983160/143‑1 (MQ=255)
                                                                         gTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGAGGGAt                                                                    >  2:936748/1‑143 (MQ=255)
                                                                                 ttCCGCTGTTGCAGGTGCTGTTTCTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCatatat                                                            <  1:85760/143‑1 (MQ=255)
                                                                                  tCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACACGGTCAATATCCTGTCCAGCAATGGCTGGACGGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGGGTGGGATGCAtatata                                                           >  2:1208398/1‑143 (MQ=255)
                                                                                  tCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCAtatata                                                           <  1:820139/143‑1 (MQ=255)
                                                                                      cTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTc                                                       <  2:459510/143‑1 (MQ=255)
                                                                                        tttGCAGGTGCTGTTTTTTACATTGGCAACACGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTCCGTTGAGTGGGATGCATATATAGTTCCt                                                     <  1:1219796/143‑1 (MQ=255)
                                                                                            cAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGGAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAAt                                                 <  1:385670/143‑1 (MQ=255)
                                                                                              ggTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGACTTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCt                                               <  1:701743/143‑1 (MQ=255)
                                                                                              ggTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCt                                               >  1:1184585/1‑143 (MQ=255)
                                                                                                tGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGt                                             <  2:789194/143‑1 (MQ=255)
                                                                                                    gTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGCGAGGCGATGCAAAACAGTAGTGAGAGCAACAGAAGAATACAACAAACAAGCTTTCGTTGAGTGGGATACATATATAGTTCCTGAATCTGTAAGt                                         >  2:1245035/1‑143 (MQ=255)
                                                                                                    gTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGCATCTGTACGt                                         >  2:413661/1‑143 (MQ=255)
                                                                                                    gTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCCCTAAgc                                         >  2:1094923/1‑142 (MQ=255)
                                                                                                            cATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTgttagt                                 >  1:387736/1‑143 (MQ=255)
                                                                                                               tGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTtat                              <  1:790490/143‑1 (MQ=255)
                                                                                                                ggCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGAAAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTAAGTTGAGTGGGCTGCATATATGGTTACGGAATCTGTGAGTCTgtgagttaga                             >  1:82787/1‑141 (MQ=255)
                                                                                                                   aacaaGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGt                          >  2:580577/1‑143 (MQ=255)
                                                                                                                          tCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATCTACAGTTCCTGACTCTGCAAGTCTGTTAGTTATAGGTTACTGTa                   >  2:999223/1‑143 (MQ=255)
                                                                                                                               cgccTGTCCAGCCATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAAtgt              >  1:257106/3‑143 (MQ=255)
                                                                                                                                gtcgtcccATCAATGGCTGAACAGGCGATGCAAAACAGTAATGACAGAAACAGTATTATTCTCCTCATCAGCTTACGTTGAGGGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAAtgtg             <  2:383388/137‑1 (MQ=255)
                                                                                                                                 ccTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGTTTTATTTTCTTCATTATCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTgg            <  1:1068989/143‑1 (MQ=255)
                                                                                                                                 ccTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGGTACTGTAAATGTgg            >  1:55774/1‑143 (MQ=255)
                                                                                                                                  cTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTGATAGGTTACTGTAAACGTGGc           >  2:176540/1‑143 (MQ=255)
                                                                                                                                     tCCAGCCATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTgg        >  1:217695/1‑143 (MQ=255)
                                                                                                                                      ccAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGt       <  1:183643/143‑1 (MQ=255)
                                                                                                                                      ccAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGACTGGt       <  1:1201720/143‑1 (MQ=255)
                                                                                                                                           aaTGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGa  <  1:587715/143‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                              
GTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGA  >  NC_000913/151380‑151661

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