Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,880,582 +G 18.9% coding (1003/1053 nt) pdxA ← 4‑hydroxythreonine‑4‑phosphate dehydrogenase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,880,5822.→G18.9% 19.2 ‑1.1 ‑0.037coding (1003/1053 nt)pdxA4‑hydroxythreonine‑4‑phosphate dehydrogenase
Reads supporting (aligned to +/- strand):  new base (4/3):  ref base (24/6):  total (28/9)
Fisher's exact test strand distribution p-value = 3.27e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.87e-01

CAAGGGCGGCAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GC.........C..CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCA  >  NC_014034/2880484‑2880673
                                                                                                                      |                                                                                            
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)
 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)
    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:5714889/101‑7 (MQ=255)
      cggcAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GC.........g..cgtcacag                                                                                     >  1:2537502/1‑95 (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:306240/1‑101 (MQ=9)
                                               ttCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGa                                                       >  1:3473087/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                                               >  1:4170815/1‑101 (MQ=12)
                                                       gCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACtg                                               >  2:3709694/1‑101 (MQ=12)
                                                        cTTCAAATACCCCGGGGGTCGGGGGGCAGAGCCCCCGGCG..........GCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACtgt                                              <  1:626368/101‑1 (MQ=12)
                                                         ttCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACtgtg                                             >  2:237304/1‑101 (MQ=12)
                                                         ttCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACtgtg                                             <  1:6870740/101‑1 (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:4542350/1‑101 (MQ=12)
                                                           cAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCg                                           >  1:340773/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                                        >  1:2411654/1‑101 (MQ=12)
                                                              atatCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGc                                        >  2:5768259/1‑101 (MQ=12)
                                                              atatCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGc                                        <  2:4602184/101‑1 (MQ=12)
                                                              atatCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GCGCGTCACAGCCGCGCCGCCCCCACCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGc                                        >  1:4268087/1‑101 (MQ=12)
                                                                                                   ctcacagcc.........c..cGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCa  <  2:6057350/93‑1 (MQ=255)
                                                                                                                      |                                                                                            
CAAGGGCGGCAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GC.........C..CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCA  >  NC_014034/2880484‑2880673

Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 27 ≤  ATCG  < 33 ≤  ATCG  < 41 ≤  ATCG