Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,007,718 2 bp→CG coding (943‑944/1023 nt) pldB → lysophospholipase L(2)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,007,7180GC100.0% 267.4 / NA 96A315R (GCC→CGC) pldBlysophospholipase L(2)
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (50/46);  total (50/46)
*NC_0009134,007,7190CG99.0% 259.1 / NA 96A315R (GCC→CGC) pldBlysophospholipase L(2)
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base G (50/45);  total (50/46)
Fisher's exact test for biased strand distribution p-value = 4.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGGCCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGT  >  NC_000913/4007651‑4007794
                                                                   ||                                                                           
tgAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTc                                                                      >  2:3919248/2‑76 (MQ=255)
cggaggaggTCGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTc                                                                      <  2:5197996/73‑1 (MQ=255)
cGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTc                                                                      <  1:661491/76‑1 (MQ=255)
cGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTc                                                                      >  1:4392383/1‑76 (MQ=255)
 gAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTcg                                                                     >  1:301898/1‑76 (MQ=255)
      aggaCGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTc                                                                >  1:5410929/1‑76 (MQ=255)
      aggaCGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTc                                                                >  1:371392/1‑76 (MQ=255)
       ggaCGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTcc                                                               >  1:1257774/1‑76 (MQ=255)
         gcGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGACAAGGACGCAATGCGCTCAGTCGCGCTCCAc                                                             <  2:1680576/75‑1 (MQ=255)
         aCGGTCTTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCAc                                                             <  1:1794393/76‑1 (MQ=255)
         aCGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCAc                                                             <  1:10018/76‑1 (MQ=255)
         aCGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCAc                                                             <  1:2264881/76‑1 (MQ=255)
          cGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTccccc                                                            >  2:703273/1‑73 (MQ=255)
          cGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACg                                                            <  2:1719764/76‑1 (MQ=255)
           ggCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGc                                                           >  2:596363/1‑76 (MQ=255)
           ggCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGc                                                           >  1:2260655/1‑76 (MQ=255)
           ggCCGTTGGGAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGGCGCGCTCCACGc                                                           >  2:648903/1‑76 (MQ=255)
            gCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCTCGcc                                                          <  1:378024/76‑1 (MQ=255)
            gCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGcc                                                          <  1:5481396/76‑1 (MQ=255)
            gCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGcc                                                          >  2:959315/1‑76 (MQ=255)
            gCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGcc                                                          >  2:2600940/1‑76 (MQ=255)
             ccGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCa                                                         >  2:1190628/1‑76 (MQ=255)
              cGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCAt                                                        <  2:2424693/76‑1 (MQ=255)
              cGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCAt                                                        <  2:5521197/76‑1 (MQ=255)
               gTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATc                                                       >  2:89918/1‑76 (MQ=255)
                ttGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCg                                                      >  2:4645499/1‑76 (MQ=255)
                 tGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTGCACGCCATCGt                                                     >  1:1281978/1‑76 (MQ=255)
                 tGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGt                                                     <  2:1137453/76‑1 (MQ=255)
                  ggTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGtt                                                    >  2:1172553/1‑76 (MQ=255)
                  ggTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGtt                                                    >  2:5108347/1‑76 (MQ=255)
                  ggTAATTAAAGGTGATTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGtt                                                    <  1:658318/76‑1 (MQ=255)
                    tAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTc                                                                      <  2:738922/56‑1 (MQ=255)
                    tAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTc                                                                      >  1:738922/1‑56 (MQ=255)
                    tAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGa                                                  >  1:5774644/1‑76 (MQ=255)
                     aaTTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGAt                                                 <  1:3163652/76‑1 (MQ=255)
                       ttAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGAttt                                               <  2:539075/76‑1 (MQ=255)
                       ttAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGAttt                                               >  1:668365/1‑76 (MQ=255)
                        tAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGCtttt                                              >  2:4691495/1‑76 (MQ=255)
                         aaaGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGAttttt                                             >  1:2553962/1‑76 (MQ=255)
                         aaaGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGAttttt                                             <  2:1106909/76‑1 (MQ=255)
                         aaaGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGAttttt                                             <  1:328551/76‑1 (MQ=255)
                         aaaGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGAttttt                                             >  2:3462193/1‑76 (MQ=255)
                          aaGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGAtttttt                                            >  2:3904893/1‑76 (MQ=255)
                          aaGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGAtttttt                                            <  2:5389705/76‑1 (MQ=255)
                              tGCTCACCATGAGATACTTTTTGAAAAGGACGCAATGCGCTCAGGCGCGCTCCACGCCCTCGCTGATTTTTTCAAc                                        >  2:3879512/1‑76 (MQ=255)
                               gCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACa                                       >  2:1672332/1‑76 (MQ=255)
                                cTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAg                                      >  2:2257741/1‑76 (MQ=255)
                                  tACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGc                                    >  1:5257996/1‑76 (MQ=255)
                                  tACCAGGAGATGCTTTTTGAAAAGGACGCAATGCCGTCAGTCGCGATCCAGGCCATCGTTTATTTTTTCAACGGGc                                    <  2:5508053/76‑1 (MQ=255)
                                   aCCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCa                                   >  2:2541457/1‑76 (MQ=255)
                                    ccATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCTCTCCACGCCATCGTTGATTTTTTCAACAGGCAt                                  <  1:1146936/76‑1 (MQ=255)
                                    ccATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGCTTTTTTCAACAGGCAt                                  >  2:1013165/1‑76 (MQ=255)
                                    ccATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGCTTTTTTCAACACGCAt                                  >  2:77726/1‑76 (MQ=255)
                                    ccATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCAt                                  >  1:3587812/1‑76 (MQ=255)
                                    ccATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCAt                                  >  1:1531715/1‑76 (MQ=255)
                                     cATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGCTTTTTTCAACACGCaaa                                 >  2:5616634/1‑74 (MQ=255)
                                     cATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATa                                 >  2:5613864/1‑76 (MQ=255)
                                     aaTGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGACGCGCTCCACGCCATCGTTGCTTTTGGCCACAGGCATa                                 >  2:784037/2‑76 (MQ=255)
                                      aTGAGATTCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATaa                                >  1:3105532/1‑76 (MQ=255)
                                      aTGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATaa                                <  2:1653888/76‑1 (MQ=255)
                                      aTGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATaa                                <  1:5781882/76‑1 (MQ=255)
                                      aTGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATaa                                >  1:4491959/1‑76 (MQ=255)
                                       tGAGCTCCTTTTTGGAAGTGACGCCATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAAc                               <  2:5792168/76‑1 (MQ=255)
                                       tGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCGACAGGCATAAc                               <  1:5578277/76‑1 (MQ=255)
                                         agaTCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTc                             <  2:668365/76‑1 (MQ=255)
                                          gaTCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTGCACGCCATCGTTGATTTTTTCAACAGGCATAACTCa                            <  2:1281978/76‑1 (MQ=255)
                                           aTCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCAc                           >  1:964578/1‑76 (MQ=255)
                                             ccTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCAccc                         <  2:5410928/76‑1 (MQ=255)
                                             ccTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCAccc                         <  1:3030054/76‑1 (MQ=255)
                                             ccTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCAccc                         <  2:505254/76‑1 (MQ=255)
                                              cTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCa                        >  2:3612559/1‑76 (MQ=255)
                                              cTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCa                        <  1:2291532/76‑1 (MQ=255)
                                              cTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCa                        <  1:3919249/76‑1 (MQ=255)
                                              cTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCa                        >  1:642227/1‑76 (MQ=255)
                                              cTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCa                        <  1:2584588/76‑1 (MQ=255)
                                                ttttGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGc                      <  1:4645500/76‑1 (MQ=255)
                                                ttttGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGc                      <  1:5222815/76‑1 (MQ=255)
                                                  tggAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCgg                    <  1:821965/74‑1 (MQ=255)
                                                   tGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGa                   <  2:5356669/76‑1 (MQ=255)
                                                    gAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGaa                  <  2:4918889/76‑1 (MQ=255)
                                                    gAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGaa                  <  1:1341296/76‑1 (MQ=255)
                                                     aaaaGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGaaa                 >  2:5025725/1‑76 (MQ=255)
                                                     aaaaGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGaaa                 >  1:4725291/1‑76 (MQ=255)
                                                     aaaaGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGaaa                 <  2:4392382/76‑1 (MQ=255)
                                                     aaaaGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGaaa                 <  1:690427/76‑1 (MQ=255)
                                                      aaaGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATATGTCACCCAGCGGAAAc                <  2:3830582/76‑1 (MQ=255)
                                                      aaaGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAAc                >  1:4540477/1‑76 (MQ=255)
                                                      aaaGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAAc                >  2:3245886/1‑76 (MQ=255)
                                                       aaGGACGCAATGCGCTGAGTCGCGCTCCACGCCGTCGTTGATTTTTTCACGAGGCATATGTCGCCCAGCGGaaccc               <  2:4371146/76‑4 (MQ=255)
                                                       aaGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAAcc               >  1:4779935/1‑76 (MQ=255)
                                                       aaGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAAca               <  2:932229/76‑2 (MQ=255)
                                                        aGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCg              <  1:2541457/76‑1 (MQ=255)
                                                          gACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGct            >  2:4377811/1‑76 (MQ=255)
                                                           aCGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGctc           >  2:2985346/1‑76 (MQ=255)
                                                           aCGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGctc           <  1:3338912/76‑1 (MQ=255)
                                                            cGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGctct          >  1:5553057/1‑76 (MQ=255)
                                                                    ccTCAGTCGCGCTCCCCGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGt  <  1:82120/76‑1 (MQ=255)
                                                                   ||                                                                           
CGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGGCCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTACAGAGGT  >  NC_000913/4007651‑4007794

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