Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,488,331 C→G 7.6% I338M (ATC→ATG lptG → lipopolysaccharide transport system protein LptG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,488,3310CG 7.6% 135.8 / 6.3 53I338M (ATC→ATGlptGlipopolysaccharide transport system protein LptG
Reads supporting (aligned to +/- strand):  ref base C (30/19);  new base G (2/2);  total (32/21)
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.27e-02

GCGCTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGA  >  NC_000913/4488191‑4488472
                                                                                                                                            |                                                                                                                                             
gcgcTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATCTTCGGCACGCTGGCGTTGGTTTATGGCATCCCGCCGatcat                                                                                                                                             >  1:312924/1‑143 (MQ=255)
 ctcTGTCGTTCATCTTCCGCACACAGCGTAGCGTACCGACGGGCGCGCGTGTGGTCACCGGTATCAGTTTCGGTTTTGCCTTCTACGTACTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGatcatc                                                                                                                                            <  1:4394/141‑1 (MQ=255)
   cTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTATCAGTTTCGGTGTTGTCTTCTACGTACTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCCCCGATCATCgg                                                                                                                                          >  2:189567/1‑143 (MQ=255)
      tCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCGGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGgcgc                                                                                                                                       >  1:281556/1‑143 (MQ=255)
      tCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATGATCGgcgg                                                                                                                                       <  1:44594/143‑2 (MQ=255)
      tCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATGATCGgcgg                                                                                                                                       <  1:269652/143‑2 (MQ=255)
         ttCATCTTTGGCCCACTGAGTAGAGTACCGATGGGCGTGCGTGTGGTCACCGGTATCAGTTTCGGTTTTGTCTTATACGTACTGGACCAGATCATCGGCCCGCTGCCGTTGGTTTATTGCATCCCGCCCATCACCgaccgcca                                                                                                                                    >  2:444413/1‑136 (MQ=255)
              cTTTGGCCCACTGCGTAGCGTTCCGATGGGCGTGCGTGTGGTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCCCTGTTGc                                                                                                                               >  2:88267/1‑143 (MQ=255)
              cTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGCTCATCGGCGCACTGTTGc                                                                                                                               >  1:109389/1‑143 (MQ=255)
                   gCCCACTGCGTATCGTACCGTTGTGTGTGCTTGTGGTCACCTGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAgc                                                                                                                          <  1:204062/143‑1 (MQ=255)
                           cgtagcgtaCCGATGGGCGTGCGTGTGGTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGCTCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGctc                                                                                                                  >  2:164094/1‑142 (MQ=255)
                            gtagcgtaCCGATGGGCGTGCGTGTGGTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGTCCATATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGcttc                                                                                                                 <  1:102603/143‑1 (MQ=255)
                                cgtaCCGATGGGCGTGCGTGTGGTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCGAGCGCCAGcttcttct                                                                                                             >  1:228834/1‑143 (MQ=255)
                                   aCCGATGGGCGTGCGTGTGGTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTaa                                                                                                          >  1:245213/1‑143 (MQ=255)
                                      gATGGGCGTGCGTGTGGTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGACCATAGGCCGACTGTTCCCAAGCGCCAGCTCCTTCTTAATcc                                                                                                       >  1:194415/1‑142 (MQ=255)
                                       aTGGGCGTGCGTGTGGTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAg                                                                                                      >  1:359898/1‑143 (MQ=255)
                                       aTGGGCGTGCGGGTGGGCACCTGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAg                                                                                                      <  2:379129/143‑1 (MQ=255)
                                             gtgcgtgTGGTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCtgtg                                                                                                >  1:348076/1‑143 (MQ=255)
                                             gggcgtgTGGTCACCGGTATTAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCtgtg                                                                                                <  2:145299/141‑1 (MQ=255)
                                                    tgGTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGCTCATCGGCGCACTGTTGCCAAGGGCCAGCTTCTTCTTAATCAGCCTGTGGCTGtga                                                                                         >  2:439327/1‑141 (MQ=255)
                                                      gTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGGTGCCAAGCGCCAGCTTCTTCTTACTCAGCCTGTGGCGGTTGAt                                                                                       >  2:36926/1‑143 (MQ=255)
                                                      gTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATCTTCGGCCCGCTGACGTTGGCTTATGGCCTCCCGCCGATCATCGGCGCACAGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTaaa                                                                                       >  2:114571/1‑142 (MQ=255)
                                                         aCCGGTATCAGTTTCGGTTTTGTCTTCTACGTCCTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATgag                                                                                    >  1:198774/1‑143 (MQ=255)
                                                              tATCAGTTTCGGTTTTGTCTTCTACGTACTTGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAAt                                                                               <  1:391783/143‑1 (MQ=255)
                                                               aTCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATc                                                                              >  2:221260/1‑143 (MQ=255)
                                                                   gTTTCGGTTTTGTCTTCTACGTACTGGACCAGATCTTCGGCCCGCTGACGTTGGCTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTGAATGAGGAAATCGTaa                                                                          >  1:412720/1‑143 (MQ=255)
                                                                     ttCGGTTTTGTCTTCTACGTACTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAAcc                                                                        >  1:284399/1‑143 (MQ=255)
                                                                        ggTTTTGTCTTCTACGTACTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATAATACTTTGGCTGCTAAGGGGAAAATCGTAACCGaa                                                                     >  2:205181/1‑143 (MQ=255)
                                                                         gtgttGTCGTCTACGTACTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAAt                                                                    <  1:21462/140‑1 (MQ=255)
                                                                         gctttGTCTTCTGCGTCCTGGACCTGAGCTTCTGCCCGCTGGCGTTGGTCTATGCCATCCCGCCGATCATCGGCCCCCGGTTGTCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAGCCGAAt                                                                    <  2:189845/141‑1 (MQ=255)
                                                                         gTTTTGTCTTCTATGTACTGGCCCGGATCTTCGGCCCGCGGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGATTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAAt                                                                    <  2:464198/143‑1 (MQ=255)
                                                                         gTTTTGTCGTCTGCGTACTGTACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAAt                                                                    <  2:362286/143‑1 (MQ=255)
                                                                            gtggcttctACGTACTGGACCAGATCTTCGGCCCGCTGACGGTGGTGTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGaa                                                                 <  1:172625/139‑1 (MQ=255)
                                                                             tgccttcgACGTACTGGACCAGATCTTCGGCCCGCGGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAAc                                                                <  1:182035/140‑1 (MQ=255)
                                                                               tctgctAGGGACGGGACCAGATGTTCGGGCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGGTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGa                                                              <  1:66270/143‑1 (MQ=255)
                                                                                 ttctACGTACTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTCTGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTGATGAGAACATCGTAACCGAATGAACGGAt                                                            >  2:78634/1‑143 (MQ=255)
                                                                                 ttctACGTACTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAAt                                                            >  2:70082/1‑143 (MQ=255)
                                                                                           tGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAGTCGTAACCGAATGAACGAATAAAGGGAGCg                                                  >  1:65446/1‑143 (MQ=255)
                                                                                           tGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATCAAGGGAGCg                                                  >  1:243326/1‑143 (MQ=255)
                                                                                           tGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACGGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGGTAATGAGAAAATCGTAAACGAATGAACGAATAAAGGGAGCg                                                  >  2:275230/1‑143 (MQ=255)
                                                                                            ggACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAAGAGTAACCGAATGAACGAATAAAGGGAGCGa                                                 >  1:90456/1‑143 (MQ=255)
                                                                                              aCCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGaaa                                               >  1:48088/1‑143 (MQ=255)
                                                                                              aCCAGACCTGCGGACCGCGCACGTTGGTTCAGGGCATCCCGCCGATCACCGGCGCACTTTTGCCAAGCGCCACCTTCGTCTTAATGAGCATGTAGCTGGTAAGGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGaaa                                               <  2:33851/143‑1 (MQ=255)
                                                                                                  gATCTTCGGCCCGCTGACGTTGGTTTATGGCAGCCCGCCGATCAGCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGATc                                           <  1:327566/143‑1 (MQ=255)
                                                                                                           cccGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGGTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCTTTCGCTCCCTttattc                                  <  1:205181/143‑2 (MQ=255)
                                                                                                             cGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGTTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTattatt                                <  2:281556/143‑1 (MQ=255)
                                                                                                                tGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGtt                             <  1:11918/143‑1 (MQ=255)
                                                                                                                   cGTTGGTTTATGGCATCCCGCCGATCATCGGCACACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCCCCCTTTATTCTCGTTAGc                          >  2:199139/1‑143 (MQ=255)
                                                                                                                       ggTTTATGGCCTCCCGCCGCTCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAg                      <  2:312924/143‑1 (MQ=255)
                                                                                                                               gCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTc              >  2:410093/1‑143 (MQ=255)
                                                                                                                                 agtgcgccgATGATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGACCGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGt            >  1:124958/5‑143 (MQ=255)
                                                                                                                                  gcccGCCGATCATAGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGtt           <  1:197816/142‑1 (MQ=255)
                                                                                                                                   tgcgccgATGATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATAATTGTTAGCAAAGGGTGCTTCGTTc          >  2:402135/3‑143 (MQ=255)
                                                                                                                                     cgccgATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAGCGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCCAAGTGTGCTTCGTTcct        >  1:56576/1‑141 (MQ=255)
                                                                                                                                     cgccgATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTCTCATTGTTCGCACTCTCTGCTTCGTTcct        >  1:56577/1‑141 (MQ=255)
                                                                                                                                           tcatcGGCGCACTGTTGCCAATCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTATCCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGa  <  2:271150/143‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                             
GCGCTGTCGTTCATCTTTGGCCCACTGCGTAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGA  >  NC_000913/4488191‑4488472

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