Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,372,084 A→C 12.9% intergenic (+36/+76) gdhB → / ← PP_2081 NAD‑specific glutamate dehydrogenase/phosphoenolpyruvate synthase regulatory protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,372,0840AC12.9% 117.8 / 8.2 53intergenic (+36/+76)gdhB/PP_2081NAD‑specific glutamate dehydrogenase/phosphoenolpyruvate synthase regulatory protein
Reads supporting (aligned to +/- strand):  ref base A (26/20);  new base C (3/4);  total (29/24)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

GACGACCTGCGCAATGCCACGGGCACTGACTATGCGATGTACGCGGTGGCCAACCGCGAGCTGGTCGACCTGGCCATGAGCGGGCAGGCAGCGGTGGTGCCGTCCTGAGCCTTGTAGTGAGGTAAAACAAAGCCCCGGCAGTGATGCCGGGGCTTTTTTTTGTATGCCTGGGCGGGCCATTTCGCGGGTGAACCCGCTCCTACAGGGGGGGCGGGGCTGTTACTTGAACCGTCTCTCCACGCCTTTTTCGACCAGGATCTTGGCGGAAATCTCTTCCACCGAAAAATGCGT  >  NC_002947/2371941‑2372231
                                                                                                                                               |                                                                                                                                                   
gacgacCTGCGCAATGCCACGGGCACTGACTATGCGATGTACGCGGTGGCCAACCGCGAGCTGGTCGACCTGGCCATGAGCGGGCAGGCAGCGGTGGTGCCGTCCTGAGCCTTGTAGTGAGGTAAAACAAAGCCCCGGCAGTGATGCCg                                                                                                                                                >  8:34681/1‑149 (MQ=255)
gacgacCTGCGCAATGCCACGGGCACTGACTATGCGATGTACGCGGTGGCCAACCGCGAGCTGGTCGACCTGGCCATGAGCGGGCAGGCAGCGGTGGTGCCGTCCTGAGCCTTGTAGTGAGGTAAAACAAAGCCCCGGCAGTGATGCCg                                                                                                                                                >  8:298612/1‑149 (MQ=255)
   gacCTGCGCAATGCCACGGGCACTGACTATGCGATGTACGCGGTGGCCAACCGCGAGCTGGTCGACCTGGCCATGAGCGGGCAGGCAGCGGTGGTGCCGTCCTGAGCCTTGTAGTGAGGTAAAACAAAGCCCCGGCAGTGATGCCgggg                                                                                                                                             >  6:2809/1‑149 (MQ=255)
         cgcAATGCCACGGGCACTGACTATGCGATGTACGCGGTGGCCAACCGCGAGCTGGTCGACCTGGCCATGAGCGGGCAGGCAGCTTTGGTGCCGTCCTGAGCCTTGTAGTGAGGTAAAACAAAGCCCCGGCATCACTGCCGGGGCttt                                                                                                                                         <  1:159242/147‑1 (MQ=255)
         cgcAATGCCACGGGCACTGACTATGCGATGTACGCGGTGGCCAACCGCGAGCTGGTCGACCTGGCCATGAGCGGGCAGGCAGCGGTGGTGCCGTCCTGAGCCTTGTAGTGAGGTAAAACAAAGCCCCGGCATCACTGCCGGGGCttt                                                                                                                                         <  7:326643/147‑1 (MQ=255)
            aaTGCCACGGGCACTGACTATGCGATGTACGCGGTGGCCAACCGCGAGCTGGTCGACCTGGCCATGAGCGGGCAGGCAGCGGTGGTGCCGTCCTGAGCCTTGTAGTGAGGTAAAACAAAGCCCCGGCAGTGATGCCGGGGCtttttttt                                                                                                                                    <  6:301047/149‑1 (MQ=255)
                 cACGGGCACTGACTATGCGATGTACGCGGTGGCCAACCGCGAGATGGTCGACCTGGCCATGAGCGGGCAGGCAGCGGTGGTGCCGTCCTGAGCCTTGTAGTGAGGTAAAACAAAGCACCGGCATCACTGCCGGGGCttg                                                                                                                                         >  2:153467/1‑138 (MQ=255)
                 cACGGGCACTGACTATGCGATGTACGCGGGGGCCACCCTCGAGCTGGTCGACCTGGCCATGAGCGGGCAGGCAGCGGTGGTGCCGTCCTGAGCCTTGTAGTGAGGGAAAACAAAGCCCCGGCATCACTGCCGGGGCttt                                                                                                                                         <  1:153467/139‑1 (MQ=255)
                  aCGGGCACTGACTATGCGATGTACGCGGTGGCCAACCGCGAGCTGGTCGACCTGGCCATGAGCGGGCAGGCAGCGGTGGTGCCGTCCTGAGCCTTGTAGTGAGGTAAAACAAAGCCCCGGCAGTGATGCCGGGGCTTTTTTTTGTATGc                                                                                                                              >  6:222329/1‑149 (MQ=255)
                      gCACTGACTATGCGATGTACGCGGTGGCCAACCGCGAGCTGGTCGACCTGGCCATGAGCGGGCAGGCAGCGGTGGTGCCGTCCTGAGCCTTGTAGTGAGGTAAAACAAAGCCCCGGCAGTGATGCCGGGGCTTTTTTTTGTATGCCTgg                                                                                                                          >  4:5135/1‑149 (MQ=255)
                       cACTGACTATGCGATGTACGCGGTGGCCAACCGCGAGCTGGTCGACCTGGCCATGAGCGGGCAGGCAGCGGTGGTGCCGTCCTGAGCCTTGTAGTGAGGTAAAACAAAGCCCCGGCAGTGATGCCGGGGCTTTTTTTTGTATGCCTggg                                                                                                                         >  2:301276/1‑149 (MQ=255)
                         cTGACTATGCGATGTACGCGGTGGCCAACCGCGAGCTGGTCGACCTGGCCATGAGCGGGCAGGCAGCGGTGGTGCCGTCCTGAGCCTTGTAGTGAGGTAAAACAAAGCCCCGGCAGTGATGCCGGGGCTTTTTTTTGTATGCCTgggcg                                                                                                                       <  4:359558/149‑1 (MQ=255)
                         cTGACTATGCGATGTACGCGGTGGCCAACCGCGAGCTGGTCGACCTGGCCATGAGCGGGCAGGCAGCGGTGGTGCCGTCCTGAGCCTTGTAGTGAGGTAAAACAAAGCCCCGGCAGTGATGCCGGGGCTTTTTTTTGTATGCCTgggcg                                                                                                                       <  5:90364/149‑1 (MQ=255)
                            aCTATGCGATGTACGCGGTGGCCAACCGCGAGCTGGTCGACCTGGCCATGCGCGGGCAGGCAGCGGTGGTGCCGTCCTGAGCCTTGTAGTGAGGTAAAACAAAGCCCCGGCAGTGATGCCGGGGCTTTTTTTTGTATGCCTgggcgg                                                                                                                      <  3:5135/147‑1 (MQ=255)
                                tGCGATGTACGCGGTGGCCAACCGCGAGCTGGTCGACCTGGCCATGAGCGGGCAGGCAGCGGTGGTGCCGTCCTGAGCCTTGTAGTGAGGTAAAACAAAGCCCCGGCAGTGATGCCGGGGCTTTTTTTTGTATGCCTGGGCGGGCCAtt                                                                                                                >  6:108007/1‑149 (MQ=255)
                                   gATGTACGCGGTGGCCAACCGCGAGCTGGTGGACCTGGCCATGAGCGGGCAGGCAGCGGTGGTGCCGTCCTGAGCCTTGTAGTGAGGTAAAACAAAGCCCCGGCAGTGATGCCGGGGCTTTTTTTTGTATGCCTGGGCGGGCCATTTcg                                                                                                             >  1:32303/1‑149 (MQ=255)
                                    cagctcgcgGTTGGCCACCGCGAGCTGGTCGACCTGGCCATGAGCGGGCAGGCAGCGGTGGTGCCGTCCTGAGCCTTGTAGTGAGGTAAAACAAAGCCCCGGCAGTGATGCCGGGGCTTTTTTTTGTATGCCTGGGCGGGCCATTTcgc                                                                                                            >  1:109191/6‑149 (MQ=255)
                                    aTGTACGCGGTGGCCAACCGCGAGCTGGTCGACCTGGCCATGAGCGGGCAGGCAGCGGTGGTGCCGTCCTGAGCCTTGTAGTGAGGTAAAACAAAGCCCCGGCAGTGATGCCGGGGCTTTTTTTTGTATGCCTGGGCGGGCCATTTcgc                                                                                                            <  7:34681/149‑1 (MQ=255)
                                       tACGCGGTGGCCAACCGCGAGCTGGTCGACCTGGCCATGAGCGGGCAGGCAGCGGTGGTGCCGTCCTGAGCCTTGTAGTGAGGTAAAACAAAGCCCCGGCAGTGATGCCGGGGCTTTTTTTTGTATGCCTGGGCGGGCCATTTCGCggg                                                                                                         >  4:79401/1‑149 (MQ=255)
                                       tACGCGGTGGCCAACCGCGAGCTGGTCGACCTGGCCATGAGCGGGCAGGCAGCGGTGGTGCCGTCCTGAGCCTTGTAGTGAGGTAAAACAAAGCCCCGGCAGTGATGCCGGGGCTTTTTTTTGTATGCCTGGGCGGGCCATTTCGCgg                                                                                                          <  5:2809/148‑1 (MQ=255)
                                                  cAACCGCGAGCTGGTCGACCTGGCCATGAGCGGGCAGGCAGCGGTGGTGCCGTCCTGAGCCTTGTAGTGAGGTAAAACAAAGCCCCGGCAGTGATGCCGGGGCTTTTTTTTGTATGCCTTGGCGGGCCATTTCGCGGGTGAACCCGCTc                                                                                              <  1:239700/149‑1 (MQ=255)
                                                    aCCGCGAGCTGGTCGACCTGGCCATGAGCGGGCAGGCAGCGGTGGTGCCGTCCTGAGCCTTGTAGTGAGGTAAAACAAAGCCCCGGCAGTGATGCCGGGGCTTTTTTTTGTATGCCTGGGCGGGCCATTTCGCGGGTGAACCCGCTCCt                                                                                            <  5:222329/149‑1 (MQ=255)
                                                             tGGTCGACCTGGCCATGAGCGGGCAGGCAGCGGTGGTGCCGTCCTGAGCCTTGTAGTGAGGTAAAACAAAGCCCCGGCAGTGATGCCGGGGCTTTTTTTTGTATGCCTGGGCGGGCCATTTCGCGGGTGAACCCGCTCCTACAgggggg                                                                                   >  5:312914/1‑149 (MQ=255)
                                                                  gACCTGGCCATGAGCGGGCAGGCAGCGGTGGTGCCGTCCTGAGCCTTGTAGTGAGGTAAAACAAAGCCCCGGCAGTGATGCCGGGGCTTTTTTTTGTATGCCTGGGCGGGCCATTTCGCGGGTGAACCCGCTCCTACAGGGggggcgg                                                                               <  3:79401/148‑1 (MQ=255)
                                                                   aCCTGGCCATGAGCGGGCAGGCAGCGGTGGTGCCGTCCTGAGCCTTGTAGTGAGGTAAAACAAAGCCCCGGCAGTGATGCCGGGGCTTTTTTTTGTATGCCTGGGCGGGCCATTTCGCGGGTGAACCCGCTCCTAc                                                                                          <  4:406447/136‑1 (MQ=255)
                                                                   aCCTGGCCATGAGCGGGCAGGCAGCGGTGGTGCCGTCCTGAGCCTTGTAGTGAGGTAAAACAAAGCCCCGGCAGTGATGCCGGGGCTTTTTTTTGTATGCCTGGGCGGGCCATTTCGCGGGTGAACCCGCTCCTAc                                                                                          >  3:406447/1‑136 (MQ=255)
                                                                          cATGAGCGGGCAGGCAGCGGTGGTGCCGTCCTGAGCCTTGTAGTGAGGTAAAACAAAGCCCCGGCAGTGATGCCGGGGCTTTTTTTTGTATGCCTGGGCGGGCCATTTCGCGGGTGACCCCGCTCCTACAGGGGGGGCGGGGCTGTTAc                                                                      <  2:109191/149‑1 (MQ=255)
                                                                          cATGAGCGGGCAGGCAGCGGGGGGGCCGTCCTGAGCCTTGTAGTGAGGTAAAACAAAGCCCCGGCAGTGATGCCGGGGCTTTTTTTTGTATGCCTGGGCGGGCCATTTCGCGGGTGAACCCGCTCCTACAGGGGGGGCGGGGCTGTTAc                                                                      <  7:112551/149‑1 (MQ=255)
                                                                            tGAGCGGGCAGGCAGCGGTGGTGCCGTCCTGAGCCTTGTAGTGAGGTAAAACAAAGCCCCGGCAGTGATGCCGGGGCTTTTTTTTGTATGCCTGGGCGGGCCATTTCGCGGGTGAACCCGCTACTACATGGGGGGCGGGGCTGTTACtt                                                                    >  6:321604/1‑149 (MQ=255)
                                                                             gAGCGGGCAGGCAGCGGTGGTGCCGTCCTGAGCCTTGTAGTGAGGTAAAACAAAGCCCCGGCAGTGATGCCGGGGCTTTTTTTTGTATGCCTGGGCGGGCCATTTCGCGGGTGAACCCGCTCCTACAGGGGGGGCGGGGCTGTTACtt                                                                    >  6:379004/1‑148 (MQ=255)
                                                                               gCGGGCAGGCAGCGGTGGTGCCGTCCTGAGCCTTGTAGTGAGGTAAAACAAAGCCCCGGCAGTGATGCCGGGGCTTTTTTTTGTATGCCTGGGCGGGCCATTTCGCGGGTGAACCCGCTCCTACAGGGGGGGCGGGGCTGTTACTTGaa                                                                 >  7:297444/1‑149 (MQ=255)
                                                                                 gggCAGGCAGCGGTGGTGCCGTCCTGAGCCTTGTAGTGAGGTAAAACAAAGCCCCGGCAGTGATGCCGGGGCTTTTTTTTGTATGCCTGGGCGGGCCATTTCGCGGGTGAACCCGCTCCTACAGGGGGGGCGGGGCTGTTACTTGAAcc                                                               <  1:149570/149‑1 (MQ=255)
                                                                                    caggcagCGGTGGTGCCGTCCTGAGCCTTGTAGTGAGGTAAAACAAAGCCCCGGCAGTGCTGCCGGGGCTGTTTTTTGTATGCCTGGGCGGGCCATTTCGCGGGTGAACCCGGTCGTACAGGGGGGGCGGGGCTGTTACTTGAACCGtg                                                            >  7:342328/1‑148 (MQ=255)
                                                                                             gtggtgCCGTCCTGAGCCTTGTAGTGAGGTAAAACAAAGCCCCGGCAGTGATGCCGGGGCTTTTTTTTGTATGCCTGGGCGGGCCATTTCGCGGGTGAACCCGCTCCTACAGGGGGGGAGGGGCTGTTACTTGAACCTTCTCTCCCCGc                                                   >  8:339273/1‑149 (MQ=255)
                                                                                                 tgCCGTCCTGAGCCTTGTAGTGAGGTAAAACAAAGCCCCGGCAGTGATGCCGGGGCTTTTTTTTGTATGCCTGGGCGGGCCATTTCGCGGGTGAACCCGCTCCTACAGGGGGGGCGGGGCTGTTACTTGAACCGTCTCTCCACGCCttt                                               <  2:32303/149‑1 (MQ=255)
                                                                                                  gCCGTCCTGCGCCTTGTAGTGCGGTCAAACCACCCCCCGGCCGGGCTGCCGGGGCTTTTTTTTGTAGGCCGGGGCGGGCCATTTCGCGGGTGAACCCGCTCCTACCGGGGGGGCGGGGCTGTTACTTGAACCGTCTCTCCACGCCtttt                                              <  7:298612/149‑1 (MQ=255)
                                                                                                   ccGTCCTGAGCCTTGTAGTGAGGTAAAACAAAGCCCCGGCAGTGATGCCGGGGCTTTTTTTTGTATGCCTGGGCGGGCCATTTCGCGGGTGAACCCGCTCCTACAGGGGGGGCGGGGCTGTTACTTGAACCGTCTCTCCACGCCttttt                                             <  8:297444/149‑1 (MQ=255)
                                                                                                       ccTGAGCCTTGTAGTGAGGTAAAACAAAGCCCCGGCAGTGATGCCGGGGCTTTTTTTTGTATGCCTGGGCGGGCCATTTCGCGGGTGAACCCGCTCCTACAGGGGGGGCGGGGCTGTTACTTGAACCGTCTCTCCACGCCTTTTTCGAc                                         >  7:380637/1‑149 (MQ=255)
                                                                                                          gAGCCTTGTAGTGAGGTAAAACAAAGCCCCGGCAGTGATGCCGGGGCTTTTTTTTGTATGCCTGGGCGGGCCATTTCGCGGGTGAACCCGCTCCTACAGGGGGGGCGGGGCTGTTACTTGAACCGTCTCTCCACGCCTTTTTCGACCAg                                      >  3:355391/1‑149 (MQ=255)
                                                                                                           aGCCTTGTAGTGAGGTAAAACAAAGCCCCGGCAGTGATGCCGGGGCTTTTTTTTGTATGCCTGGGCGGGCCATTTCGCGGGTGAACCCGCTCCTACAGGGGGGGCGGGGCTGTTACTTGAACCGTCTCTCCACGCCTTTTTCGACCAg                                      >  2:204897/1‑148 (MQ=255)
                                                                                                           aGCCTTGTAGGGAGGTAAAACAAAGCCCCGGCTGTGATGCCGGGGCTTTTTTTTGTATGCCTGGGCGGGCCATTTCGCGGGTGAACCCGCTCCTACCGGGGGGGCGGGGCTGTTACTTGAACCGTCTCTCCACGCCTTTTTCGACCAg                                      <  5:108007/148‑1 (MQ=255)
                                                                                                               ttGTAGTGAGGTAAAACAAAGCCCCGGCAGTGATGCCGGGGGTTTTTTTTGTATGCCTGGGCGGGCCATTTCGCGGGTGAACCCGCTCCTACAGGGGGGGCGGGGCTGTTACTTGAACCGTGTGTCCCCGCCTTTTTCGACCAGGATCt                                 >  7:186625/1‑149 (MQ=255)
                                                                                                                    gTGAGGTAAAACAAAGCCCCGGCAGTGATGCCGGGGCTTTTTTTTGTATGCCTGGGCGGGCCATTTCGCGGGTGAACCCGCTCCTACAGGGGGGGAGGGGCTGTTACTTGAACCGTCTATCCACGCCTTTTTCGACCAGCATCTTggag                            >  4:253395/1‑147 (MQ=255)
                                                                                                                     tGAGGTAAAACAAAGCCCCGGCAGTGATGCCGGGGCTTTTTTTTGTATGCCTGGGCGGGCCATTTCGCGGGTGAACCCGCTCCTACAGGGGGGGCGGGGCTGTTACTTGAACCGTCTCTCCAGGCCTTTTTGGACCAGGATCTTGGCg                            >  3:313892/1‑148 (MQ=255)
                                                                                                                     tGAGGTAAAACAAAGCCCCGGCAATGATGCCGGGGCTTTTTTTTGTATGCCTGGGCGGGCCATTTCGCGGGTGAACCCGCTCCTACAGGGGGGGCGGGGCTGTTAATTGAACCGTCTCTCCACGCCTTTTTCGACCAGGATCTTGGCg                            >  7:259589/1‑148 (MQ=255)
                                                                                                                               cAACGCCCCGGCAGGGATGCCGGGGCTTTTTTTTGTATGCGTGGGCGGGCCATTGCGCGGGTGACCCCGCTCCTACCGGGGGGGCGGGGCTGTTACTTGAACCGTCTCTCCACGCCTTTTTCGACCAGGATCTTGGCGGAAATCTCTTc                 <  7:339273/149‑1 (MQ=255)
                                                                                                                                aaaGCCCCGGCAGTGATGCCGGGGCTTTTTTTTGTATGCCTGGGCGGGCCATTTCGCGGGTGAACCCGCTCCTACAGGGGGGGCGGGGCTGTTACTTGaa                                                                 >  7:174309/1‑100 (MQ=255)
                                                                                                                                aaaGCCCCGGCAGTGATGCCGGGGCTTTTTTTTGTATGCCTGGGCGGGCCATTTCGCGGGTGAACCCGCTCCTACAGGGGGGGCGGGGCTGTTACTTGaa                                                                 <  8:174309/100‑1 (MQ=255)
                                                                                                                                   gCCCCGGCAGTGATGCCGGGGCTTTTTTTTGTATGCCTGGGCGGGCCATTTCGCGGGTGAACCCGCTCCTCCAGGGGGGGCGGGGCTGTTACTTGAACCTTCTCTCCACGCCTTTTTCGACCAGGATCTTGGCGGAAATCTCTTCCAcc             <  1:301276/149‑1 (MQ=255)
                                                                                                                                      ccGGCAGTGATGCCGGGGCTTTTTTTTGTATGCCTGGGCGGGCCATTTCGCGGGTGAACCCGCTCCTACAGGGGGGGCGGGGCTGTTACTTGAACCGTCTCTCCACGCCTTTTTCGACCAGGATCTTGGCGGAAATCTCTTCCACCGa           >  5:368370/1‑148 (MQ=255)
                                                                                                                                      ccGGCAGCACTGCCGGGGCTTTTTTTTGTATGCCTGGGCGGGCCATTTCGCGGGTGAACCCGCTCCTACAGGGGGGGCGGGGCTGTTAATTGAACCGTCTCTCCACGCCTTTTTCGCCCAGGATCTTGGAGGAAATCTCTTCCACCGa           >  7:252478/1‑148 (MQ=39)
                                                                                                                                             tGATGCCGGGGCTTTTTTTTGTATGCCTGGGCGGGCCATTTCGCGGGTGAACCCGCTCCTACAGGGGGGGCGGGGCTGTTACTTGAACCGTCTCTCCACGCCTTTTTCGACCa                                       >  5:237629/1‑113 (MQ=255)
                                                                                                                                             tGATGCCGGGGCTTTTTTTTGTATGCCTGGGCGGGCCATTTCGCGGGTGAACCCGCTCCTACAGGGGGGGCGGGGCTGTTACTTGAACCGTCTCTCCACGCCTTTTTCGACCa                                       <  6:237629/113‑1 (MQ=255)
                                                                                                                                               aTGCCGGGGCTTTTTTTTGTATGCCTGCGCGGGCCATTCCCCGGCGGAACCCGCTCCTACAGGGGGGGCGGGGCTGTTACTTGAACCGTCTCTCCACGCCTTTTTCGACCAGGATCTTGGCGGAAATCTCTTCCACCGAAAAATGCGt  <  4:313892/148‑1 (MQ=255)
                                                                                                                                               aTGCCGGGGCTGTTTTTTGTATGCCTGGCCGGCCCATTTCGCGGCTGAACCCGCTCCTACCGGGGTGGCGGGGCTGTTACTTGAACCGTCTCTCCACGCCTTTTTCGACCAGGATCTTGGCGGAAATCTCTTCCACCGAAAAATGCGt  <  6:368370/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
GACGACCTGCGCAATGCCACGGGCACTGACTATGCGATGTACGCGGTGGCCAACCGCGAGCTGGTCGACCTGGCCATGAGCGGGCAGGCAGCGGTGGTGCCGTCCTGAGCCTTGTAGTGAGGTAAAACAAAGCCCCGGCAGTGATGCCGGGGCTTTTTTTTGTATGCCTGGGCGGGCCATTTCGCGGGTGAACCCGCTCCTACAGGGGGGGCGGGGCTGTTACTTGAACCGTCTCTCCACGCCTTTTTCGACCAGGATCTTGGCGGAAATCTCTTCCACCGAAAAATGCGT  >  NC_002947/2371941‑2372231

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑