Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913_3_pae_tpiA 2,106,741 A→T 100% S162S (TCT→TCA wbbH ← O‑antigen polymerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_tpiA2,106,7410AT98.3% 179.8 / ‑3.7 58S162S (TCT→TCAwbbHO‑antigen polymerase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base T (29/28);  total (30/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.38e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ACTCTATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTAGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAACC  >  NC_000913_3_pae_tpiA/2106595‑2106883
                                                                                                                                                  |                                                                                                                                              
gctctATTAACACATACGATGAATGCATAAGAGATTATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAAATGCAAAGATCAATACAATAAAAGCAAGTAAAGTAAATGTTTTACTTTCTTTTGTATTAGTAAATTTTTTTGa                                                                                                                                              <  2:265712/148‑1 (MQ=255)
aCTCTATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATCAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGa                                                                                                                                              <  1:477630/149‑1 (MQ=255)
 ctctATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGAc                                                                                                                                             >  2:422633/1‑149 (MQ=255)
  tctATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGAcc                                                                                                                                            >  1:289393/1‑149 (MQ=255)
  tctATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGAcc                                                                                                                                            >  1:248117/1‑149 (MQ=255)
   ctATTAACACCTGCGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCa                                                                                                                                           <  1:572409/149‑1 (MQ=255)
      ttAACACCGACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCaaat                                                                                                                                        <  2:570507/149‑1 (MQ=255)
        aaCACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAAAGTTTTACTTACCTTTGTATAAGTAAATTTTTTTGACCaaatac                                                                                                                                      >  2:366089/1‑148 (MQ=255)
             cTACGATGAATGCATAAGGGATGATAATCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAAtaaataaa                                                                                                                                 <  2:627984/149‑1 (MQ=255)
                  aTGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCaag                                                                                                                            >  1:365302/1‑148 (MQ=255)
                    gAAGGCATAAGAGATGATAACCATAAAGAAAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAag                                                                                                                          <  2:257904/149‑1 (MQ=255)
                      aTGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAGGGTGAATACAATAAAAACAATTAAAGTAAATGTTTTACTTACATTTGTATTAGTAAATGTTTGTGACCAAATAAATAAAGCAAAagta                                                                                                                        >  1:725595/1‑149 (MQ=255)
                       tGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAagtag                                                                                                                       <  2:650365/149‑1 (MQ=255)
                       tGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAaatag                                                                                                                       >  1:243607/1‑149 (MQ=255)
                        gCAAAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAATACAACGAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAagtagt                                                                                                                      <  2:447014/149‑1 (MQ=255)
                           tAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATAAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAg                                                                                                                   >  1:490063/1‑149 (MQ=255)
                                    atAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTTATTATTTATATAAAGCAAAAGTAGTTAGAATGATTgg                                                                                                          >  2:55952/1‑149 (MQ=255)
                                          aTAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCAt                                                                                                    >  1:447808/1‑149 (MQ=255)
                                              aGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATAATTGGCTGCATGTAt                                                                                                >  2:641961/1‑149 (MQ=255)
                                                 agagttGCGTACTAGTACTCAGTATAATTGCAAAGATGAATACAATAAAAAAAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTGAGAATGATTGGCTGTATTTATGCt                                                                                             <  2:316842/145‑1 (MQ=255)
                                                  aaTTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTg                                                                                            >  1:27314/1‑149 (MQ=255)
                                                      tGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTTTTTTTATTTTTTAAGCAAAAGTAGTTATAATGATTTTATGCATGTATGCTGAGaa                                                                                        >  2:563633/1‑149 (MQ=255)
                                                        cTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAAt                                                                                      >  2:271315/1‑149 (MQ=255)
                                                         ttACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTAGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAAtt                                                                                     >  1:466660/1‑149 (MQ=255)
                                                             cAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTTTTGTTTTTATTAAAGGAAAAGTAGTTAGAATAATTGGCTGGATGTTTGCTGAGAAATTTCtt                                                                                 >  2:678143/1‑149 (MQ=255)
                                                             cAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGATCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGAATGTATGCTGAGAAATTTAtt                                                                                 >  2:5057/1‑149 (MQ=255)
                                                                   tCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGAtgtg                                                                           >  1:252843/1‑149 (MQ=255)
                                                                    cAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGTAAAAGTAGTTAGAATGATTGGCTGCCTGTATGCTGAGAAATTTCTTGATgtgt                                                                          <  2:662844/149‑1 (MQ=255)
                                                                      gTATAATTGCAAAGATGAATATAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCt                                                                        <  2:232514/149‑1 (MQ=255)
                                                                      gTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCt                                                                        >  2:610471/1‑149 (MQ=255)
                                                                       tataATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCtt                                                                       >  1:215881/1‑149 (MQ=255)
                                                                        aaaaTTTCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTc                                                                      <  1:82277/147‑1 (MQ=255)
                                                                               cAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCa                                                               <  1:3793/149‑1 (MQ=255)
                                                                                aaaGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGAACAAATAAATAAAGCAAAAGTAGTTAAAATGATTGGCTGTATGTATGCTGCAAAATTTCTTGATTTGTCTCCTACATcgt                                                              >  2:8799/1‑147 (MQ=255)
                                                                                aaaGATGAATAAAATAAAAACAATAAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATCCACCATGAAAAATTTCTTGATGTGTCTTCAACATCAg                                                              <  1:545172/149‑1 (MQ=255)
                                                                                         tACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGTTGTGTCTTTAACATAAGCATCTCTTa                                                     >  1:113057/1‑149 (MQ=255)
                                                                                         tACAATAAAAACAAGAAAAATGAATGTTTTAATTATCTTTGTCTTAGTAAATTTTTTTGATCAAATAAATAAAGCAAACGTAGTTAGAATGATTGGTTGAACGTATGCTGAGATATGTCTTGATGTGTCTTCAACATCAGCATCTCTTa                                                     <  1:62893/149‑1 (MQ=255)
                                                                                          aCAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAATAAAAGACAATGTCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTAt                                                    >  1:538206/1‑149 (MQ=255)
                                                                                            aaTAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTAATAGTAAATTTTTTTGACGAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTCTGCTCACAAATTTCTTAATATGTTTTCAAAATAAGCTTCTCTTTTCa                                                  >  2:452716/1‑149 (MQ=255)
                                                                                              tAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCaaa                                                >  1:560365/1‑149 (MQ=255)
                                                                                               aaaaaCAAGTAAAGTAAAGTGTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAAt                                               <  2:480363/149‑1 (MQ=255)
                                                                                                  aaCAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCa                                            <  2:548480/149‑1 (MQ=255)
                                                                                                  aaCAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACAAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCa                                            <  1:378249/149‑1 (MQ=255)
                                                                                                   aCAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATGTTTTTGATTAAATAAATAAACCAAAAGTGGTTAGAATGAATGGCTGGATGTATGCTTAGAAATTTCTTGATGTTTCTTCAACATCACCTTATATTATCTTATTCat                                           >  2:553937/1‑149 (MQ=255)
                                                                                                    cAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCata                                          <  1:375847/149‑1 (MQ=255)
                                                                                                    cAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAAGAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTGTATGTGTCTTCAACCTCAGCATCTCTTATCAAATTCata                                          <  2:7390/149‑1 (MQ=255)
                                                                                                       gtaaagtcaaTGTTTTACTTACCTGGGTATTAGTAAATTTTTTTTTCTTATTATATAAAGAAAAAGGAGTTAGAATGATTGGCTGGATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAg                                       >  2:619356/1‑149 (MQ=255)
                                                                                                                   ttttACTTACCTTTGTATTAGTAACATTTTTTGACCAAATAAATAAAGCAAACGTAGTTAGAATGATTGGGTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCGAATTTATATAGCTAAGTAAGCTa                           <  1:678667/149‑1 (MQ=255)
                                                                                                                        cttacCTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTccc                      <  2:611817/149‑1 (MQ=255)
                                                                                                                             cTTTGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCGTGTTGGGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACt                 >  1:723761/1‑149 (MQ=255)
                                                                                                                              tttGTATTAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTg                <  2:358760/149‑1 (MQ=255)
                                                                                                                                tGTATTAGTAAATTTTTTTGACCAAAGAAATAAACCAAAAGTAGTTAGAAAGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGt              <  2:484365/149‑1 (MQ=255)
                                                                                                                                    ttAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAAGTATTGACTACATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGAATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGtt          >  1:726921/1‑149 (MQ=255)
                                                                                                                                     tAGTAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTa         <  1:288356/149‑1 (MQ=255)
                                                                                                                                        tacaTTTGTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACtt      <  1:194346/146‑1 (MQ=255)
                                                                                                                                        tAAATTTTTTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACtt      >  1:216176/1‑149 (MQ=255)
                                                                                                                                         aaaTTGTCTTGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTa     <  2:442193/149‑1 (MQ=255)
                                                                                                                                            ttttttTGACCAAGAAAGTAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGGTAGTCCAGAACTGGTAGTTACTTAAcc  <  1:715700/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
ACTCTATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTAGACCAAATAAATAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAACC  >  NC_000913_3_pae_tpiA/2106595‑2106883

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: