Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,043,495 G→T 7.5% P239H (CCC→CAC)  etk ← tyrosine‑protein kinase, role in O‑antigen capsule formation

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,043,4950GT 7.5% 178.2 / 5.5 53P239H (CCC→CAC) etktyrosine‑protein kinase, role in O‑antigen capsule formation
Reads supporting (aligned to +/- strand):  ref base G (24/25);  new base T (2/2);  total (26/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.92e-01

TTCTTCCGCTTGGTCCAGCTCGCTGCGCACTTCAGGTAACTGGCGCTGTAAGAATTCAAGGCTTTGTGAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCTGGCTCAGGACAAACTGTGTTCCTGGTTTGGCCTTAATGTCCGC  >  NC_000913/1043347‑1043643
                                                                                                                                                    |                                                                                                                                                    
ttcttcCGCTTGGTCCAGCTCGCTGCGCACTTCAGGTAACTGGCGCTGTAAGAATTCAAGGCTTTGTGAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTggg                                                                                                                                                      <  1:356294/149‑1 (MQ=255)
ttcttcCGCTTGGTCCAGCTCGCTGCGCACTTCAGGTAACTGGCGCTGTAAGAATTCAAGGCTTTGTGAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTggg                                                                                                                                                      <  1:357615/149‑1 (MQ=255)
  cttcCGCTTGGTCCAGCTCGCTGCGCACTTCAGGTAACTGGCGCTGTAAGAATTCAAGGCTTTGTGAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGa                                                                                                                                                    >  1:412646/1‑149 (MQ=255)
   ttcCGCTTGGTCCAGCTCGCTGCGCACTTCAGGTAACTGGCGCTGTAAGAATTCAAGGCTTTGTGAATCCTGCGTCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGat                                                                                                                                                   >  4:243621/1‑149 (MQ=255)
   ttcCGCTTGGTCCAGCTCGCTGCGCACTTCAGGTAACTGGCGCTGTAAGAATTCAAGGCTTTGTGAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGat                                                                                                                                                   <  3:420675/149‑1 (MQ=255)
    tcCGCTTGGTCCAGCTCGCTGCGCACTTCAGGTAACTGGCGCTGTAAGAATTCAAGGCTTTGTGAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGatc                                                                                                                                                  <  3:105964/149‑1 (MQ=255)
        cTTGGTCCAGCTCGCTGCGCACTTCAGGTAACTGGCGCTGTAAGAATTCAAGGCTTTGTGAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGAtcatca                                                                                                                                              >  4:91997/1‑149 (MQ=255)
                aGCTCGCTGCGCACTTCAGGTAACTGGCGCTGTAAGAATTCAAGGCTTTGTGAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCAt                                                                                                                                      <  2:265544/149‑1 (MQ=255)
                aGCTCGCTGCGCACTTCAGGTAACTGGCGCTGTAAGAATTCAAGGCTTTGTGAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCAt                                                                                                                                      >  1:454680/1‑149 (MQ=255)
                 gCTCGCTGCGCACTTCAGGTAACTGGCGCTGTAAGAATTCAAGGCTTTGTGAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATg                                                                                                                                     <  2:133595/149‑1 (MQ=255)
                   tCGCTGCGCACTTCAGGTAACTGGCGCTGTAAGAATTCAAGGCTTTGTGAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGt                                                                                                                                   >  4:172723/1‑149 (MQ=255)
                   tCGCTGCGCACTTCAGGTAACTGGCGCTGTAAGAATTCAAGGCTTTGTGAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGt                                                                                                                                   >  2:19915/1‑149 (MQ=255)
                                       cTGGCGCTGTAAGAATTCAAGGCTTTGTGAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGcttt                                                                                                               <  3:91997/149‑1 (MQ=255)
                                            gcTGTAAGAATTCAAGGCTTTGTGAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTa                                                                                                          <  3:45238/149‑1 (MQ=255)
                                            gcTGTAAGAATTCAAGGCTTTGTGAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTa                                                                                                          >  3:37451/1‑149 (MQ=255)
                                                tAAGAATTCAAGGCTTTGTGAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGc                                                                                                      >  1:115371/1‑149 (MQ=255)
                                                  aGAATTCAAGGCTTTGTGAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGt                                                                                                    >  2:373346/1‑149 (MQ=255)
                                                    aaTTCAAGGCTTTGTGAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTc                                                                                                  <  4:161434/149‑1 (MQ=255)
                                                                tgtgAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCGACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTaaa                                                                                      <  2:187365/149‑1 (MQ=255)
                                                                tgtgAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTaaa                                                                                      <  3:388840/149‑1 (MQ=255)
                                                                 gtgAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTc                                                                                                  <  4:220673/136‑1 (MQ=255)
                                                                 gtgAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTc                                                                                                  >  3:220673/1‑136 (MQ=255)
                                                                   gAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGt                                                                                                                               <  2:401975/105‑1 (MQ=255)
                                                                   gAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGt                                                                                                                               >  1:401975/1‑105 (MQ=255)
                                                                     aTCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGttt                                                                                 >  1:363279/1‑149 (MQ=255)
                                                                      tCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTc                                                                                >  2:308100/1‑149 (MQ=255)
                                                                           cgccgcCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCa                                                                           >  2:432165/1‑149 (MQ=255)
                                                                            gccgccTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCaa                                                                          <  2:115371/149‑1 (MQ=255)
                                                                               gccTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGTGATGATCCCCAGtt                                                                                                                                        >  2:34856/1‑83 (MQ=255)
                                                                               gccTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGTGATGATCCCCAGtt                                                                                                                                        <  1:34856/84‑2 (MQ=255)
                                                                                  tGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGtt                                                                    >  3:438364/1‑149 (MQ=255)
                                                                                     cgagcgaTATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAAt                                                                 <  4:443797/149‑1 (MQ=255)
                                                                                                  gTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCg                                                    <  1:19915/149‑1 (MQ=255)
                                                                                                   ttGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCgg                                                   <  2:454680/149‑1 (MQ=255)
                                                                                                         aTAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCt                                             >  2:413326/1‑149 (MQ=255)
                                                                                                                ttAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCTGGCTCAg                                      >  1:114832/1‑149 (MQ=255)
                                                                                                                  aGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCTGGCTCAGGa                                    <  1:308100/149‑1 (MQ=255)
                                                                                                                     gATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCTGGCTCAGGACaa                                 >  1:342394/1‑149 (MQ=255)
                                                                                                                       tGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCTGGCTCAGGACAAAc                               <  2:443736/149‑1 (MQ=255)
                                                                                                                        gCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCTGGCTCAGGACAAACt                              >  3:437222/1‑149 (MQ=255)
                                                                                                                             tCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGct                                                                                                                 >  1:365720/1‑61 (MQ=255)
                                                                                                                             tCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGct                                                                                                                 <  2:365720/61‑1 (MQ=255)
                                                                                                                                gAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCTGGCTCAGGACAAACTGTGTTCCt                      <  4:92410/149‑1 (MQ=255)
                                                                                                                                 aaTACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCTGGCTCAGGACAAACTGTGTTCCTg                     >  2:16843/1‑149 (MQ=255)
                                                                                                                                  aTACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCTGGCTCAGGACAAACTGTGTTCCTgg                    <  2:342394/149‑1 (MQ=255)
                                                                                                                                     cGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCTGGCTCAGGACAAACTGTGTTCCTGGttt                 <  1:432165/149‑1 (MQ=255)
                                                                                                                                      gAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTc                                                           <  4:123926/106‑1 (MQ=255)
                                                                                                                                      gAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTc                                                           >  3:123926/1‑106 (MQ=255)
                                                                                                                                        gTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCTGGCTCAGGACAAACTGTGTTCCTGGTTTGGc              <  3:243621/149‑1 (MQ=255)
                                                                                                                                          ccatgaCTGGTGATGATCCCCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCTGGCTCAGGACAAACTGTGTTCCTGGTTTGGCCt            >  4:411779/6‑149 (MQ=255)
                                                                                                                                           catgaCTGGTGATGATCCCCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCTGGCTCAGGACAAACTGTGTTCCTGGTTTGGCCtt           <  3:411779/145‑1 (MQ=255)
                                                                                                                                            tCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCTGGCTCAGGACAAACTGTGTTCCTGGTTTGGCCTTa          >  4:474527/1‑149 (MQ=255)
                                                                                                                                             cAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCTGGCTCAGGACAAACTGTGTTCCTGGTTTGGCCTTaa         <  2:114832/149‑1 (MQ=255)
                                                                                                                                              aaCTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCTGGCTCAGGACAAACTGTGTTCCTGGTTTGGCCTTAAt        >  4:334490/1‑149 (MQ=255)
                                                                                                                                              aaCTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCTGGCTCAGGACAAACTGTGTTCCTGGTTTGGCCTTAAt        <  4:438364/149‑1 (MQ=255)
                                                                                                                                                  ggggATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCTGGCTCAGGACAAACTGTGTTCCTGGTTTGGCCTTAATGTcc    >  1:390795/1‑149 (MQ=255)
                                                                                                                                                   gggATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCt                                                                                               <  1:275689/57‑1 (MQ=255)
                                                                                                                                                   gggATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCt                                                                                               >  2:275689/1‑57 (MQ=255)
                                                                                                                                                    ggATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCTGGCTCAGGACAAACTGTGTTCCTGGTTTGGCCTTAATGTCcgc  <  4:37451/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTCTTCCGCTTGGTCCAGCTCGCTGCGCACTTCAGGTAACTGGCGCTGTAAGAATTCAAGGCTTTGTGAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCTGGCTCAGGACAAACTGTGTTCCTGGTTTGGCCTTAATGTCCGC  >  NC_000913/1043347‑1043643

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑