Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 354,592 G→A *629* (TAG→TAA prpE → propionate‑‑CoA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913354,5920GA96.9% 198.8 / ‑6.7 64*629* (TAG→TAAprpEpropionate‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (34/28);  minor base C (0/2);  total (34/30)
Fisher's exact test for biased strand distribution p-value = 2.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGGAAAAATGCTGCGCCGCACGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACC  >  NC_000913/354474‑354723
                                                                                                                      |                                                                                                                                   
ccGGAAAAATGCTGCGCCGCACGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTcgc                                                                                                                 <  2:129756/139‑1 (MQ=255)
            tgcGCCGCACGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCg                                                                                                     <  2:20893/139‑1 (MQ=255)
                 cgcACGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTACGTTATTGTAGGATGCGACGCGCGGTGCATCCGGCACt                                                                                                <  1:268682/139‑1 (MQ=255)
                                  ttGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGGGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCa                                                                                                  >  2:883038/1‑120 (MQ=255)
                                  ttGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTcgcgcg                                                                                                              <  1:790526/108‑1 (MQ=255)
                                  ttGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTcgcgcg                                                                                                              >  2:790526/1‑108 (MQ=255)
                                  ttGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCa                                                                                                  <  1:1009539/120‑1 (MQ=255)
                                  ttGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCa                                                                                                  >  2:1009539/1‑120 (MQ=255)
                                  ttGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCa                                                                                                  <  1:883038/120‑1 (MQ=255)
                                  ttGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGAt                                                                                      <  1:253025/132‑1 (MQ=255)
                                  ttGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGAt                                                                                      >  2:253025/1‑132 (MQ=255)
                                  ttGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGAt                                                                               >  1:924396/1‑139 (MQ=255)
                                  ttGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGAt                                                                               >  2:579726/1‑139 (MQ=255)
                                  ttGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGAt                                                                               >  2:437692/1‑139 (MQ=255)
                                  ttGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGAt                                                                               >  2:1052646/1‑139 (MQ=255)
                                  ttGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGAt                                                                               >  1:811456/1‑139 (MQ=255)
                                                    tGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGtttt                                                             >  2:410130/1‑139 (MQ=255)
                                                         aTCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCAt                                                        <  1:579726/139‑1 (MQ=255)
                                                           cTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGt                                                                                                            >  1:1204897/1‑85 (MQ=255)
                                                           cTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGt                                                                                                            <  2:1204897/85‑1 (MQ=255)
                                                           cTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGt                                                                >  1:917468/1‑129 (MQ=255)
                                                           cTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGt                                                                <  2:917468/129‑1 (MQ=255)
                                                           cTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATc                                                          <  1:809973/135‑1 (MQ=255)
                                                           cTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATc                                                          >  2:809973/1‑135 (MQ=255)
                                                           cTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGc                                                      >  1:15320/1‑139 (MQ=255)
                                                           cTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGc                                                      >  1:416469/1‑139 (MQ=255)
                                                           cTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGc                                                      >  2:881284/1‑139 (MQ=255)
                                                           cTGACGACCATTGATGATCCTGCGTCGTTGGAGCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTATGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGc                                                      >  1:725365/1‑139 (MQ=255)
                                                              acgacCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTATGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTa                                                   <  2:725365/139‑1 (MQ=255)
                                                                     ttGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTACGTTATTGTCGGATGCGTCGCGCGGTGCATCCg                                                                                                     <  2:236690/82‑1 (MQ=255)
                                                                     ttGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCg                                                                                                     >  1:236690/1‑82 (MQ=255)
                                                                     ttGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTa                                                   >  1:905374/1‑132 (MQ=255)
                                                                     ttGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTa                                                   <  2:905374/132‑1 (MQ=255)
                                                                       gatgatCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGAc                                                                    >  1:943911/1‑113 (MQ=255)
                                                                       gatgatCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGAc                                                                    <  2:943911/113‑1 (MQ=255)
                                                                            tCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGt                                     <  1:437692/139‑1 (MQ=255)
                                                                               tGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAgg                                  <  1:184330/139‑1 (MQ=255)
                                                                                gCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGt                                 >  1:819257/1‑139 (MQ=255)
                                                                                 cgtcgtTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTc                                >  2:688169/1‑139 (MQ=255)
                                                                                 cgtcgtTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGAAGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTc                                >  1:338866/1‑139 (MQ=255)
                                                                                 cgtcgtTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGAGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTc                                >  1:871331/1‑139 (MQ=255)
                                                                                      ttGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATa                           >  2:981101/1‑139 (MQ=255)
                                                                                      ttGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATa                           >  1:1043680/1‑139 (MQ=255)
                                                                                      ttGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATa                           >  2:1045613/1‑139 (MQ=255)
                                                                                           tCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGa                             >  2:1150150/1‑132 (MQ=255)
                                                                                           tCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGa                             <  1:1150150/132‑1 (MQ=255)
                                                                                              gATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTc                   >  2:1126067/1‑139 (MQ=255)
                                                                                                tCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTcgc                 >  2:1199058/1‑139 (MQ=255)
                                                                                                 ccGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTATGTCGGATAAGGCGCTcgcg                <  2:811456/139‑1 (MQ=255)
                                                                                                   gCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTc              <  1:981101/139‑1 (MQ=255)
                                                                                                    ccAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTa                                                   >  1:1014990/1‑101 (MQ=255)
                                                                                                    ccAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTa                                                   <  2:1014990/101‑1 (MQ=255)
                                                                                                     cAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGc            <  1:410130/139‑1 (MQ=255)
                                                                                                        gCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGc                                                                                                          <  1:858252/42‑1 (MQ=255)
                                                                                                        gCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGc                                                                                                          >  2:858252/1‑42 (MQ=255)
                                                                                                         cGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGt                                                                                                            >  1:451013/1‑39 (MQ=255)
                                                                                                         cGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGt                                                                                                            <  2:451013/39‑1 (MQ=255)
                                                                                                           aTGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGa      >  2:473777/1‑139 (MQ=255)
                                                                                                             ggAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGaca    <  1:688169/139‑1 (MQ=255)
                                                                                                             ggAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGaca    <  1:1199058/139‑1 (MQ=255)
                                                                                                             ggAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGaca    <  1:1126067/139‑1 (MQ=255)
                                                                                                              gAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGacac   <  2:924396/139‑1 (MQ=255)
                                                                                                               aaGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACAcc  <  2:15320/139‑1 (MQ=255)
                                                                                                               aaGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGAAGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACAcc  <  2:338866/139‑1 (MQ=255)
                                                                                                                      |                                                                                                                                   
CCGGAAAAATGCTGCGCCGCACGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACC  >  NC_000913/354474‑354723

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: