Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,049,176 A→T 7.0% D576V (GAC→GTC)  quiA → quinate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,049,1760AT 7.0% 166.2 / 3.1 57D576V (GAC→GTC) quiAquinate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (24/28);  new base T (2/2);  total (27/30)
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.95e-01

TACCGGTCAAGCCAAGCAACATCCCGAACGAGCCGTACTCCCCAACCCAGCCCAAGTCGGTGGGCATGCCACAGATCGGCGCGCAGACCCTGACCGAATCGGACATGTGGGGCGCCACCCCGTATGACCAGCTGCTATGCCGCATCGACTTCAAGAAGATGCGCTACGACGGCCTGTACACCGCGCCGGGTACCGACCTGTCGCTGAGCTTCCCGGGTTCGCTGGGTGGCATGAACTGGGGCAGCATTTCTACCGACCCGGTACATGGCTTCATCTTCGTCAACGACATGCGTCT  >  NC_002947/4049029‑4049323
                                                                                                                                                   |                                                                                                                                                   
tACCGGTCAAGCCAAGCAACATCCCGAACGAGCCGTACTCCCCAACCCAGCCCAAGTCGGTGGGCATGCCACAGATCGGCGCGCAGACCCTGACCGAATCGGACATGTGGGGCGCCACCCCGTATGACCAGCTGCTATGCCGCATCGAc                                                                                                                                                    <  5:46807/149‑1 (MQ=255)
     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                                                                                                                                              >  6:302526/1‑124 (MQ=255)
                               gCCGTACTCCCCAACCCAGCCCAAGTCGGTGGGCATGCCACAGATCGGCGCGCAGACCCTGACCGAATCGGACATGTGGGGCGCCACCCCGTATGACCAGCTGCTATGCCGCATCGACTTCaag                                                                                                                                              <  5:302526/124‑1 (MQ=255)
                                ccGTACTCCCCAACCCAGCCCAAGTCGGTGGGCATGCCACAGATCGGCGCGCAGACCCTGCCCGAATCGGACATGTGGGGCGCCACCCCGTATGACCAGCTGCTATTCAGCATAGCCTTCAAGAAGATGCCCTACGACGGACTGAacac                                                                                                                    >  7:248580/1‑149 (MQ=255)
                                  gTACTCCCCAACCCAGCCCAAGTCGGTGGGCATGCCAAAGATCGGCGCGCAGACCCTGACCGAATCGGACATGTGGGGCGCCACCCCGTATGACCAGCTGCTATGCCGCATCGACTTCAAGAAGATGCGCTACGACGGCCTGTACACcg                                                                                                                  >  3:223631/1‑149 (MQ=255)
                                        cccAACCCAGCCCAAGTCGGTGGGCATGCCACAGATCGGCGCGCAGACCCTGACCGAATCGGACATGTGGGGCGCCACCCCGTATGACCAGCTGCTATGCCGCATCGACtt                                                                                                                                                  <  8:161374/111‑1 (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                                                                                                    <  8:50849/139‑1 (MQ=255)
                                                          ggTGGGCATGCCACAGATCGGCGCGCAGACCCTGACCGAATCGGACATGTGGGGCGCCACCCCGTATGACCAGCTGCTATGCCGCATCGACTTCAAGAAGATGCGCTACGACGGCCTGTACACCGCGCCGGGTACCGAc                                                                                                    >  7:50849/1‑139 (MQ=255)
                                                                    ccACAGATCGGCGCGCAGACCCTGACCGAATCGGACATGTGGGGCGCCACCCCGTATGACCAGCTGCTATGCCGCATCGACTTCAAGAAGATGCGCTACGACGGCCTGTACACCGCGCCGGGTAcc                                                                                                       <  6:185812/126‑1 (MQ=255)
                                                                    ccACAGATCGGCGCGCAGACCCTGACCGAATCGGACATGTGGGGCGCCACCCCGTATGACCAGCTGCTATGCCGCATCGACTTCAAGAAGATGCGCTACGACGGCCTGTACACCGCGCCGGGTAcc                                                                                                       >  5:185812/1‑126 (MQ=255)
                                                                            cGGCGCGCAGACCCTGACCGAATCGGACATGTGGGGCGCCACCCCGTATGACCAGCTGCTATGCCGCATCGACTTCAAGAAGATGCGCTACGACGGCCTGTACACCGCGCCGGGTACCGACCTGTCGCTGAGCTTCCCGGGtt                                                                              <  5:198628/143‑1 (MQ=255)
                                                                            cGGCGCGCAGACCCTGACCGAATCGGACATGTGGGGCGCCACCCCGTATGACCAGCTGCTATGCCGCATCGACTTCAAGAAGATGCGCTACGACGGCCTGTACACCGCGCCGGGTACCGACCTGTCGCTGAGCTTCCCGGGtt                                                                              >  6:198628/1‑143 (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                                                            >  4:296224/1‑149 (MQ=255)
                                                                                        ccTGACCGAATCGGACATGTGGGGCGCCACCCCGTATGACCAGCTGCTATGCCGCATCGACTTCAAGAAGATGCGCTACGACGGCCTGTACACCGCGCCGGGTACCGACCTGTCGCTGAGCTTCCCGGGTTCGCTGGGTGGCATGAACt                                                            <  5:234824/149‑1 (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                                                         <  2:224681/104‑1 (MQ=255)
                                                                                                                                        aTGCCGCATCGACTTCAAGAAGATGCGCTACGACGGCCTGTACACCGCGCCGGGTACCGACCTGTCGCTGAGCTTCCCGGGTTCGCTGGGTGGCATGAACTggg                                                         >  1:224681/1‑104 (MQ=255)
                                                                                                                                         tGCCGCATCGACTTCAAGAAGATGCGCTACGACGGCCTGTACACCGCGCCGGGTACCGACCTGTCGCTGAGCTTCCCGGGTTCGCTGGGTGGCATGAACTGGGGCAGCATTTCTACCGACCCGGTACATGGCTTCATCTTCGTCAACGa           <  4:41748/149‑1 (MQ=255)
                                                                                                                                         tGCCGCATCGACTTCAAGAAGATGCGCTACGACGGCATGTACACCGCGCCGGGTACCGACCTGTCGCTGAGCTTCCCGGGTTCTCTGTGTGGCATTAACTGGTGCAGCATTTCTACCGACCCGGTACATGGCTTCATCTTCGTCAAc             <  7:290207/147‑1 (MQ=255)
                                                                                                                                           gcGCATCTTCTTGAAGTCGATGCGCTACGACGGCCTGTACACCGCGCCGGGTACCGACCTGTCGCTGAGCTTCCCGGGTTCGCTGGGTGGCATGAACTGGGGCAGCATTTCTACCGACCCGGTACATGGCTTCATCTTCGTCAACGACa         >  5:17526/2‑149 (MQ=255)
                                                                                                                                           ccGCATCGACTTCAAGAAGATGCGCTACGACGGCCTGTACACCGCGCCGGGTACCGACCTGTCGCTGAGCTTCCCGGGTTCGCTGGGTGGCATGAACTGGGGCAGCATTTCTACCGACCCGGTACATGGCTTCATCTTCGTCAACGACa         >  6:158169/1‑149 (MQ=255)
                                                                                                                                           ccGCATCGACTTCAAGAAGATGCGCTACGACGGCCTGTACACCGCGCCGGGTACCGACCTGTCGCTGAGCTTCCCGGGTTCGCTGGGTGGCATGAACTGGGGCAGCATTTCTACCGACCCGGTACATGGCTTCATCTTCGTCAACGACa         >  4:219630/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)
                                                                                                                                                   |                                                                                                                                                   
TACCGGTCAAGCCAAGCAACATCCCGAACGAGCCGTACTCCCCAACCCAGCCCAAGTCGGTGGGCATGCCACAGATCGGCGCGCAGACCCTGACCGAATCGGACATGTGGGGCGCCACCCCGTATGACCAGCTGCTATGCCGCATCGACTTCAAGAAGATGCGCTACGACGGCCTGTACACCGCGCCGGGTACCGACCTGTCGCTGAGCTTCCCGGGTTCGCTGGGTGGCATGAACTGGGGCAGCATTTCTACCGACCCGGTACATGGCTTCATCTTCGTCAACGACATGCGTCT  >  NC_002947/4049029‑4049323

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