Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,774,279 C→T 7.4% A164V (GCC→GTC)  ydiB → quinate/shikimate 5‑dehydrogenase, NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,774,2790CT 7.4% 165.1 / 2.8 54A164V (GCC→GTC) ydiBquinate/shikimate 5‑dehydrogenase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base C (26/24);  new base T (2/2);  total (28/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.74e-01

TTAAAGAGAGCGGTTTTGATATCAAAGGCAAAACGATGGTGCTGTTAGGGGCCGGTGGTGCCTCAACGGCAATTGGCGCGCAGGGGGCAATTGAAGGTTTAAAAGAAATTAAACTCTTTAACCGTCGGGATGAGTTCTTCGATAAAGCCCTCGCCTTCGCGCAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTT  >  NC_000913/1774132‑1774421
                                                                                                                                                   |                                                                                                                                              
ttAAAGAGAGCGGTTTTGATATCAAAGGCAAAACGATGGTGCTGTTAGGGGCCGGTGGTGCCTCAACGGCAATTGGCGCGCAGGGGGCAATTGAAGGTTTAAAAGAAATTAAACTCTTTAACCGTCGGGATGAGTTCTTCGATAAAGcc                                                                                                                                               >  3:101016/1‑149 (MQ=255)
      agagCGGTTTTGATATCAAAGGCAAAACGATGGTGCTGTTAGGGGCCGGTGGTGCCTCAACGGCAATTGGCGCGCAGGGGGCAATTGAAGGTTTAAAAGAAATTAAACTCTTTAACCGTCGGGATGAGTTCTTCGATAAAGCCCTCGcc                                                                                                                                         <  7:45723/149‑1 (MQ=255)
       gagCGGTTTTGATATCAAAGGCAAAACGATGGTGCTGTTAGGGGCCGGTGGTGCCTCAACGGCAATTGGCGCGCAGGGGGCAATTGAAGGTTTAAAAGAAATTAAACTCTTTAACCGTCGGGATGAGTTCTTCGATAAAGCCCTCGCCt                                                                                                                                        >  6:148106/1‑149 (MQ=255)
         gCGGTTTTGATATCAAAGGCAAAACTATGGTGCTGTTAGGGGCCGGTGGTGCCTCAACGGCAATTGGCGCGCAGGGGGCAATTGAAGGTTTAAAAGAAATTAAACTCTTTAACAGTCGGGATGAGTTCTTCGATAAAGCCCTCGCCTTc                                                                                                                                      <  3:139832/149‑1 (MQ=255)
          cGGTTTTGATATCAAAGGCAAAACGATGGTGCTGTTAGGGGCCGGTGGTGCCTCAACGGCAATTGGCGCGCAGGGGGCAATTGAAGGTTTAAAAGAAATTAAACTCTTTAACCGTCGGGATGAGTTCTTCGATAAAGCCCTCGCCTTcg                                                                                                                                     >  8:76367/1‑149 (MQ=255)
                   tatCAAAGGCAAAACGATGGTGCTGTTAGGGGCCGGTGGTGCCTCAACGGCAATTGGCGCGCAGGGGGCAATTGAAGGTTTAAAAGAAATTAAACTCTTTAACCGTCGGGATGAGTTCTTCGATAAAGCCCTCGCCTTCGCGCAGCggg                                                                                                                            >  8:139452/1‑149 (MQ=255)
                   tatCAAAGGCAAAACGATGGTGCTGTTAGGGGCCGGTGGTGCCTCAACGGCAATTGGCGCGCAGGGGGCAATTGAAGGTTTAAAAGAAATTAAACTCTTTAACCGTCGGGATGAGTTCTTCGATAAAGCCCTCGCCTTCGCGCAGCggg                                                                                                                            >  4:230347/1‑149 (MQ=255)
                    atCAAAGGCAAAACGATGGTGCTGTTAGGGGCCGGTGGTGCCTCAACGGCAATTTGCGCGCAGGGGGCAATTGAAGGTTTAAAAGAAATTAAACTCTTTAACCGTCGGGATGAGTTCTTCGATAAAGCCCTCGCCTTCGCGCAGCGGGt                                                                                                                           >  4:95696/1‑149 (MQ=255)
                     tCAAAGGCAAAACGATGGTGCTGTTAGGGGCCGGTGGTGCCTCAACGGCAATTGGCGCGCAGGGGGCAATTGAAGGTTTAAAAGAAATTAAACTCTTTAACCGTCGGGATGAGTTCTTCGATAAAGCCCTCGCCTTCGCGCAGc                                                                                                                               <  6:136143/144‑1 (MQ=255)
                     tCAAAGGCAAAACGATGGTGCTGTTAGGGGCCGGTGGTGCCTCAACGGCAATTGGCGCGCAGGGGGCAATTGAAGGTTTAAAAGAAATTAAACTCTTTAACCGTCGGGATGAGTTCTTCGATAAAGCCCTCGCCTTCGCGCAGc                                                                                                                               >  5:136143/1‑144 (MQ=255)
                                       tGCTGTTAGGGGCCGGTGGTGCCTCAACGGCAATTGGCGCGCAGGGGGCAATTGAAGGTTTAAAAGAAATTAAACTCTTTAACCGTCGGGATGAGTTCTTCGATAAAGCCCTCGCCTTCGCGCAGCGGGTTAATGAAAACACCGATtgt                                                                                                        <  5:153905/149‑1 (MQ=255)
                                       tGCTGTTAGGGGCCGGTGGTGCCTCAACGGCAATTGGCGCGCAGGGGGCAATTGAAGGTTTAAAAGAAATTAAACTCTTTAACCGTCGGGATGAGTTCTTCGATAAAGCCCTCGCCTTCGCGCAGCGGGTTAATGAAAACACCGATtgt                                                                                                        <  8:215890/149‑1 (MQ=255)
                                              aGGGGCCGGTGGTGCCTCAACGGCAATTGGCGCGCAGGGGGCAATTGAAGGTTTAAAAGAAATTAAACTCTTTAACCGTCGGGATGAGTTCTTCGATAAAGCCCTCGCCTTCGCGCAGCGGGTTAATGAAAACACCGATTGTGTCGTCa                                                                                                 >  7:37973/1‑149 (MQ=255)
                                                       tggtgCCTCAACGGCAATTGGCGCGCAGGGGGCAATTGAAGGTTTAAAAGAAATTAAACTCTTTAACCGTCGGGATGAGTTCTTCGATAAAGCCCTCGCCTTCGCGCAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCAcc                                                                                         >  1:110648/1‑148 (MQ=255)
                                                             cTCAACGGCAATTGGCGCGCAGGGGGCAATTGAAGGTTTAAAAGAAATTAAACTCTTTAACCGTCGGGATGAGTTCTTCGATAAAGCCCTCGCCTTCGCGCAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCg                                                                                  >  3:104846/1‑149 (MQ=255)
                                                                  cGGCAATTGGCGCGCAGGGGGCAATTGAAGGTTTAAAAGAAATTAAACTCTTTAACCGTCGGGATGAGTTCTTCGATAAAGCCCTCGCCTTCGCGCAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGAt                                                                             <  4:4106/149‑1 (MQ=255)
                                                                   ggCAATTGGCGCGCAGGGGGCAATTGAAGGTTTAAAAGAAATTAAACTCTTTAACCGTCGGGATGAGTTCTTCGATAAAGCCCTCGCCTTCGCGCAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGAt                                                                             >  5:88560/1‑148 (MQ=255)
                                                                        ttGGCGCGCAGGGGGCAATTGAAGGTTTAAAAGAAATTAAACTCTTTAACCGTCGGGATGAGTTCTTCGATAAAGCCCTCGCCTTCGCGCAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCaa                                                                       <  8:37973/149‑1 (MQ=255)
                                                                          ggCGCGCAGGGGGCAATTGAAGGTTTAAAAGAAATTAAACTCTTTAACCGTCGGGATGAGTTCTTCGATAAAGCCCTCGCCTTCGCGCAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAg                                                                      <  4:104846/148‑1 (MQ=255)
                                                                             gcgcAGGGGGCAATTGAAGGTTTAAAAGAAATTAAACTCTTTAACCGTCGGGATGAGTTCTTCGATAAAGCCCTCGCCTTCGCGCAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCtt                                                                  <  7:139452/149‑1 (MQ=255)
                                                                               gcAGGGGGCAATTGAAGGTTTAAAAGAAATTAAACTCTTTAACCGTCGGGATGAGTTCTTCGATAAAGCCCTCGCCTTCGCGCAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTg                                                                <  3:95696/149‑1 (MQ=255)
                                                                               gcAGGGGGCAATTGAAGGTTTAAAAGAAATTAAACTCTTTAACCGTCGGGATGAGTTCTTCGATAAAGCCCTCGCCTTCGCGCAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTg                                                                <  5:26505/149‑1 (MQ=255)
                                                                                 aGGGGGCAATTGAAGGTTTAAAAGAAATTAAACTCTTTAACCGTCGGGATGAGTTCTTCGATAAAGCCCTCGCCTTCGCGCAGCGGGTTAATGAAAACCCCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCt                                                              >  1:26002/1‑149 (MQ=255)
                                                                                      gCAATTGAAGGTTTAAAAGAAATTAAACTCTTTAACCGTCGGGATGAGTTCTTCGATAAAGCCCTCGCCTTCGCGCAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCggc                                                                                <  5:205572/126‑3 (MQ=255)
                                                                                      gCAATTGAAGGTTTAAAAGAAATTAAACTCTTTAACCGTCGGGATGAGTTCTTCGATAAAGCCCTCGCCTTCGCGCAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGcc                                                                                >  6:205572/1‑126 (MQ=255)
                                                                                        aaTTGAAGGTTTAAAAGAAATTAAACTCTTTAACCGTCGGGATGAGTTCTTCGATAAAGTCCTCGCCTTCGCGCAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCAGATCAGCAAGCCTTTTCTGAAGccc                                                       >  6:182700/1‑149 (MQ=255)
                                                                                        aaTTGAAGGTTTAAAAGAAATTAAACTCTTTAACCGTCGGGATGAGTTCTTCGATAAAGCCCTCGCCTTCGCGCAGCGGGGTAATGAAACCACCGATTGTGTCGTCACGGTCACCGCTCTCGCCGATCAGCAAGCCTTTGCTGAAGccc                                                       >  4:132988/1‑149 (MQ=255)
                                                                                         aTTGAAGGTTTAAAAGAAATTAAACTCTTTAACCGTCGGGATGAGTTCTTCGATAAAGCCCTCGCCTTCGCGCAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCt                                                      <  5:148106/149‑1 (MQ=255)
                                                                                         aTTGAAGGTTTAAAAGAAATTAAACTCTTTAACCGTCGGGATGAGTTCTTCGATAAAGCCCTCGCCTTCGCGCAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCt                                                      <  3:102988/149‑1 (MQ=255)
                                                                                          ttGAAGGTTTAAAAGAAATTAAACTCTTTAACCGTCGGGATGAGTTCTTCGATAAAGCCCTCGCCTTCGCGCAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTg                                                     >  4:190023/1‑149 (MQ=255)
                                                                                           tgtacggttAAAACCAAATAAACTCTTTAACCGTCGGGATGAGTTCTTCGATAAAGCCCTCGCCTTCGCGCAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTgg                                                    <  4:210885/142‑1 (MQ=255)
                                                                                             aaGGTTTAAAAGAAATTAAACTCTTTAACCGTCGGGATGAGTTCTTCGATAAAGCCCTCGCCTTCGCGCAGCGGGTTAATGAAAACACCGAtt                                                                                                          <  7:29449/93‑1 (MQ=255)
                                                                                             aaGGTTTAAAAGAAATTAAACTCTTTAACCGTCGGGATGAGTTCTTCGATAAAGCCCTCGCCTTCGCGCAGCGGGTTAATGAAAACACCGAtt                                                                                                          >  8:29449/1‑93 (MQ=255)
                                                                                             aaGGTTTAAAAGAAATTAAACTCTTTAACCGTCGGGATGAGTTCTTCGATAAAGCCCTCGCCTTCGCGCAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCt                                                  >  4:222466/1‑149 (MQ=255)
                                                                                                  ttAAAAGAAATTAAACTCTTTAACCGTCGGGATGAGTTCTTCGATAAAGCCCTCGCCTTCGCGCAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTccgc                                             >  4:17037/1‑149 (MQ=255)
                                                                                                      aaGAAATTAAACTCTTTAACCGTCGGGCTGAGTTCTTCGATAAAGCCCTCGCCTTCGCGCAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCAGGGCTTcagcaaa                                          <  2:66944/148‑4 (MQ=255)
                                                                                                      aaGAAATTAAACTCTTTAACCGTCGGGCTGAGTTCTTCGATAAAGCCCTCGCCTTCGCGCAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCAGGGCTTcagca                                            >  1:66944/1‑145 (MQ=255)
                                                                                                       aGAAATTAAACTCTTTAACCGTCGGGATGAGTTCTTCGATAAAGCCCTCGCCTTCGCGCAGCGGGTTAATGAAAACACCGAtt                                                                                                          <  8:184904/83‑1 (MQ=255)
                                                                                                       aGAAATTAAACTCTTTAACCGTCGGGATGAGTTCTTCGATAAAGCCCTCGCCTTCGCGCAGCGGGTTAATGAAAACACCGAtt                                                                                                          >  7:184904/1‑83 (MQ=255)
                                                                                                            ttAAACTCTTTAACCGTCGGGATGAGTTCTTCGATAAAGCCCTCGCCTTCGCGCAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACAtttt                                    >  5:27792/1‑148 (MQ=255)
                                                                                                                aCTCTTTAACCGTCGGGATGAGTTCTTCGATAAAGCCCTCGCCTTCGCGCAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGaa                                                           <  7:116913/121‑1 (MQ=255)
                                                                                                                aCTCTTTAACCGTCGGGATGAGTTCTTCGATAAAGCCCTCGCCTTCGCGCAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGaa                                                           >  8:116913/1‑121 (MQ=255)
                                                                                                                 ctctTTAACCGTCGGGATGAGTTCTTCGATAAAGCCCTCGCCTTCGCGCAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCaa                              <  4:101016/149‑1 (MQ=255)
                                                                                                                  tctTTAACCGTCGGGATGAGTTCTTCGATAAAGCCCTCGCCTTCGCGCAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAAt                             >  8:171597/1‑149 (MQ=255)
                                                                                                                    tttAACCGTCGGGATGAGTTCTTCGATAAAGCCCTCGCCTTCGCGCAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAAt                             <  6:88560/147‑1 (MQ=255)
                                                                                                                        aCCGTCGGGATGAGTTCTTCGATAAAGCCCTCGCCTTCGCGCAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGcac                        >  3:166913/1‑148 (MQ=255)
                                                                                                                          cGTCGGGATGAGTTCTTCGATAAAGCCCTCGCCTTCGCGCAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACaa                      >  2:166308/1‑148 (MQ=255)
                                                                                                                               ggATGAGTTCTTCGATAAAGTCCTCGCCTTCGCGCAGCGGGTTAATGAAAACACCGATTGTGTCGTCACTGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTgg                <  5:182700/149‑1 (MQ=255)
                                                                                                                                      ttcttcGATAAAGTCCTCGCCTTCGCGCAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGa          >  1:77068/1‑148 (MQ=255)
                                                                                                                                        cttcGATAAAGCCCTCGCCTTCGCGCAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGaaa        <  7:171597/148‑1 (MQ=255)
                                                                                                                                           cGATAAAGTCCTCGCCTTCGCGCAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAAccc     <  2:77068/148‑1 (MQ=255)
                                                                                                                                            gATAAAGCCCTCGCCTTCGCGCAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAAt                             >  7:205268/1‑123 (MQ=255)
                                                                                                                                            gATAAAGCCCTCGCCTTCGCGCAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAAt                             <  8:205268/123‑1 (MQ=255)
                                                                                                                                             aTAAAGCCCTCGCCTTCGCGCAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCtt  <  3:190023/149‑1 (MQ=255)
                                                                                                                                               aaaGCCCTCGCCCTCGCGCAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCtt  >  1:116341/1‑147 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
TTAAAGAGAGCGGTTTTGATATCAAAGGCAAAACGATGGTGCTGTTAGGGGCCGGTGGTGCCTCAACGGCAATTGGCGCGCAGGGGGCAATTGAAGGTTTAAAAGAAATTAAACTCTTTAACCGTCGGGATGAGTTCTTCGATAAAGCCCTCGCCTTCGCGCAGCGGGTTAATGAAAACACCGATTGTGTCGTCACGGTCACCGATCTCGCCGATCAGCAAGCCTTTGCTGAAGCCCTGGCTTCCGCCGACATTTTAACCAATGGCACAAAAGTGGGTATGAAACCCCTT  >  NC_000913/1774132‑1774421

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