Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,008,075 A→G 100% D189G (GAC→GGC)  ygeX → 2,3‑diaminopropionate ammonia lyase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,008,0750AG98.7% 275.6 / ‑4.7 78D189G (GAC→GGC) ygeX2,3‑diaminopropionate ammonia lyase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (42/35);  total (43/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.45e-01

TCGGTGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGACACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAATGG  >  NC_000913/3007979‑3008170
                                                                                                |                                                                                               
tCGGTGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGgcac                                                                                              >  2:2195985/1‑100 (MQ=255)
tCGGTGCCGAGTGCATCGGCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGgcag                                                                                              <  2:717420/100‑2 (MQ=255)
      ccGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATg                                                                                        <  1:250745/100‑1 (MQ=255)
       cGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCGGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATgg                                                                                       <  2:1258913/100‑1 (MQ=255)
            gCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAgg                                                                                  >  1:1739673/1‑100 (MQ=255)
            gCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAgg                                                                                  <  2:2062272/100‑1 (MQ=255)
             cATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGt                                                                                 <  1:1501335/100‑1 (MQ=255)
              aTCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGtt                                                                                >  1:44309/1‑100 (MQ=255)
              aTCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGtt                                                                                >  2:821915/1‑100 (MQ=255)
                cGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTac                                                                              <  2:355509/100‑1 (MQ=255)
                 gTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTaca                                                                             >  2:559414/1‑100 (MQ=255)
                   cACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGCAGTGGTACAGGGCACGGACTGGGAAGGTTACAcc                                                                           >  2:467363/1‑100 (MQ=255)
                    aCGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCa                                                                          >  2:240858/1‑100 (MQ=255)
                     cGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACAGCTGGGAAGTGGTCAAGGGCACGGCATGGGAAGGTTACATCa                                                                          >  2:833368/1‑99 (MQ=255)
                       gATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCaaca                                                                       >  1:1217423/1‑98 (MQ=255)
                       gATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCaaaa                                                                       <  2:1944384/100‑1 (MQ=255)
                       gATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCaaaa                                                                       >  1:993855/1‑100 (MQ=255)
                         tatGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATc                                                                     >  2:1816930/1‑100 (MQ=255)
                           tGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATccc                                                                   >  2:83652/1‑100 (MQ=255)
                           tGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATccc                                                                   >  1:965956/1‑100 (MQ=255)
                           tGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATccc                                                                   >  2:993228/1‑100 (MQ=255)
                           tGAACTATCAAGATACCGTTCGCCAGACCATGCAACACGCGCACCAGCACGGATGGGAAGTAGTACAGGGCACTGCATGGGAAGGTTACACCAAAATccc                                                                   <  2:311065/100‑1 (MQ=255)
                            gAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCa                                                                  <  1:2237432/100‑1 (MQ=255)
                            gAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCa                                                                  >  1:2170470/1‑100 (MQ=255)
                             aaCTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCaa                                                                 <  2:1067920/100‑1 (MQ=255)
                             aaCTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCaa                                                                 <  2:635438/100‑1 (MQ=255)
                              aCTAGGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAAc                                                                <  2:996186/100‑1 (MQ=255)
                                tATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCt                                                              <  1:723745/100‑1 (MQ=255)
                                 aTGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTg                                                             <  2:779802/100‑1 (MQ=255)
                                  tGACGATACCGTTCGCCTGATCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTgg                                                            >  1:2214395/1‑100 (MQ=255)
                                  tGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTgg                                                            <  2:616021/100‑1 (MQ=255)
                                   ggcGATACCGTTCGCCTGACCATGCAACACGGGCAGCAGCACAGCTGGGAAGTGGTACAGGACACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGa                                                           >  1:1588366/3‑100 (MQ=255)
                                    aCGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGAt                                                          >  2:1758593/1‑100 (MQ=255)
                                    aCGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGAt                                                          <  2:1954226/100‑1 (MQ=255)
                                     cGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATc                                                         >  1:2206/1‑100 (MQ=255)
                                          ccGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGCTCATGCa                                                    >  2:1387184/1‑100 (MQ=255)
                                            gTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAg                                                  >  1:295643/1‑100 (MQ=255)
                                            gTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAg                                                  <  2:1745190/100‑1 (MQ=255)
                                            gTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAg                                                  <  2:937428/100‑1 (MQ=255)
                                             ttCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAgg                                                 <  2:857261/100‑1 (MQ=255)
                                              tCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGc                                                <  1:1401971/100‑1 (MQ=255)
                                              tCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGc                                                >  1:1008581/1‑100 (MQ=255)
                                                  cTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACg                                            <  1:84403/100‑1 (MQ=255)
                                                   tGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGc                                           <  2:2021015/100‑1 (MQ=255)
                                                   tGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGc                                           <  2:413705/100‑1 (MQ=255)
                                                    gACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTGCACCAAAATCCCAACCTGGATCATGCAAGGCTACGCa                                          <  2:1252371/100‑1 (MQ=255)
                                                    gACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCa                                          >  2:142291/1‑100 (MQ=255)
                                                    gACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCa                                          >  1:2347477/1‑100 (MQ=255)
                                                       cATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAAcc                                       >  2:2068720/1‑100 (MQ=255)
                                                            aaCACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTgg                                   >  2:1636690/1‑99 (MQ=255)
                                                            aaCACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGc                                  >  1:34247/1‑100 (MQ=255)
                                                            aaCACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGc                                  <  1:1857800/100‑1 (MQ=255)
                                                              cacGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAg                                >  2:1628980/1‑100 (MQ=255)
                                                              cacGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAg                                <  2:2347455/100‑1 (MQ=255)
                                                               acGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGa                               <  2:775423/100‑1 (MQ=255)
                                                                 gcgcAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATg                             <  2:1556545/100‑1 (MQ=255)
                                                                 gcgcAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATg                             <  1:2108441/100‑1 (MQ=255)
                                                                     agcagcaCGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGc                         <  2:44309/100‑1 (MQ=255)
                                                                      gcagcaCGGCTGGGAGGTGGTACAGGGCACGGCATGCGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGc                         <  2:518085/99‑1 (MQ=255)
                                                                      gcagcaCGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAATGCTACGCAACCCTGGCAGATGAAGcc                        <  1:1431781/100‑1 (MQ=255)
                                                                      gaagcaCGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGcc                        >  1:1595013/3‑100 (MQ=255)
                                                                        agcaCGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCg                       >  2:585482/1‑99 (MQ=255)
                                                                          caCGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCTTCg                    >  2:1054738/1‑100 (MQ=255)
                                                                            cGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAg                  <  2:1784723/100‑1 (MQ=255)
                                                                            cGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAg                  >  2:2082781/1‑100 (MQ=255)
                                                                            cGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTAAGCAACCCTGGCAGATGAAGCCGTAGa                   >  2:1433160/1‑99 (MQ=255)
                                                                             ggCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGc                 <  2:34057/100‑1 (MQ=255)
                                                                             ggCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGc                 <  1:1922397/100‑1 (MQ=255)
                                                                               cTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCaa               >  1:930191/1‑100 (MQ=255)
                                                                                tGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCaaa              >  1:1109068/1‑100 (MQ=255)
                                                                                  ggAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATg            >  2:1617153/1‑100 (MQ=255)
                                                                                    aaGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCg          <  2:1124454/100‑1 (MQ=255)
                                                                                     aGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGCAGCCGTCGAGCATATGCGt         >  2:653903/1‑100 (MQ=255)
                                                                                        ggTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGaa      >  2:1788502/1‑100 (MQ=255)
                                                                                        ggTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGaa      >  1:2327707/1‑100 (MQ=255)
                                                                                         gTACAGGGCCCGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGaaa     >  2:1801499/1‑100 (MQ=255)
                                                                                         gTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGaaa     >  1:2277689/1‑100 (MQ=255)
                                                                                          tACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAAt    >  1:1547102/1‑100 (MQ=255)
                                                                                           aCAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAATg   >  1:1142392/1‑100 (MQ=255)
                                                                                            cAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAATgg  >  1:656312/1‑100 (MQ=255)
                                                                                                |                                                                                               
TCGGTGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGACACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAATGG  >  NC_000913/3007979‑3008170

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