Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,338,493 A→T 9.7% intergenic (+113/+29) PP_3812 → / ← PP_5632 lipase/allo‑isoleucine metabolism protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,338,4930AT 9.7% 113.2 / 2.4 51intergenic (+113/+29)PP_3812/PP_5632lipase/allo‑isoleucine metabolism protein
Reads supporting (aligned to +/- strand):  ref base A (20/26);  new base T (3/2);  total (23/28)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTT‑CGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGAAT  >  NC_002947/4338351‑4338637
                                                                                                                                               |                                                                                                                                                
ggCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTT‑CGCGGGTAAACCcg                                                                                                                                            <  5:4188/149‑1 (MQ=255)
     cGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTT‑CGCGGGTTTACCCGCGaaa                                                                                                                                       <  2:292001/149‑1 (MQ=255)
     cGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTT‑CGCGGGTAAACCCGCGaaa                                                                                                                                       <  5:368472/149‑1 (MQ=255)
     cGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTT‑CGCGGGTAAACCCGCGaaa                                                                                                                                       <  8:274723/149‑1 (MQ=255)
      gAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTT‑CGCGGGTTTACCCGCGAAAg                                                                                                                                      <  3:133898/149‑1 (MQ=255)
        aCTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTT‑CGCGGGTAAACCCGCGAAAgg                                                                                                                                     >  4:219389/1‑148 (MQ=255)
        aCTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTT‑CGCGGGTAAACCCGCGAAAgg                                                                                                                                     <  3:219389/148‑1 (MQ=255)
          ttAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTT‑CGCGGGTAAACCCGCGAAAGGGcc                                                                                                                                  >  7:67512/1‑149 (MQ=255)
          ttAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTT‑CGCGGGGAAGCCCGCGACAggggc                                                                                                                                  >  4:92752/1‑147 (MQ=255)
           tAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTACCCGGGAAAACCGACAAAAGGccc                                                                                                                                  >  7:300901/1‑149 (MQ=255)
            aaaGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTT‑CGCGGGTAAACCCGCGAAAGGGCCgg                                                                                                                                >  5:287541/1‑149 (MQ=255)
               gTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTT‑CGCGGGTAAACCCGCGAAAGGGCCGGCa                                                                                                                              >  5:291384/1‑148 (MQ=255)
                        cAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTT‑CGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTg                                                                                                                       <  8:339875/146‑1 (MQ=255)
                           tGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTT‑CGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAgg                                                                                                                 <  8:316235/149‑1 (MQ=255)
                            gAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCGTT‑GGCGGGGTATCCCGCGCAAGGGGCGGCAAATCCTGGTCAggg                                                                                                                >  5:208393/1‑149 (MQ=255)
                             aGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTT‑CGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGc                                                                                                               >  7:143205/1‑149 (MQ=255)
                               gCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGTATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTT‑CGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCt                                                                                                              >  2:224101/1‑148 (MQ=255)
                                    aaGCAAAAGACCCGCCCAGGCGGGTCTATTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTT‑CGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAAcc                                                                                                        >  3:119470/2‑149 (MQ=255)
                                        aaTAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTT‑AGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTACAACCCtta                                                                                                    >  4:57273/1‑148 (MQ=255)
                                         atcgACCCGCCGGGGCGGGTCTTTTGCTTCTATTTGCGCAAACAGCCCGCAGGCTTCAATCAATCTGCCCGCGGGGGGGCAGGGCCGGCCCTTT‑CGCGGGTAAACCCGCTAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTcc                                                                                                   <  6:208393/146‑1 (MQ=255)
                                             aCCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTT‑CGCGGGTAAACCCCCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACaa                                                                                               <  1:253670/149‑1 (MQ=255)
                                               ccGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTT‑CGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAAt                                                                                              >  2:66880/1‑148 (MQ=255)
                                                       gCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTT‑CGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTc                                                                                     <  1:224101/149‑1 (MQ=255)
                                                        cGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTT‑CGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCa                                                                                    <  6:291384/149‑1 (MQ=255)
                                                                 tGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTT‑CGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGc                                                                            <  1:66880/148‑1 (MQ=255)
                                                                    ttCTATTTGCGGAAACAGCCCGCAGGCTTCTATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTT‑CGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCATCGCCaag                                                                        >  1:28920/1‑147 (MQ=255)
                                                                    ttCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTT‑CGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCAcc                                                                        >  5:322620/1‑149 (MQ=255)
                                                                    ttCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTT‑CGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCATCGCCaag                                                                        >  1:261620/1‑147 (MQ=255)
                                                                             cGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTT‑CGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGa                                                               <  3:57273/149‑1 (MQ=255)
                                                                              ggAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTT‑CGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACTATGATCACCTCAGCCTTCGCCACCCCTTGGCGAt                                                              <  1:48587/149‑1 (MQ=255)
                                                                              ggAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTT‑CGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGAt                                                              >  7:158062/1‑149 (MQ=255)
                                                                              ggAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTT‑CGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGAt                                                              <  8:67512/149‑1 (MQ=255)
                                                                              ggAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTT‑CGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGAt                                                              <  4:119470/149‑1 (MQ=255)
                                                                                aaaCAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTT‑CGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATgg                                                            <  6:322620/149‑1 (MQ=255)
                                                                                       ccGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTT‑CGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCAc                                                                         >  6:171232/1‑129 (MQ=255)
                                                                                       ccGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTT‑CGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCAc                                                                         <  5:171232/129‑1 (MQ=255)
                                                                                         gCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTT‑AGCGGGTTAACCCGCGCAAGGGCCGGCCAATACTGCTCAGGGCTGAAACCCATGCACAATGAAAACCTCAGCCATGGCCCCCCCGTGGCGATGGCtgcctgcc                                                   >  6:314992/1‑149 (MQ=255)
                                                                                              cTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTT‑CGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTg                                              <  5:314992/149‑1 (MQ=255)
                                                                                                            ccAGCGGTGCGGCAGGTCCGGCCCTTT‑AGCCGGTAAACCCCCGAAAGGGGCGGCAAATACTGCTAAGGGCTCAAACCCTACCACAATGATCACCACAGCCTTAGCCCCCCCATGGCCATGGCTGCATGGCTCTTTGTACAGCGGGGAGc                                >  8:144104/1‑149 (MQ=255)
                                                                                                             caggtttGCGGCAGGTACGGCGCTTT‑GGCTGTTAAACCGGGGAAAGGGCCGGCAAATCGTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGGCTTAGCCACCCATTGGCGATGGCTGCATGCCTATTGGTACAGCGCGGAGCg                               <  8:300901/143‑1 (MQ=255)
                                                                                                               gCGGTGCGGCAGGTCCGGCCCTTT‑CGCGGTTAAACCCGCTAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGt                             <  3:164430/149‑1 (MQ=255)
                                                                                                               gCGGTGCGGCAGGTCCGGCCCTTT‑CGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGt                             <  8:158062/149‑1 (MQ=255)
                                                                                                                 ggTGCGGCAGGTCCGGCCCTTT‑CGCGGGTAACCCCGCTAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCAccc                                                                       <  1:158951/105‑1 (MQ=255)
                                                                                                                    gCGGCAGGTCCGGCCCTTT‑CGCGGGTAACCCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAg                        <  7:144104/149‑1 (MQ=255)
                                                                                                                           gTCCGGCCCTTT‑CGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGa                                  >  5:306928/1‑132 (MQ=255)
                                                                                                                           gTCCGGCCCTTT‑CGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGa                                  <  6:306928/132‑1 (MQ=255)
                                                                                                                            tCCGGCCCTTT‑CGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGc                <  1:339770/149‑1 (MQ=255)
                                                                                                                                  ccTTT‑CGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTgcgc          >  7:274497/1‑149 (MQ=255)
                                                                                                                                   ctat‑CGCGGGTTTACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGcgcg         >  7:92445/4‑149 (MQ=255)
                                                                                                                                   cTTT‑CGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGcgcg         >  3:65782/1‑149 (MQ=255)
                                                                                                                                           gggTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGAAt  <  8:143205/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                
GGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTT‑CGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGAAT  >  NC_002947/4338351‑4338637

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