Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,919,006 A→T 8.0% T469S (ACC→TCC)  PP_4328 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,919,0060AT 8.0% 131.3 / 3.3 50T469S (ACC→TCC) PP_4328hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (23/22);  new base T (2/2);  total (26/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 = 7.93e-01

CCAACATGAAGCCCAAGACCCTCGACAAGCCCTGTGGCGCATCGAGCTGGCCTTCGACCTGTCCCCGCTCGGCCCCCTGCAAGTGCAGGCACAACTGAGCGA‑GGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACACCCAGCTTGGCGCGCTGCGCGCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGCCACCCGGGCATCCCGCCGCAGGGGCCACGCACACGGATGGAACAACGCTGGGTGGACGAAACCG  >  NC_002947/4918864‑4919153
                                                                                                                                               |                                                                                                                                                   
ccAACATGAAGCCCAAGACCCTCGACAAGCCCTGTGGCGCATCGAGCTTGCCTTCGACCTGTCCCCGCTCGGCCCCCTGCAAGTGCAGGCACAACTGAGCGA‑GGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACcgcgc                                                                                                                                               <  2:92185/149‑1 (MQ=255)
 cAACATGAAGCCCAAGACCCTCGACAAGCCCTGTGGCGCATCGAGCTGGCCTTCGACCTGTCCCCGCTCGGCCCCCTGCAAGTGCAGGCACAACTGAGCGA‑GGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGTCcgcgcg                                                                                                                                              >  7:58694/1‑149 (MQ=255)
 cAACATGAAGCCCAAGACCCTCGACAAGCCCTGTGGCGCATCGAGCTGGCCTTCGACCTGTCCCCGCTCGGCCCCCTGCAAGTGCAGGCACAACTGAGCGA‑GGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGTCcgcgcg                                                                                                                                              >  7:54472/1‑149 (MQ=255)
      tGAAGCCCAAGACCCTCGACAAGCCCTGTGGCGCATCGAGCTGGCCTTCGACCTGTCCCCGCTCGGCCCCCTGCAAGTGCAGGCACAACTGAGCGA‑GGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGa                                                                                                                                         >  4:149465/1‑149 (MQ=255)
             cAAGACCCTCGACAAGCCCTGTGGCGCATCGAGCTGGCCTTCGACCTGTCCCCGCTCGGCCCCCTGCAAGTGCAGGCACAACTGAGCGA‑GGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGacac                                                                                                                                  >  2:294773/1‑149 (MQ=255)
              aaGACCCTCGACAAGCCCTGTGGCGCATCGAGCTGGCCTTCGACCTGTCCCCGCTCGGCCCCCTGCAAGTGCAGGCACAACTGAGCGA‑GGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACAcc                                                                                                                                 <  7:15668/149‑1 (MQ=255)
                gACCCTCGACAAGCCCTGTGGCGCATCGAGCTGGCCTTCGACCTGTCCCCGCTCGGCCCCCTGCAAGTGCAGGCACAACTGAGCGA‑GGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACACCCa                                                                                                                               <  1:256025/149‑1 (MQ=255)
                    cTCGACAAGCCCTGTGGCGCATCGAGCTGGCCTTCGACCTGTCCCCGCTCGGCCCCCTGCAAGTGCAGGCACAACTGAGCGA‑GGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACACCCAGCtt                                                                                                                           >  7:213091/1‑149 (MQ=255)
                    cTCGACAAGCCCTGTGGCGCATCGAGCTGGCCTTCGACCTGTCCCCGCTCGGCCCCCTGCAAGTGCAGGCACAACTGAGCGA‑GGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACACCCAGCtt                                                                                                                           <  3:82992/149‑1 (MQ=255)
                     tCGACAAGCCCTGTGGCGCATCGAGCTGGCCTTCGACCTGTCCCCGCTCGGCCCCCTGCAAGTGCAGGCACAACTGAGCGA‑GGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACACCCAGCTTg                                                                                                                          >  2:264138/1‑149 (MQ=255)
                            gCCCTGTGGCGCATCGAGCTGGCCTTCGACCTGTCCCCGCTCGGCCCCCTGCAAGTGCAGGCACAACTGAGCGA‑GGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGTCCGCGCGGTTGATCGACACCCAGCTTGGCGCGCt                                                                                                                   <  8:58694/149‑1 (MQ=255)
                            gCCCTGTGGCGCATCGAGCTGGCCTTCGACCTGTCCCCGCTCGGCCCCCTGCAAGTGCAGGCACAACTGAGCGA‑GGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGTCCGCGCGGTTGATCGACACCCAGCTTGGCGCGCt                                                                                                                   <  8:54472/149‑1 (MQ=255)
                            gCCCTGTGGCGCATCGAGCTGGCCTTCGACCTGTCCCCGCTCGGCCCCCTGCAAGTGCAGGCACAACTGAGCGA‑GGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACACCCAGCTTGGCGCGCt                                                                                                                   >  2:232602/1‑149 (MQ=255)
                                           gAGCTGGCCTTCGACCTGTCCCCGCTCGGCCCCCTGCAAGTGCAGGCACAACTGAGCGA‑GGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACACCCAGCTTGGCGCGCTGCGCGCGCGCctgct                                                                                                    <  1:286433/149‑1 (MQ=255)
                                                ggCCTTCGACCTGTCCCCGCTCGGCCCCCTGCAAGTGCAGGCACAACTGAGCGA‑GGGCCGCCTGAGCGGGCAGCGATGGGCGGAGCGCGAGCAGCCCGCGCGGTGGCTCGACACCCAGCTTGGCGCGCTGCGCGCGCGCCTGCTGGccc                                                                                               >  5:101245/1‑149 (MQ=255)
                                                      cGACCTGTCCCCGCTCGGCCCCCTGCAAGTGCAGGCACAACTGAGCGA‑GGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACACCCAGCTTGGCGCGCTGCGCGCGCGCCTGCTGGCCCgggg                                                                                           <  1:232602/147‑1 (MQ=255)
                                                      cGACCTGTCCCCGCTCGGCCCCCTGCAAGTGCAGGCACAACTGAGCGA‑GGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACACCCAGCTTGGCGCGCTGCGCGCGCGCCTGCTGGCCCGGGGc                                                                                          <  2:103819/148‑1 (MQ=255)
                                                        aCCTGTCCCCGCTCGGCCCCCTGCAAGTGCAGGCACAACTGAGCGA‑GGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACACCCAGCTTGGCGCGCTGCGCGCGCGCCTGCTGGCCCGGGGCCTg                                                                                       <  3:149465/149‑1 (MQ=255)
                                                               cccGCTCGGCCCCCTGCAAGTGCAGGCACAACTGAGCGAGGGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACACCCAGCTTGGCGCGCTGCGCGCGCGCCTGCTGGCCCGGGGCCTGGAAGTc                                                                                 >  1:214097/1‑149 (MQ=255)
                                                               cccGCTCGGCCCCCTGCAAGTGCAGGCACAACTGAGCGA‑GGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACACCCAGCTTGGCGCGCTGCGCGCGCGCCTGCTGGCCCGGGGCCTGGAAGTCg                                                                                <  5:290084/149‑1 (MQ=255)
                                                                       gCCCCCTGCAAGTGCAGGCACAACTGAGCGA‑GGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACACCCAGCTTGGCGCGCTGCGCGCGCGCCTGCTGGCCCGGGGcc                                                                                         <  6:267328/132‑1 (MQ=255)
                                                                       gCCCCCTGCAAGTGCAGGCACAACTGAGCGA‑GGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACACCCAGCTTGGCGCGCTGCGCGCGCGCCTGCTGGCCCGGGGcc                                                                                         >  5:267328/1‑132 (MQ=255)
                                                                           ccTGCAAGTGCAGGCACAACTGAGCGA‑GGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACACCCAGCTTGGCGCGCTGCGCGCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGt                                                                    <  7:29582/149‑1 (MQ=255)
                                                                            cTGCAAGTGCAGGCACAACTGAGCGA‑GGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACACCCAGCTTGGCGCGCTGCGCGCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTg                                                                   >  3:89168/1‑149 (MQ=255)
                                                                                    gCAGGCACAACTGAGCGA‑GGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACACCCAGCTTGGCGCGCTGCGCGCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGCCACCCgg                                                           >  7:205620/1‑149 (MQ=255)
                                                                                      aGGCACAACTGAGCGA‑GGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACACCCAGCTTGGCGCGCTGCGCGCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGCCACCCGGGc                                                         <  1:264138/149‑1 (MQ=255)
                                                                                            aaCTGAGCGA‑GGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACACCCAGCTTGGCGCGCTGCGCGCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGCCACCCGGGCATCccg                                                   >  1:259168/1‑149 (MQ=255)
                                                                                            aaCTGAGCGA‑GGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACACCCAGCTTGGCGCGCTGCGCGCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGCCACCCGGGCATCccg                                                   >  1:78779/1‑149 (MQ=255)
                                                                                             aCTGAGCGA‑GGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACACCCAGCTTGGCGCGCTGCGCGCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGCCACCCGGGCATCccgc                                                  >  6:21228/1‑149 (MQ=255)
                                                                                              cTGAGCGA‑GGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACACCCAGCTTGGCGCGCTGCGCGCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGCCACCCGGGCAt                                                       >  1:291983/1‑143 (MQ=255)
                                                                                              cTGAGCGA‑GGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACACCCAGCTTGGCGCGCTGCGCGCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGCCACCCGGGCAt                                                       <  2:291983/143‑1 (MQ=255)
                                                                                              cTGAGCGA‑GGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACACCCAGCTTGGCGCGCTGCGCGCGCGACTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGCCACCCGGGCATCccgcc                                                 >  4:38297/1‑149 (MQ=255)
                                                                                                          ccgccTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACACCCAGCTTGGCGCGCTGCGCGCGCGCCTGCTAGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGCCACCCgg                                                           <  1:281072/128‑1 (MQ=255)
                                                                                                          ccgccTGAGCGGGCAGCTATGGGCGGAGCGCGAGCACACCGCGCGGTTGATCGACACCCAGCTTGGCGCGCTGCCCCCGCGCCAGCTAGCCAGGGGCCTCGAAGTCGGAGACCTGCAGAGCCAacccg                                                           >  2:281072/1‑126 (MQ=255)
                                                                                                                    gggCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACACCCAGCGTGGCGCGCTGCGCGCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGCCACCCGGGCATCCCGCCGCAGGGGCCACGCACACGGAt                            >  3:30801/1‑149 (MQ=255)
                                                                                                                     ggCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTGGATCGACACCCAGCTTGGCGCGCTGCGCGCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGCCACCCGGGCATCCCGCCGCAGGGGCCACGCACACGGATg                           <  4:19822/149‑1 (MQ=255)
                                                                                                                      gCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACACCCAGCTTGGCGCGCTGCGCGCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGCCACCCGGGCATCCCGCCGCAGGGGCCACGCACACGGATgg                          <  6:101245/149‑1 (MQ=255)
                                                                                                                        aGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACACCCAGCTTGGCGCGCTGCGCGCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGCCACCCGGGCATCCCGCCGCAGGGGCCACGCACACGGATGGaa                        <  2:259168/149‑1 (MQ=255)
                                                                                                                          cTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACACCCAGCTTGGCGCGCTGCGCGCGCGCCAGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGCCACCCGGGCATCCCGCCGCAGGGGCCACGCACACGGATGGaac                       >  8:46045/1‑148 (MQ=255)
                                                                                                                           tATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACACCCAGCTTGGCGCGCTGCGCGCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGCCACCCGGGCATCCCGCCGCAGGGGCCACGCACACGGATGGaacca                     >  8:265878/1‑147 (MQ=255)
                                                                                                                             tGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACACCCAGCTTGGCGCGCTGCGCGCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGCCACCCGGGCATCCCGCCGCAGGGGCCACGCACACGGATGGaacaac                    <  2:78779/148‑1 (MQ=255)
                                                                                                                             tGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACACCCAGCTTGGCGCGCTGCGCGCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGCCACCCGGGCATCCCGCCGCAGGGGCCACGCACACGGATGGAACAACg                   >  5:87872/1‑149 (MQ=255)
                                                                                                                                    aGCGCGAGCAGACCGCGCGGTTGATCGACACCCAGCTTGGCGCGCTGCGCGCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGCCACCCGGGCATCCCGCCGCAGGGGCCACGCACACGGATGGAACAACGCTGGGTg            <  6:87872/149‑1 (MQ=255)
                                                                                                                                      cgcgAGCAGACCGCGCGGTTGATCGACACCCAGCTTGGCGCGCTGCGCGCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGCCACCCGGGCATCCCGCCGCAGGGGCCACGGACACGGATGGAACAACGCTGGGTGGa          <  2:214097/149‑1 (MQ=255)
                                                                                                                                       gcgAGCAGACCGCGCGGTTGATCGACACCCCGCTTGGCGCGCTGCGCGCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGCCACCCGGGCATCCCGCCGCAGGGGCCACGCACACGGATGGAACAACGCTg                <  8:222218/142‑1 (MQ=255)
                                                                                                                                       gcgAGCAGACCGCGCGGTTGATCGACACCCAGCTTGGCGCGCTGCGCGCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGCCACCCGGGCATCCCGCCGCAGGGGCCACGCACACGGATGGAACAACGCTg                >  7:222218/1‑142 (MQ=255)
                                                                                                                                       gcgAGCAGACCGCGCGGTTGATCGACACCCAGCTTGGCGCGCTGCGCGCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGCCACCCGGGCATCCCGCCGCAGGGGCCACGCACACGGATGGAACAACGCTGGGTGGAc         <  4:89168/149‑1 (MQ=255)
                                                                                                                                             aGACCGCGCGGTTGATCGACACCCAGCTTGGCGCGCTGCGCGCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGCCACCCGGGCATCCCGCCGCAGGGGCCACGCACACGGATGGAACAACGCTGGGTGGACGAAAcc   >  8:32847/1‑149 (MQ=255)
                                                                                                                                             aGACCGCGCGGTTGATCGACACCCAGCTTGGCGCGCTGCGCGCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGCCACCCGGGCATCCCGCCGCAGGGGCCACGCACACGGATCGAACAACGCTGCGTGCACGAAAcc   >  4:30259/1‑149 (MQ=255)
                                                                                                                                              gACCGCGCGGTTGATCGACACCCAGCTTGGCGCGCTGCGCGCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGCCACCCGGGCATCCCGCCGCAGGGGCCACGCACACGGATGGAACAACGCTGGGTGGACGAAACCg  >  7:133647/1‑149 (MQ=255)
                                                                                                                                               aCCGCGCGGTTGATCGACACCCAGCTTGGCGCGCTGCGCGCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGCCACCCGGGCATCCCGCCGCAGGGGCCACGCACACGGATGGaacaa                     >  4:18986/1‑129 (MQ=255)
                                                                                                                                               aCCGCGCGGTTGATCGACACCCAGCTTGGCGCGCTGCGCGCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGCCACCCGGGCATCCCGCCGCAGGGGCCACGCACACGGATGGaacaa                     <  3:18986/129‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
CCAACATGAAGCCCAAGACCCTCGACAAGCCCTGTGGCGCATCGAGCTGGCCTTCGACCTGTCCCCGCTCGGCCCCCTGCAAGTGCAGGCACAACTGAGCGA‑GGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACACCCAGCTTGGCGCGCTGCGCGCGCGCCTGCTGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTGCCACCCGGGCATCCCGCCGCAGGGGCCACGCACACGGATGGAACAACGCTGGGTGGACGAAACCG  >  NC_002947/4918864‑4919153

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