Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,206,285 A→T 7.8% Y530F (TAC→TTC)  PP_2811 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,206,2850AT 7.8% 148.5 / 3.0 51Y530F (TAC→TTC) PP_2811hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (24/23);  new base T (2/2);  total (26/25)
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.70e-01

GCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGC‑CACCTTCACCGAAGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTCCCTCAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGCCTGCCGTGAACGGTTCTGTAAAGCTATCTGAAGCCCCGGCACGTGGCCTGTTGTCACGCGTCGAGGGTGTG  >  NC_002947/3206138‑3206417
                                                                                                                                                    |                                                                                                                                    
gctcgctTGGCGGGCGAACGCCGAACCTGAGCGCACCCCTTCACCGAAGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTCCCTCAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTa                                                                                                                                      >  1:223914/1‑149 (MQ=255)
    gctTGGCGGGCGAACGCCGAACCTGAGCGC‑CACCTTCACCGAAGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTCCCTCAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGAc                                                                                                                                  >  1:152958/1‑148 (MQ=255)
        ggcgggcgAACGCCGAACCTGAGCGC‑CACCTTCACCGAAGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTCCCTCAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGc                                                                                                                              <  7:259241/148‑1 (MQ=255)
             gcgAACGCCGAACCTGAGCGC‑CACCTTCACCGAAGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTCCCTCAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATc                                                                                                                        <  4:333967/149‑1 (MQ=255)
                      gAACCTGAGCGC‑CACCTTCACCGAAGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATCTCCCTCAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACcgc                                                                                                                >  5:13993/1‑148 (MQ=255)
                           tGAGCGC‑CACCTTCACCGAAGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTCCCTCAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTAt                                                                                                          >  6:146277/1‑149 (MQ=255)
                           tGAGCGC‑CACCTTCACCGAAGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTCCCTCAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTAt                                                                                                          >  5:184135/1‑149 (MQ=255)
                           tGAGCGC‑CACCTTCACCGAAGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTCCCTCAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTAt                                                                                                          >  5:106975/1‑149 (MQ=255)
                           tGAGCGC‑CACCTTCACCGAAGACGCGAGCAGAATGAACCTGTACCTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTCCCTCAACTACGCGAAGTTCATGGGCGGCGAGACCCCCAACGACCAGCTTTAACGTGACCGAGAatat                                                                                                          >  4:245134/1‑149 (MQ=255)
                              gcgc‑CACCTTCACCGAAGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTCCCTCAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTc                                                                                                       <  2:223914/149‑1 (MQ=255)
                                             gAAGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTCCCTCAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTc                                                                                         >  7:81682/1‑149 (MQ=255)
                                               aGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATCTCCCTCAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTccc                                                                                       <  6:13993/149‑1 (MQ=255)
                                                       gcagcaTGAACCTGTACCTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTCCCTCAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTtg                                                                               >  1:309674/1‑149 (MQ=255)
                                                          gcaTGAACCTGTACCTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTCCCTCAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCGTCATCGCGATTTCCCGCACCTTGTGa                                                                            <  2:309674/149‑1 (MQ=255)
                                                              gAACCTGTACCTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTCCCTCAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGcc                                                                        >  5:328778/1‑149 (MQ=255)
                                                                 ccTGTACCTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTCCCTCAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGCCTGc                                                                     <  8:81682/149‑1 (MQ=255)
                                                                  cTGTACCTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTCCCTCAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATc                                                                                                                        <  4:95382/97‑1 (MQ=255)
                                                                  cTGTACCTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTCCCTCAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATc                                                                                                                        >  3:95382/1‑97 (MQ=255)
                                                                     tACCTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTCCCTCAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGCCTGCCg                                                                   <  6:106975/147‑1 (MQ=255)
                                                                       ccTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTCCCTCAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGCCTGCCGTGaa                                                               >  5:102763/1‑149 (MQ=255)
                                                                        cTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTCCCTCAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTc                                                                                                       <  7:313027/108‑1 (MQ=255)
                                                                        cTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTCCCTCAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTc                                                                                                       >  8:313027/1‑108 (MQ=255)
                                                                        cTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTCCCTCAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGCCTGCCGTGAAc                                                              <  2:152958/149‑1 (MQ=255)
                                                                          gAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTCCCTCAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGCCTGCCGTGAACgg                                                            >  6:192680/1‑149 (MQ=255)
                                                                                    cGAAACCCGGCCAGTTGGCAGATTTCCCTCAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGCCTGCCGTGAACGGTTCTGTAAAg                                                  <  6:184135/149‑1 (MQ=255)
                                                                                          ccGGCCAGTTGGCAGATTTCCCTCAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGCCTGCCGTGAACGGTTCTGTAAAGCTATCt                                            >  4:324266/1‑149 (MQ=255)
                                                                                           cGGCCAGTTGGCAGATTTCCCTCAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGCCTGCCGTGAACGGTTCTGTAAAGCTATCTg                                           <  6:328778/149‑1 (MQ=255)
                                                                                             gCCAGTTGGCAGATTTCCGTCAACTACGCGAAGTTCATGGGCTGCGAAACCCCATACGTCCAGCTTTATCGTGACCGCGATTATGTCGGCATCTCGATTTCCCGCACCTTGTGAGGCCTGCCGTGAACGGTTCTGTAAAGCTATCTGaa                                         <  3:245134/149‑1 (MQ=255)
                                                                                                    ggCAGATTTCCCTCAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGCCTGCCGTGAACGGTTCTGTAAAGCTATCTGAAGCCCCg                                   <  5:146277/148‑1 (MQ=255)
                                                                                                              ccTCAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGCCTGCCGTGAACGGTTCTGTAAAGCTATCTGAAGCCCCGGCACGTGGCCt                        >  3:321102/1‑149 (MQ=255)
                                                                                                               cTCAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCgg                                                                                                     >  5:190502/1‑71 (MQ=255)
                                                                                                               cTCAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCgg                                                                                                     <  6:190502/71‑1 (MQ=255)
                                                                                                                tCAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTTCGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGCCTGCCGTGAACGATTCTGTAAAGCTATCTGAAGCCCCGGCACGTGGCCtgt                      >  5:105988/1‑149 (MQ=255)
                                                                                                                tCAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTTCGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGCCTGCCGTGAACGATTCTGTAAAGCTATCTGAAGCCCCGGCACGTGGCCtgt                      >  5:101015/1‑149 (MQ=255)
                                                                                                                tCAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGCCTGCCGTGAACGGTTCTGTAAAGCTATCTGAAGCCCCGGCACGTGGCCtgt                      >  8:295064/1‑149 (MQ=255)
                                                                                                                tCAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGCCTGCCGTGAACGGTTCTGTAAAGCTATCTGAAGCCCCGGCACGTGACCtgt                      >  5:121585/1‑149 (MQ=255)
                                                                                                                             gTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGCCTGCCGTGAACGGTTCTGTAAAGCTATCTGAAGCCCCGGCACGTGGCCTGTTGTCACGCGTCGa         <  3:324266/149‑1 (MQ=255)
                                                                                                                              ttCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGCCTGCCGTGAACGGTTCTGTAAAGCTATCTGAAGCCCCGGCACGTGGCCTGTTGTCACGCGTCGAg        <  7:295064/149‑1 (MQ=255)
                                                                                                                               tCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGCCTGCCGTGAACGGTTCTGTAAAGCAATCTGAAGCCCCGGCACGTGGCCTGTTGTCACGCGTCGAgg       >  5:305405/1‑149 (MQ=255)
                                                                                                                                cATGGGCGGCGAAACCCCCTTCGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGCCTGCCGTGAACGATTCTGTAAAGCTATCTGAAGCCCCGGCACGTGGCCTGTTGTCACGCGTCGAgg       <  6:105988/148‑1 (MQ=255)
                                                                                                                                cATGGGCGGCGAAACCCCCTTCGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGCCTGCCGTGAACGATTCTGTAAAGCTATCTGAAGCCCCGGCACGTGGCCTGTTGTCACGCGTCGAgg       <  6:101015/148‑1 (MQ=255)
                                                                                                                                cATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGCCTGCCGTGAACGGTTCTGTAAAGCTATCTGAAGCCCCGGCACGTGGCCTGTTGTCACGCGTCGAgg       <  6:305405/148‑1 (MQ=255)
                                                                                                                                 aTGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGCCTGCCGTGAACGGTTCTGTAAAGCTATCTGAAGCCCCGGCACGTGGCCtg                       <  6:274913/131‑1 (MQ=255)
                                                                                                                                 aTGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGCCTGCCGTGAACGGTTCTGTAAAGCTATCTGAAGCCCCGGCACGTGGCCtg                       >  5:274913/1‑131 (MQ=255)
                                                                                                                                   gggCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGCCTGCCGTGAACGGTTCTGTAAAGCTATCTGAAGCCCCGGCACGTGGCCTGTTGTCACGCGTCGAGGgtgt   >  1:169/1‑149 (MQ=255)
                                                                                                                                    ggcggcGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGCCTGCCGTGAACGGTTCTGTAAAGCTATCTGAAGCCCCGGCACGTGGCCTGTTGTCACGCGTCGAGGgtgt   <  5:192680/148‑1 (MQ=255)
                                                                                                                                       ggcgAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGCCTGCCGTGAACGGTTCTGTAAAGCTATCTGAAGCCCCGGCACGTGGCCTGTTGTCAc                <  1:158181/132‑1 (MQ=255)
                                                                                                                                       ggcgAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGCCTGCCGTGAACGGTTCTGTAAAGCTATCTGAAGCCCCGGCACGTGGCCTGTTGTCAc                >  2:158181/1‑132 (MQ=255)
                                                                                                                                               ccccTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGCCTGCCGTGAACGGTTCTGTAAAGCTATCTGAAGCCCCGGCACGTGGCCTGTTGTCACGCGTCGAGGgtgtg  >  3:81675/1‑138 (MQ=255)
                                                                                                                                               ccccTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGCCTGCCGTGAACGGTTCTGTAAAGCTATCTGAAGCCCCGGCACGTGGCCTGTTGTCACGCGTCGAGGgtgtg  <  4:81675/138‑1 (MQ=255)
                                                                                                                                                  cTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGCCTGCCGTGAACGGTTCTGTAAAGCTATCTGAAGCCCCGGCACGTGGCCTGTTGTCACGCGTc           <  6:203346/126‑1 (MQ=255)
                                                                                                                                                  cTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGCCTGCCGTGAACGGTTCTGTAAAGCTATCTGAAGCCCCGGCACGTGGCCTGTTGTCACGCGTc           >  5:203346/1‑126 (MQ=255)
                                                                                                                                                    |                                                                                                                                    
GCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGC‑CACCTTCACCGAAGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTCCCTCAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGCCTGCCGTGAACGGTTCTGTAAAGCTATCTGAAGCCCCGGCACGTGGCCTGTTGTCACGCGTCGAGGGTGTG  >  NC_002947/3206138‑3206417

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