Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,880,579 +A 61.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.→A61.8% 42.4 ‑2.0 ‑1.834coding (1006/1053 nt)pdxA4‑hydroxythreonine‑4‑phosphate dehydrogenase
Reads supporting (aligned to +/- strand):  new base (21/0):  ref base (11/2):  total (32/2)
Fisher's exact test strand distribution p-value = 1.39e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.75e-02

GCAAGGGCGGCAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GC.........C..CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGC  >  NC_014034/2880483‑2880635
                                                                                                 |                                                                            
gcAAGGGCGGCAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCggcg..........gcgc                                                                 >  2:464558/1‑99 (MQ=255)
 cAAGGGCGGCAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GC.........g..cg                                                     >  2:834944/1‑101 (MQ=255)
  aaGGGCGGCAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGNGGTCCGGGGGCAGAGCCCCCggcg..........gcggcc                                                               >  1:137584/1‑97 (MQ=255)
  aaGGGCGGCAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GG.........c..cgt                                                    >  1:1830145/1‑100 (MQ=255)
  aaGGGCGGCAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GC.........g..cgt                                                    >  1:3203415/1‑100 (MQ=255)
  aaGGGCGGCAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GC.........g..cgt                                                    >  2:3162036/1‑100 (MQ=255)
  aaGGGCGGCAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GC.........g..cgt                                                    >  1:4598412/1‑100 (MQ=255)
   aGGGCGGCAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GC.........g..cgtc                                                   >  2:3820948/1‑99 (MQ=255)
   aGGGCGGCAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GC.........g..cggc                                                   >  1:188331/1‑99 (MQ=255)
      gcggcAGCCCGTCAATCGCACCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATACCCCGGGGGTCCGGGGGCAGAGCCCCGGGCG..........GC.........g..cgtcaca                                                <  2:3169527/101‑6 (MQ=25)
        ggcAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GC.........g..cgtcacagc                                              >  2:3253654/1‑94 (MQ=255)
         gcAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GC.........g..cgtcacagcc                                             >  2:422996/1‑93 (MQ=255)
         ccAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GC.........g..cgtcacagcc                                             <  2:2748722/100‑9 (MQ=25)
                                               tttCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCcgcg                  >  1:767228/1‑101 (MQ=2)
                                               tttCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCcgcg                  >  1:445722/1‑101 (MQ=2)
                                                ttCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGa                 >  1:4653693/1‑101 (MQ=9)
                                                ttCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGa                 >  2:706911/1‑101 (MQ=9)
                                                 tCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGAt                >  2:2608944/1‑101 (MQ=9)
                                                  cATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATc               >  2:2880881/1‑101 (MQ=9)
                                                   aTCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCa              >  1:2415107/1‑101 (MQ=9)
                                                     cTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGa            >  1:1316668/1‑101 (MQ=9)
                                                     cTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGa            >  1:2202890/1‑101 (MQ=9)
                                                        gCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACtg         >  1:454003/1‑101 (MQ=12)
                                                         cTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACtgt        >  2:1226658/1‑101 (MQ=12)
                                                         cTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACtgt        >  1:2895720/1‑101 (MQ=12)
                                                          ttCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACtgtg       >  1:3252792/1‑101 (MQ=12)
                                                          ttCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACtgtg       >  2:418576/1‑101 (MQ=12)
                                                           tCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGAcggaaa      >  1:55676/1‑96 (MQ=9)
                                                            cAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCg     >  1:527154/1‑101 (MQ=12)
                                                            cAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCg     >  2:3867631/1‑101 (MQ=12)
                                                             aaaTATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACCGC.........CG.CGCCGCCCCCCTCTTCGCCGCCCTCGTCAACGCCGCGATCAGACTGTGCgg    >  1:56261/1‑101 (MQ=9)
                                                             aaaTATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCgg    >  1:4808524/1‑101 (MQ=12)
                                                             aaaTATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCgg    >  1:3107660/1‑101 (MQ=12)
                                                             aaaTATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCgg    >  2:2291770/1‑101 (MQ=12)
                                                              aaTATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCggg   >  2:70316/1‑101 (MQ=12)
                                                               atatCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGc  >  1:2314399/1‑101 (MQ=12)
                                                               atatCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGc  >  2:4340487/1‑101 (MQ=12)
                                                               atatCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGc  >  2:3464172/1‑101 (MQ=12)
                                                               atatCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGACCAGACTGTGCGGGc  >  2:4115729/1‑101 (MQ=12)
                                                                                                 |                                                                            
GCAAGGGCGGCAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GC.........C..CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGC  >  NC_014034/2880483‑2880635

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