Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,880,579 +A 58.8% coding (1006/1053 nt) pdxA ← 4‑hydroxythreonine‑4‑phosphate dehydrogenase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,880,57910.→A58.8% 73.4 ‑1.1 ‑0.951coding (1006/1053 nt)pdxA4‑hydroxythreonine‑4‑phosphate dehydrogenase
Reads supporting (aligned to +/- strand):  new base (24/6):  ref base (14/7):  total (38/13)
Fisher's exact test strand distribution p-value = 3.38e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.24e-01

CGCGCAAGGGCGGCAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GC.........C..CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGC  >  NC_014034/2880480‑2880635
                                                                                                    |                                                                            
cgcgCAAGGGCGGCAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCcggcg..........g                                                                    >  2:4131015/1‑101 (MQ=255)
 gcgcAAGGGCGGCAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCggcg..........gc                                                                   <  2:5307972/101‑1 (MQ=255)
  cgcAAGGGCGGCAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCggcg..........gcg                                                                  <  1:6510548/101‑2 (MQ=255)
   gcAAGGGCGGCAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCggcg..........gcgc                                                                 >  1:1229102/1‑99 (MQ=255)
   gcAAGGGCGGCAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGGCCGAGGGGAGAGCCCCCggcg..........gcgc                                                                 >  2:5356768/1‑99 (MQ=255)
    cAAGGGCGGCCGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCACCTATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GC.........g..cg                                                     >  1:430535/1‑101 (MQ=37)
    cAAGGGCGGCAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GC.........g..cg                                                     >  1:1093948/1‑101 (MQ=255)
    cAAGGGCGGCAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GC.........g..cg                                                     >  1:5937049/1‑101 (MQ=255)
     aaGGGCGGCAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGTGAGGTGGGGGGCGGAGGCCCcggcg..........gggggt                                                               >  1:6938373/1‑96 (MQ=255)
     aaGGGCGGCAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GC.........g..cgt                                                    >  2:312885/1‑100 (MQ=255)
     aaGGGCGGCAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GC.........g..cgt                                                    >  2:5037859/1‑100 (MQ=255)
      aGGGCGGCAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GC.........g..cgtc                                                   >  1:358983/1‑99 (MQ=255)
       gggCGGCAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATACCGCGGGTCCGGGGGCAGGGCCCCcgggg..........gagggtgc                                                             >  1:3335745/1‑94 (MQ=255)
        ggcggcAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGTGGGTCCGGGGGCAGAGCCCCCGGCG..........GC.........g..cgtcac                                                 <  2:3699044/101‑5 (MQ=255)
        ggcggcAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGTGGGCAGAGCCCCCGGCG..........GC.........g..cgtcac                                                 >  1:4190755/1‑97 (MQ=255)
        ggcggcAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GC.........g..cgtcac                                                 >  1:45703/1‑97 (MQ=255)
        ggcggcAGCCCGTCAATCGCAGCCAAAGCCCGCCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GC.........g..cgtcac                                                 <  1:5577049/101‑5 (MQ=255)
        aatgacAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GC.........g..cgtcac                                                 >  2:4473962/6‑97 (MQ=25)
          cggcAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GC.........g..cgtcacag                                               >  1:2537502/1‑95 (MQ=255)
          cggcAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GC.........g..cgtcacag                                               <  1:5714889/101‑7 (MQ=255)
          cggcAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GC.........g..cgtcacag                                               >  2:6570197/1‑95 (MQ=255)
            gcAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GC.........g..cgtcacagcc                                             >  1:1917141/1‑93 (MQ=255)
             cAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GC.........g..cgtcacagccg                                            <  1:980955/101‑10 (MQ=25)
                                                 ttttCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCcgc                   >  2:706357/1‑101 (MQ=2)
                                                  tttCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCcgcg                  >  2:53308/1‑101 (MQ=2)
                                                  tttCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCcgcg                  >  2:2984580/1‑101 (MQ=2)
                                                   ttCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGa                 >  1:3473087/1‑101 (MQ=9)
                                                   ttCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGa                 >  1:306240/1‑101 (MQ=9)
                                                    tCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGAt                <  1:4539919/101‑1 (MQ=9)
                                                    tCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGAGAt                >  1:1067246/1‑101 (MQ=2)
                                                     cATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGCTc               >  2:4414547/1‑101 (MQ=2)
                                                     cATCTTGCTCCAAATACCCCGGGGGTCCGGGGGCAGGCCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATc               <  1:1847620/101‑1 (MQ=9)
                                                       tCTTGCTTCAAATATCCCGGGGGGCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAg             <  2:5616254/101‑1 (MQ=9)
                                                        cTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGGGATCAGa            >  1:1040235/1‑101 (MQ=9)
                                                         ttGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGAc           >  1:2662480/1‑101 (MQ=9)
                                                          tGCTTCAAATATCCGGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACt          <  2:3975010/101‑1 (MQ=9)
                                                          tGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACt          >  1:2039672/1‑101 (MQ=9)
                                                          tGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACt          >  1:2658808/1‑101 (MQ=9)
                                                          tGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACt          >  1:6016920/1‑101 (MQ=9)
                                                          tGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACt          >  1:6544169/1‑101 (MQ=9)
                                                          tGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACt          >  1:6171083/1‑101 (MQ=9)
                                                           gCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACtg         >  2:3709694/1‑101 (MQ=12)
                                                           gCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACtg         >  1:4170815/1‑101 (MQ=12)
                                                            cTTCAAATACCCCGGGGGTCGGGGGGCAGAGCCCCCGGCG..........GCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACtgt        <  1:626368/101‑1 (MQ=12)
                                                             ttCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACtgtg       <  1:6870740/101‑1 (MQ=12)
                                                             ttCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACtgtg       >  2:237304/1‑101 (MQ=12)
                                                             ttCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGAGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACtgtg       >  2:4898503/1‑101 (MQ=12)
                                                              tCAAATATCCCGGGGGTCCGGGGGGAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGc      <  2:6938000/101‑1 (MQ=12)
                                                              tCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGc      >  1:1849678/1‑101 (MQ=12)
                                                               cAAATATCCGGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCg     <  2:1229054/101‑1 (MQ=12)
                                                               cAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCg     >  1:340773/1‑101 (MQ=12)
                                                               cAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCg     >  1:4542350/1‑101 (MQ=12)
                                                                aaaTATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCgg    >  2:4943680/1‑101 (MQ=12)
                                                                 aaTATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCggg   >  1:5317028/1‑101 (MQ=12)
                                                                 aaTATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCggg   >  1:1133669/1‑101 (MQ=12)
                                                                  atatCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCCAGGCCGCGATCAGACTGTGCGGGc  >  2:4748683/1‑101 (MQ=12)
                                                                  atatCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGc  <  2:4602184/101‑1 (MQ=12)
                                                                  atatCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGc  >  1:2411654/1‑101 (MQ=12)
                                                                  atatCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGc  >  2:5768259/1‑101 (MQ=12)
                                                                  atatCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GCGCGTCACAGCCGCGCCGCCCCCACCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGc  >  1:4268087/1‑101 (MQ=12)
                                                                                                    |                                                                            
CGCGCAAGGGCGGCAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GC.........C..CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGC  >  NC_014034/2880480‑2880635

Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 26 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG