Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,008,075 A→G 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.0% 170.1 / ‑4.1 53D189G (GAC→GGC) ygeX2,3‑diaminopropionate ammonia lyase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (29/21);  total (30/23)
Fisher's exact test for biased strand distribution p-value = 4.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.84e-01

GCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGACACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAATG  >  NC_000913/3007984‑3008169
                                                                                           |                                                                                              
gCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAACAGCACGGCTGGGAAGTGGTACAGGGCACGGCAt                                                                                        >  2:1708/1‑100 (MQ=255)
 ccGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATg                                                                                       <  1:223271/100‑1 (MQ=255)
  cGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATgg                                                                                      >  1:1505131/1‑100 (MQ=255)
     gTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGaa                                                                                   <  1:455603/100‑1 (MQ=255)
                  gATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATAGGAAGGTTACACCaaaa                                                                      >  2:992916/1‑100 (MQ=255)
                    tatGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATc                                                                    >  1:111093/1‑100 (MQ=255)
                     atGAACTATTACGATACTGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATcc                                                                   <  2:761032/100‑1 (MQ=255)
                      tGAACTATGACGATACCGTTAGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATccc                                                                  >  1:1111319/1‑100 (MQ=255)
                       gAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCa                                                                 <  1:1429534/100‑1 (MQ=255)
                         aCTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCaa                                                                >  2:1174906/1‑99 (MQ=255)
                          cTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGACTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAAACCCAAcc                                                              >  2:875841/1‑100 (MQ=255)
                           tATGACGATACCGTTCGCCTTACCATGCAACTCGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGGACGGCATGGGAAGGTTACACCAAAATCCCAAc                                                               <  2:99192/98‑1 (MQ=255)
                           tATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCt                                                             >  1:1389332/1‑100 (MQ=255)
                             tGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTgg                                                           >  2:384602/1‑100 (MQ=255)
                             tGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGCCATGGGAAGGTTACACCAAAATCCCAACCTgg                                                           >  1:578941/1‑100 (MQ=255)
                             tGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGACACGGCATGGGAAGGTTACACCAAAATCCCAACCTgg                                                           <  2:1152080/100‑1 (MQ=255)
                              gACAATAACGTTCGACTGACCATGCAACACGCGCAGCAGCACGGCAGGACAGAGGTAAAGGCCGCGGCATGGGAAGGTTACACCAAAATCCCAAc                                                               >  1:1081848/1‑95 (MQ=255)
                                cGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATc                                                        <  2:1603201/100‑1 (MQ=255)
                                cGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATc                                                        >  2:1530787/1‑100 (MQ=255)
                                 gATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCa                                                       <  1:1463563/100‑1 (MQ=255)
                                     ccGTTCGCCTGACCATGCAACACGCGCAGCAGCCCGGCTGGGAAGTGGTACAGGGCACGGCATGGGAATATAACACAAAAATCCAAACATGGATCATTCa                                                   >  2:573645/1‑100 (MQ=255)
                                     ccGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCa                                                   >  2:893601/1‑100 (MQ=255)
                                             cTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACg                                           <  2:1389342/100‑1 (MQ=255)
                                              tGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGc                                          <  1:424465/100‑1 (MQ=255)
                                              tGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGc                                          <  2:1137516/100‑1 (MQ=255)
                                                   aTGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAAccc                                     <  1:64782/100‑1 (MQ=255)
                                                    tGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCt                                    >  2:1474824/1‑100 (MQ=255)
                                                         cacGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAg                               >  2:34012/1‑100 (MQ=255)
                                                            gcgcAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATg                            <  1:1104295/100‑1 (MQ=255)
                                                               cagctgcaCGGATGGGAAGTGGTACAGGTCACGGTATGGGAAGGTTACACCAAAATCTCAACCTTGTTCATGCAAGGCTACGCAACCCTGGCAGATGAAg                         <  1:385466/100‑1 (MQ=255)
                                                                 gcagcaCGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGcc                       >  2:112544/1‑100 (MQ=255)
                                                                   agcaCGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGt                     <  1:1621521/100‑1 (MQ=255)
                                                                    gcaCGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCg                      >  2:225471/1‑98 (MQ=255)
                                                                    gcaCGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTc                    <  2:651977/100‑1 (MQ=255)
                                                                       cGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAg                 >  1:1185556/1‑100 (MQ=255)
                                                                        ggaggggAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCGAACCTGGAGCATGCAAGGCTACGCAACCCTGGCAGATGAATCCGTCGa                  <  1:77393/94‑1 (MQ=255)
                                                                        ggCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGc                <  1:1737326/100‑1 (MQ=255)
                                                                        ggCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGc                >  1:1816858/1‑100 (MQ=255)
                                                                         gCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCa               >  2:290428/1‑100 (MQ=255)
                                                                         gCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCa               <  1:1010031/100‑1 (MQ=255)
                                                                            gggAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCaaa             >  2:1586943/1‑99 (MQ=255)
                                                                            gggAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAAt            >  2:1439099/1‑100 (MQ=255)
                                                                            gggAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAAt            <  1:1708/100‑1 (MQ=255)
                                                                            gggAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAAt            >  2:935088/1‑100 (MQ=255)
                                                                              gAAGTGGTACAGGGCACGGCATAGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATg           >  1:719849/1‑99 (MQ=255)
                                                                               aaGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCg         >  1:502163/1‑100 (MQ=255)
                                                                                 gTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGTAAATGCGTg       <  1:1505721/100‑1 (MQ=255)
                                                                                 gTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTg       <  2:73370/100‑1 (MQ=255)
                                                                                   ggTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGaa     >  1:1368938/1‑100 (MQ=255)
                                                                                   ggTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGaa     >  1:1359715/1‑100 (MQ=255)
                                                                                     tACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGAGAAAt   >  2:352575/1‑100 (MQ=255)
                                                                                      aCAGGGCACGGCATGGGAAGGTTAAACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAATg  <  1:1711086/100‑1 (MQ=255)
                                                                                        aGGGCACGGCATGGGAAGGTTAAACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGaa     >  1:6739/1‑95 (MQ=255)
                                                                                           |                                                                                              
GCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGACACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAATG  >  NC_000913/3007984‑3008169

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