Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,591,521 C→T 8.0% A146V (GCG→GTG)  PP_2269 → N‑acetylmuramoyl‑L‑alanine amidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,591,5210CT 8.0% 141.6 / 2.6 50A146V (GCG→GTG) PP_2269N‑acetylmuramoyl‑L‑alanine amidase
Reads supporting (aligned to +/- strand):  ref base C (22/24);  new base T (2/2);  total (24/26)
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.06e-01

AACTTCACGCCCGAGCAGTTTGCAAGTCTCAAGCACCTGCTGGGCGAGCTGAAAGAGAAATACCCCAGCGCGACAATCCAAGGCCACCGGGACTTCCCGAAAGTAGCCAAGGCTTGCCCGAGCTTCGACGTTAAACCGTGGGTAGCGGCCAACTTATAAGAGTCGTGGGCCCTCAGTTGGATATGGCAGCCGCCACGTCTATCACTCAACTCTGAGAAAGGGCTGGCCTTCGGGCTGGCCCTCTTTATTGCCCATGAAAAATCGAAAGCATGAAACAGTGAGTGAGTTC  >  NC_002947/2591376‑2591664
                                                                                                                                                 |                                                                                                                                               
aacTTCACGCCCGAGCAGTTTGCAAGTCTCAAGCACCTGCTGGGCGAGCTGAAAGAGAAATACCCCAGCGCGACAATCCAAGGCCACCGGGACTTCCCGAAAGTAGCCAAGGCTTGCCCGAGCTTCGACGTTAAACCGTGGGTAGCGGc                                                                                                                                              >  2:44395/1‑149 (MQ=255)
      aCGCCCGAGCAGTTTGCAAGTCTCAAGCACCTGCTGGGCGAGCTGAAAGAGAAATACCCCAGCGCGACAATCCAAGGCCACCGGGACTTCCCGAAAGTAGCCAAGGCTTGCTCGAGCTTCGACGTTAAACCGTGGGTAGCGGCCAACtt                                                                                                                                        <  4:361229/149‑1 (MQ=255)
      aCGCCCGAGCAGTTTGCAAGTCTCAAGCACCTGCTGGGCGAGCTGAAAGAGAAATACCCCAGCGCGACAATCCAAGGCCACCGGGACTTCCCGAAAGTAGCCAAGGCTTGCCCGAGCTTCGACGTTAAACCGTGGGTAGCGGCCAACt                                                                                                                                         >  3:220433/1‑148 (MQ=255)
        gCCCGAGCAGTTTGCAAGTCTCAAGCACCTGCTGGGCGAGCTGAAAGAGAAATACCCCAGCGCGACAATCCAAGGCCACCGGGACTTCCCGAAAGTAGCCAAGGCTTGCCCGAGCTTCGACGTTAAACCGTGGGTAGCGGCCAACTtat                                                                                                                                      <  3:260297/149‑1 (MQ=255)
        gCCCGAGCAGTTTGCAAGTCTCAAGCACCTGCTGGGCGAGCTGAAAGAGAAATACCCCAGCGCGACAATCCAAGGCCACCGGGACTTCCCGAAAGTAGCCAAGGCTTGCCCGAGCTTCGACGTTAAACCGTGGGTAGCGGCCAACTta                                                                                                                                       >  1:283910/1‑148 (MQ=255)
          ccGAGCAGTTTGCAAGTCTCAAGCACCTGCTGGGCGAGCTGAAAGAGAAATACCCCAGCGCGACAATCCAAGGCCACCGGGACTTCCCGAAAGTAGCCAAGGCTTGCCCGAGCTTCGACGTTAAACCGTGGGTAGCGGCCAACTTATaa                                                                                                                                    <  3:105111/149‑1 (MQ=255)
          ccGAGCAGTTTGCAAGTCTCAAGCACCTGCTGGGCGAGCTGAAAGAGAAATACCCCAGCGCGACAATCCAAGGCCACCGGGACTTCCCGAAAGTAGCCAAGGCTTGCCCGAGCTTCGACGTTAAACCGTGGGTAGCGGCCAACTTATaa                                                                                                                                    <  6:254923/149‑1 (MQ=255)
                 gTTTGCAAGTCTCAAGCACCTGCTGGGCGAGCTGAAAGAGAAATACCCCAGCGCGACAATCCAAGGCCACCGGGACTTCCCGAAAGTAGCCAAGGCTTGCCCGAGCTTCGACGTTAAACCGTGGGTAGCGGCCAACTTATAAGAGTCGt                                                                                                                             <  1:45991/149‑1 (MQ=255)
                  tttGCAAGTCTCAAGCACCTGCTGGGCGAGCTGAAAGAGAAATACCCCAGCGCGACAATCCAAGGCCACCGGGACTTCCCGAAAGTAGCCAAGGCTTGCCCGAGCTTCGACGTTAAACCGAGGGTAGCGGCCAACTTATAAGAGTCGTg                                                                                                                            >  7:268544/1‑149 (MQ=255)
                            tcAAGCACCTGCTGGGCGAGCTGAAAGAGAAATACCCCAGCGCGACAATCCAAGGCCACCGGGACTTCCCGAAAGTAGCCAAGGCTTGCCCGAGCTTCGACGTTAAACCGTGGGTAGCGGCCAACTTATAAGAGTCGTGGGCCCTCAGt                                                                                                                  >  5:340346/1‑149 (MQ=255)
                             cAAGCACCTGCTGGGCGAGCTGAAAGAGAAATACCCCAGCGCGACAATCCAAGGCCACCGGGACTTCCCGAAAGTAGCCAAGGCTTGCCCGAGCTTCGACGTTAAACCGTGGGTAGCGGCCAACTTATAAGAGTCGTGGGCCCTCAGtt                                                                                                                 <  4:220433/149‑1 (MQ=255)
                              aaGCACCTGCTGGGCGAGCTGAAAGAGAAATACCCCAGCGCGACAATCCAAGGCCACCGGGACTTCCCGAAAGTAGCCAAGGCTTGCCCGAGCTTCGACGTTAAACCGTGGGTAGCGGCCAACTTATAAGAGTCGTGGGCCCTCAg                                                                                                                   <  8:277028/146‑1 (MQ=255)
                              aaGCACCTGCTGGGCGAGCTGAAAGAGAAATACCCCAGCGCGACAATCCAAGGCCACCGGGACTTCCCGAAAGTAGCCAAGGCTTGCCCGAGCTTCGACGTTAAACCGTGGGTAGCGGCCAACTTATAAGAGTCGTGGGCCCTCAGtt                                                                                                                 <  2:283910/148‑1 (MQ=255)
                                gCACCTGCTGGGCGAGCTGAAAGAGAAATACCCCAGCGCGACAATCCAAGGCCACCGGGACTTCCCGAAAGTAGCCAAGGCTTGCCCGAGCTTCGACGTTAAACCGTGGGTAGCGGCCAACTTATAAGAGTCGTGGGCCCTCAGTTGGa                                                                                                              >  5:117517/1‑149 (MQ=255)
                                gCACCTGCTGGGCGAGCTGAAAGAGAAATACCCCAGCGCGACAATCCAAGGCCACCGGGACTTCCCGAAAGTAGCCAAGGCTTGCCCGAGCTTCGACGTTAAACCGTGGGTAGCGGCCAACTTATAAGAGTCGTGGGCCCTCAGTTGGa                                                                                                              >  3:69595/1‑149 (MQ=255)
                                    ctgctgGGCGAGCTGAAAGAGAAATACCCCAGCGCGACAATCCAAGGCCACCGGGACTTCCCGAAAGTAGCCAAGGCTTGCCCGAGCTTCGACGTTAAACCGTGGGTAGCGGCCAACTta                                                                                                                                       <  3:176339/120‑1 (MQ=255)
                                    ctgctgGGCGAGCTGAAAGAGAAATACCCCAGCGCGACAATCCAAGGCCACCGGGACTTCCCGAAAGTAGCCAAGGCTTGCCCGAGCTTCGACGTTAAACCGTGGGTAGCGGCCAACTta                                                                                                                                       >  4:176339/1‑120 (MQ=255)
                                     tgctgGGCGAGCTGAAAGAGAAATACCCCAGCGCGACAATCCAAGGCCACCGGGACTTCCCGAAAGTAGCCAAGGCTTGCCCGAGCTTCGACGTTAAACCGTGGGTAGCGGCCAACTTATAAGAGTCGTGGGCCCTCAGTTGGATATgg                                                                                                         >  3:327694/1‑149 (MQ=255)
                                                 tGAAAGAGAAATACCCCAGCGCGACAATCCAAGGCCACCGGGACTTCCCGAAAGTAGCCAAGGCTTGCCCGAGCTTCGACGTTAAACCGTGGGTAGCGGCCAACTTATAAGAGTCGTGGGCCCGCAGTTGGATATGGCAGCCGCCACg                                                                                              >  6:263599/1‑148 (MQ=255)
                                                   aaaGTGAAATACCCCAGCGGGACAATCGAAGGCCACCGGGACTTCCTTATAGTAGCCAAGGCGTGCCGGAGCTTCGACGTTAAACCGTGGGTAGTGGCCAACTTATAAGAGTCGTGGGCCCTCAGTTTGATATGGCAGCCGCCACGt                                                                                             <  8:383851/147‑1 (MQ=255)
                                                   aaaGAGAAATACCCCAGCGCGACAATCCAAGGCCACCGGGACTTCCCGAAAGTAGCCAAGGCTTGCCCGAGCTTCGACGTTAAACCGTGGGTAGTGGCCAACTTATAAGAGTCGTCGGCCCTCAGTTGGATATGGCAgcagcccagg                                                                                             >  7:383851/1‑143 (MQ=255)
                                                               cccAGCGCGACAATCCAAGGCCACCGGGACTTCCCGAAAGTAGCCAAGGCTTGCCCGAGCTTCGACGTTAAACCGTGGGTAGCGGCCAACTTATAAGAGTCGTGGGCCCTCAGTTGGATATGGCAGCCGCCACGTCTATCACTCAActc                                                                               >  7:377025/1‑149 (MQ=255)
                                                                 cAGCGCGACAATCCAAGGCCACCGGGACTTCCCGAAAGTAGCCAAGGCTTGCCCGAGCTTCGACGTTAAACCGTGGGTAGCGGCCAACTTATAAGAGTCGTGGGCCCTCAGTTGGATATGGCAGCCGCCACGTCTATCACTCAActct                                                                              >  6:388958/1‑148 (MQ=255)
                                                                   gcgcGACAATCCAAGGCCACCGGGACTTCCCGAAAGTAGCCAAGGCTTGCCCGAGCTTCGACGTTAAACCGTGGGTAGCGGCCAACTTATAAGAGTCGTGGGCCCTCAGTTGGATATGGCAGCCGCCACGTCTATCACTCAACTCTgag                                                                           <  6:117517/149‑1 (MQ=255)
                                                                       gACAATCCAAGGCCACCGGGACTTCCCGAAAGTAGCCAAGGCTTGCCCGAGCTTCGACGTTAAACCGTGGGTAGTGGCCAACTTATAAGAGTCGTGGGCCCTCAGTTGGATATGGCAGCCGCCACGTCTATCACTCAACTCTGAGAAAg                                                                       >  6:155304/1‑149 (MQ=255)
                                                                          aaTCCAAGGCCACCGGGACTTCCCGAAAGTAGCCAAGGCTTGCCCGAGCTTCGACGTTAAACCGTGGGTAGCGGCCAACTTATAAGCGTCGTGGGCCCTCAGTTGGATATGGCAGCCGCCACGTCTATCACTCAACTCTGAGAAAGGGc                                                                    <  4:278715/149‑1 (MQ=255)
                                                                          aaTCCAAGGCCACCGGGACTTCCCGAAAGTAGCCAAGGCTTGCCCGAGCTTCGACGTTAAACCGTGGGTAGCGGCCAACTTATAAGAGTCGTGGGCCCTCAGTTGGATATGGCAGCCGCCACGTCTATCACTCAACTCTGAGAAAGGGc                                                                    >  1:359560/1‑149 (MQ=255)
                                                                             ccAAGGCCACCGGGACTTCCCGAAAGTAGCCAAGGCTTGCCCGAGCTTCGACGTTAAACCGTGGGTAGCGGCCAACTTATAAGAGTCGTGGGCCCTCAGTTGGATATGGCAGCCGCCACGTCTATCACTCAACTCTGAGAAAGGGCTgg                                                                 >  3:173300/1‑149 (MQ=255)
                                                                                   ccACCGGGACTTCCCGAAAGTAGCCAAGGCTTGCCCGAGCTTCGACGTTAAACCGTGGGTAGCGGCCAACTTATAAGAGTCGTGGGCCCTCAGTTGGATATGGCAGCCGCCACGTCTATCACTCAACTCTGAGAAAGGGCTGGCCTTCg                                                           <  4:69595/149‑1 (MQ=255)
                                                                                   ccACCGGGACTTCCCGAAAGTAGCCAAGGCTTGCCCGAGCTTCGACGTTAAACCGTGGGTAGCGGCCAACTTATAAGAGTCGTGGGCCCTCAGTTGGATATGGCAGCCGCCACGTCTATCACTCAACTCTGAGAAAGGGCTGGCCTTCg                                                           <  7:90511/149‑1 (MQ=255)
                                                                                       cGGGACTTCCCGAAAGTAGCCAACGCTTGCCCGAGCTTCGACGTTAAACCGTGGGTAGCGGCCAACTTATAAGAGTCGTGGGCCCTCAGGTGGATATGGCAGCCGCCACGTCTATCACTCAACTCTGAGAAAGGGCTGGCCTTCGGGCt                                                       <  6:222964/149‑1 (MQ=255)
                                                                                         ggACTTCCCGAAAGTAGCCAAGGCTTGCCCGAGCTTCGACGTTAAACCGTGGGTAGCGGCCAACTTATAAGAGTCGTGGGCCCTCAGTTGGATATGGCAGCCGCCACGTCTATCACTCAACTCTGAGAAAGGGCTGGCCTTCGGGCTgg                                                     <  1:44395/149‑1 (MQ=255)
                                                                                               cccGAAAGTAGCCAAGGCTTGCCCGAGCTTCGACGTTAAACCGTGGGTAGCGGCCAACTTATAAGAGTCGTGGGCCCTCAGTTGGATATGGCAgccgcc                                                                                                 <  4:110034/99‑1 (MQ=255)
                                                                                               cccGAAAGTAGCCAAGGCTTGCCCGAGCTTCGACGTTAAACCGTGGGTAGCGGCCAACTTATAAGAGTCGTGGGCCCTCAGTTGGATATGGCAgccgcc                                                                                                 >  3:110034/1‑99 (MQ=255)
                                                                                                    aaGTAGCCAAGGCTTGCCCGAGCTTCGACGTTAAACCGTGGGTAGCGGCCAACTTATAAGAGTCGTGGGCCCTCACTTGGATATGGCAGCCGCGACGTCTATCACTCAACAGTGAGAAAGGGCTGGCCTTCGGGCTGGCCCTCTTTAtt                                          >  2:284681/1‑149 (MQ=255)
                                                                                                          ccAAGGCTTGCCCGAGCTTCGACGTTAAACCGTGAGTAGCGGCCAACTTATAAGAGTCGTGGGCCCTCAGTTGGATATGGCAGCCGCCACGTCTATCACTCAACTCTGAGAAAGGGCTGGCCTTCGGGCTGGCCCTCTTTATTGCCCAt                                    >  1:53035/1‑149 (MQ=255)
                                                                                                            aaGGCTTGCCCGAGCTTCGACGTTAAACCGTGGGTAGCGGCCAACTTATAAGAGTCGTGGGCCCTCAGTTGGATATGGCAGCCGCCACGTCTATCACTCAACTCTGAGAAAGGGCTGGCCTTCGGGCTGGCCCTCTTTATTGCCCATGa                                  <  5:388958/149‑1 (MQ=255)
                                                                                                                     ccGAGCTTCGACGTTAAACCGTGGGTAGCGGCCAACTTATAAGAGTCGTGGGCCCTCAGTTGGATATGGCAGCCGCCACGTCTATCACTCAACTCTGAGAAAGGGCTGGCCTTCGGGCTGGCCCTCTTTATTGCCCATGAAAAATCGa                          >  4:343928/1‑148 (MQ=255)
                                                                                                                      cGAGCTTCGACGTTAAACCGTGGGTAGCGGCCAACTTATAAGAGTCGTGGGCCCTCAGTTGGATATGGCAGCCGCCACGTCTATCACTCAACTCTGAGAAAGGGCTGGCCTTCGGGCTGGCCCTCTTTATTGCCCATGAAAAATCGaaa                        <  2:359560/149‑1 (MQ=255)
                                                                                                                       gAGCTTCGACGTTAAACCGTGGGTAGCGGCCAACTTATAAGAGTCGTGGGCCCTCAGTTGGATATGGCAGCCGCCACGTCTATCACTCAACTCTGAGAAAGGGCTGGCCTTCGGGCTGGCCCTCTTTATTGCCCATGAAAAATCGAAAg                       <  3:343928/149‑1 (MQ=255)
                                                                                                                        aGCTTCGACGTTAAACCGTGGGTAGCGGCCAACTTATAAGAGTCGTGGGCCCTCAGTTGGATATGGCAGCCGCCACGTCTATCACTCAACTCTGAGAAAGGGCTGGCCTTCGGGCTGGCCCTCTTTATTGCCCATGAAAAATCGAAAGc                      <  6:340346/149‑1 (MQ=255)
                                                                                                                           ttCGACGTTAAACCGTGGGTAGCGGCCAACTTATAAGAGTCGTGGGCCCTCAGTTTGATATGGCAGCCGCCACGTCTATCACTCAACTATGAGAAAGGGCTGGCCTTCGGGCTGGCCCTCTTTATTGccc                                      >  5:3301/1‑130 (MQ=255)
                                                                                                                           ttCGACGTTAAACCGTGGGTAGCGGCCAACTTATAAGAGTCGTGGGCCCTCAGTTGGATATGGCAGCCGCCACGTCTATCACTCAACTCTGAGAAAGGGCTGGCCTTCGGGCTGGCCCTCTTTATTGccc                                      <  6:3301/130‑1 (MQ=255)
                                                                                                                              gACGTTAAACCGTGGGTAGTGGCCAACTTATAAGAGTCGTGGGCCCTCAGTTGGATATGGCAGCCGCCACGTCTATCACTCAACTCTGAGAAAGGGCTGGCCTTCGGGCTGGCCCTCTTTATTGCCCATGAAAAATCGAAAGCATGaaa                <  5:155304/149‑1 (MQ=255)
                                                                                                                              gACGTTAAACCGTGGGTAGCGGCCAACTTATAAGAGTCGTGGGCCCTCAGTTGGATATGGCAGCCGCCACGTCTATCACTCAACTCTGAGAAAGGGCTGGCCTTCGGGCTGGCCCTCTTTATTGCCCATGAAAAATCGAAAGCATGaaa                <  5:263599/149‑1 (MQ=255)
                                                                                                                               aCGTTTAACCGTGGGTAGCGGCCAACTTATAAGAGTCGTGGGCCCTCAGTTGGATATGGCAGCCGCCACGTCTATCACTCAACTCTGAGAAAGGGCTGGCCTTCGGGCTGGCCCTCTTTATTGCCCATGAAAAATCGAAAGCATGaaa                >  6:289893/1‑148 (MQ=255)
                                                                                                                                 gTTAAACCGTGGGTAGCGGCCAACTTATAAGAGTCGTGGGCCCTCAGTTGGATATGGCAGCCGCCACGTCTATCACTCAACTCTGAGAAAGGGCTGGCCTTCGGGCAGGCCCTCTTGATTGCCCATGAAAAATCGAAAGCAGGAAACag             >  5:362367/1‑149 (MQ=255)
                                                                                                                                  ttAAACCGTGGGTAGCGGCCAACTTATAAGAGTCGTGGGCCCTCAGTTGGATATGGCAGCCGCCACGTCTATCACTCAACTCTGAGAAAGGGCTGGCCTTCGGGCTGGCCCTCTTTATTGCCCATGAAAAATCGAAAGCATGAAACagt            <  2:53035/149‑1 (MQ=255)
                                                                                                                                     aaCCGTGGGTAGCGGCCAACTTATAAGAGTCGTGGGCCCTCAGTTGGATATGGCAGCCGCCACGTCTATCACTCAACTCTGAGAAAGGGCTGGCCTTCGGGCTGGCCCTCTTTATTGCCCATGAAAAATCGAAAGCATGAAACagtga          <  8:377025/148‑1 (MQ=255)
                                                                                                                                            ggTAGCGGCCAACTTATAAGAGTCGTGGGCCCTCAGTTGGATATGGCAGCCGCCACGTCTATCACTCAACTCTGAGAAAGGGCTGGCCTTCGGGCTGGCCCTCTTTATTGCCCATGAAAAATCGAAAGCATGAAACAGTGAGTGAGTTc  >  7:298654/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
AACTTCACGCCCGAGCAGTTTGCAAGTCTCAAGCACCTGCTGGGCGAGCTGAAAGAGAAATACCCCAGCGCGACAATCCAAGGCCACCGGGACTTCCCGAAAGTAGCCAAGGCTTGCCCGAGCTTCGACGTTAAACCGTGGGTAGCGGCCAACTTATAAGAGTCGTGGGCCCTCAGTTGGATATGGCAGCCGCCACGTCTATCACTCAACTCTGAGAAAGGGCTGGCCTTCGGGCTGGCCCTCTTTATTGCCCATGAAAAATCGAAAGCATGAAACAGTGAGTGAGTTC  >  NC_002947/2591376‑2591664

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