Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 591,152 G→A 8.6% S133S (TCG→TCA)  PP_0501 → NAD‑dependent epimerase/dehydratase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947591,1520GA 8.6% 130.9 / 2.8 47S133S (TCG→TCA) PP_0501NAD‑dependent epimerase/dehydratase family protein
Reads supporting (aligned to +/- strand):  ref base G (21/21);  new base A (2/2);  total (23/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.40e-01

GGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCAGTACCTGGACTTCTACCGTCGCCAGCATGGTCTGGAGCCGGTGGTGTTTCGCTTCTTCAACATCTTCGGGCCACGGCAGGAC  >  NC_002947/591005‑591299
                                                                                                                                                   |                                                                                                                                                   
ggCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAACCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGa                                                                                                                                                    <  2:233184/149‑1 (MQ=255)
 gCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGAt                                                                                                                                                   >  5:173897/1‑149 (MQ=255)
       gTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCg                                                                                                                                              >  2:117251/1‑148 (MQ=255)
            aGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGac                                                                                                                                         <  8:50573/148‑1 (MQ=255)
               cccGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGc                                                                                                                                     >  4:197845/1‑149 (MQ=255)
                ccGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGc                                                                                                                                     >  4:205126/1‑148 (MQ=255)
                      cGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGgc                                                                                                                              >  8:212715/1‑149 (MQ=255)
                          cccACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCg                                                                                                                          <  6:277227/149‑1 (MQ=255)
                                                   ccTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGa                                                                                                                                          >  4:77416/1‑108 (MQ=255)
                                                   ccTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGa                                                                                                                                          <  3:77416/108‑1 (MQ=255)
                                                   ccTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCAATTTCCGAAGACACGCCAAAGGCGCCGtt                                                                                                                        >  6:24237/1‑126 (MQ=255)
                                                   ccTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCAATTTCCGAAGACACGCCAAAGGCGCCGtt                                                                                                                        >  6:230770/1‑126 (MQ=255)
                                                   ccTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCAATTTCCGAAGACACGCCAAAGGCGCCGtt                                                                                                                        <  5:24237/126‑1 (MQ=255)
                                                   ccTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCAATTTCCGAAGACACGCCAAAGGCGCCGtt                                                                                                                        <  5:230770/126‑1 (MQ=255)
                                                      cAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGa                                                                                                                      >  3:199302/1‑125 (MQ=255)
                                                      cAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGa                                                                                                                      <  4:199302/125‑1 (MQ=255)
                                                         cGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCa                                                                                           <  3:197845/149‑1 (MQ=255)
                                                         cGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCATGGACAAGCTGGCCa                                                                                           <  3:205126/149‑1 (MQ=255)
                                                              gCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCt                                                                                                                 >  2:257481/1‑122 (MQ=255)
                                                              gCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCt                                                                                                                 <  1:257481/122‑1 (MQ=255)
                                                              gCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGa                                                                                       <  7:212715/148‑1 (MQ=255)
                                                               cGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCTCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAg                                                                                      >  7:182346/1‑148 (MQ=255)
                                                               cGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAg                                                                                      >  5:207303/1‑148 (MQ=255)
                                                               cGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAg                                                                                      >  7:58396/1‑148 (MQ=255)
                                                                gAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTg                                                                                               >  2:92907/1‑138 (MQ=255)
                                                                 acgCCATGCACGCGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCAg                                                                                   >  2:58854/3‑149 (MQ=255)
                                                                            gTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCAGTACCTGGACt                                                                         >  3:204699/1‑148 (MQ=255)
                                                                                      tGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCAGTACCTGGACTTCTACCGTCGc                                                              <  5:73892/149‑1 (MQ=255)
                                                                                           cgcgTGCTATTCGCCTCCATCGCAGCTGTGTACGTCAACAATGTTGAATGGGATTCTTTTTCCGAAGAAACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCAGTACCTTGACTTCTACCTTCGCCAgaa                                                         <  8:182346/146‑3 (MQ=255)
                                                                                           cGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCAGTACCTGGACTTCTACCGTCGCCAGCa                                                         >  1:146549/1‑149 (MQ=255)
                                                                                                 ctctTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCAGTACCTGGACTTCTACCGTCGCCAGCATGGTc                                                    <  1:58854/148‑1 (MQ=255)
                                                                                                 ctctTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCAGTACCTGGACTTCTACCGTCGCCAGCATGGTCt                                                   <  6:207303/149‑1 (MQ=255)
                                                                                                            cagcgcagcgGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCAGTACCTGGACTTCTACCGTCGCCAGCATGGTCTGGAGCCggtgg                                        <  1:117251/149‑1 (MQ=255)
                                                                                                            cagcgcagcgGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCAGTACCTGGACTTCTACCGTCGCCAGCATGGTCTGGAGCCggtg                                         <  6:173897/148‑1 (MQ=255)
                                                                                                               cgcagcgGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCAGTACCTGGACTTCTACCGTCGCCAGCATGGTCTg                                                  >  4:73430/1‑136 (MQ=255)
                                                                                                               cgcagcgGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCAGTACCTGGACTTCTACCGTCGCCAGCATGGTCTg                                                  <  3:73430/136‑1 (MQ=255)
                                                                                                                 cagcgGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCAGTACCTGGACTTCTACCGTCGCCAGCATGGTCTGGAGCCGGTGGTGttt                                   <  2:146549/149‑1 (MQ=255)
                                                                                                                    cgGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCAGTACCTGGACTTCTACCGTCGCCAGCATGGTCTGGAGCCGGTGGTGTTTCGc                                >  7:165646/1‑149 (MQ=255)
                                                                                                                           cGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCAGTACCTGGACTTCTACCGTCGCCAGCATGGTCTGGAGCCGGTGGTGTTTCGCTTCTTCa                         <  8:58396/149‑1 (MQ=255)
                                                                                                                             gCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCAGTACCTGGACTTCTACCGTCGCCAGCATGGTCTGGAGCCGGTGGTGTTTCGCTTCTTCAAc                       <  4:119920/149‑1 (MQ=255)
                                                                                                                              caacaaTGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCAGTACCTGGACTTCTACCGTCGCCAGCATGGTCTGGAGCCGGTGGTGTTTCGCTTCTTCAACa                      >  3:154978/1‑149 (MQ=255)
                                                                                                                                  aaTGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTgg                                                                                              >  2:4919/1‑73 (MQ=255)
                                                                                                                                  aaTGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTgg                                                                                              <  1:4919/73‑1 (MQ=255)
                                                                                                                                  aaTGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGAGCCGTTAACCCCTTATGCCGTGGACAACCTGGCCAGTGAGCAGTACCTGGACTTCTACAGTCGCCAGCATGGTCTGGAGCCCCTGGTCTTTCGcttc                             >  2:92414/1‑138 (MQ=255)
                                                                                                                                       tGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCAGTACCTGGACTTCTACCGTCGCCAGCATGGTCTGGAGCCGGTGGTGTTTCGCTTCTTCAACATCTTCGGGc             <  4:154978/149‑1 (MQ=255)
                                                                                                                                                 tCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCAGTACCTGGACTTCTACCGTCGCCAGCATGGTCTGGAGCCGGTGGTGttt                                   >  2:176800/1‑117 (MQ=255)
                                                                                                                                                 tCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCAGTACCTGGACTTCTACCGTCGCCAGCATGGTCTGGAGCCGGTGGTGttt                                   <  1:176800/117‑1 (MQ=255)
                                                                                                                                                   gATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCAGTACCTGGACTTCTACCGTCGCCAGCATGGTCTGGAGCCGGTGGTGTTTCGCTTCTTCAACATCTTCGGGCCACGGCAGGAc  <  8:165646/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCAGTACCTGGACTTCTACCGTCGCCAGCATGGTCTGGAGCCGGTGGTGTTTCGCTTCTTCAACATCTTCGGGCCACGGCAGGAC  >  NC_002947/591005‑591299

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: ‑