Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,398,903 C→G 7.8% intergenic (+123/‑119) aroQ‑III → / → PP_3004 type II 3‑dehydroquinate dehydratase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,398,9030CG 7.8% 160.0 / 4.6 51intergenic (+123/‑119)aroQ‑III/PP_3004type II 3‑dehydroquinate dehydratase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (23/24);  new base G (2/2);  total (25/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 = 4.26e-01

CGTGCGGCGATGAGCATGGTCTGAGCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTACAACGATCAGTCAAC  >  NC_002947/3398757‑3399049
                                                                                                                                                  |                                                                                                                                                  
cgtgcgGCGATGAGCATGGTCTGAGCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGc                                                                                                                                                    <  5:204496/147‑1 (MQ=255)
     ggcgATGAGCATGGTCTGAGCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCatggatg                                                                                                                                             >  6:110060/1‑149 (MQ=255)
         aTGAGCATGGTCTGAGCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCATGctggatgggta                                                                                                                                         <  4:211084/149‑4 (MQ=255)
         aTGAGCATGGTCTGAGCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGaaa                                                                                                                                         <  8:12803/149‑1 (MQ=255)
          tGAGCATGGTCTGAGCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAg                                                                                                                                        <  3:256764/149‑1 (MQ=255)
            aGCATGGTCTGAGCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGc                                                                                                                                       <  4:118077/148‑1 (MQ=255)
                  gTCTGAGCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGc                                                                                                                                >  4:143221/1‑149 (MQ=255)
                       aGCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACtt                                                                                                                           >  1:133677/1‑149 (MQ=255)
                             gcgcAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGa                                                                                                                     <  8:254687/149‑1 (MQ=255)
                                        ggATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTa                                                                                                           <  1:240141/148‑1 (MQ=255)
                                             tcCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTcc                                                                                                     >  1:166939/1‑149 (MQ=255)
                                             tcCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTcc                                                                                                     >  1:172106/1‑149 (MQ=255)
                                               cGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACt                                                                                                                   <  8:97888/133‑1 (MQ=255)
                                               cGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACt                                                                                                                   >  7:97888/1‑133 (MQ=255)
                                                gAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTg                                                                                                                  >  7:202786/1‑133 (MQ=255)
                                                gAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTg                                                                                                                  <  8:202786/133‑1 (MQ=255)
                                                  gcgcGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCtt                                                                                                <  2:8236/149‑1 (MQ=255)
                                                    gcgTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTc                                                                                               >  5:142793/1‑148 (MQ=255)
                                                      gTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGt                                                                                                       <  7:15426/138‑1 (MQ=255)
                                                      gTCTACCAGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGt                                                                                                       >  8:15426/1‑138 (MQ=255)
                                                       tCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTcaac                                                                                            >  1:242897/1‑148 (MQ=255)
                                                            ccGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCa                                                                                      >  6:30537/1‑149 (MQ=255)
                                                             cgcgAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGGGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCAc                                                                                     <  8:66135/149‑1 (MQ=255)
                                                             cgcgAAGAAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGGGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCAc                                                                                     <  8:61760/149‑1 (MQ=255)
                                                               cgAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCAccc                                                                                   <  5:110060/149‑1 (MQ=255)
                                                                gAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCg                                                                                  >  3:68633/1‑149 (MQ=255)
                                                                      ggTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCATGctggatgggt                                                                                                                                          <  4:58317/87‑3 (MQ=255)
                                                                      ggTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCATGctggatgggt                                                                                                                                          >  3:58317/1‑85 (MQ=255)
                                                                          gcgcGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACg                                                                                                         >  3:91690/1‑116 (MQ=255)
                                                                          gcgcGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACg                                                                                                         <  4:91690/116‑1 (MQ=255)
                                                                                            cGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCgg                                                      <  4:68633/149‑1 (MQ=255)
                                                                                             ggCGCCTGAAGTGACCGTAGGTGTGGTGGAAGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCAGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCggg                                                     >  6:170732/1‑149 (MQ=255)
                                                                                              gcgcCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCgggg                                                    <  2:133677/149‑1 (MQ=255)
                                                                                               cgcCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCAc                                                             >  3:101289/1‑139 (MQ=255)
                                                                                               cgcCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCAc                                                             <  4:101289/139‑1 (MQ=255)
                                                                                                     aaGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCtctc                                               <  6:142793/147‑1 (MQ=255)
                                                                                                       gTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGa                                           >  2:263795/1‑149 (MQ=255)
                                                                                                                        ggATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAgg                                  >  1:238714/1‑141 (MQ=255)
                                                                                                                        ggATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAgg                                  <  2:238714/141‑1 (MQ=255)
                                                                                                                          aTGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGAc                        <  5:30537/149‑1 (MQ=255)
                                                                                                                                     ccTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAg                                   >  1:180444/1‑127 (MQ=255)
                                                                                                                                     ccTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAg                                   <  2:180444/127‑1 (MQ=255)
                                                                                                                                       tACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGa                                    >  4:14711/1‑124 (MQ=255)
                                                                                                                                       tACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGa                                    <  3:14711/124‑1 (MQ=255)
                                                                                                                                        ttccatccatGCTGGATGAAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTACAACGAt          >  4:15209/3‑149 (MQ=255)
                                                                                                                                         cccATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCg                               >  8:125504/1‑127 (MQ=255)
                                                                                                                                         cccATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCg                               <  7:125504/127‑1 (MQ=255)
                                                                                                                                         cccATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTACAACGATc         <  2:172106/149‑1 (MQ=255)
                                                                                                                                         cccATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTACAACGATc         <  2:166939/149‑1 (MQ=255)
                                                                                                                                          ccATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTACAACGAt          >  1:92400/1‑147 (MQ=255)
                                                                                                                                                 gCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGCCCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTACAACGATCAGTcaac  >  4:214412/1‑148 (MQ=255)
                                                                                                                                                 gCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTACAACGATCAGTcaac  >  4:110745/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CGTGCGGCGATGAGCATGGTCTGAGCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTACAACGATCAGTCAAC  >  NC_002947/3398757‑3399049

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