Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,717,162 C→T 5.8% V148I (GTT→ATT)  pdxY ← pyridoxal kinase 2/pyridoxine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,717,1620CT 5.8% 228.7 / 3.0 69V148I (GTT→ATT) pdxYpyridoxal kinase 2/pyridoxine kinase
Reads supporting (aligned to +/- strand):  ref base C (32/33);  new base T (2/2);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

AGTTCGCGCGCTGCCAGAACCGCTTCTTCGACGTTATTTACCGCATGCTCACAGAGTATTTCCAGCTCAACCAGATTTGGCGCAATGATATCGCTGGCAGGCAAACCGTGCCGCACATGAAACTCTGCGACACCCGGTGC  >  W3110S.gb/1717092‑1717231
                                                                      |                                                                     
aGTTCGCGCGCTGCCAGAACCGCTTCTTCGACGTTATTTACCGCATGCTCACAGAGTATTTCCAGCTCAAc                                                                       >  1:4579390/1‑71 (MQ=255)
    cgcgcgCTGCCAGAACCGCTTCTTCGACGTTATTTACCGCATGCTCACAGAGTATTTCCAGCTCAAc                                                                       >  1:2836740/1‑67 (MQ=255)
    cgcgcgCTGCCAGAACCGCTTCTTCGACGTTATTTACCGCATGCTCACAGAGTATTTCCAGCTCAAc                                                                       >  1:397620/1‑67 (MQ=255)
    cgcgcgCTGCCAGAACCGCTTCTTCGACGTTATTTACCGCATGCTCACAGAGTATTTCCAGCTCAAc                                                                       >  1:2259243/1‑67 (MQ=255)
    cgcgcgCTGCCAGAACCGCTTCTTCGACGTTATTTACCGCATGCTCACAGAGTATTTCCAGCTCAAc                                                                       >  1:1961541/1‑67 (MQ=255)
    cgcgcgCTGCCAGAACCGCTTCTTCGACGTTATTTACCGCATGCTCACAGAGTATTTCCAGCTCAAc                                                                       >  1:3763684/1‑67 (MQ=255)
    cgcgcgCTGCCAGAACCGCTTCTTCGACGTTATTTACCGCATGCTCACAGAGTATTTCCAGCTCAACCAg                                                                    <  1:3506093/70‑1 (MQ=255)
    cgcgcgCTGCCAGAACCGCTTCTTCGACGTTATTTACCGCATGCTCACAGAGTATTTCCAGCTCAACCAg                                                                    <  1:2029183/70‑1 (MQ=255)
     gcgcgcTGCCAGAACCGCTTCTTCGACGTTATTTACCGCATGCTCACAGAGTATTTCCAGCTCAACCAGa                                                                   <  1:3767618/70‑1 (MQ=255)
     gcgcgcTGCCAGAACCGCTTCTTCGACGTTATTTACCGCATGCTCACAGAGTATTTCCAGCTCAACCAGa                                                                   <  1:1969361/70‑1 (MQ=255)
      cgcgcTGCCAGAACCGCTTCTTCGACGTTATTTACCGCATGCTCACAGAGTATTTCCAGCTCAACCAGAtt                                                                 >  1:3424370/1‑71 (MQ=255)
       gcgcTGCCAGAACCGCTTCTTCGACGTTATTTACCGCATGCTCACAGAGTATTTCCAGCTCAACCAGAttt                                                                >  1:530225/1‑71 (MQ=255)
       gcgcTGCCAGAACCGCTTCTTCGACGTTATTTACCGCATGCTCACAGAGTATTTCCAGCTCAACCAGAttt                                                                <  1:2219437/71‑1 (MQ=255)
         gcTGCCAGAACCGCTTCTTCGACGTTATTTACCGCATGCTCACAGAGTATTTCCAGCTCAACCAGATTTgg                                                              >  1:3348385/1‑71 (MQ=255)
         gcTGCCAGAACCGCTTCTTCGACGTTATTTACCGCATGCTCACAGAGTATTTCCAGCTCAACCAGATTTg                                                               >  1:3989178/1‑70 (MQ=255)
            gCCAGAACCGCTTCTTCGACGTTATTTACCGCATGCTCACAGAGTATTTCCAGCTCAACCAGa                                                                   <  1:1990640/63‑1 (MQ=255)
            gCCAGAACCGCTTCTTCGACGTTATTTACCGCATGCTCACAGAGTATTTCCAGCTCAACCAGATTTGgcg                                                            >  1:3504686/1‑70 (MQ=255)
            gCCAGAACCGCTTCTTCGACGTTATTTACCGCATGCTCACAGAGTATTTCCAGCTCAACCAGATTTGgcg                                                            >  1:2212074/1‑70 (MQ=255)
               aGAACCGCTTCTTCGACGTTATTTACCGCATGCTCACAGAGTATTTCCAGCTCAACCAGATTTGGCGCAAt                                                        <  1:1565845/71‑1 (MQ=255)
                 aaCCGCTTCTTCGACGTTATTTACCGCATGCTCACAGAGTATTTCCAGCTCAACCAGATTTGGCGCAATGa                                                      <  1:4013855/71‑1 (MQ=255)
                   ccGCTTCTTCGACGTTATTTACCGCATGCTCACAGAGTATTTCCAGCTCAACCAGATTTGGCGCAATGa                                                      >  1:3212140/1‑69 (MQ=255)
                     gCTTCTTCGACGTTATTTACCGCATGCTCACAGAGTATTTCCAGCTCAACCAGAttt                                                                >  1:1914644/1‑57 (MQ=255)
                     gCTTCTTCGACGTTATTTACCGCATGCTCACAGAGTATTTCCAGCTCAACCAGAttt                                                                >  1:1065450/1‑57 (MQ=255)
                     gCTTCTTCGACGTTATTTACCGCATGCTCACAGAGTATTTCCAGCTCAACCAGAttt                                                                >  1:3779565/1‑57 (MQ=255)
                     gCTTCTTCGACGTTATTTACCGCATGCTCACAGAGTATTTCCAGCTCAACCAGATTTGGCGCAATGATATc                                                  <  1:3944165/71‑1 (MQ=255)
                      cttcttCGAGGTTATTTTCCGCATGCTCACAGAGTATTTCCAGCTCTACCAGATTTGGCGCa                                                          <  1:2419440/62‑1 (MQ=255)
                        tcttcGACGTTATTTACCGCATGCTCACAGAGTATTTCCAGCTCAACCAGATTTGGCGCAATGATATCGCt                                               >  1:4780178/1‑71 (MQ=255)
                          ttcGACGTTATTTACCGCATGCTCACAGAGTATTTCCAGCTCAACCAGATTTGGCGCAATGATATc                                                  >  1:3552382/1‑66 (MQ=255)
                           tcGACGTTATTTACCGCATGCTCACAGAGTATTTCCAGCTCAACCAGATTTGGCGCAATGATATCGCTgg                                             >  1:4055866/1‑70 (MQ=255)
                            cGACGTTATTTACCGCATGCTCACAGAGTATTTCCAGCTCAACCAGATTTGGCGCAATGata                                                    >  1:3335188/1‑62 (MQ=255)
                            cGACGTTATTTACCGCATGCTCACAGAGTATTTCCAGCTCAACCAGATTTGGCGCAATGATATCGCTgg                                             >  1:3662643/1‑69 (MQ=255)
                            cGACGTTATTTACCGCATGCTCACAGAGTATTTCCAGCTCAACCAGATTTGGCGCAATGATATCGCTgg                                             >  1:1436471/1‑69 (MQ=255)
                                gTTATTTACCGCATGCTCACAGAGTATTTCCAGCTCAACCAGATTTGGCGCAATGATATCGCTggcaggca                                       <  1:1392406/71‑1 (MQ=255)
                                 ttATTTACCGCATGCTCACAGAGTATTTCCAGCTCAACCAGATTTTGCGCAATGATATCGCTGGCAGGCaa                                      <  1:3978969/71‑1 (MQ=255)
                                 ttATTTACCGCATGCTCACAGAGTATTTCCAGCTCAACCAGATTTGGCGCaa                                                         <  1:3982527/52‑1 (MQ=255)
                                 ttATTTACCGCATGCTCACAGAGTATTTCCAGCTCAACCAGATTTGGCGCAAt                                                        <  1:3412798/53‑1 (MQ=255)
                                 ttATTTACCGCATGCTCACAGAGTATTTCCAGCTCAACCAGATTTGGCGCAAt                                                        <  1:1385300/53‑1 (MQ=255)
                                 ttATTTACCGCATGCTCACAGAGTATTTCCAGCTCAACCAGATTTGGCGCAATGATATCGCTGGCAGGCaa                                      <  1:1002229/71‑1 (MQ=255)
                                 ttATTTACCGCATGCTCACAGAGTATTTCCAGCTCAACCAGATTTGGCGCAATGATATCGCTGGCAGGCaa                                      <  1:1490134/71‑1 (MQ=255)
                                   aTTTACCGCATGCTCACAGAGTATTTCCAGCTCAACCAGATTTGGCGCAATGATATCGCTGGCAGGCAAAc                                    >  1:4334920/1‑71 (MQ=255)
                                   aTTTACCGCATGCTCACAGAGTATTTCCAGCTCAACCAGATTTGGCGCAATGATATCGCTGGCAGGCAAAc                                    >  1:1278629/1‑71 (MQ=255)
                                    tttACCGCATGCTCACAGAGTATTTCCAGCTCAACCAGATTTGGCGCAATGATATCGCTGGCAGGCAAAcc                                   <  1:2573541/71‑1 (MQ=255)
                                    tttACCGCATGCTCACAGAGTATTTCCAGCTCAACCAGATTTGGCGCAATGATATCGCTGGCAGGCAAAcc                                   <  1:2319943/71‑1 (MQ=255)
                                    tttACCGCATGCTCACAGAGTATTTCCAGCTCAACCAGATTTGGCGCAATGATATCGCTGGCAGGCAAAcc                                   >  1:474974/1‑71 (MQ=255)
                                      tACCGCATGCTCACAGAGTATTTCCAGCTCAACCAGATTTGGCGCAATGATATCGCTGGCAGGCAAACCg                                  <  1:261673/70‑1 (MQ=255)
                                       aCCGCATGCTCACAGAGTATTTCCAGCTCAACCAGATTTGGCTCAATGATATCGCTGGCAGGCAAACCGt                                 <  1:1781302/70‑1 (MQ=255)
                                            aTGCTCACAGAGTATTTCCAGCTCAACCAGATTTGGCGCAATGATATCGCTGGCAGGCAAACCGTGCCGca                           >  1:748577/1‑71 (MQ=255)
                                            aTGCTCACAGAGTATTTCCAGCTCAACCAGATTTGGCGCAATGATATCGCTGGCAGGCAAACCGTGCCGca                           >  1:222529/1‑71 (MQ=255)
                                             tGCTCACAGAGTATTTCCAGCTCAACCAGATTTGGCGCAATGATATCGCTgg                                             <  1:910483/52‑1 (MQ=255)
                                                tCACAGAGTATTTCCAGCTCAACCAGATTTGGCGCAAt                                                        <  1:3828639/38‑1 (MQ=255)
                                                tCACAGAGTATTTCCAGCTCAACCAGATTTGGCGCAATGATATCGCTGGAGGCAAACCGTg                                >  1:4687688/1‑61 (MQ=255)
                                                  acaGAGTATTTCCAGCTCAACCAGATTTGGCGCAATGATATCGCTGGCAGGCAAACCGTGCCGCTCATGa                      <  1:3515697/70‑1 (MQ=255)
                                                  acaGAGTATTTCCAGCTCAACCAGATTTGGCGCAATGATATCGCTGGCAGGCAAACCGTGCCGCACATGa                      <  1:4863880/70‑1 (MQ=255)
                                                   caGAGTATTTCCAGCTCAACCAGATTTGGCGCAATGATATCGCTGGCAGGCAAACCGt                                 <  1:2745793/58‑1 (MQ=255)
                                                    agagTATTTCCAGCTCAACCAGATTTGGCGCAATGATATCGc                                                >  1:3255957/1‑42 (MQ=255)
                                                    agagTATTTCCAGCTCAACCAGATTTGGCGCAATGATATCGCTGGCAGGCAAACCGTGCCGCACATGaaa                    >  1:3909117/1‑70 (MQ=255)
                                                    agagTATTTCCAGCTCAACCAGATTTGGCGCAATGATATCGCTGGCAGGCAAACCGTGCCGCACATGaaa                    >  1:2868878/1‑70 (MQ=255)
                                                    agagTATTTCCAGCTCAACCAGATTTGGCGCAATGATATCGCTGGCAGGCAAACCGTGCCGCACATGAAAc                   >  1:2111558/1‑71 (MQ=255)
                                                    agagTATTTCCAGCTCAACCAGATTTGGCGCAAAGata                                                    >  1:3465503/1‑38 (MQ=255)
                                                    agagTATTTCCAGCTAAACCAGATTTGGCGCAATGATATCGCTGGCAGGAAAACCGTGCCGCACATGaaa                    >  1:461354/1‑70 (MQ=255)
                                                     gagTATTTCCAGCTCAACCAGATTTGGCGCAATGATATCGCTGGCAGGCAAACCGTGCCGCACATGAAAct                  >  1:3170211/1‑71 (MQ=255)
                                                      agTATTTCCAGCTCAATCAGATTTGGCGCAATGat                                                     <  1:1018209/35‑1 (MQ=255)
                                                      agTATTTCCAGCTCAATCAGATTTGGCGCAATGat                                                     <  1:1635533/35‑1 (MQ=255)
                                                        tATTTCCAGCTCAATCAGATTTGGCGCAATGatat                                                   >  1:2047708/1‑35 (MQ=255)
                                                        tATTTCCAGCTCAATCAGATTTGGCGCAATGatat                                                   >  1:2998348/1‑35 (MQ=255)
                                                          tttCCAGCTCAACCAGATTTGGCGCAATGATATCGCTggc                                            >  1:1448550/1‑40 (MQ=255)
                                                          tttCCAGCTCAACCAGATTTGGCGCAATGATATCGCTGGCAGGCAAACCGTGCCGCACATGAAACTCTGc              <  1:4056438/70‑1 (MQ=255)
                                                          tttCCAGCTCAACCAGATTTGGCGCAATGATATCGCTGGCAGGCAAACCGTGCCGCACATGAAACTCTGCg             <  1:3655110/71‑1 (MQ=255)
                                                          tttCCAGCTCAACCAGATTTGGCGCAATGATATCGCTGGCAGGCAAACCGTGCCGCACATGAAACTCTGCg             <  1:3197445/71‑1 (MQ=255)
                                                            tCCAGCTCAACCAGATTTGGCGCAATGATATCGCTGGCAGGCAAACCGTGCCGCACATGAAACTCTGCGa            >  1:4219769/1‑70 (MQ=255)
                                                               aGCTCAACCAGATTTGGCGCAATGATATCGCTGGCAGGCAAACCGTGCCGCACATGAAACTCTGCGACAcc        <  1:4646084/71‑1 (MQ=255)
                                                                gCTCAACCAGATTTGGCGCAATGATATCGCTGGCAGGCAAACCGTGCCGCACATGAAACTCTGCGACAccc       >  1:4526693/1‑71 (MQ=255)
                                                                gCTCAACCAGATTTGGCGCAATGATATCGCTGGCAGGCAAACCGTGCCGCACATGAAACTCTGCGACAccc       <  1:2755693/71‑1 (MQ=255)
                                                                  tCAACCAGATTTGGCGCAATGATATCGCTGGCAGGCAAACCGTGCCGCACATGAAActct                <  1:2901950/60‑1 (MQ=255)
                                                                     aCCAGATTTGGCGCAATGATATCGCTGGCAGGCAAACCGTGCCGCACATGAAACTCTGCGACACCCGGTGc  <  1:2913817/71‑1 (MQ=255)
                                                                      ccAGATTTGGCGCAATGATATCGCTGGCAGGCAAACCGTGc                               >  1:466050/1‑41 (MQ=255)
                                                                      ccAGATTTGGCGCAATGATATCGCTGGCAGGCAAACCGTGc                               >  1:4519659/1‑41 (MQ=255)
                                                                      ccAGATTTGGCGCAATGATATCGCTGGCAGGCAAACCGTGc                               >  1:3542490/1‑41 (MQ=255)
                                                                      |                                                                     
AGTTCGCGCGCTGCCAGAACCGCTTCTTCGACGTTATTTACCGCATGCTCACAGAGTATTTCCAGCTCAACCAGATTTGGCGCAATGATATCGCTGGCAGGCAAACCGTGCCGCACATGAAACTCTGCGACACCCGGTGC  >  W3110S.gb/1717092‑1717231

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: