Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 855,327 G→A 8.0% A60T (GCC→ACC)  murI → glutamate racemase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947855,3270GA 8.0% 142.4 / 2.2 51A60T (GCC→ACC) murIglutamate racemase
Reads supporting (aligned to +/- strand):  ref base G (22/24);  new base A (2/2);  total (24/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.06e-01

GACTCGGGGGTTGGCGGATTGTCGGTGCTGGCCGAGATCCAGCGCCTGCTGCCCAACGAAACGTTGCTGTATGTGGGCGATTGCGGGCATATACCCTACGGTGAGAAGTCGCCGGACTATATCCGCGAGCGCTGCCGACGCATTGCCGGGTTCTTCCATGAACAAGGCGCAAAGGCCATGGTCCTGGCCTGCAACACTGCCACGGTGGCGGCAGTGGCCGACCTGCGTGAGTTGTATCCGACCTGGCCGTTGGTGGGCATGGAACCGGCCGTCAAGCCTGCAGCGGCAG  >  NC_002947/855183‑855471
                                                                                                                                                |                                                                                                                                                
gACTCGGGGGTTGGCGGATTGTCGGTGCTGGCCGAGATCCAGCGCCTGCTGCCCAACGAAACGTTGCTGTATGTGGGCGATTGCGGGCATATACCCTACGGTGAGAAGTCGCCGGACTATATCCGCGAGCGCTGCCGACGCATTGCCgg                                                                                                                                              <  6:239874/149‑1 (MQ=255)
gACTCGGGGGTTGGCGGATTGTCGGTGCTGGCCGAGATCCAGCGCCTGCTGCCCAACGAAACGTTGCTGTATGTGGGCGATTGCGGGCATATACCCTACGGTGAGAAGTCGCCGGACTATATCCGCGAGCGCTGCCGACGCATTGCCgg                                                                                                                                              <  2:153587/149‑1 (MQ=255)
        ggTTGGCGGATTGTCGGTGCTGGCCGAGATCCAGCGCCTGCTGCCCAACGAAACGTTGCTGTATGTGGGCGATTGCGGGCATATACCCTACGGTGAGAAGTCGCCGGACTATATCCGCGAGCGCTGCCGACGCATTGCCGGGTTCTTcc                                                                                                                                      <  8:45562/149‑1 (MQ=255)
         gTTGGCGGATTGTCGGTGCTGGCCGAGATCCAGCGCCTGCTGCCCAACGAAACGTTGCTGTATGTGGGCGATTGCGGGCATATACCCTACGGTGAGAAGTCGCCGGACTATATCCGCGAGCGCTGCCGACGCATTGCTGGGTTCTTCCa                                                                                                                                     <  1:197700/149‑1 (MQ=255)
         gTTGGCGGATTGTCGGTGCTGGCCGAGATCCAGCGCCTGCTGCCCAACGAAACGTTGCTGTATGTGGGCGATTGCGGGCATATACCCTACGGTGAGAAGTCGCCGGACTATATCCGCGAGCGCTGCCGACGCATTGCCGGGTTCTTCCa                                                                                                                                     >  3:64566/1‑149 (MQ=255)
               ggATTGTCGGTGCTGGCCGAGATCCAGCGCCTGCTGCCCAACGAAACGTTGCTGTATGTGGGCGATTGCGGGCATATACCCTACGGTGAGAAGTCGCCGGACTATATCCGCGAGCGCTGCCGACGCATTGCCGGGTTCTTCCAt                                                                                                                                    <  1:175825/144‑1 (MQ=255)
               ggATTGTCGGTGCTGGCCGAGATCCAGCGCCTGCTGCCCAACGAAACGTTGCTGTATGTGGGCGATTGCGGGCATATACCCTACGGTGAGAAGTCGCCGGACTATATCCGCGAGCGCTGCCGACGCATTGCCGGGTTCTTCCAt                                                                                                                                    >  2:175825/1‑144 (MQ=255)
                  ttGTCGGTGCTGGCCGAGATCCAGCGCCTGCTGCCCAACGAAACGTTGCTGTATGTGGGCGATTGCGGGCATATACCCTACGGTGAGAAGTCGCCGGACTATATCCGCGAGCGCTGCCGACGCATTGCCGGGTTCTTCCATGAACAAgg                                                                                                                            <  1:63018/149‑1 (MQ=255)
                     tCGGTGCTGGCCGAGATCCAGCGCCTGCTGCCCAACGAAACGTTGCTGTATGTGGGCGATTGCGGGCATATACCCTACGGTGAGAAGTCGCCGGACTATATCCGCGAGCGCTGCCGACGCATTGCCGGGTTCTTCCATGAACAAGgcg                                                                                                                          >  4:227730/1‑148 (MQ=255)
                               cctcgaTCCAGAGCCTGCTGCCCCACGAAACGTTGCTGTATGTGGGCGATTGCGGGCATATACCCTACGGTGAGAAGTCGCCGGACTATATCCGCGAGCGCTGCCGACGCATTGCCGGGTTCTTCCATGAACAAGGCGCAAAGGCCATg                                                                                                               <  6:233764/145‑1 (MQ=255)
                               ccGAGATCCAGCGCCTGCTGCCCAACGAAACGTTGCTGTATGTGGGCGATTGCGGGCATATACCCTACGGTGAGAAGTCGCCGGACTATATCCGCGAGCGCTGCCGACGCATTGCCGGGTTCTTCCATGAACAAGGCGCAAAGGCCATg                                                                                                               >  4:176876/1‑149 (MQ=255)
                                cGAGATCCAGCGCCTGCTGCCCAACGAAACGTTGCTGTATGTGGGCGATTGCGGGCATATACCCTACGGTGAGAAGTCGCCGGACTATATCCGCGAGCGCTGCCGACGCATTGCCGGGTTCTTCCATGAACAAGGCGCAAAGGCCATgg                                                                                                              >  8:219414/1‑149 (MQ=255)
                                cGAGATCCAGCGCCTGCTGCCCAACGAAACGTTGCTGTATGTGGGCGATTGCGGGCATATACCCTACGGTGAGAAGTCGCCGGACTATATCCGCGAGCGCTGCCGACGCATTGCCGGGTTCTTCCATGAACAAGGCGCAAAGGCCATgg                                                                                                              >  4:256149/1‑149 (MQ=255)
                                cGAGATCCAGCGCCTGCTGCCCAACGAAACGTTGCTGTATGTGGGCGATTGCCGGCATATACCCTACGGTGAGAAGTCGCCGGACTATATCCGCGAGCGCTGCCGACGCATTGCCGGGTTCTTCCATGAACAAGGCGCAAAGGCCATgg                                                                                                              >  1:205007/1‑149 (MQ=255)
                                      ccAGCGCCTGCTGCCCAACGAAACGTTGCTGTATGTGGGCGATTGCGGGCATATACCCTACGGTGAGAAGTCGCAGGACTATATCCGCGAGCGCTGCCGACGCATTGCCGGGTTCTTCCATGAACAAGGCGCAAAGGCCATGGTCCTgg                                                                                                        >  3:182307/1‑149 (MQ=255)
                                       cAGCGCCTGCTGCCCAACGAAACGTTGCTGTATGTGGGCGATTGCGGGCATATACCCTACGGTGAGAAGTCGCCGGACTATATCCGCGAGCGCTGCCGACGCATTGCCGGGTTCTTCCATGAACAAGGCGCAAAGGCCATGGTCCTGGc                                                                                                       <  7:50125/149‑1 (MQ=255)
                                             ctgctgCCCAACGAAACGTTGCTGTATGTGGGCGATTGCGGGCATATACCCTACGGTGAGAAGTCGCCGGACTATATCCGCGAGCGCTGCCGACGCATTGCCGGGTTCTTCCATGAACAAGGCGCAAAGGCCATGGTCCTg                                                                                                         >  6:126186/1‑141 (MQ=255)
                                             ctgctgCCCAACGAAACGTTGCTGTATGTGGGCGATTGCGGGCATATACCCTACGGTGAGAAGTCGCCGGACTATATCCGCGAGCGCTGCCGACGCATTGCCGGGTTCTTCCATGAACAAGGCGCAAAGGCCATGGTCCTg                                                                                                         <  5:126186/141‑1 (MQ=255)
                                                     cAACGAAACGTTGCTGTATGTGGGCGATTGCGGGCATATACCCTACGGTGAGAAGTCGCCGGACTATATCCGCGAGCGCTGCCGACGCATTACCGGGTTCTTCCATGAACAAGGCGCAAAGGCCATGGTCCTGGCCTGCAACACTGCCa                                                                                         >  1:238924/1‑149 (MQ=255)
                                                     cAACGAAACGTTGCTGTATGTGGGCGATTGCGGGCATATACCCTACGGTGAGAAGTCGCCGGACTATATCCGAGAGCGCTGCCGACGCATTACCGGGTTCTTCCATGAACAAGGCGCAAAGGCCATGGTCCTGGCCTGCAACACTGCCa                                                                                         >  1:234934/1‑149 (MQ=255)
                                                            aCGTTGCTGTATGTGGGCGATTGCGGGCATATACCCTACGGTGAGAAGTCGCCGGACTATATCCGCGAGCGCTGCCGACGCATTGCCGGGTTCTTCCATGAACAAGGCGCAAAGGCCATGGTCCTGGCCTGCAACACTGCCACGGTggc                                                                                  >  6:137355/1‑149 (MQ=255)
                                                                    gTATGTGGGCGATTGCGGGCATATACCCTACGGTGAGAAGTCGCCGGACTATATCCGCGAGCGCTGCCGACGCATTGCCGGGTTCTTCCATGAACAAGGCGCAAAGGCCATGGTCCTGGCCTGCAACACTGCCACGGTGGCGGCAGTgg                                                                          >  6:204485/1‑149 (MQ=255)
                                                                       tgtgGGCGATTGCGGGCATATACCCTACGGTGAGAAGTCGCCGGACTATATCCGCGAGCGCTGCCGACGCATTGCCGGGTTCTTCCATGAACAAGGCGCAAAGGCCATGGTCCTGGCCTGCAACACTGCCACGGTGGCGGCAGTGGcc                                                                        >  1:136470/1‑148 (MQ=255)
                                                                                  gCGGGCATATACCCTACGGTGAGAAGTCGCCGGACTATATCCGCGAGCGCTGCCGACGCATTGCCGGGTTCTTCCATGAACAAGGCGCAAAGGCCATGGTCCTGGCCTGCAACACTGCCACGGTGGCGGCAGTGGCCGACCTGCGTGa                                                             <  7:219414/148‑1 (MQ=255)
                                                                                  gCGGGCATATACCCTACGGTGAGAAGTCGCCGGACTATATCCGCGAGCGCTGCCGACGCATTGCCGGGTTCTTCCATGAACAAGGCGCAAAGGCCATGGTCCTGGCCTGCAACACTGCCACGGTGGCGGCAGTGGCCGACCTGCGTGa                                                             <  3:256149/148‑1 (MQ=255)
                                                                                    gggCATATACCCTACGGTGAGAAGTCGCCGGACTATATCCGCGAGCGCTGCCGACGCATTGCCGGGTTCTTCCATGAACAAGGCGCCAAGGCCATGGTCCTGGCCTGCAACACTGCCACGGTGGCGGCAGTGGCCGACCTGCGTGAGtt                                                          <  4:64566/149‑1 (MQ=255)
                                                                                    gggCATATACCCTACGGTGAGAAGTCGCCGGACTATATCCGCGAGCGCTGCCGACGCATTGCCGGGTTCTTCCATGAACAAGGCGCAAAGGCCATGGTCCTGGCCTGCAACACTGCCACGGTGGCGGCAGTGGCCGACCTGCGTGAGtt                                                          <  5:204485/149‑1 (MQ=255)
                                                                                          ataCCCTACGGTGAGAAGTCGCCGGACTATATCCGCGAGCGCTGCCGACGCATTGCCGGGTTCTTCCATGAACAAGGCGCAAAGGCCATGGTCCTGGCCTGCAACACTGCCACGGTGGCGGCAGTGGCCGACCTGCGTGAGTTGTATcc                                                    >  3:92940/1‑149 (MQ=255)
                                                                                            aCCCTACGGTGAGAAGTCGCCGGACTATATCCGCGAGCGCTGCCGACGCATTGCCGGGTTCTTCCATGAACAAGGCGCAAAGGCCATGGTCCTGGCCTGCAACACTGCCACGGTGGCGGCAGTGGCCGACCTGCGTGAGTTGTATCCGa                                                  >  1:206840/1‑149 (MQ=255)
                                                                                              ccTACGGTGAGAAGTCGCCGGACTATATCCGCGAGCGCTGCCGACGCATTGCCGGGTTCTTCCATGAACAAGGCGCAAAGGCCATGGTCCTGGCCTGCAACACTGCCACGGTGGCGGCAGTGGCCGACCTGCGt                                                               <  7:95544/134‑1 (MQ=255)
                                                                                              ccTACGGTGAGAAGTCGCCGGACTATATCCGCGAGCGCTGCCGACGCATTGCCGGGTTCTTCCATGAACAAGGCGCAAAGGCCATGGTCCTGGCCTGCAACACTGCCACGGTGGCGGCAGTGGCCGACCTGCGt                                                               >  8:95544/1‑134 (MQ=255)
                                                                                               cTACGGTGAGAATTCGCCGGACTATATCCGCGAGCGCTGCCGACGCATTACCGGGTTCTTCCATGAACAAGGCGCAAAGGCCATGGTCCTGGCCTGCAACACTGCCACGGTGGCGGCAGTGGCCGACCTGCGTGAGTTGTATCCGACCt                                               <  2:238924/149‑1 (MQ=255)
                                                                                               cTACGGTGAGAAGTCGCCGGACTATATCCGCGAGCGCTGCCGACGCATTACCGGGTTCTTCCATGAACAAGGCGCAAAGGCCATGGTCCTGGCCTGCAACACTGCCACGGTGGCGGCAGTGGCCGACCTGCGTGAGTTGTATCCGACCt                                               <  2:234934/149‑1 (MQ=255)
                                                                                                tACGTTGAGAAGTCGCCGGACTATATCCGCGAGCGCTGCCGACGCATTGCCGGGTTCTTCCATGAACAAGGCGCAAAGGCCATGGTc                                                                                                            <  2:81061/87‑1 (MQ=255)
                                                                                                       agaAGTCGCCGGACTATATCCGCGAGCGCTGCCGACGCATTGCCGGGTTCTTCCATGAACAAGGCGCAAAGGCCATGGTCCTGGCCTGCAACACTGCCACGGTGGCGGCAGTGGCCGACCTGCGTGAGTTGTATCCGACCTGGCCGTtg                                       >  7:50321/1‑149 (MQ=255)
                                                                                                       agaAGTCGCCGGACTATATCCGCGAGCGCTGCCGACGCATTGCCGGGTTCTTCCATGAACAAGGCGCAAAGGCCATGGTCCTGGCCTGCAACACTGCCACGGTGGCGGCAGTGGCCGACCTGCGTGAGTTGTATCCGACCTGGCCGTtg                                       <  6:41797/149‑1 (MQ=255)
                                                                                                        gaAGTCGCCGGACTATATCCGCGAGCGCTGCCGACGCATTTCCGGGTTCTTCCATGAACAAGGCGCAAAGGCCATGGTCCTGGCCTGCAACACTGCCACGGTGGCGGCAGTGGCCGACCTGCGTGAGTTGTATCCGACCTGGCCGTtg                                       <  4:182307/148‑1 (MQ=255)
                                                                                                         aaGTCGCCGGACTATATCCGCGAGCGCTGCCGACGCATTGCCGGGTTCTTCCATGAAAAAGGCGCAAAGGCCATGGTCCTGGCCTGCAACACTGCCACGGTGGCGGCAGTGGCCGACCTGCGTGAGTTGTATCCGACCTGGCCGTtggt                                     >  6:13809/1‑149 (MQ=255)
                                                                                                               ccGGACTATATCCGCGAGCGCTGCCGACGCATTGCCGGGTTCTTCGATGAACAAGGCGCAAAGGCCATGGTCCTGGCCTGCAACACTGGCACGGTGGCGGCAGTGGCCGACCTGCGTGAGTTGTATCCGACCTGGCCGTTGGTGGGCAt                               <  2:205007/149‑1 (MQ=255)
                                                                                                                cGGACTATATCCGCGAGCGCTGCCGACGCATTGCCGGGTTCTTCCATGAACAAGGCGCAAAGGCCATGGTCCTGGCCTGCAACACTGCCACGGTGGCGGCAGTGGCCGACCTGCGTGAGTTGTATCCGACCTGGCCGTTGGTGGGCATg                              <  3:227730/149‑1 (MQ=255)
                                                                                                                   aCTATATCCGCGAGCGCTGCCGACGCATTGCCGGGTTCTTCCATGAACAAGGCGCAAAGGCCATGGTCCTGGCCTGCAACACTGCCACGGTGGCGGCAGt                                                                            >  1:197036/1‑100 (MQ=255)
                                                                                                                   aCTATATCCGCGAGCGCTGCCGACGCATTGCCGGGTTCTTCCATGAACAAGGCGCAAAGGCCATGGTCCTGGCCTGCAACACTGCCACGGTGGCGGCAGt                                                                            <  2:197036/100‑1 (MQ=255)
                                                                                                                       tatCCGCGAGCGCTGCCGACGCATTGCCGGGTTCTTCCATGAACAAGGCGCAAAGGCCATGGTCCTGGCCTGCAACACTGCCACGGTGGCGGCAGTGGCCGACCTGCGTGAGTTGTATCCGACCTGGCCGTTGGTGGGCATGGAACCgg                       >  7:208394/1‑149 (MQ=255)
                                                                                                                              gAGCGCTGCCGACGCATTGCCGGGTTCTTCCATGAACAAGGCGCAAAGGCCATGGTCCTGGCCTGCAACACTGCCACGGTGGCGGCAGTGGCCGACCTGCGTGAGTTGTATCCGACCTGGCCGTTGGTGGGCATGGAACCGGCCGTCaa                <  2:136470/149‑1 (MQ=255)
                                                                                                                                gcgcTGCCGACGCATTGCCGGGTTCTTCCATGAACAAGGCGCAAAGGCCATGGTCCTGGCCTGCAACACTGCCACGGTGGCGGCAGTGGCCGACCTGCGTGAGTTGTATCCGACCTGGCCGTTGGTGGGCATGGAACCGGCCGTCAAGc              <  8:50321/149‑1 (MQ=255)
                                                                                                                                gcgcTGCCGACGCATTGCCGGGTTCTTCCATGAACAAGGCGCAAAGGCCATGGTCCTGGCCTGCAACACTGCCACGGTGGCGGCAGTGGCCGACCTGCGTGAGTTGTATCCGACCTGGCCGTTGGTGGGCATGGAACCGGCCGTCAAGc              <  5:137355/149‑1 (MQ=255)
                                                                                                                                    tGCCGACGCATTGCCGGGTTCTTCCATGAACAAGGCGCAAAGGCCATGGTCCTGGCCTGCAACACTGCCACGGTGGCGGCAGTGGCCGACCTGCGTGAGTTGTATCCGACCTGGCCGTTGGTGGGCATGGAACCGGCCGTCAAGCCTg           >  4:243856/1‑148 (MQ=255)
                                                                                                                                          cGCATTGCCGGGTTCTTCCATGAACAAGGCGCAAAGGCCATGGTCCTGGCCTGCAACACTGCCACGGTGGCGGCAGTGGCCGACCTGCGTGAGTTGTATCCGACCTGGCCGTTGGTGGGCATGGAACCGGCCGTCAAGCCTGCAGCGGc    <  2:206840/149‑1 (MQ=255)
                                                                                                                                            cATTGCCGGGTTCTTCCATGAACAAGGCGCAAAGGCCATGGTCCTGGCCTGCAACACTGCCACGGTGGCGGCAGTGGCCGACCTGCGTGAGTTGTATCCGACCTGGCCGTTGGTGGGCATGGAACCGGCCGTCAAGCCTGCAGCGGCAg  >  3:280962/1‑149 (MQ=255)
                                                                                                                                               tGCCGGGTTCTTCCATGAACAAGGCGCAAAGGCCATGGTCCTGGCCTGCAACACTGCCACGGTGGCGGCAGTGGCCGACCTGCGTGAGTTGTATCCGACCTGGCCGTTGGTGGGCATGGAACCgg                       <  6:180996/125‑1 (MQ=255)
                                                                                                                                               tGCCGGGTTCTTCCATGAACAAGGCGCAAAGGCCATGGTCCTGGCCTGCAACACTGCCACGGTGGCGGCAGTGGCCGACCTGCGTGAGTTGTATCCGACCTGGCCGTTGGTGGGCATGGAACCgg                       >  5:180996/1‑125 (MQ=255)
                                                                                                                                                |                                                                                                                                                
GACTCGGGGGTTGGCGGATTGTCGGTGCTGGCCGAGATCCAGCGCCTGCTGCCCAACGAAACGTTGCTGTATGTGGGCGATTGCGGGCATATACCCTACGGTGAGAAGTCGCCGGACTATATCCGCGAGCGCTGCCGACGCATTGCCGGGTTCTTCCATGAACAAGGCGCAAAGGCCATGGTCCTGGCCTGCAACACTGCCACGGTGGCGGCAGTGGCCGACCTGCGTGAGTTGTATCCGACCTGGCCGTTGGTGGGCATGGAACCGGCCGTCAAGCCTGCAGCGGCAG  >  NC_002947/855183‑855471

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