Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,887,056 T→G 7.4% intergenic (‑41/‑89) PP_5161 ← / → PP_5162 lysophospholipase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,887,0560TG 7.4% 169.2 / 4.1 54intergenic (‑41/‑89)PP_5161/PP_5162lysophospholipase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (23/27);  new base G (2/2);  total (25/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GGGGTGAAAGGCAGGGGGCATCGGGTACATCCACATTGATGCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCTGTGCACGGTCATCGCCTTGTTGGTCGGCGCC  >  NC_002947/5886920‑5887201
                                                                                                                                        |                                                                                                                                                 
ggggTGAAAGGCAGCGGGCATCGGGGACATCCACATTGATGCAAGTATTCATCTGGAATCTTCACCTCTGGGCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGTAGCGGCCGTTCGCGACGCAAgg                                                                                                                                       <  4:179755/149‑1 (MQ=255)
  ggTGTAAGGCAGGGGGCTGCGGGGACATCCACATTGATGCCAGTATTGCTCGGGAATCTTCAGCTCTGGCCGGCTCCATGGCAAGCTTGCTCACCTCGAACAACCAGCCTCGAAGTAATTGCGGGAGCGGCCTTGCGGCGCGAAAGGcc                                                                                                                                     <  4:6057/149‑1 (MQ=255)
         ggCAGGGGGCATCGGGTACATCCACATTGATGCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCca                                                                                                                              <  3:61882/149‑1 (MQ=255)
         ggCAGGGGGCAGCGGGTACATCCACATTGATGCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCca                                                                                                                              <  7:144611/149‑1 (MQ=255)
           cAGGGGGCAGCGGGTACATCCACATTGATGCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCcaca                                                                                                                            <  5:219654/149‑1 (MQ=255)
            aGGGGGCATCGGGTACATCCACATTGATGCAAGTATTGATCTGGAAGCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAg                                                                                                                           <  5:179298/149‑1 (MQ=255)
             gggggCATCGGGTACATCCACATTGATGCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAgg                                                                                                                          <  5:255080/149‑1 (MQ=255)
                  cATCGGGTACATCCACATTGATGCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCg                                                                                                                     <  8:224923/149‑1 (MQ=255)
                        gTACATCCACATTGATGCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTCCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCa                                                                                                               <  3:69998/149‑1 (MQ=255)
                             tCCACATTGATGCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCttt                                                                                                          >  6:142883/1‑149 (MQ=255)
                                acaTTGATGCTAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAg                                                                                                       <  5:77926/149‑1 (MQ=255)
                                   ttGATGCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGaaa                                                                                                    <  2:114707/149‑1 (MQ=255)
                                      aTGCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAAtt                                                                                                  <  6:4854/148‑1 (MQ=255)
                                      aTGCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAAtt                                                                                                  >  3:21880/1‑148 (MQ=255)
                                        gCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGaaa                                                                                                    >  7:236362/1‑144 (MQ=255)
                                        gCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGaaa                                                                                                    <  8:236362/144‑1 (MQ=255)
                                        gCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAg                                                                                               >  1:245443/1‑149 (MQ=255)
                                          aaGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGc                                                                                              <  3:26695/148‑1 (MQ=255)
                                          aaGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACAAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCTGCTTTTAGAAATTGAGc                                                                                              >  1:137416/1‑148 (MQ=255)
                                           aGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCaa                                                                                            <  4:21880/149‑1 (MQ=255)
                                            gTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGc                                                                                              >  5:258003/1‑146 (MQ=255)
                                              aTTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGa                                                                                          >  6:148944/1‑148 (MQ=255)
                                                tGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGgc                                                                                                                   >  1:68909/1‑121 (MQ=255)
                                                tGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGgc                                                                                                                   <  2:68909/121‑1 (MQ=255)
                                                 gATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAg                                                                                       <  1:160121/148‑1 (MQ=255)
                                                     tGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCt                                                                                  >  3:116451/1‑149 (MQ=255)
                                                      ggAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCtt                                                                                 >  4:142039/1‑149 (MQ=255)
                                                       gAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAAGTGAGCAAGACAGTTGCtgc                                                                                >  4:98421/1‑147 (MQ=255)
                                                             tCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAg                                                                          >  4:96567/1‑149 (MQ=255)
                                                               aGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAgg                                                                         <  3:142039/148‑1 (MQ=255)
                                                               aGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAgg                                                                         <  2:245443/148‑1 (MQ=255)
                                                                 ctctGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTccc                                                                                                                               <  7:127368/92‑1 (MQ=255)
                                                                 ctctGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTccc                                                                                                                               >  8:127368/1‑92 (MQ=255)
                                                                    tGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCAcc                                                                   >  3:203515/1‑149 (MQ=255)
                                                                    tGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCAcc                                                                   >  1:11160/1‑149 (MQ=255)
                                                                           cgcCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCtt                                                            <  3:32767/149‑1 (MQ=255)
                                                                            gcCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTc                                                           <  5:148944/149‑1 (MQ=255)
                                                                            gcCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTc                                                           <  5:142883/149‑1 (MQ=255)
                                                                                                   gAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGa                                                                                                      <  2:270741/83‑1 (MQ=255)
                                                                                                   gAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGa                                                                                                      >  1:270741/1‑83 (MQ=255)
                                                                                                   gAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAAc                                     <  4:116451/148‑1 (MQ=255)
                                                                                                    aacaacCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAAc                                     <  1:115675/147‑1 (MQ=255)
                                                                                                    aacaacCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAAc                                     >  2:115675/1‑147 (MQ=255)
                                                                                                          cAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCtgtg                              >  6:257217/1‑148 (MQ=255)
                                                                                                                cGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTATTAGAAATTGCGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGATTAGCCGAATCTTCAAACCGCTGTGCACGGt                        >  8:142189/1‑148 (MQ=255)
                                                                                                                  aaCTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTc                                                                                <  7:142072/90‑1 (MQ=255)
                                                                                                                  aaCTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTc                                                                                >  8:142072/1‑90 (MQ=255)
                                                                                                                            gggAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCTGTGCACGGTCATCGCCttgttg           <  1:113950/149‑1 (MQ=255)
                                                                                                                              gAGCGGCCTTGCGTCGCGAAAGGCCGCTCCCGCAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCTGTGCACGGTCATCGCCttgttg           >  6:47984/1‑147 (MQ=38)
                                                                                                                              gAGCGGCCTTGCGTCGCGAAAGGCCGCTCCCGCAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCTGTGCACGGTCATCGCCttgttg           <  5:47984/147‑1 (MQ=38)
                                                                                                                                 cGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCTGTGCACGGTCATCGCCTTGTTGGTCgg      >  3:174577/1‑149 (MQ=255)
                                                                                                                                 cGGCCTTGCGTCGCGAAAGGCCGCTCCCGCAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCTGTGCACGGTCATCGCCTTGTTGGTCgg      >  3:51230/1‑149 (MQ=38)
                                                                                                                                  ggCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCTGTGCACGGTCATCGCCTTGTTGGTCGgc     >  7:230467/1‑149 (MQ=255)
                                                                                                                                     cTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCTGTGCACGGTCATCGCCTTGTTGGTCGGCGcc  >  4:103779/1‑149 (MQ=255)
                                                                                                                                        |                                                                                                                                                 
GGGGTGAAAGGCAGGGGGCATCGGGTACATCCACATTGATGCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCTGTGCACGGTCATCGCCTTGTTGGTCGGCGCC  >  NC_002947/5886920‑5887201

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