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,0880AG96.0% 255.8 / ‑0.5 77G31G (GGT→GGCargDbifunctional acetylornithine aminotransferase and succinyldiaminopimelate aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base G (36/37);  total (38/39)
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.73e-01

CCAGCGCAGGATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCA  >  NC_000913/3489991‑3490184
                                                                                                 |                                                                                                
ccAGCGCAGGATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGcc                                                                                                >  2:490908/1‑100 (MQ=255)
  aGCGCAGGATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCtt                                                                                              >  1:1712759/1‑100 (MQ=255)
  aGCGCAGGATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCtt                                                                                              <  2:435075/100‑1 (MQ=255)
   gcgcAGGATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCttt                                                                                             >  1:1534197/1‑100 (MQ=255)
   gcgcAGGATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCttt                                                                                             <  1:1240230/100‑1 (MQ=255)
   gcgcAGGATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCttt                                                                                             <  2:1284861/100‑1 (MQ=255)
    cgcAGGATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCtttt                                                                                            >  2:1175530/1‑100 (MQ=255)
      cAGGATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTAc                                                                                          <  1:1283008/100‑1 (MQ=255)
      cAGGATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTAc                                                                                          >  2:135750/1‑100 (MQ=255)
        gcaTGGCACTGGCACACCGCCGTAACTGCAATGCCTCCCGCGAAATCGACATACACCTTGCCTAGCTGATCCCAGATGCGGCTGCCCTGTCCTTTTACCg                                                                                        >  2:1908442/3‑100 (MQ=255)
         gATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCgg                                                                                       <  1:212857/100‑1 (MQ=255)
         gATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCgg                                                                                       <  1:1911026/100‑1 (MQ=255)
         gATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCgg                                                                                       <  1:845280/100‑1 (MQ=255)
         gATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCgg                                                                                       >  1:553580/1‑100 (MQ=255)
           tGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGaa                                                                                     <  1:1491786/100‑1 (MQ=255)
           tGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGaa                                                                                     <  1:122131/100‑1 (MQ=255)
            ggCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAAt                                                                                    >  1:1460732/1‑100 (MQ=255)
             gCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATa                                                                                   >  2:1681393/1‑100 (MQ=255)
                aTGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAAc                                                                                <  2:1103826/100‑1 (MQ=255)
                aTGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAAc                                                                                <  1:1897282/100‑1 (MQ=255)
                   gCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAActct                                                                             <  2:277474/100‑1 (MQ=255)
                    cccAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAActctn                                                                            <  2:1414354/100‑2 (MQ=255)
                    cccAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTg                                                                            >  2:1530025/1‑100 (MQ=255)
                    cccAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTg                                                                            >  2:102448/1‑100 (MQ=255)
                          gCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAg                                                                      <  1:961826/100‑1 (MQ=255)
                            cGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCa                                                                    <  1:1006319/100‑1 (MQ=255)
                              tAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATa                                                                  <  2:1794804/100‑1 (MQ=255)
                                 cTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAAt                                                               >  1:72761/1‑100 (MQ=255)
                                  tGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATc                                                              <  2:1174011/100‑1 (MQ=255)
                                    cAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCgg                                                            <  1:191769/100‑1 (MQ=255)
                                    cAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCgg                                                            <  1:1198036/100‑1 (MQ=255)
                                      aTGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCa                                                          <  1:314841/100‑1 (MQ=255)
                                        gCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAgg                                                        >  2:1566126/1‑100 (MQ=255)
                                        gCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAgg                                                        <  2:72127/100‑1 (MQ=255)
                                        gCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTGTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAgg                                                        <  2:1255086/100‑1 (MQ=255)
                                         ccACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGa                                                       >  2:452582/1‑100 (MQ=255)
                                           aCCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATc                                                     <  1:1356642/100‑1 (MQ=255)
                                              cgcgAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACt                                                  >  1:1755552/1‑100 (MQ=255)
                                              cgcgAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACt                                                  >  2:1284696/1‑100 (MQ=255)
                                               gcgAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACtt                                                 <  1:1369295/100‑1 (MQ=255)
                                                cgAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTc                                                <  2:1755538/100‑1 (MQ=255)
                                                   aaTCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATc                                             <  1:102449/100‑1 (MQ=255)
                                                       gACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGaaa                                         <  1:336332/100‑1 (MQ=255)
                                                         cATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGt                                       <  1:1530034/100‑1 (MQ=255)
                                                           tACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTcg                                     <  1:1566136/100‑1 (MQ=255)
                                                           tACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTcg                                     >  2:216161/1‑100 (MQ=255)
                                                           tACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTcg                                     >  2:140155/1‑100 (MQ=255)
                                                           tACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCCTGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTACATCGAAAGTc                                      >  2:1659106/1‑99 (MQ=255)
                                                            aCTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTcgc                                    >  2:626376/1‑100 (MQ=255)
                                                            aCTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTcgc                                    <  2:1534188/100‑1 (MQ=255)
                                                               ccttgccttgCTGGTCCCAGATTCGGCTGCCCTGGCCTATTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTcgcgc                                  >  2:1259622/1‑99 (MQ=255)
                                                               ccttgccttgCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGCACAAACTCTGCCGGAGCATCAATCGGCAGGATCACTTCATCGAAAGTcgcgc                                  >  1:1709130/1‑99 (MQ=255)
                                                               ccttgccttgCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTcgcgcg                                 >  2:971840/1‑100 (MQ=255)
                                                               ccttgccttgCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTcgcgcg                                 >  2:1847725/1‑100 (MQ=255)
                                                               ccttgccttgCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTcgcgcg                                 >  1:993427/1‑100 (MQ=255)
                                                               ccttgccttgCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTcgcgcg                                 >  1:514286/1‑100 (MQ=255)
                                                                cttgccttgcTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCgt                                >  2:990866/1‑100 (MQ=255)
                                                                cttgccttgcTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCgt                                <  1:1681406/100‑1 (MQ=255)
                                                                  tgccttgcTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCgtgt                              >  2:592192/1‑100 (MQ=255)
                                                                     cttgcTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAAt                           <  2:277447/100‑1 (MQ=255)
                                                                     cttgcTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAAt                           <  2:1111641/100‑1 (MQ=255)
                                                                      ttgcTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAAtt                          <  1:1175536/100‑1 (MQ=255)
                                                                       tgcTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTg                         >  1:1750905/1‑100 (MQ=255)
                                                                       tgcTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTg                         >  2:413614/1‑100 (MQ=255)
                                                                           gATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCtgt                     >  1:1508102/1‑100 (MQ=255)
                                                                            aTCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCtgtt                    <  1:971844/100‑1 (MQ=255)
                                                                            aTCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCtgtt                    >  2:785894/1‑100 (MQ=255)
                                                                              cccAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCtgtttg                  <  1:140156/100‑1 (MQ=255)
                                                                               ccAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCtgtttgt                 >  2:1532961/1‑100 (MQ=255)
                                                                                 aGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTc               <  2:514285/100‑1 (MQ=255)
                                                                                 aGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTc               <  1:135751/100‑1 (MQ=255)
                                                                                       ggCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGc         >  1:1681711/1‑100 (MQ=255)
                                                                                         cTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCa       <  2:1810912/100‑1 (MQ=255)
                                                                                           gCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCAtc     >  1:642647/1‑100 (MQ=255)
                                                                                           gCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCAtc     >  2:35332/1‑100 (MQ=255)
                                                                                           gCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCAtc     >  1:1670290/1‑100 (MQ=255)
                                                                                              cTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCa  >  1:87162/1‑100 (MQ=255)
                                                                                                 |                                                                                                
CCAGCGCAGGATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCA  >  NC_000913/3489991‑3490184

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