Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,185,405 T→G 7.3% D352A (GAT→GCT) ‡ aceK ← isocitrate dehydrogenase kinase/phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,185,4050TG 7.3% 169.7 / 4.9 55D352A (GAT→GCT) ‡aceKisocitrate dehydrogenase kinase/phosphatase
Reads supporting (aligned to +/- strand):  ref base T (27/24);  new base G (2/2);  total (29/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 = 9.55e-01

GCGCCACCTCCAGCAACTCGGCCAGGCAATCCGGCTCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCATGACAAAGCGGTCATCACTGCTGGCCAGGTGGTTGATCAGTGC  >  NC_002947/5185259‑5185553
                                                                                                                                                  |                                                                                                                                                    
gcgcCACCTCCATCAACTCGGCCAGGCAATCCGGCTCTAAATTGCTGCGCGTAAAGCTAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTGTCGATCACCGTGGCGCGATCCa                                                                                                                                                   <  7:52449/149‑1 (MQ=255)
       cTCCAGCAACTCGGCCAGGCAATCCGGCTCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGATCGCGCCACGGtga                                                                                                                                             <  3:233710/149‑3 (MQ=255)
        tCCAGCAACTCGGCCAGGCAATCCGGCTCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGtttt                                                                                                                                            >  1:350626/1‑149 (MQ=255)
           aGCAACTCGGCCAGGCAATCCGGCTCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGa                                                                                                                                         <  8:344960/149‑1 (MQ=255)
            gCAACTCGGCCAGGCAATCCGGCTCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGAt                                                                                                                                        <  7:65034/149‑1 (MQ=255)
              aaCTCGGCCAGGCAATCCGGCTCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATgg                                                                                                                                      >  4:91177/1‑149 (MQ=255)
                cTCGGCCAGGCAATCCGGCTCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGc                                                                                                                                     >  1:133729/1‑148 (MQ=255)
                cTCGGCCAGGCAATCCGGCTCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCg                                                                                                                                    <  3:91177/149‑1 (MQ=255)
                cTCGGCCAGGCAATCCGGCTCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGATCGCGCCACGGtga                                                                                                                                             >  8:313226/1‑138 (MQ=255)
                cTCGGCCAGGCAATCCGGCTCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGATCGCGCCACGGtga                                                                                                                                             <  7:313226/140‑3 (MQ=255)
                     ccAGGCAATCCGGCTCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAg                                                                                                                               <  2:59200/149‑1 (MQ=255)
                             tCCGGCTCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCt                                                                                                                        <  3:136139/148‑1 (MQ=255)
                             tCCGGCTCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCt                                                                                                                        >  4:136139/1‑148 (MQ=255)
                               cGGCTCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTtg                                                                                                                      <  7:23187/148‑1 (MQ=255)
                                   tCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTgatgat                                                                                                                 >  5:322031/1‑149 (MQ=255)
                                     gAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCt                                                                                                               >  1:237182/1‑149 (MQ=255)
                                         ttGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGAt                                                                                                                                        >  1:235508/1‑120 (MQ=255)
                                         ttGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGAt                                                                                                                                        <  2:235508/120‑1 (MQ=255)
                                              gctcgGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTAGACACTCTTCCACAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTAGATGATCTTGAACACCg                                                                                                      >  5:259148/4‑149 (MQ=255)
                                                      aGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAg                                                                                              <  2:350626/149‑1 (MQ=255)
                                                          aaaaTCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCgg                                                                                          <  2:237182/149‑1 (MQ=255)
                                                                ggCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCGACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAg                                                                                     <  2:133729/148‑1 (MQ=255)
                                                                   aaaCTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACCTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTg                                                                                  >  4:309051/1‑148 (MQ=255)
                                                                     aCTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCg                                                                                           <  6:152202/137‑1 (MQ=255)
                                                                     aCTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCg                                                                                           >  5:152202/1‑137 (MQ=255)
                                                                       tCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACAGTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGt                                                                                   <  3:144295/143‑1 (MQ=255)
                                                                            ggTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCa                                                                        <  5:353562/149‑1 (MQ=255)
                                                                                tCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGAc                                                                    >  3:208564/1‑149 (MQ=255)
                                                                                  ggCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCa                                                                  >  6:248104/1‑149 (MQ=255)
                                                                                          gCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCg                                                           >  6:238721/1‑148 (MQ=255)
                                                                                          gCCCCACCAGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCg                                                           >  7:173568/1‑148 (MQ=255)
                                                                                            cccACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGcccc                                                            >  5:36539/1‑145 (MQ=255)
                                                                                            cccACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGcccc                                                            <  6:36539/145‑1 (MQ=255)
                                                                                             ccACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAg                                                                                              >  8:313221/1‑110 (MQ=255)
                                                                                             ccACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAg                                                                                              <  7:313221/110‑1 (MQ=255)
                                                                                             ccACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGccc                                                             >  5:129870/1‑143 (MQ=255)
                                                                                             ccACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGccc                                                             <  6:129870/143‑1 (MQ=255)
                                                                                                   ggTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGt                                                 <  5:248104/149‑1 (MQ=255)
                                                                                                       gACACTCTTCACCAGCCGATACATGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCACGCACCCCGGGTGcc                                              >  7:236929/1‑148 (MQ=255)
                                                                                                         cacTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCAt                                            <  4:208564/148‑1 (MQ=255)
                                                                                                                 aCCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCATGACAAAGc                                    <  8:173568/148‑1 (MQ=255)
                                                                                                                   cAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACCCGCTCATGACCATGCCCCGCACCACTGGTGCCATGACAAAGCGGt                                 >  8:328395/1‑149 (MQ=255)
                                                                                                                   cAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCATGACAAAGCGGt                                 >  2:24654/1‑149 (MQ=255)
                                                                                                                    aGCCGATACCTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCTGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCATGACAAAGCGGtc                                <  3:309051/149‑1 (MQ=255)
                                                                                                                     gCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCAt                                            >  6:373914/1‑136 (MQ=255)
                                                                                                                     gCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCAt                                            <  5:373914/136‑1 (MQ=255)
                                                                                                                          tACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCATGACAAAGCGGTCATCAct                          >  1:231783/1‑149 (MQ=255)
                                                                                                                             ttGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGCACCCGCTCATGACCATGCCCCGCAACCCTGGTGCCCTGACAACGCGGGCATCCctgct                       >  8:201040/1‑149 (MQ=255)
                                                                                                                             ttGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCATGACAAAGCGGTCATCActgct                       >  3:369011/1‑149 (MQ=255)
                                                                                                                               gTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGc                                                                           <  7:125397/95‑1 (MQ=255)
                                                                                                                               gTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGc                                                                           >  8:125397/1‑95 (MQ=255)
                                                                                                                                     aaCCGTGGATCGCGCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCATGACAAAGCGGTCATCACTGCTGGCCAggt               >  6:171017/2‑149 (MQ=255)
                                                                                                                                       ccGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCATGACAAAGCGGTCATCACTGCTGGCCAggtgg             >  6:193477/1‑149 (MQ=255)
                                                                                                                                              gcgATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCATGACAAAGCg                                   >  4:206934/1‑120 (MQ=255)
                                                                                                                                              gcgATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCATGACAAAGCg                                   <  3:206934/120‑1 (MQ=255)
                                                                                                                                                  tCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCATGACAAAGCGGTCATCACTGCTGGCCAGGTGGTTGATCAGTgc  <  8:236929/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GCGCCACCTCCAGCAACTCGGCCAGGCAATCCGGCTCGAACTTGCTGCGCGGAAAGCGAAAATCGGCAAACTCCTGGGTGTCGGCCATACGCCCCACCCGGTCGACACTCTTCACCAGCCGATACTTGTCGATCACCGTGGCGCGATCCACGGTTTTCGATGGCGAGAAGCGGTCCTTGATGATCTTGAACACCGTATTGAAGCCGGGCAAGGTGAACACGCTCATGACCATGCCCCGCACCCCTGGTGCCATGACAAAGCGGTCATCACTGCTGGCCAGGTGGTTGATCAGTGC  >  NC_002947/5185259‑5185553

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: