Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,148,681 C→T G529S (GGC→AGC)  eptC ← LPS heptose I phosphoethanolamine transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,148,6810CT98.5% 208.1 / NA 65G529S (GGC→AGC) eptCLPS heptose I phosphoethanolamine transferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (31/33);  minor base A (1/0);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCCCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTAC  >  NC_000913/4148548‑4148808
                                                                                                                                     |                                                                                                                               
tGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCa                                                                                                                            >  1:277975/1‑139 (MQ=255)
  aTCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCaaa                                                                                                                          >  2:427141/1‑139 (MQ=255)
        aTAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGAc                                                                                                                    <  2:597677/139‑1 (MQ=255)
            ggCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAgg                                                                                                                >  2:545196/1‑139 (MQ=255)
                 tgtgtCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTgg                                                                                                           >  1:670054/1‑139 (MQ=255)
                 tgtgtCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTgg                                                                                                           >  2:88272/1‑139 (MQ=255)
                 tgtgtCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTgg                                                                                                           >  1:616970/1‑139 (MQ=255)
                 tgtgtCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTgg                                                                                                           >  2:716512/1‑139 (MQ=255)
                 tgtgtCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTgg                                                                                                           >  1:766050/1‑139 (MQ=255)
                   tgtCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGAt                                                                                                         >  2:58450/1‑139 (MQ=255)
                   tgtCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCACGCCAAATCTGACCAGGTGTGGAt                                                                                                         >  1:430917/1‑139 (MQ=255)
                           tCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCg                                                                                                 <  2:430917/139‑1 (MQ=255)
                            cGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGc                                                                                                <  1:716512/139‑1 (MQ=255)
                               tCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAg                                                                                             >  1:372913/1‑139 (MQ=255)
                                cAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAgg                                                                                            <  2:670054/139‑1 (MQ=255)
                                cAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAgg                                                                                            <  1:749441/139‑1 (MQ=255)
                                cAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAgg                                                                                            <  2:766050/139‑1 (MQ=255)
                                 aGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGc                                                                                           <  1:595904/139‑1 (MQ=255)
                                 aGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGc                                                                                           <  2:277975/139‑1 (MQ=255)
                                 aGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGc                                                                                           <  2:616970/139‑1 (MQ=255)
                                    gCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGAcc                                                                                                                   <  2:730728/112‑1 (MQ=255)
                                    gCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGAcc                                                                                                                   >  1:730728/1‑112 (MQ=255)
                                          ttCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTAc                                                                                  <  1:156919/139‑1 (MQ=255)
                                             ttATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGc                                                                                                                              >  1:26818/1‑92 (MQ=255)
                                             ttATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGc                                                                                                                              <  2:26818/92‑1 (MQ=255)
                                                  cGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTAc                                                                                  <  2:205305/131‑1 (MQ=255)
                                                  cGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTAc                                                                                  >  1:205305/1‑131 (MQ=255)
                                                     gTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTAt                                                                                      <  2:134024/124‑1 (MQ=255)
                                                     gTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTAt                                                                                      >  1:134024/1‑124 (MQ=255)
                                                     gTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATcc                                                                       >  2:170966/1‑139 (MQ=255)
                                                       tACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTg                                                                     >  2:525267/1‑139 (MQ=255)
                                                        aCCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGc                                                                    <  1:58450/139‑1 (MQ=255)
                                                               cAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAAc                                                                           >  1:424858/1‑125 (MQ=255)
                                                               cAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAAc                                                                           <  2:424858/125‑1 (MQ=255)
                                                                  cGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCAc                                                          >  1:555156/1‑139 (MQ=255)
                                                                   gggTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACg                                                         >  2:137009/1‑139 (MQ=255)
                                                                    ggTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACgg                                                        >  1:78019/1‑139 (MQ=255)
                                                                       aGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGa                                                     >  2:293358/1‑139 (MQ=255)
                                                                           tCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAg                                                                                             <  2:169040/95‑1 (MQ=255)
                                                                           tCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAg                                                                                             >  1:169040/1‑95 (MQ=255)
                                                                                      cGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGAcc                                                                                                                   >  1:353055/1‑62 (MQ=255)
                                                                                      cGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGAcc                                                                                                                   <  2:353055/62‑1 (MQ=255)
                                                                                      cGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCAt                                      >  2:392550/1‑139 (MQ=255)
                                                                                                   gAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAg                                                                                             >  1:218569/1‑71 (MQ=255)
                                                                                                   gAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAg                                                                                             <  2:218569/71‑1 (MQ=255)
                                                                                                    aaCGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTAc                                                                                  <  2:230396/81‑1 (MQ=255)
                                                                                                    aaCGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTAc                                                                                  >  1:230396/1‑81 (MQ=255)
                                                                                                          ttGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCAcagcag                  <  2:372913/139‑1 (MQ=255)
                                                                                                            ggATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGa                <  1:137009/139‑1 (MQ=255)
                                                                                                            ggATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGa                >  1:236106/1‑139 (MQ=255)
                                                                                                             gATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGaa               <  1:293358/139‑1 (MQ=255)
                                                                                                               tCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACg             <  2:555156/139‑1 (MQ=255)
                                                                                                                 gTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGa           <  1:88272/139‑1 (MQ=255)
                                                                                                                 gTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGa           <  1:211768/139‑1 (MQ=255)
                                                                                                                  tAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGAt                                                                                                         >  2:255558/1‑44 (MQ=255)
                                                                                                                  tAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGAt                                                                                                         <  1:255558/44‑1 (MQ=255)
                                                                                                                  tAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGaa          <  1:427141/139‑1 (MQ=255)
                                                                                                                  tAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGaa          <  2:78019/139‑1 (MQ=255)
                                                                                                                  tAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGaa          <  1:170966/139‑1 (MQ=255)
                                                                                                                   aaCCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAAt         >  2:37890/1‑139 (MQ=255)
                                                                                                                     ccGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATgg       <  1:706586/139‑1 (MQ=255)
                                                                                                                       gtcgtAAGATAAGCTCGCCAAATCTGACCATGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGgtg     <  1:525267/139‑1 (MQ=255)
                                                                                                                          gtAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTAc  >  2:189128/1‑139 (MQ=255)
                                                                                                                             aGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGaa          >  2:27069/1‑128 (MQ=255)
                                                                                                                             aGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGaa          <  1:27069/128‑1 (MQ=255)
                                                                                                                                     |                                                                                                                               
TGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCCCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTAC  >  NC_000913/4148548‑4148808

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