Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,049,180 C→G 7.1% F577L (TTC→TTG)  quiA → quinate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,049,1800CG 7.1% 172.4 / 4.2 56F577L (TTC→TTG) quiAquinate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (26/26);  new base G (2/2);  total (28/28)
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.39e-01

GTCAAGCCAAGCAACATCCCGAACGAGCCGTACTCCCCAACCCAGCCCAAGTCGGTGGGCATGCCACAGATCGGCGCGCAGACCCTGACCGAATCGGACATGTGGGGCGCCACCCCGTATGACCAGCTGCTATGCCGCATCGACTTCAAGAAGATGCGCTACGACGGCCTGTACACCGCGCCGGGTACCGACCTGTCGCTGAGCTTCCCGGGTTCGCTGGGTGGCATGAACTGGGGCAGCATTTCTACCGACCCGGTACATGGCTTCATCTTCGTCAACGACATGCGTCTGGG  >  NC_002947/4049034‑4049326
                                                                                                                                                  |                                                                                                                                                  
gTCAAGCCAAGCAACATCCCGAACGAGCCGTACTCCCCAACCCAGCCCAAGTCGGTGGGCATGCCACAGATCGGCGCGCAGACCCTGACCGAATCGGACATGTGGGGCGCCACCCCGTATGACCAGCTGCTATGCCGCATCGACTTCaa                                                                                                                                                  <  3:94789/149‑1 (MQ=255)
           cAACATCCCGAACGAGCCGTACTCCCCAACCCAGCCCAAGTCGGTGGGCATGCCACAGATCGGCGCGCAGACCCTGACCGAATCGGACATGTGGGGCGCCACCCCGTATGACCAGCTGCTATGCCGCATCTTCTTGAAGTCGATgcggc                                                                                                                                       <  3:60494/149‑3 (MQ=255)
                  ccGAACGAGCCGTACTCCCCAACCCAGCCCAAGTCGGTGGGCATGCCACAGATCGGCGCGCAGACCCTGACCGAATCGGACATGTGGGGCGCCACCCCGTATGACCAGCTGCTATGCCGCATCGACTTCAAGAAGATGCGCTACGACgg                                                                                                                                >  6:234824/1‑149 (MQ=255)
                     aaCGAGCCGTACTCCCCAACCCAGCCCAAGTCGGTGGGCATGCCACAGATCGGCGCGCAGACCCTGGCCGAATCGGACATGTGGGGCGCCACCCCGTATGACCAGCTGCTATGCCGCATCGACTTCAAGAAGATGCGCTACGACGGCCt                                                                                                                             <  1:84882/149‑1 (MQ=255)
                     aaCGAGCCGTACTCCCCAACCCAGCCCAAGTCGGTGGGCATGCCACAGATCGGCGCGCAGACCCTGACCGAATCGGACATGTGGGGCGCCACCCCGTATGACCAGCTGCTATGCCGCATCGACTTCAAGAAGATGCGCTACGACGGCCt                                                                                                                             >  3:292602/1‑149 (MQ=255)
                       cGAGCCGTACTCCCCAACCCAGCCCAAGTCGGTGGGCATGCCACAGATCGGCGCGCAGACCCTGACCGAATCGGACATGTGGGGCGCCACCCCGTATGACCAGCTGCTATGCCGCATCGACTTCAAGAAGATGCGCTACGACGGCCTGt                                                                                                                           <  6:68874/149‑1 (MQ=255)
                          gCCGTACTCCCCAACCCAGCCCAAGTCGGTGGGCATGCCACAGATCGGCGCGCAGACCCTGACCGAATCGGACATGTGGGGCGCCACCCCGTATGACCAGCTGCTATGCCGCATCGACTTCaag                                                                                                                                                 <  5:302526/124‑1 (MQ=255)
                          gCCGTACTCCCCAACCCAGCCCAAGTCGGTGGGCATGCCACAGATCGGCGCGCAGACCCTGACCGAATCGGACATGTGGGGCGCCACCCCGTATGACCAGCTGCTATGCCGCATCGACTTCaag                                                                                                                                                 >  6:302526/1‑124 (MQ=255)
                           ccGTACTCCCCAACCCAGCCCAAGTCGGTGGGCATGCCACAGATCGGCGCGCAGACCCTGCCCGAATCGGACATGTGGGGCGCCACCCCGTATGACCAGCTGCTATTCAGCATAGCCTTCAAGAAGATGCCCTACGACGGACTGAacac                                                                                                                       >  7:248580/1‑149 (MQ=255)
                             gTACTCCCCAACCCAGCCCAAGTCGGTGGGCATGCCAAAGATCGGCGCGCAGACCCTGACCGAATCGGACATGTGGGGCGCCACCCCGTATGACCAGCTGCTATGCCGCATCGACTTCAAGAAGATGCGCTACGACGGCCTGTACACcg                                                                                                                     >  3:223631/1‑149 (MQ=255)
                                    ccaaaacAGCACTAGTCGGTGGGAATGCCACAGTTAGGCGCGCATACCCTGACCGAATAGTTCATGTGGGGCGACACACGGTATGACCAGCTGCTATGCCGCATCGACTTCAATAAGATGCGCTACGACGGCCTGTACACCGCGCAggg                                                                                                              <  8:248580/143‑1 (MQ=255)
                                    ccAACCCAGCCCAAGTCGGTGGGCATGCCACAGATCGGCGCGCAGACCCTGACCGAATCGGACATGTGGGGCGCCACCCCGTATGACCAGCTGCTATGCCGCATCTTCTTGAAGTCGATgcgg                                                                                                                                        <  8:96577/123‑2 (MQ=255)
                                    ccAACCCAGCCCAAGTCGGTGGGCATGCCACAGATCGGCGCGCAGACCCTGACCGAATCGGACATGTGGGGCGCCACCCCGTATGACCAGCTGCTATGCCGCATCTTCTTGAAGTCGATgcgg                                                                                                                                        >  7:96577/1‑122 (MQ=255)
                                      aaCCCAGCCCAAGTCGGTGGGCATGCCACAGATCGGCGCGCAGACCCTGACCTAATCGGACATGTGGGGCGCCACCCCGTATGACCAGCTGCTATGCCTCATCGACTTCAAGAAGATTCGCTACGACGGCCTGTACACCGCGCCGGGTa                                                                                                            <  2:25103/149‑1 (MQ=255)
                                         ccAGCCCAAGTCGGTGGGCATGCCACAGATCGGCGCGCAGACCCTGACCGAATCGGACATGTGGGGCGCCACCCCGTATGACCAGCTGCTATGCCGCATCGACTTCAAGAAGATGCGCTACGACGGCCTGTACACCGCGCCGGGTACCg                                                                                                         <  7:1017/149‑1 (MQ=255)
                                                     ggTGGGCATGCCACAGATCGGCGCGCAGACCCTGACCGAATCGGACATGTGGGGCGCCACCCCGTATGACCAGCTGCTATGCCGCATCGACTTCAAGAAGATGCGCTACGACGGCCTGTACACCGCGCCGGGTACCGAc                                                                                                       >  7:50849/1‑139 (MQ=255)
                                                     ggTGGGCATGCCACAGATCGGCGCGCAGACCCTGACCGAATCGGACATGTGGGGCGCCACCCCGTATGACCAGCTGCTATGCCGCATCGACTTCAAGAAGATGCGCTACGACGGCCTGTACACCGCGCCGGGTACCGAc                                                                                                       <  8:50849/139‑1 (MQ=255)
                                                               ccACAGATCGGCGCGCAGACCCTGACCGAATCGGACATGTGGGGCGCCACCCCGTATGACCAGCTGCTATGCCGCATCGACTTCAAGAAGATGCGCTACGACGGCCTGTACACCGCGCCGGGTAcc                                                                                                          <  6:185812/126‑1 (MQ=255)
                                                               ccACAGATCGGCGCGCAGACCCTGACCGAATCGGACATGTGGGGCGCCACCCCGTATGACCAGCTGCTATGCCGCATCGACTTCAAGAAGATGCGCTACGACGGCCTGTACACCGCGCCGGGTAcc                                                                                                          >  5:185812/1‑126 (MQ=255)
                                                                       cGGCGCGCAGACCCTGACCGAATCGGACATGTGGGGCGCCACCCCGTATGACCAGCTGCTATGCCGCATCGACTTCAAGAAGATGCGCTACGACGGCCTGTACACCGCGCCGGGTACCGACCTGTCGCTGAGCTTCCCGGGtt                                                                                 >  6:198628/1‑143 (MQ=255)
                                                                       cGGCGCGCAGACCCTGACCGAATCGGACATGTGGGGCGCCACCCCGTATGACCAGCTGCTATGCCGCATCGACTTCAAGAAGATGCGCTACGACGGCCTGTACACCGCGCCGGGTACCGACCTGTCGCTGAGCTTCCCGGGtt                                                                                 <  5:198628/143‑1 (MQ=255)
                                                                        ggCGCGCAGACCCTGACCGAATCGGACATGTGGGGCGCCACCCCGTATGACCAGCTGCTATGCCGCATCGACTTCAAGAAGATGCGCTACGACGGCCTGTACACCGCGCCGGGTACCGACCTGTCGCTGAGCTTCCCGGGTTCGCTgg                                                                           <  4:292602/148‑1 (MQ=255)
                                                                         gcgcgcAGACCCTGACCGAATCGGACATGTGGGGCGCCACCCCGTATGACCAGCTGCTATGCCGCATCGACTTCAAGAAGATGCGCTACGACGGCCTGTACACCGCGCCGGGTACCGACCTGTCGCTGAGCTTCCCGGGTTCGCTGGGt                                                                         >  8:24255/1‑149 (MQ=255)
                                                                          cgcgcAGACCCTGACCGAATCGGACATGTGGGGCGCCACCCCGTATGACCAGCTGCTATGCCGCATCGACTTCAAGAAGATGCGCTACGACGGCCTGTACACCGCGCCGGGTACCGACCTGTCGCTGAGCTTCCCGGGTTCGATGGGTg                                                                        >  7:104203/1‑149 (MQ=255)
                                                                            cgcAGACCCTGACCGAATCGGACATGTGGGGCGCCACCCCGTATGACCAGCTGCTATGCCGCATCGACTTCAAGAAGATGCGCTACGACGGCCTGTACACCGCGCCGGGTACCGACCTGTCGCTGAGCTTCCCGGGTTCGCTGGGTGGc                                                                      <  7:24255/149‑1 (MQ=255)
                                                                                   ccTGACCGAATCGGACATGTGGGGCGCCACCCCGTATGACCAGCTGCTATGCCGCATCGACTTCAAGAAGATGCGCTACGACGGCCTGTACACCGCGCCGGGTACCGACCTGTCGCTGAGCTTCCCGGGTTCGCTGGGTGGCATGAACt                                                               <  5:234824/149‑1 (MQ=255)
                                                                                   ccTGACCGAATCGGACATGTGGGGCGCCACCCCGTATGACCAGCTGCTATGCCGCATCGACTTCAAGAAGATGCGCTACGACGGCCTGTACACCGCGCCGGGTACCGACCTGTCGCTGAGCTTCCCGGGTTCGCTGGGTGGCATGAACt                                                               >  4:296224/1‑149 (MQ=255)
                                                                                   ccTGACCGAATCGGACATGTGGGGCGCCACCCCGTATGACCAGCTGCTATGCCGCATCGACTTCAAGAAGATGCGCTACGACGGCCTGTACACCGCGCCGGGTACCGACCGGTCGCTGAGCTTCGCGGGTTCGCTGGGTGGCATGAACt                                                               >  8:290207/1‑149 (MQ=255)
                                                                                       aCCGAATCGGACATGTGGGGCGCCACCCCGTATGACCAGCTGCTATGCCGCATCGACTTCAAGAAGATGCGCTACGACGGCCTGTACACCGCGCCGGGTACCGACCTGTCGCTGAGCTTCCCGGGTTCGCTGGGTGGCATGAACTgggg                                                           <  3:280194/149‑1 (MQ=255)
                                                                                       aCCGAATCGGACATGTGGGGCGCCACCCCGTATGACCAGCTGCTATGCCGCATCGACTTCAAGAAGATGCGCTACGACGGCCTGTACACCGCGCCGGGTACCAACCTGTCGCTGAGCTTCCCGGGTTCGCTGGGTGGCATGAACTgggg                                                           <  4:116140/149‑1 (MQ=255)
                                                                                         cGAATCGGACATGTGGGGCGCCACCCCGTATGACCAGCTGCTATGCCGCATCGACTTCAAGAAGATGCGCTACGACGGCCTGTACACCGCGCCGGGTACCg                                                                                                         <  5:12281/101‑1 (MQ=255)
                                                                                         cGAATCGGACATGTGGGGCGCCACCCCGTATGACCAGCTGCTATGCCGCATCGACTTCAAGAAGATGCGCTACGACGGCCTGTACACCGCGCCGGGTACCg                                                                                                         >  6:12281/1‑101 (MQ=255)
                                                                                         cGAATCGGACATGTGGGGCGCCACCCCGTATGACCAGCTGCTATGCCGCATCGACTTCAAGAAGATGCGCTACGACGGCCTGTACACCGCGCCGGGTACCGACCTGTCGCTGAGCTTCCCGGGTTCGCTGGGTGGCATGAACTgggg                                                           >  5:45917/1‑147 (MQ=255)
                                                                                             tCGGACATGTGGGGCGCCACCCCGTATGACCAGCTGCTATGCCGCATCGACTTCAAGAAGATGCGCTACGACGGCCTGTACACCGCGCCGGGTACCGACCTGTCGCTGAGCTTCCCGGGTTCGCTGGGTGGCATGAACTGGGGCAGCAt                                                     <  4:140405/149‑1 (MQ=255)
                                                                                                      tgGGGCGCCACCCCGTATGACCAGCTGCTATGCCGCATCGACTTCAAGAAGATGCGCTACGACGGCCTGTACACCGCGCCGGGTACCGACCTGTCGCTGAGCTTCCCGGGTTCGCTGGGTGGCATGAACTGGGGCAGCATTTCTACCGa                                            >  4:11726/1‑149 (MQ=255)
                                                                                                       ggggCGCCACCCCGTATGACCAGCTGCTATGCCGCATCGACTTCAAGAAGATGCGCTACGACGGCCTGTACACCGCGCCGGGTACCGACCTGTCGCTGAGCTTCCCGGGTTCGCTGGGTGGCATGAACTGGGGCAGCATTTCTACCGAc                                           >  3:41748/1‑149 (MQ=255)
                                                                                                        gggCGCCACCCCGTATGACCAGCTGCTATGCCGCATCGACTTCAAGAAGATGCGCTACGACGGCCTGTACACCGCGCCGGGTACCGACCTGTCGCTGAGCTTCCCGGGTTCGCTGGGTGGCATGAACTGGGGCAGCATTTCTACCGAcc                                          >  2:244998/1‑149 (MQ=255)
                                                                                                        gggCGCCACCCCGTATGACCAGCTGCTATGCCGCATCGACTTCAAGAAGATGCGCTACGACGGCCTGTACACCGCGCCGGGTACCGACCTGTCGCTGAGCTTCCCGGGTTCGCTGGGTGGCATGAACTGGGGCAGCATTTCTACCGAcc                                          <  6:45917/149‑1 (MQ=255)
                                                                                                             ccACCCCGTATGTCCAGCTGCTATGCCGCATCGACTTCAAGAAGATGCGCTACGACGGCCTGTACACCGCGCCGGGTACCGACCTGTCGCTGAGCTTCCCGGGTTCGCTGGGTGGCATGAACTGGGGCAGCATTTCTACCGACCCGGTa                                     >  5:28193/1‑149 (MQ=255)
                                                                                                                       tGACCAGCTGCTATGCCGCATCGACTTCAAGAAGATGCGCTACGACGGCCTGTACACCGCTCCGGGTACCGACCTGTCGCTGAGCTTCCCGGGTTCGCTGGGTGGCATGAACTGGGGCAGCATTTCTACCGACCCGGTACATGGCTTc                            <  1:244998/148‑1 (MQ=255)
                                                                                                                            aGCTGCTATGCCGCATCGACTTCAAGAAGATGCGCTACGACGGCCTGTACACCGCGCCGGGTACCGACCTGTCGCTGAGCTTCCCGGGTTCGCTGGGTGGCATGAACTGGGGCAGCATTTCTACCGACCCGGTACATGGCTTCATCTTc                      <  4:223631/149‑1 (MQ=255)
                                                                                                                            aGCTGCTATGCCGCATCGACTTCAAGAAGATGCGCTACGACGGCCTGTACACCGCGCCGGGTACCGACCTGTCGCTGAGCTTCCCGGGTTCGCTGGGTGGCATGAACTGGGGCAGCATTTCTACCGACCCGGTACATGGCTTCATCTTc                      >  8:150521/1‑149 (MQ=255)
                                                                                                                             gctgctATGCCGCATCGACTTCAAGAAGATGCGCTACGACGGCCTGTACACCGCGCCGGGTACCGACCTGTCGCTGAGCTTCCCGGGTTCGCTGGGTGGCATGAACTGGGGCAGCATTTCTACCGACCCGGTACATGGCTTCATCTTCg                     <  7:150521/149‑1 (MQ=255)
                                                                                                                                gctATGCCGCATCGACTTCAAGAAGATGCGCTACGACGGCCTGTACACCGCGCCGGGTACCGACCTGTCGCTGAGCTTCCCGGGTTCGCTGGGTGGCATGAACTGGGGCAGCATTTCTACCGACCCGGTACATGGCTTCATCTTc                      <  2:45176/145‑1 (MQ=255)
                                                                                                                                gctATGCCGCATCGACTTCAAGAAGATGCGCTACGACGGCCTGTACACCGCGCCGGGTACCGACCTGTCGCTGAGCTTCCCGGGTTCGCTGGGTGGCATGAACTGGGGCAGCATTTCTACCGACCCGGTACATGGCTTCATCTTc                      >  1:45176/1‑145 (MQ=255)
                                                                                                                                   aTGCCGCATCGACTTCAAGAAGATGCGCTACGACGGCCTGTACACCGCGCCGGGTACCGACCTGTCGCTGAGCTTCCCGGGTTCGCTGGGTGGCATGAACTggg                                                            >  1:224681/1‑104 (MQ=255)
                                                                                                                                   aTGCCGCATCGACTTCAAGAAGATGCGCTACGACGGCCTGTACACCGCGCCGGGTACCGACCTGTCGCTGAGCTTCCCGGGTTCGCTGGGTGGCATGAACTggg                                                            <  2:224681/104‑1 (MQ=255)
                                                                                                                                    tGCCGCATCGACTTCAAGAAGATGCGCTACGACGGCCTGTACACCGCGCCGGGTACCGACCTGTCGCTGAGCTTCCCGGGTTCGCTGGGTGGCATGAACTGGGGCAGCATTTCTACCGACCCGGTACATGGCTTCATCTTCGTCAACGa              <  4:41748/149‑1 (MQ=255)
                                                                                                                                    tGCCGCATCGACTTCAAGAAGATGCGCTACGACGGCATGTACACCGCGCCGGGTACCGACCTGTCGCTGAGCTTCCCGGGTTCTCTGTGTGGCATTAACTGGTGCAGCATTTCTACCGACCCGGTACATGGCTTCATCTTCGTCAAc                <  7:290207/147‑1 (MQ=255)
                                                                                                                                      gcGCATCTTCTTGAAGTCGATGCGCTACGACGGCCTGTACACCGCGCCGGGTACCGACCTGTCGCTGAGCTTCCCGGGTTCGCTGGGTGGCATGAACTGGGGCAGCATTTCTACCGACCCGGTACATGGCTTCATCTTCGTCAACGACa            >  5:17526/2‑149 (MQ=255)
                                                                                                                                      ccGCATCGACTTCAAGAAGATGCGCTACGACGGCCTGTACACCGCGCCGGGTACCGACCTGTCGCTGAGCTTCCCGGGTTCGCTGGGTGGCATGAACTGGGGCAGCATTTCTACCGACCCGGTACATGGCTTCATCTTCGTCAACGACa            >  4:219630/1‑149 (MQ=255)
                                                                                                                                      ccGCATCGACTTCAAGAAGATGCGCTACGACGGCCTGTACACCGCGCCGGGTACCGACCTGTCGCTGAGCTTCCCGGGTTCGCTGGGTGGCATGAACTGGGGCAGCATTTCTACCGACCCGGTACATGGCTTCATCTTCGTCAACGACa            >  6:158169/1‑149 (MQ=255)
                                                                                                                                       cGCATCGACTTCAAGAAGATGCGCTACGACGGCCTGTACACCGCGCCGGGTACCGACCTGTCGCTGAGCTTCCCGGGTTCGCTGGGTGGCATGAACTGGGGCAGCATTTCAACCGACCCGGTACATGGCTTCATCTTCGTCAACGACAt           >  5:280085/1‑149 (MQ=255)
                                                                                                                                             gACTTCAAGAAGATGCGCTACGACGGCCTGTACACCGCGCCGGGTACCGACCTGTCGCTGAGCTTCCCGGGTTCGCTGGGTGGCATGAACTGGGGCAGCATTTCTACCGACCCGGTACATGGCTTCATCTTCGTCAACGACATGCGTc      <  6:28193/148‑1 (MQ=255)
                                                                                                                                             gACTTCAAGAAGATGCGCTACGACGGACTGTACACCGCGCCGGGTACCGACCTGTCGCTGAGCTTCCCGGGTTCGCTGTGTGGCATGAACTGGGGCAGCATTTCTACCGACCAGGTACATGGCTTCATCTTCGTCAACGACATGCTTCt     >  2:55240/1‑149 (MQ=255)
                                                                                                                                                ttCAAGAAGATGCGCTACGACGGCCTGTACACCGCGCCGGGTACCGACCTGTCGCTGAGCTTCCCGGGTTCGCTGGGTGGCATGAACTGGGGCAGCATTTCTACCGACCCGGTACATGGCTTCATCTTCGTCAACGACATGCGTCTggg  >  5:131153/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GTCAAGCCAAGCAACATCCCGAACGAGCCGTACTCCCCAACCCAGCCCAAGTCGGTGGGCATGCCACAGATCGGCGCGCAGACCCTGACCGAATCGGACATGTGGGGCGCCACCCCGTATGACCAGCTGCTATGCCGCATCGACTTCAAGAAGATGCGCTACGACGGCCTGTACACCGCGCCGGGTACCGACCTGTCGCTGAGCTTCCCGGGTTCGCTGGGTGGCATGAACTGGGGCAGCATTTCTACCGACCCGGTACATGGCTTCATCTTCGTCAACGACATGCGTCTGGG  >  NC_002947/4049034‑4049326

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