Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 53,012 A→G G135G (GGT→GGC)  pdxA ← 4‑hydroxy‑L‑threonine phosphate dehydrogenase, NAD‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091353,0120AG93.7% 197.1 / 2.5 63G135G (GGT→GGC) pdxA4‑hydroxy‑L‑threonine phosphate dehydrogenase, NAD‑dependent
Reads supporting (aligned to +/- strand):  ref base A (2/2);  new base G (22/37);  total (24/39)
Fisher's exact test for biased strand distribution p-value = 6.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

GTAAATGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCTGTGATCAGCGCGGCAAATTCGCCGTTCAGACAACCATCGCACGCT  >  NC_000913/52915‑53108
                                                                                                 |                                                                                                
gTAAATGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGGcc                                                                                                <  1:1133851/100‑1 (MQ=255)
  aaaTGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGGCCgg                                                                                              <  1:270460/100‑1 (MQ=255)
     tGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGGCCGGTaa                                                                                           >  1:108250/1‑100 (MQ=255)
       cgtcgtTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGGCCGGTaaaa                                                                                         >  1:1533552/1‑100 (MQ=255)
       cgtcgtTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGGCCGGTaaaa                                                                                         >  2:581434/1‑100 (MQ=255)
        gtcgtTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGGCCGGTAAAAg                                                                                        <  1:1203002/100‑1 (MQ=255)
        gtcgtTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGGCCGGTAAAAg                                                                                        <  1:1220800/100‑1 (MQ=255)
        gtcgtTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGGCCGGTAAAAg                                                                                        <  1:851021/100‑1 (MQ=255)
              gCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAATCTCGGTATGGCCGGTAAAAGGAATGc                                                                                  <  2:1427359/100‑1 (MQ=255)
              gCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGc                                                                                  >  1:242419/1‑100 (MQ=255)
              gCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTCCGAACGCTCTTCGAAAAACTCGGTATGGCCGGTAAAAGGAATGc                                                                                  >  2:1126579/1‑100 (MQ=255)
               ccAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGGCCGGTAAAAGGAATGcc                                                                                 <  2:108250/100‑1 (MQ=255)
                cAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGGCCGGTAAAAGGAATGCCa                                                                                >  1:730261/1‑100 (MQ=255)
                  gcgcCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGGCCGGTAAAAGGAATGCCAGc                                                                              <  1:201324/100‑1 (MQ=255)
                     ccACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGGCCGGTAAAAGGAATGCCAGcgtc                                                                           <  1:1426576/100‑1 (MQ=255)
                      cACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGGCCGGTAAAAGGAATGCCAGcgtcg                                                                          >  1:202065/1‑100 (MQ=255)
                        cgcgAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGGCCGGTAAAAGGAATGCCAGCGTCGtt                                                                        <  2:1226441/100‑1 (MQ=255)
                        cgcgAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGGCCGGTAAAAGGAATGCCAGCGTCGtt                                                                        <  2:38300/100‑1 (MQ=255)
                        cgcgAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGGCCGGTAAAAGGAATGCCAGCGTCGtt                                                                        <  1:1539964/100‑1 (MQ=255)
                        cgcgAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGGCCGGTAAAAGGAATGCCAGCGTCGtt                                                                        >  1:193532/1‑100 (MQ=255)
                           gAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGGCCGGTAAAAGGAATGCCAGCGTCGTtaat                                                                     <  1:397459/100‑1 (MQ=255)
                             aGTTCGTCGGGCGCCAGCATCATCACCACCTGTTTCGCCGGCGAACGCTCTTCGAAAAATTCGGTATGGCCGGTAAAAGGAATGCCAGCGTCGTtaataa                                                                   <  2:991340/100‑1 (MQ=255)
                              gTTCTGCGGTCGCGAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAATCTCGGTATGGCCGGTAAAAGGAATGCCAGCGTCGTTAATAAc                                                                  <  2:382875/100‑1 (MQ=255)
                               ttcttcGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGGCCGGTAAAAGGAATGCCAGCGTCGTTAATAACg                                                                 <  2:630746/100‑1 (MQ=255)
                                 cttcGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGGCCGGTAAAAGGAATGCCAGCGTCGTTAATAACGcc                                                               <  2:730259/100‑1 (MQ=255)
                                     gctCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGGCCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTa                                                           <  2:961521/98‑1 (MQ=255)
                                       tCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGGCCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATg                                                         <  1:1010711/100‑1 (MQ=255)
                                        cGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGGCCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGc                                                        <  2:790199/100‑1 (MQ=255)
                                         gCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGGCCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCa                                                       >  1:677592/1‑100 (MQ=255)
                                          ccAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGGCCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCAc                                                      <  2:52520/100‑1 (MQ=255)
                                          ccAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAAATCGGTATGGCCGGTAAAAGGCATGCCAGCGTCGTTAATAACGCCTTTATGCAc                                                      >  1:207987/1‑100 (MQ=255)
                                             gCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGGCCGGTAAAAGGAATGCCAGCGACGTTAATCACGCCTTTATGCAcc                                                     >  2:1536893/1‑98 (MQ=255)
                                             gCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCgg                                                   <  1:1386048/100‑1 (MQ=255)
                                                tcatcaCCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGGCCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGAcc                                                <  1:937785/100‑1 (MQ=255)
                                                 catcaCCACCTTTTTCGCCTGCGAACGCTCTGCGAAAAACTCGGTATGACCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGAcc                                                <  2:329547/99‑1 (MQ=255)
                                                   tcaCCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGGCCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCtgt                                             >  1:903503/1‑100 (MQ=255)
                                                    caccacCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGGCCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCtgtg                                            >  1:1361323/1‑100 (MQ=255)
                                                      ccaccTTTTTCTCTTTCGTTCGCTCTTCGAAAAACTCGGTATGGCCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCTGTGAt                                          >  2:191341/1‑100 (MQ=255)
                                                      ccaccTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGGCCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCTGTGAt                                          <  1:1286493/100‑1 (MQ=255)
                                                       caccTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGGCCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCTGTGATc                                         <  1:1025073/100‑1 (MQ=255)
                                                          cTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGGCCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCTGTGATCAgc                                      <  2:1161913/100‑1 (MQ=255)
                                                          cTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCTGTGATCAgc                                      >  1:1307964/1‑100 (MQ=255)
                                                               tCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGGCCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCTGTGATCAGCGCGGc                                 <  1:161111/100‑1 (MQ=255)
                                                                  ccTGCGAACGCTCTTCGAAAAACTCGGTATGGCCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCTGTGATCAGCGCGGCaaa                              >  1:131511/1‑100 (MQ=255)
                                                                      cGAACGCTCTTCGAAAAACTCGGTATGGCCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCTGTGATCAGCGCGGCAAATTCg                          <  1:1481518/100‑1 (MQ=255)
                                                                      cGAACGCTCTTCGAAAAACTCGGTATGGCCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCTGTGATCAGCGCGGCAAATTCg                          <  1:851736/100‑1 (MQ=255)
                                                                            ctctTCGAAAAACTCGGTATGGCCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCTGTGATCAGCGCGGCAAATTCGCCGTTc                    <  1:960476/100‑1 (MQ=255)
                                                                            ctctTCGAAAAACGCGGTATGGCCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCATCGGACCTGTGATCAGCGCGGCAAATTCGCCGTTc                    <  1:776106/100‑1 (MQ=255)
                                                                              ctTCGAAAAACTCGGTATGGCCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCTGTGATCAGCGCGGCAAATTCGCCGTTCAg                  <  2:1552055/100‑1 (MQ=255)
                                                                              ctTCGAAAAACTCGGTATGGCCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCTGTGATCAGCGCGGCAAATTCGCCGTTCAg                  <  2:1009018/100‑1 (MQ=255)
                                                                               ttCGAAAAACTCGGTATGGCCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCTGTGATCAGCGCGGCAAATTCGCCGTTCAGa                 <  1:281955/100‑1 (MQ=255)
                                                                                tCGAAAAACTCGGTATGGCCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCTGTGATCAGCGCGGCAAATTCGCCGTTCAGAc                >  2:984195/1‑100 (MQ=255)
                                                                                tCGAAAAACTCGGTATGGCCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCTGTGATCAGCGCGGCAAATTCGCCGTTCAGAc                <  1:1254954/100‑1 (MQ=255)
                                                                                 cGAAAAACTCGGTATGGCCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCTGTGATCAGCGCGGCAAATTCGCCGTTCAGACa               <  1:1126584/100‑1 (MQ=255)
                                                                                 cGAAAAACTCGGTATGGCCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCTGTGATCAGCGCGGCAAATTCGCCGTTCAGACa               >  2:1177992/1‑100 (MQ=255)
                                                                                  gAAAAACTCGGTATGGCCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCTGTGATCAGCGCGGCAAATTCGCCGTTCAGACaa              >  2:397061/1‑100 (MQ=255)
                                                                                  gAAAAACTCGGTATGGCCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCTGTGATCAGCGCGGCAAATTCGCCGTTCAGACaa              >  2:1089977/1‑100 (MQ=255)
                                                                                     aaaCTCGGTATGGCCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCTGTGATCAGCGCGGCAAATTCGCCGTTCAGACAACCa           >  1:1099534/1‑100 (MQ=255)
                                                                                        cTCGGTATGGCCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCTGTGATCAGCGCGGCAAATTCGCCGTTCAGACAACCATCg        <  2:344535/100‑1 (MQ=255)
                                                                                          cGGTATGGCCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCTGTGATCAGCGCGGCAAATTCGCCGTTCAGACAACCATCGCa      <  2:1578279/100‑1 (MQ=255)
                                                                                           ggTATGGCCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCTGTGATCAGCGCGGCAAATTCGCCGTTCAGACAACCATCGCAc     >  2:109570/1‑100 (MQ=255)
                                                                                            gTATGGCCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCTGTGATCAGCGCGGCAAATTCGCCGTTCAGACAACCATCGCACg    >  1:1160448/1‑100 (MQ=255)
                                                                                              aTGGCCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCTGTGATCAGCGCGGCAAATTCGCCGTTCAGACAACCATCGCACGCt  >  1:1357733/1‑100 (MQ=255)
                                                                                                 |                                                                                                
GTAAATGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCTGTGATCAGCGCGGCAAATTCGCCGTTCAGACAACCATCGCACGCT  >  NC_000913/52915‑53108

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