Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,579,784 C→G 8.7% V361L (GTG→CTG) ‡ serB ← phosphoserine phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,579,7840CG 8.7% 138.0 / 4.7 46V361L (GTG→CTG) ‡serBphosphoserine phosphatase
Reads supporting (aligned to +/- strand):  ref base C (21/21);  new base G (2/2);  total (23/23)
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.41e-01

ACGCAGGCCCAGCAGGTACAGCACCCCGTCCAGCCCCAAGGTGGAAATCGCCTGCTTGGCCGACTGGCGCACCAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATCGACAATCGGCTCAACAGCCACGCCCGTCACCTTGCCATCGACCACTTCAAGTTCGTTGGCGAACACATAGTC  >  NC_002947/5579638‑5579931
                                                                                                                                                  |                                                                                                                                                   
aCGCAGGCCCAGCAGGTACAGCACCCCGTCCAGCCCCAAGGTGGAAATCGCCTGCTTGGCCGACTGGCGCACCAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACgg                                                                                                                                                   <  5:229460/149‑1 (MQ=255)
     ggCCCAGCAGGTTCAGCACCCCGTCCAGCCCCAAGGTGGAAATCGCCTGCTTGGCCGACTGGCGCACCAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATg                                                                                                                                              >  2:198273/1‑149 (MQ=255)
     ggCCCAGCAGGTACAGCACCCCGTCCAGCCCCAAGGTGGAAATCGCCTGCTTGGCCGACTGGCGCACCAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATg                                                                                                                                              <  2:194455/149‑1 (MQ=255)
       cccAGCAGGTACAGCACCCCGTCCAGCCCCAAGGTGGAAATCGCCTGCTTGGCCGACTGGCGCACCAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCGTGGGCGATggc                                                                                                                                            <  8:317431/149‑2 (MQ=255)
        ccAGCAGGTACAGCACCCCGTCCAGCCCCAAGGTGGAAATCGCCTGCTTGGCCGACTGGCGCACCAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCGTGGGCGATggc                                                                                                                                            >  4:200090/1‑147 (MQ=255)
        ccAGCAGGTACAGCACCCCGTCCAGCCCCAAGGTGGAAATCGCCTGCTTGGCCGACTGGCGCACCAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCGTGGGCGATggc                                                                                                                                            <  3:200090/148‑2 (MQ=255)
               gTACAGCACCCCGTCCAGCCCCAAGGTGGAAATCGCCTGCTTGGCCGACTGGCGCACCAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCa                                                                                                                                    <  7:17376/149‑1 (MQ=255)
               gTACAGCACCCCGTCCAGCCCCAAGGTGGAAATCGCCTGCTTGGCCGACTGGCGCACCAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCa                                                                                                                                    <  1:111017/149‑1 (MQ=255)
                  cAGCACCCCGTCCAGCCCCAAGGTGGAAATCGCCTGCTTGGCCGACTGGCGCACCAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGtt                                                                                                                                 <  3:266692/149‑1 (MQ=255)
                    gCACCCCGTCCAGCCCCAAGGTGGAAATCGCCTGCTTGGCCGACTGGCGCACCAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGc                                                                                                                                        <  8:114326/140‑1 (MQ=255)
                    gCACCCCGTCCAGCCCCAAGGTGGAAATCGCCTGCTTGGCCGACTGGCGCACCAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGc                                                                                                                                        >  7:114326/1‑140 (MQ=255)
                        cccGTCCAGCCCCAAGGTGGAAATCGCCTGCTTGGCCGACTGGCGCACCAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCa                                                                                                                              >  7:5787/1‑146 (MQ=255)
                        cccGTCCAGCCCCAAGGTGGAAATCGCCTGCTTGGCCGACTGGCGCACCAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCa                                                                                                                              <  8:5787/146‑1 (MQ=255)
                        cccGTCCAGCCCCAAGGTGGAAATCGCCTGCTTGGCCGACTGGCGCACCAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGc                                                                                                                           >  5:230495/1‑149 (MQ=255)
                             ccAGCCCCAAGGTGGAAATCGCCTGCTTGGCCGACTGGCGCACCAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTtcg                                                                                                                      >  5:102279/1‑149 (MQ=255)
                                    cAAGGTGGAAATCGCCTGCTTGGCCGACTGGCGCACCAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTgg                                                                                                               >  8:110026/1‑149 (MQ=255)
                                    cAAGGTGGAAATCGCCTGCTTGGCCGACTGGCGCAACAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGCCGATGGTATGCTCCAGTTGCAAGCCTTCGTCGCTgg                                                                                                               >  8:99151/1‑149 (MQ=255)
                                           gAAATCGCCTGCTTGGCCGACTGGCGCACCAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGct                                                                                                         >  2:261266/1‑148 (MQ=255)
                                              aTCGCCTGCTTGGCCGACTGGCGCACCAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGcttctt                                                                                                     <  6:102279/149‑1 (MQ=255)
                                                  cgtggttGGCGGTCTGGCGCATGAGCGTCTTGGCACGTAACGCCAGGCGCAACCCAGCCAGCGACAGCATCGTCATGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTGTTcaac                                                                                                 <  6:264438/144‑1 (MQ=255)
                                                      cTTGGCCGACTGGCGCACCAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAAt                                                                                             >  2:31250/1‑149 (MQ=255)
                                                      cTTGGCCGACTGGCGCACCAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAAt                                                                                             <  6:230495/149‑1 (MQ=255)
                                                           ccGACTGGCGCACCAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGcc                                                                                        <  1:295518/149‑1 (MQ=255)
                                                              aCTGGCGCACCAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTg                                                                                     <  1:261266/149‑1 (MQ=255)
                                                                   cgcACCAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTg                                                                                <  7:134588/149‑1 (MQ=255)
                                                                     cACCAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTggg                                                                              >  3:176241/1‑149 (MQ=255)
                                                                      aCCAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTggg                                                                              >  3:124904/1‑148 (MQ=255)
                                                                         aGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATc                                                                          <  1:198273/149‑1 (MQ=255)
                                                                           cGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATCGa                                                                        <  1:237060/149‑1 (MQ=255)
                                                                                  gCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACATATCGGCCTTGCGCTGGGCATCGACAATCgg                                                                 <  1:31250/149‑1 (MQ=255)
                                                                                   cACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATCGACAATCGGc                                                                >  6:105171/1‑149 (MQ=255)
                                                                                    aCGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATCGACAATCGGCt                                                               >  5:78989/1‑149 (MQ=255)
                                                                                        aaCGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATCGACAATCGGCTCAAc                                                           <  7:110026/149‑1 (MQ=255)
                                                                                         acgccacgccCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATCGACAATCGGCTCAACa                                                          <  4:176241/149‑1 (MQ=255)
                                                                                              acgccCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATCGACAATCGGCTCAACAGCCAc                                                     <  4:124904/149‑1 (MQ=255)
                                                                                                                gACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATCGACAATCGGCTCAACAGCCACGCCCGTCACCTTGCCATc                                   <  5:105171/149‑1 (MQ=255)
                                                                                                                     cATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCAGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATCGACAATCGGCTCAACAGCCACGCCCGTCACCTTGCCATCGACCa                              >  8:120334/1‑149 (MQ=255)
                                                                                                                       tCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATCGACAATCGGCTCAACAGCCACGCCCGTCACCTTGCCATCGACCACt                            >  4:301334/1‑149 (MQ=255)
                                                                                                                          gCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATCGACAATCGGCTCAACAGCCACGCCCGTCACCTTGCCATCGACCACTTc                          >  7:316408/1‑148 (MQ=255)
                                                                                                                             ggTCATAGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGCCCAGCTTCTTCAACAAATAGGCCTTGCGCTGGGCATCGACAATCGGCTAAAAAGCCACGCCCGTCACCTTGCCATCGCCAACTACAAGt                      >  8:87447/1‑149 (MQ=255)
                                                                                                                               tCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATCGACAATCGGCTCAACAGCCACGCCCGTCACCTTGCCATCGACCACTTCAAGTTc                    >  3:63485/1‑149 (MQ=255)
                                                                                                                                 aTTGGCGACATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATCGACAATCGGCTCAACAGCCACGCCCGTCACCATCCCATCGACCACTTCAAGTTcct                  >  6:69162/1‑147 (MQ=255)
                                                                                                                                       accATCGCCGTGGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATCGACAATCGGCTCAACAGCCACGCCCGTCACCTTGCCATCGACCACTTCAAGTTCGTTGGCg             >  5:88026/2‑148 (MQ=255)
                                                                                                                                             gCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATCGACAATCGGCTCAACAGCCACGCCCGTCACCTTGCCATCGACCACTTCAAGTTCGTTGGCGAAcaca       <  3:301334/148‑1 (MQ=255)
                                                                                                                                               ccACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATCGACAATCGGCTCAACAGCCACGCCCGTCACCTTGCCATCGACCACTTCAAGTTCGTTGGCGAACACATAg    >  7:97622/1‑149 (MQ=255)
                                                                                                                                                 aCGGCGATGGTCTGCTCAAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATCGACAATCGGCTCAACAGCCACGCCCGTCACCTTGCCATCGACCACTTCAAGTTCGTTGGCGAACACATAGTc  >  1:6214/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
ACGCAGGCCCAGCAGGTACAGCACCCCGTCCAGCCCCAAGGTGGAAATCGCCTGCTTGGCCGACTGGCGCACCAGCGGCTTGGCACGAAACGCCACGCCCAACCCAGCCAGCGACAGCATCGGCAGGTCATTGGCGCCATCGCCCACGGCGATGGTCTGCTCCAGTTGCAAGCCTTCGTCGCTGGCCAGCTTCTTCAACAAATCGGCCTTGCGCTGGGCATCGACAATCGGCTCAACAGCCACGCCCGTCACCTTGCCATCGACCACTTCAAGTTCGTTGGCGAACACATAGTC  >  NC_002947/5579638‑5579931

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