Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,177,554 A→T 7.3% K146M (AAG→ATG) ‡ guaA → glutamine‑hydrolyzing GMP synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,177,5540AT 7.3% 163.7 / 3.3 55K146M (AAG→ATG) ‡guaAglutamine‑hydrolyzing GMP synthase
Reads supporting (aligned to +/- strand):  ref base A (26/25);  new base T (2/2);  total (28/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 = 8.33e-01

CAAGGTAGAAGGTTCCGACCTGCGTGAGTTCGGCTATGCCCGTGTGGACGTTGTCGGCAA‑GAGCCGCCTGCTCGACGGCATCGAAGACCACGTTGATGACGACGGCGTGCTGGGCCTGGACGTGTGGATGAGCCACGGTGACAAGGTCACCCAGATGCCTGGCAACTTCAACATTCTGGCCAGCACCCCAAGCT‑GCCCGATCGCCGGCATGTACGACGACGTCCGCGGCTACTACGGCGTGCAGTTCCACCCCGAAGTGACCCACACCAAGCAAGGTGGTCGTATCCTGTC  >  NC_002947/1177411‑1177701
                                                                                                                                                |                                                                                                                                                    
cAAGGTAGAAGGTTCCGACCTGCGTGAGTTCGGCTATGCCCGTGTGGACGTTGTCGGCAA‑GAGCCGCCTGCTCGACGGCATCGAAGACCACGTTGATGACGACGGCGTGCTGGGCCTGGACGTGTGGATGAGCCACGGTGACAAGGTCa                                                                                                                                                 <  8:276151/149‑1 (MQ=255)
   ggTAGAAGGTTCCGACCTGCGTGAGTTCGGCTATCCCCGTGTGGACGTTGTCGGCAA‑GAGCCGCCTGCTCGACGGCATCGAAGACCACGTTGATGACGACGGCGTGCTGGGCCAGGACGTGTGCATGAGCCACGGTGACAAGGTCAcc                                                                                                                                               >  5:257888/1‑148 (MQ=255)
                 aCCTGCGTGAGTTCGGCTATGCCCGTGTGGACGTTGTCGGCAA‑GAGCCGCCTGCTCGACGGCATCGAAGACCACGTTGATGACGACGGCGTGCTGGGCCTGGACGTGTGGATGAGCCACGGTGACAAGGTCACCCAGATGCCTGGCAAc                                                                                                                                <  5:249353/149‑1 (MQ=255)
                   cTGCGTGAGTTCGGCTATGCCCGTGTGGACGTTGTCGGCAAGGAGCCGCCTGCTCGACGGCATCGAAGACCACGTTGATGACGACGGCGTGCTGGGCCTGGACGTGTGGATGAGCCACGGTGACAAGGTCACCCAGATGCCTGGCAACt                                                                                                                               >  1:243100/1‑149 (MQ=255)
                   cTGCGTGAGTTCGGCTATGCCCGTGTGGACGTTGTCGGCAA‑GAGCCGCCTGCTCGACGGCATCGAAGACCACGTTGATGACGACGGCGTGCTGGGCCTGGACGTGTGGATGAGCCACGGTGACAAGGTCACCCAGATGCCTGGCAACtt                                                                                                                              <  2:106645/149‑1 (MQ=255)
                    tGCGTGAGTTCGGCTATGCCCGTGTGGACGTTGTCGGCAA‑GAGCCGACTCCTCCACGGAATCGAAGACCACTTAGATGACGACGGCGTGCTGGGCCTATACGTGTGTATGAGACAAGGTGACAAGGTCACCAAGATCCCTGGCACCta                                                                                                                              >  8:82974/1‑147 (MQ=255)
                     gCGTGAGTTCGGCTATGCCCGTGTGGACGTTGTCGGCAA‑GAGCCGCCTGCTCGACGGCATCGAAGACCACGTTGATGACGACGGCGTGCTGGGCCTGGACGTGTGGATGAGCCACGGTGACAAGGTCACCCAGATGCCTGGCAACTTc                                                                                                                             >  5:309592/1‑148 (MQ=255)
                      cGTGAGTTCGGCTATGCCCGTGTGGACGTTGTCGGCAA‑GAGCCGCCTGCTCGACGGCATCGAAGACCACGTTGATGACGACGGCGTGCTGGGCCTGGACGTGTGGATGAGCCACGGTGACAAGGTCACCCAGATGCCTGGCAACTTCaa                                                                                                                           <  6:309592/149‑1 (MQ=255)
                       gTGAGTTCGGCTATGCCCGTGTGGACGTTGTCGGCAA‑GAGCCGCCTGCTCGACGGCATCGAAGACCACGTTGATGACGACGGCGTGCTGGGCCTGGACGTGTGGATGAGCCACGGTGACAAGGTCACCCAGATGCCTGGCAACTTCAAc                                                                                                                          <  2:243100/149‑1 (MQ=255)
                            ttCGGCTATGCCCGTGTGGACGTTGTCGGCAA‑GAGCCGCCTGCTCGACGGCATCGAAGACCACGTTGATGACGACGGCGTGCTGGGCCTGGACGTGTGGATGAGCCACGGTGACCTTGTCAccgtg                                                                                                                                            >  7:113633/1‑123 (MQ=255)
                            ttCGGCTATGCCCGTGTGGACGTTGTCGGCAA‑GAGCCGCCTGCTCGACGGCATCGAAGACCACGTTGATGACGACGGCGTGCTGGGCCTGGACGTGTGGATGAGCCACGGTGACCTTGTCAccgtg                                                                                                                                            <  8:113633/126‑4 (MQ=255)
                                         gtgtGGACGTTGTCGGCAA‑GAGCCGCCTGCTCGACGGCATCGAAGACCACGTTGATGACGACGGCGTGCTGGGCCTGGACGTGTGGATGAGCCACGGTGACAAGGTCACCCAGATGCCTGGCAACTTCAACATTCTGGCCAGCACCCCa                                                                                                        <  7:202666/149‑1 (MQ=255)
                                                  ttGTCGGCAA‑GAGCCGCCTGCTCGACGGCATCGAAGACCACGTTGATGACGACGGCGTGCCGGGCCTGGACGTGTGGATGAGCCACGGTGACAAGGTCACCCAGATGCCTGGCAACTTCAACATTCTGGCCAGCACCCCAAGCT‑Gccc                                                                                               >  4:101470/1‑148 (MQ=255)
                                                        gCAA‑GAGCCGCCTGCTCGACGGCATCGAAGACCACGTTGATGACGACGGCGTGCTGGGCCTGGACGTGTGGATGAGCCACGGTGACCTTGTCAccgt                                                                                                                                             >  5:14739/1‑95 (MQ=255)
                                                        gCAA‑GAGCCGCCTGCTCGACGGCATCGAAGACCACGTTGATGACGACGGCGTGCTGGGCCTGGACGTGTGGATGAGCCACGGTGACCTTGTCAccgt                                                                                                                                             <  6:14739/97‑3 (MQ=255)
                                                        gCAA‑GAGCCGCCTGCTCGACGGCATCGAAGACCACGTTGATGACGACGGCGTGCTGGGCCTGGACGTGTGGATGAGCCACGGTGACAAGGTCACCCAGATGCCTGGCAACTTCAACa                                                                                                                         >  7:161397/1‑117 (MQ=255)
                                                        gCAA‑GAGCCGCCTGCTCGACGGCATCGAAGACCACGTTGATGACGACGGCGTGCTGGGCCTGGACGTGTGGATGAGCCACGGTGACAAGGTCACCCAGATGCCTGGCAACTTCAACa                                                                                                                         <  8:161397/117‑1 (MQ=255)
                                                        gCAA‑GAGCCGCCTGCTCGACGGCATCGAAGACCACGTTGATGACGACGGCGTGCTGGGCCTGGACGTGTGGATGAGCCACGGTGACAAGGTCACCCAGATGCCTGGCAACTTCAACATTCTGGCCAGCACCCCAAGCT‑GCCCGATCg                                                                                          <  3:12745/147‑1 (MQ=255)
                                                        gCAA‑GAGCCGCCTGCTCGACGGCATCGAAGACCACGTTGATGACGACGGCGTGCTGGGCCTGGACGTGTGGATGAGCCACGGTGACAAGGTCACCCAGATGCCTGGCAACTTCAACATTCTGGCCAGCACCCCAAGCT‑GCCCGATCg                                                                                          >  4:12745/1‑147 (MQ=255)
                                                          aa‑GAGCCGCCTGCTCGACGGCATCGAAGACCACGTTGATGACGACGGCGTGCCGGGCCTGGACGTGTGGATGAGCCACGGTGACAAGGTCACCCAGATGCCTGGCAACTTCAACATTCTGGCCAGCACCCCAAGCT‑GCCCGATCGCCgg                                                                                      <  3:101470/149‑1 (MQ=255)
                                                                        tCGACGGCATCGAAGACCACGTTGATGACGACGGCGTGCTGGGCCTGGACGTGTGGATGAGCCACGGTGACAAGGTCACCCAGATGCCt                                                                                                                                      >  1:75287/1‑89 (MQ=255)
                                                                        tCGACGGCATCGAAGACCACGTTGATGACGACGGCGTGCTGGGCCTGGACGTGTGGATGAGCCACGGTGACAAGGTCACCCAGATGCCt                                                                                                                                      <  2:75287/89‑1 (MQ=255)
                                                                        tCGACGGCATCGAAGACCACGTTGATGACGACGGCGTGCTGGGCCTGGACGTGTGGATGAGCCACGGTGACAAGGTCACCCAGATGCCTGGCAACTTCAACATTCTGGc                                                                                                                  >  6:110327/1‑109 (MQ=255)
                                                                        tCGACGGCATCGAAGACCACGTTGATGACGACGGCGTGCTGGGCCTGGACGTGTGGATGAGCCACGGTGACAAGGTCACCCAGATGCCTGGCAACTTCAACATTCTGGCCAGCACCCCAAGCT‑GCCCGATCGCCGGCATGTacgacgac                                                                         >  6:237654/1‑149 (MQ=255)
                                                                         cGACGGCATCGAAGACCACGTTGATGACGACGGCGTGCTGGGCCTGGACGTGTGGATGAGCCACGGTGACAAGGTCACCCAGATGCCTGGCAACTTCAACATTCTGGCCAGCACCCCAAGCT‑GCCCGATCGCCGGCATGTacgacgacg                                                                        >  5:88867/1‑149 (MQ=255)
                                                                          gACGGCATCGAAGACCACGTTGATGACGACGGCGTGCTGGGCCTGGACGTGTGGATGAGCCACGGTGACAAGGTCACCCAGATGCCTGGCAACTTCAACATTCTGGCCAGCACACCAAGCT‑GCCCGATCGCCGGCATGTACGACGACGt                                                                       >  2:257144/1‑149 (MQ=255)
                                                                           aCGGCATCGAAGTCCACGTTGATGACGACGTCGTGCTGGGCCTGGACGTGTGGATGAGCCACGGTGACAAGGTCACCCAGATGCCTGGCAACTTCAACATTCTGGCCAGCACCCCAAGCT‑GCCCGATCGCCGGCATGTACGACGACGTc                                                                      <  7:82974/149‑1 (MQ=255)
                                                                               cATCGAAGACCACGTTGATGACGACGGCGTGCTGGGCCTGGACGTGTGGATGAGCCACGGTGACAAGGTCACCCAGATGCCTGGCAACTTCAACATTCTGGCCAGCACCCCAAGCT‑GCCCGATCGCCGGCATGTACGACGACGTCcgcg                                                                  >  7:291604/1‑149 (MQ=255)
                                                                                       aCCACGTTGATGACGACGGCGTGCTGGGCCTGGACGTGTGGATGAGCCACGGTGACAAGGTCACCCAGATGCCTGGCAACTTCAACATTCTg                                                                                                                    <  8:245562/92‑1 (MQ=255)
                                                                                       aCCACGTTGATGACGACGGCGTGCTGGGCCTGGACGTGTGGATGAGCCACGGTGACAAGGTCACCCAGATGCCTGGCAACTTCAACATTCTg                                                                                                                    >  7:245562/1‑92 (MQ=255)
                                                                                          aCGTTGATGACGACGGCGTGCTGGGCCTGGACGTGTGGATGAGCCACGGTGACAAGGTCACCCAGATGCCTGGCAACTTCAACATTCTGGCCAGCACCCCAAGCT‑GCCCGATCGCCg                                                                                       <  4:262169/117‑1 (MQ=255)
                                                                                          aCGTTGATGACGACGGCGTGCTGGGCCTGGACGTGTGGATGAGCCACGGTGACAAGGTCACCCAGATGCCTGGCAACTTCAACATTCTGGCCAGCACCCCAAGCT‑GCCCGATCGCCg                                                                                       >  3:262169/1‑117 (MQ=255)
                                                                                                 tgaCGACGGCGTGCTGGGCCTGGACGTGTGGATTAGCCACGGTGACAAGGTCACCCAGATGCCTGGCAACTTCAACATTCTGGCCAGCACCCCAAGCTCGCCCGATCGCCGGCATGTACGACGACGTCCGCGGCTACTACGGCGTGCa                                                  <  1:257144/148‑1 (MQ=255)
                                                                                                 tgaCGACGGCGTGCTGGGCCTGGACGTGTGGATGAGCCACGGTGACAAGGTCACCCAGATGCCTGGCAACTTCAACATTCTGGCCAGCACCCCAAGCT‑GCCCGATCGCCGGCATGTACGACGACGTCCGCGGctac                                                             >  1:254865/1‑136 (MQ=255)
                                                                                                 tgaCGACGGCGTGCTGGGCCTGGACGTGTGGATGAGCCACGGTGACAAGGTCACCCAGATGCCTGGCAACTTCAACATTCTGGCCAGCACCCCAAGCT‑GCCCGATCGCCGGCATGTACGACGACGTCCGCGGctac                                                             <  2:254865/136‑1 (MQ=255)
                                                                                                 tgaCGACGGCGTGCTGGGCCTGGACGTGTGGATGAGCCACGGTGACAAGGTCACCCAGATGCCTGGCAACTTCAACATTCTGGCCAGCACCCCAAGCT‑GCCCGATCGCCGGCATGTACGACGACGTCCGCGGCTACTACGGCGTGCAGt                                                <  8:291604/149‑1 (MQ=255)
                                                                                                 tgaCGACGGCGTGCTGGGCCTGGACGTGTGGATGAGCCACGGTGACAAGGTCACCCAGATGCCTGGCAACTTCAACATTCTGGCCAGCACCCCAAGCT‑GCCCGATCGCCGGCATGTACGACGACGTCCGCGGCTACTACGGCGTGCAGt                                                >  3:121384/1‑149 (MQ=255)
                                                                                                   acgacgGCGTGCTGGGCCTGGACGTGTGGATGAGCCACGGTGACAAGGTCACCCAGATGCCTGGCAACTTCAACATTCTGGCCAGCACCCCAAGCT‑GCCCGATCGCCGGCATGTACGACGACGTCCGCGGCTACTACGGCGTGCAGTTc                                              <  2:215958/149‑1 (MQ=255)
                                                                                                     gacgGCGTGCTGGGCCTGGACGTGTGGATGAGCCACGGTGACAAGGTCACCCAGATGCCTGGCAACTTCAACATTCTGGCCAGCACCCCAAGCT‑GCCCGATCGCCGGCATGTACGACGACGTCcgc                                                                   <  4:259531/126‑1 (MQ=255)
                                                                                                     gacgGCGTGCTGGGCCTGGACGTGTGGATGAGCCACGGTGACAAGGTCACCCAGATGCCTGGCAACTTCAACATTCTGGCCAGCACCCCAAGCT‑GCCCGATCGCCGGCATGTACGACGACGTCcgc                                                                   >  3:259531/1‑126 (MQ=255)
                                                                                                            tGCTGGGCCTGGACGTGTGGATGAGCCACGGTGACAAGGTCACCCAGATGCCTGGCAACTTCAACATTCTGGCCAGCACCCCAAGCT‑GCCCGATCGCCGGCATGTACGACGACGTCCGCGGCTACTACGGCGTGCAGTTCCACCCCGaa                                     <  4:290822/149‑1 (MQ=255)
                                                                                                            tGCTGGGCCTGGACGTGTGGATGAGCCACGGTGACAAGGTCACCCAGATGCCTGGCAACTTCAACATTCTGGCCAGCACCCCAAGCT‑GCCCGATCGCCGGCATGTACGACGACGTCCGCGGCTACTACGGCGTGCAGTTCCACCCCGaa                                     <  6:88867/149‑1 (MQ=255)
                                                                                                               tGGGCCTGGACGTGTGGATGAGCCACGGTGACAAGGTCACCCAGATGCCTGGCAACTTCAACATTCTGGCCAGCACCCCAAGCT‑GCCCGATCGCCGGCATGTACGACGACGTCCGCGGCTACTACGGCGTGCAGTTCCACCCCGAAGt                                   >  7:209017/1‑148 (MQ=255)
                                                                                                                 ggCCTGGACGTGTGGATGAGCCACGGTGACAAGGTCACCCAGATGCCTGGCAACTTCAACATTCTGGCCAGCACCCCAAGCT‑GCCCGATCGCCGGCATGTACGACGACGTCCGCGGCTACTACGGCGTGCAGTTCCACCCCGAAGTGa                                 <  5:237654/148‑1 (MQ=255)
                                                                                                                 ggCCTGGACGTGTGGATGAGCCACGGTGACAAGGTCACCCAGATGCCTGGCAACTTCAACATTCTGGCCAGCACCCCAAGCT‑GCCCGATCGCCGGCATGTACGACGACGTCCGCGGCTACTACGGCGTGCAGTTCCACCCCGAAGTGAc                                >  1:54871/1‑149 (MQ=255)
                                                                                                                       gACGTGTGGATGAGCCACGGTGACAAGGTCACCCAGATGCCTGGCAACTTCAACATTCTGGCCAGCACCCCAAGCT‑GCCCGATCGCCGGCATGTACGACGACGTCAGCGGCTACTACGGCGTGCAGTTCCACCCCGAAGTGACCcac                            <  8:209017/147‑1 (MQ=255)
                                                                                                                             tgGATGAGCCACGGTGACAAGGTCACCCAGATGCCTGGCAACTTCAACATTCTGGCCAGCACCCAAAGCT‑GCCCGATCGCCGGCATGTACGACGACGTCCGCGGCTACTACGGCGTGCAGTTCCAc                                           <  1:322537/126‑1 (MQ=255)
                                                                                                                             tgGATGAGCCACGGTGACAAGGTCACCCAGATGCCTGGCAACTTCAACATTCTGGCCAGCACCCAAAGCT‑GCCCGATCGCCGGCATGTACGACGACGTCCGCGGCTACTACGGCGTGCAGTTCCAc                                           >  2:322537/1‑126 (MQ=255)
                                                                                                                                    gCCACGGTGACAAGGTCACCCAGATGCCTGGCAACTTCAACATTCTGGCCAGCACCCCAAGCT‑GCCCGATCGCCGGCATGTACGACGACGTCCGCGGCTACTACGGCGTGCAGTTCCACCCCGAAGTGACCCACACCAAGCAAggtggt             >  7:163849/1‑149 (MQ=255)
                                                                                                                                      cACGGTGACAAGGTCACCCAGATGCCTGGCAACTTCAACATTCTGGCCAGCACCCCAAGCT‑GCCCGATCGCCGGCATGTACGACGACGTCCGCGGCTACTACGGCGTGCAGTTCCACCCCGAAGTGACCCACACCAAGCAAGGTGGTCg           >  5:118165/1‑149 (MQ=255)
                                                                                                                                         ggTGACAAGGTCACCCAGATGCCTGGCAACTTCAACATTCTGGCCAGCACCCCAAGCT‑GCCCGATCGCCGGCATGTACGACGACGTCCGCGGCTACTACGGCGTGCAGTTCCACCCCGaa                                     >  3:188163/1‑120 (MQ=255)
                                                                                                                                         ggTGACAAGGTCACCCAGATGCCTGGCAACTTCAACATTCTGGCCAGCACCCCAAGCT‑GCCCGATCGCCGGCATGTACGACGACGTCCGCGGCTACTACGGCGTGCAGTTCCACCCCGaa                                     <  4:188163/120‑1 (MQ=255)
                                                                                                                                             aCAAGGTCACCCAGATGCCTGGCAACTTCAACATTCTGGCCAGCACCCCAAGCT‑GCCCGATCGCCGGCATGTACGACGACGTCCGCGGCTACTACGGCGTGCAGTTCCACCCCGAAGTGa                                 <  3:304192/120‑1 (MQ=255)
                                                                                                                                             aCAAGGTCACCCAGATGCCTGGCAACTTCAACATTCTGGCCAGCACCCCAAGCT‑GCCCGATCGCCGGCATGTACGACGACGTCCGCGGCTACTACGGCGTGCAGTTCCACCCCGAAGTGa                                 >  4:304192/1‑120 (MQ=255)
                                                                                                                                              cAAGGTCACCCAGATGCCTGGCAACTTCAACATTCTGGCCAGCACCCCAAGCT‑GCCCGATCGCCGGCATGTACGACGACGTCCGCGGCTACTACGGCGTGCAGTTCCACCCCGAAGTGACCCACACCAAGCAAGGTGGTCGTATCCTGt   >  8:267209/1‑149 (MQ=255)
                                                                                                                                               aaGGTCACCCAGATGCCTGGCAACTTCAACATTCTGGCCAGCACCCCAAGCT‑GCCCGATCGCCGGCATGTACGACGACGTCCGCGGCTACTACGGCGTGCAGTTCCACCCCGAAGTGACCCACACCAAGCAAGGTGGTCGTATCCTGTc  <  4:121384/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
CAAGGTAGAAGGTTCCGACCTGCGTGAGTTCGGCTATGCCCGTGTGGACGTTGTCGGCAA‑GAGCCGCCTGCTCGACGGCATCGAAGACCACGTTGATGACGACGGCGTGCTGGGCCTGGACGTGTGGATGAGCCACGGTGACAAGGTCACCCAGATGCCTGGCAACTTCAACATTCTGGCCAGCACCCCAAGCT‑GCCCGATCGCCGGCATGTACGACGACGTCCGCGGCTACTACGGCGTGCAGTTCCACCCCGAAGTGACCCACACCAAGCAAGGTGGTCGTATCCTGTC  >  NC_002947/1177411‑1177701

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