Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,542,321 C→T 15.5% P134P (CCG→CCA)  galE ← UDP‑glucose 4‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,542,3210CT15.5% 102.9 / 10.4 45P134P (CCG→CCA) galEUDP‑glucose 4‑epimerase
Reads supporting (aligned to +/- strand):  ref base C (17/21);  new base T (3/4);  total (20/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.01e-01

TGGGTTTCCGGCTCGTGACGCTCGCCTAGCTGGC‑CTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTGCAGCAGGGAAAGGGTGGCTGCGAC  >  NC_002947/3542175‑3542461
                                                                                                                                                   |                                                                                                                                            
tGGGTTTCCGGCTCGTGACGCTCGCCTAGCTGGCGCTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACg                                                                                                                                             >  1:197151/1‑149 (MQ=255)
 ggggTTCCGGCTCGTGACGCTCGCCTAGCTGGC‑CTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACgg                                                                                                                                            <  2:144699/148‑1 (MQ=255)
  ggTTTCCGGCTCGTGACGCTCGCCTAGCTGGC‑CTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACggg                                                                                                                                           <  6:99115/148‑1 (MQ=255)
  ggTTTCCGGCTCGTGACGCTCGCCTAGCTGGC‑CTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGt                                                                                                                                          <  6:22809/149‑1 (MQ=255)
        cGGCTCGTGACGCTCGCCTAGCTGGC‑CTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCGTATGGGTTGagc                                                                                                                                     <  1:285213/148‑3 (MQ=255)
         ggCTCGTGACGCTCGCCTAGCTGGC‑CTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGATCAACCCGTATGGGTTGagca                                                                                                                                    >  7:30567/1‑145 (MQ=255)
          gCTCGTGACGCTCGCCTAGCTGGC‑CTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGc                                                                                                                                  <  2:197151/149‑1 (MQ=255)
           cTCGTGACGCTCGCCTAGCTGGC‑CTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGATCAACCCGTATGGGTTGagcaaat                                                                                                                                 <  8:30567/149‑7 (MQ=255)
                gACGCTCGCCTAGCTGGC‑CTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTg                                                                                                                            >  4:182850/1‑149 (MQ=255)
                       gCCTAGCTGGC‑CTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGgtgtgt                                                                                                                      <  4:261528/148‑1 (MQ=255)
                        ccTAGCTGGC‑CTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGc                                                                                                                    <  8:48182/149‑1 (MQ=255)
                              tGGC‑CTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCatc                                                                                                              >  6:161342/1‑149 (MQ=255)
                                   cTTCCGGGTCAGCGCCCGCTGCATTGAAATAACGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGgtgt                                                                                                                        >  5:149116/1‑135 (MQ=255)
                                   cTTCCGGGTCAGCGCCCGCTGCATTGAAATAACGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGgtgt                                                                                                                        <  6:149116/135‑1 (MQ=255)
                                    ttCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGAt                                                                                                                                      <  5:155538/120‑1 (MQ=255)
                                    ttCCGGGTCAGCGCCCGCAGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGAt                                                                                                                                      >  6:155538/1‑120 (MQ=255)
                                       ctggTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACCTCAAACCATTGGCGCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCTATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGTTCGGc                                                                                                       <  1:237936/146‑1 (MQ=255)
                                                         aTTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCg                                                                                             >  6:183181/1‑140 (MQ=255)
                                                         aTTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCg                                                                                             <  5:183181/140‑1 (MQ=255)
                                                         aTTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGt                                                                                    <  5:122900/149‑1 (MQ=255)
                                                             aaaTAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATAc                                                                                >  1:150726/1‑149 (MQ=255)
                                                                tAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATAcggc                                                                             <  2:150726/149‑1 (MQ=255)
                                                                 aGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATAcggcg                                                                            <  5:15753/149‑1 (MQ=255)
                                                                     cAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCAc                                                                        >  4:50462/1‑149 (MQ=255)
                                                                        gCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTc                                                                      >  3:214609/1‑148 (MQ=255)
                                                                                    cAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACAcc                                                          <  3:182850/148‑1 (MQ=255)
                                                                                    cAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACAcc                                                          <  5:161342/148‑1 (MQ=255)
                                                                                        ccATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTg                                                     >  1:211008/1‑149 (MQ=255)
                                                                                         caaaGGGCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGc                                                    <  4:222943/146‑1 (MQ=255)
                                                                                                    tCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCgg                                          >  6:244390/1‑148 (MQ=255)
                                                                                                      gAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGcccc                                                                                                                               >  6:287283/1‑61 (MQ=255)
                                                                                                      gAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGcccc                                                                                                                               <  5:287283/61‑1 (MQ=255)
                                                                                                               aaGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATgg                              <  5:244390/149‑1 (MQ=255)
                                                                                                                       gCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTgcag                       >  8:173523/1‑148 (MQ=255)
                                                                                                                             ccatccatTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCt                           >  6:151617/1‑138 (MQ=255)
                                                                                                                             ccatccatTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCt                           <  5:151617/138‑1 (MQ=255)
                                                                                                                                  catTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTGCAGCAGGGAAAGGGt           <  4:214609/149‑1 (MQ=255)
                                                                                                                                   atTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTGCAGCAGGGAAAGGGTg          >  2:112481/1‑149 (MQ=255)
                                                                                                                                     ccgatCAACCCGTATGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCAtt                                      <  3:167941/115‑1 (MQ=255)
                                                                                                                                     ccgatCAACCCGTATGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCAtt                                      >  4:167941/5‑119 (MQ=255)
                                                                                                                                        atCAACCCGTATGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTGCAGCAGGGAAAGGGTgg         <  7:14537/144‑1 (MQ=255)
                                                                                                                                        atCAACCCGTATGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTGCAGCAGGGAAAGGGTgg         >  8:14537/2‑145 (MQ=255)
                                                                                                                                         tCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTGCAGCAGGGAAAGGGTGGCTGc     >  4:114373/1‑148 (MQ=255)
                                                                                                                                           aaCCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTGCAGCAGGGAAAGGGTGGCTGCGa   <  1:112481/148‑1 (MQ=255)
                                                                                                                                           aaCCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTGCAGCAGGGAAAGGGTGGCTGCGAc  >  2:289801/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                            
TGGGTTTCCGGCTCGTGACGCTCGCCTAGCTGGC‑CTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTGCAGCAGGGAAAGGGTGGCTGCGAC  >  NC_002947/3542175‑3542461

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