Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,490,088 A→G G31G (GGT→GGC argD ← bifunctional acetylornithine aminotransferase and succinyldiaminopimelate aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,490,0880AG97.3% 257.0 / ‑2.5 74G31G (GGT→GGCargDbifunctional acetylornithine aminotransferase and succinyldiaminopimelate aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (36/36);  total (37/37)
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.97e-01

AGCGCAGGATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTC  >  NC_000913/3489993‑3490183
                                                                                               |                                                                                               
aGCGCAGGATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCtt                                                                                             <  1:1156487/100‑1 (MQ=255)
aGCGCAGGATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCtt                                                                                             <  1:1156603/100‑1 (MQ=255)
aGCGCAGGATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCtt                                                                                             <  1:199340/100‑1 (MQ=255)
   gcAGGATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTa                                                                                          >  1:931328/1‑100 (MQ=255)
    cAGGATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTAc                                                                                         >  2:1879758/1‑100 (MQ=255)
    cAGGATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTAc                                                                                         <  1:1767381/100‑1 (MQ=255)
    cAGGATGGCAATGGCCCAACGCCATAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTAc                                                                                         >  1:1982028/1‑100 (MQ=255)
     aGGATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTAcc                                                                                        <  2:2165666/100‑1 (MQ=255)
      ggATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCg                                                                                       >  1:1062564/1‑100 (MQ=255)
      ggATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCg                                                                                       >  2:2330129/1‑100 (MQ=255)
       gATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCgg                                                                                      <  1:2056155/100‑1 (MQ=255)
       gATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCgg                                                                                      <  2:1062562/100‑1 (MQ=255)
       gATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCgg                                                                                      <  1:1602758/100‑1 (MQ=255)
        aTGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGa                                                                                     >  1:574531/1‑100 (MQ=255)
          ggCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAAt                                                                                   <  1:1448505/100‑1 (MQ=255)
          ggCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAAt                                                                                   >  2:2277021/1‑100 (MQ=255)
          ggCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAAt                                                                                   <  2:316667/100‑1 (MQ=255)
           gCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATa                                                                                  >  2:1178973/1‑100 (MQ=255)
            cAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATaa                                                                                 >  2:1282859/1‑100 (MQ=255)
              aTGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAAc                                                                               <  2:2101521/100‑1 (MQ=255)
                ngCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAActc                                                                             >  1:1124299/2‑100 (MQ=255)
                ggCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAActc                                                                             >  1:594637/1‑100 (MQ=255)
                ggCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAActc                                                                             >  1:1632013/1‑100 (MQ=255)
                     aaCGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCg                                                                        >  1:703470/1‑100 (MQ=255)
                     aaCGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCg                                                                        <  1:1179838/100‑1 (MQ=255)
                     aaCGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCg                                                                        <  2:1234776/100‑1 (MQ=255)
                        gCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAg                                                                     <  2:594636/100‑1 (MQ=255)
                         ccGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGc                                                                    <  1:846106/100‑1 (MQ=255)
                          cGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCa                                                                   >  2:1244960/1‑100 (MQ=255)
                              aCTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATaaa                                                               <  1:1222748/100‑1 (MQ=255)
                                tGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATc                                                             <  1:1178976/100‑1 (MQ=255)
                                tGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATc                                                             >  1:414027/1‑100 (MQ=255)
                                 gCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCg                                                            >  1:1768004/1‑100 (MQ=255)
                                   aaTGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGc                                                          <  1:559218/100‑1 (MQ=255)
                                    aTGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCa                                                         <  1:1288971/100‑1 (MQ=255)
                                      gCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAgg                                                       >  1:402648/1‑100 (MQ=255)
                                        cACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGAt                                                     >  1:960012/1‑100 (MQ=255)
                                            cgcgAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACt                                                 <  1:1959652/100‑1 (MQ=255)
                                            cgcgAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACt                                                 >  2:693861/1‑100 (MQ=255)
                                                     gACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGaaa                                        <  2:531655/100‑1 (MQ=255)
                                                     gACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGaaa                                        >  2:1846206/1‑100 (MQ=255)
                                                     gACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGaaa                                        >  1:1491057/1‑100 (MQ=255)
                                                     gACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGaaa                                        >  2:1765988/1‑100 (MQ=255)
                                                     gACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGaaa                                        <  1:1201022/100‑1 (MQ=255)
                                                      tcttACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAg                                       <  2:353189/97‑1 (MQ=255)
                                                        aTACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTc                                     >  1:2238517/1‑100 (MQ=255)
                                                         tACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAgaga                                    >  1:353191/1‑97 (MQ=255)
                                                            tCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTcgcgc                                 <  1:2165629/100‑1 (MQ=255)
                                                            tCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGAATCACTTCATCGAAAGTcgcgc                                 <  2:1632008/100‑1 (MQ=255)
                                                                tgccttgcTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCgtgt                             >  1:1481160/1‑100 (MQ=255)
                                                                   cttgcTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAAt                          <  1:1567881/100‑1 (MQ=255)
                                                                   cttgcTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAAt                          >  1:1371177/1‑100 (MQ=255)
                                                                   cttccgatctccAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAAt                          <  2:1755031/90‑1 (MQ=255)
                                                                     tgcTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTg                        <  2:931326/100‑1 (MQ=255)
                                                                     tgcTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTg                        >  1:2098180/1‑100 (MQ=255)
                                                                     tgcTGATCCCAGATTCAGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTg                        >  1:2012439/1‑100 (MQ=255)
                                                                      gcTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGc                       >  1:1989937/1‑100 (MQ=255)
                                                                       cTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCt                      <  2:277409/100‑1 (MQ=255)
                                                                         gATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCtgt                    >  1:1095156/1‑100 (MQ=255)
                                                                         gATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCtgt                    >  1:1249116/1‑100 (MQ=255)
                                                                          aTCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCtgtt                   >  2:2217054/1‑100 (MQ=255)
                                                                          aTCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCtgtt                   <  2:414025/100‑1 (MQ=255)
                                                                            cccAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCTGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCtgtttg                 <  2:1491053/100‑1 (MQ=255)
                                                                            cccAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCtgtttg                 >  2:1381486/1‑100 (MQ=255)
                                                                             ccAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAAtagatcggaag                >  1:1755036/1‑90 (MQ=255)
                                                                             ccAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCtgtttgt                <  1:748526/100‑1 (MQ=255)
                                                                             ccAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCtgtttgt                >  1:1723324/1‑100 (MQ=255)
                                                                             ccAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCtgtttgt                >  1:1681984/1‑100 (MQ=255)
                                                                                    cGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTg         <  2:1590016/100‑1 (MQ=255)
                                                                                     ggCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAAttn         <  2:1124297/99‑2 (MQ=255)
                                                                                       cTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCa      <  2:1766936/100‑1 (MQ=255)
                                                                                        tGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCAt     <  2:2083567/100‑1 (MQ=255)
                                                                                         gCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCAtc    >  1:2252938/1‑100 (MQ=255)
                                                                                           ccTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCAtctc  <  2:703469/100‑1 (MQ=255)
                                                                                               |                                                                                               
AGCGCAGGATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTC  >  NC_000913/3489993‑3490183

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