Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,049,184 A→T 7.5% K579* (AAG→TAG) ‡ quiA → quinate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,049,1840AT 7.5% 155.1 / 3.2 53K579* (AAG→TAG) ‡quiAquinate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (25/24);  new base T (2/2);  total (27/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 = 7.79e-01

CAACATCCCGAACGAGCCGTACTCCCCAACCCAGCCCAAGTCGGTGGGCATGCCACAGATCGGCGCGCAGACCCTGACCGAATCGGACATGTGGGGCGCCACCCCGTATGACCAGCTGCTATGCCGCATCGACTTCAAGAAGATGCGCTACGACGGCCTGTACACCGCGCCGGGTACCGACCTGTCGCTGAGCTTCCCGGGTTCGCTGGGTGGCATGAACTGGGGCAGCATTTCTACCGACCCGGTACATGGCTTCATCTTCGTCAACGACATGCGTCTGGGCCTGT  >  NC_002947/4049045‑4049331
                                                                                                                                           |                                                                                                                                                   
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)
                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                                                                                      <  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                                                                    <  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                           >  1:45176/1‑145 (MQ=255)
                                                                                                                     gctATGCCGCATCGACTTCAAGAAGATGCGCTACGACGGCCTGTACACCGCGCCGGGTACCGACCTGTCGCTGAGCTTCCCGGGTTCGCTGGGTGGCATGAACTGGGGCAGCATTTCTACCGACCCGGTACATGGCTTCATCTTc                           <  2:45176/145‑1 (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)
                                                                                                                                          gaagaTGCTCTACGACGGCCTGTACACCGCGCCGGGTACCGACCTGTCGCTGAGCTTCCCGGGTTCGCTGGGTGGCATGAACTGGGGCATCATTTCTACCGAACCGGTACATGGCTTCATCTTCGTCAACGACATGCGTCTGGGCCtgt  <  1:55240/149‑1 (MQ=255)
                                                                                                                                          gaagaTGCGCTACGACGGCCTGTACACCGCGCCGGGTACCGACCTGTCGCTGAGCTTCCCGGGTTCGCTGGGTGGCATGAACTGGGGCAGCATTTCAACCGACCCGGTACATGGCTTCATCTTCGTCAACGACATGCGTCTGGGCCtgt  <  6:280085/149‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                   
CAACATCCCGAACGAGCCGTACTCCCCAACCCAGCCCAAGTCGGTGGGCATGCCACAGATCGGCGCGCAGACCCTGACCGAATCGGACATGTGGGGCGCCACCCCGTATGACCAGCTGCTATGCCGCATCGACTTCAAGAAGATGCGCTACGACGGCCTGTACACCGCGCCGGGTACCGACCTGTCGCTGAGCTTCCCGGGTTCGCTGGGTGGCATGAACTGGGGCAGCATTTCTACCGACCCGGTACATGGCTTCATCTTCGTCAACGACATGCGTCTGGGCCTGT  >  NC_002947/4049045‑4049331

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