Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 125,387 C→T Q791* (CAG→TAG)  aceE → pyruvate dehydrogenase E1 component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913125,3870CT95.9% 131.5 / NA 49Q791* (CAG→TAG) aceEpyruvate dehydrogenase E1 component
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (25/22);  minor base G (0/2);  total (25/24)
Fisher's exact test for biased strand distribution p-value = 2.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGAGAA  >  NC_000913/125245‑125526
                                                                                                                                              |                                                                                                                                           
tGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGGAACGCCGCGCGTTCCGTATAGCGCTc                                                                                                                                             >  2:811622/1‑143 (MQ=255)
     tCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGtg                                                                                                                                        <  1:812775/143‑1 (MQ=255)
       cTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGAAGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGtgat                                                                                                                                      >  1:112631/1‑143 (MQ=255)
        tGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGtgatg                                                                                                                                     >  1:686206/1‑143 (MQ=255)
         ggCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGAAACATGCTGCACCCGCTGGACACTCCGCGCGTTCGGTATATCGCTTCGGtgataa                                                                                                                                    >  1:313401/1‑141 (MQ=255)
          gCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGaa                                                                                                                                   >  1:643663/1‑143 (MQ=255)
               aGATTAGGGCGTGGGGTCTGACTGTTATAGGGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAacgacg                                                                                                                              <  1:1041604/143‑1 (MQ=255)
               aGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCCCGTAGGTGATGAacgacg                                                                                                                              >  1:566377/1‑143 (MQ=255)
                  ttccGGCGTATGTTCTGACGTTTATTGTGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTAAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAATTCCGCGCTTTCCGTATATCGCTTAGGTGATGAACGACGCTc                                                                                                                           <  1:794437/140‑1 (MQ=255)
                  ttACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTc                                                                                                                           >  1:612050/1‑143 (MQ=255)
                   tACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTcc                                                                                                                          <  1:757410/143‑1 (MQ=255)
                       gCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCAGCa                                                                                                                      >  1:32386/1‑143 (MQ=255)
                         gTAGGTTCTGCCGTTTCTAGCGTGTCCTCCTTCTCCGTGCTGGCGTGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGt                                                                                                                    <  2:105793/143‑1 (MQ=255)
                           aagtacTGACGTTTATAGCGCGACCTCCTTCACCGAGCTGGCGCGTTATGGTCACGATTGTGAACGCTGGAACATGCCGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTgg                                                                                                                  <  1:467939/138‑1 (MQ=255)
                               tCTGTCTTTTATTGCGTGATCTCCTTCTCCGGGCTGGCGGGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATc                                                                                                              <  2:904985/143‑1 (MQ=255)
                                  gACGTTTATAGCGGGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCACCACCAAAGGTTCTTTGGTGATGAACGACGCTCCGGCAGTGGCATCTAc                                                                                                           <  1:795208/143‑1 (MQ=255)
                                   aCGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTAcc                                                                                                          >  1:1028488/1‑143 (MQ=255)
                                   aCGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGGCAGGATGGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGCCGCTCCGCCAGACGCATCTAcc                                                                                                          >  2:157718/1‑143 (MQ=255)
                                         ataGCGTGACCTCCTTCACCGAGCTGGCGCGTCATGGTAATGACTCTTAGCGCTGGCCCCTATGACACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCGTAGGGGATTAACGACGCTCCGGCAGTGGCATCTACCGACtat                                                                                                    <  1:1031819/143‑1 (MQ=255)
                                         ataGCGGGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACtat                                                                                                    <  1:91206/143‑1 (MQ=255)
                                             cGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGCTTGTGAACGCTGGAACCTGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCGGCGGAATCTACCGACGATAGGa                                                                                                >  1:394260/1‑143 (MQ=255)
                                             cGTGACCTCCTGGATCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGa                                                                                                <  2:739899/143‑1 (MQ=255)
                                               tgcggtccgTCAACGAGCTGGCGCGTGATGGTCAGGATTGTGAGCGCTGGAAGATGCTGCACCCGCGGGAAACTCCGCGCGGTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGaaa                                                                                              <  2:232511/138‑1 (MQ=255)
                                               tgcccccctTGACCCAGCTGGCGCGTGATGGGCAGGATTGGGAAGGGTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGaaa                                                                                              <  1:1128622/137‑1 (MQ=255)
                                                  cctcctTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTg                                                                                           >  1:198241/1‑143 (MQ=255)
                                                    tcctTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGtt                                                                                         >  1:169532/1‑143 (MQ=255)
                                                          aCCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGa                                                                                   <  1:508056/143‑1 (MQ=255)
                                                           ccGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGCTGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTGCGCTgcg                                                                                  >  2:325041/1‑141 (MQ=255)
                                                              ggctgtggcgTGAGGGTCAGGGTGGTGAACGCTGGAACATGCTGGAACGGCGGTGAGCTTCGCGCCTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAg                                                                               <  2:513058/136‑1 (MQ=255)
                                                                cTGGCGCGTGATGCCAACGCTAGTGTCCGCTGCAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTTTATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGt                                                                             <  2:1059000/143‑1 (MQ=255)
                                                                 tGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATCTCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTc                                                                            >  2:834324/1‑143 (MQ=255)
                                                                 cggCGCGTGAGGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTc                                                                            <  2:391707/142‑1 (MQ=255)
                                                                       gTGATGGTCAGGATTGTGAACGCGTGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTGAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACt                                                                      <  2:784661/143‑1 (MQ=255)
                                                                         gATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTa                                                                    <  2:686206/143‑1 (MQ=255)
                                                                             gTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTCATTACGTa                                                                >  1:150995/1‑143 (MQ=255)
                                                                               cAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGGACGACGCTCCGGCAGTGGCATCTACCGACTATATGCAACTGTTCGCTTAGCAGGTCCGTTCTTAcctccc                                                              >  2:709296/1‑138 (MQ=255)
                                                                                  gATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGGTCCGTATATCGCTTGGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCGGGTCCGTACTCACGTACCGGc                                                           >  1:21598/1‑143 (MQ=255)
                                                                                         aaCGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGCTGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGCGCAGGTCCGTACTTACGTCCCGGCTgacgac                                                    >  1:232798/1‑143 (MQ=255)
                                                                                                      gctgcACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGCCGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCGGTCCTTACGTACCGGCTGACGAGTACCGCCTACTgg                                       >  1:1028821/1‑143 (MQ=255)
                                                                                                      gctgcACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTGAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTgg                                       <  2:198063/143‑1 (MQ=255)
                                                                                                       cggctcgcGCTGGACATGCCGTGGGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTACTCGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTggg                                      <  2:392340/136‑1 (MQ=255)
                                                                                                         gcACCCGCTGGAAACTCCGCGCGTCCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGAATATAAGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTa                                    <  2:1071554/143‑1 (MQ=255)
                                                                                                                            cgcgTTCCGTATATCGCTTAGGTGATGAACGACGCCCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTc                 <  1:171052/143‑1 (MQ=255)
                                                                                                                                ttCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTGCGTTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTGCTGGGTACTGATGGCTTCGGTCGTTGCggc             >  1:1037575/1‑141 (MQ=255)
                                                                                                                                 tCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTGCGGCa            >  1:51733/1‑143 (MQ=255)
                                                                                                                                      atatCGCCTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGCGTATGTGAAACTGTTCGCTTAGCGGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGt       <  1:1095933/143‑1 (MQ=255)
                                                                                                                                        atCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAACAGGTCCGTACTTACGTACCGGCTGACGCATACCGCGTACTGGGTAATGATGGCATCGGTCGTTCCGACAGCCCTga     >  2:110264/1‑143 (MQ=255)
                                                                                                                                         tCGCTTGGGTGATGAACGAAGCTCCGGCAGTGGCAACTACCGACTCTATGAAACTGTTCGCTGAGCGGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGATACTGATGGCATCGGTCGTTCTGACAGCCATgag    >  2:108743/1‑143 (MQ=255)
                                                                                                                                         tCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCGGAGAGCCGTgag    >  1:539193/1‑143 (MQ=255)
                                                                                                                                         tCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACCGCCGTgcg    >  2:406882/1‑141 (MQ=255)
                                                                                                                                           gcgcAGGGGACGAACGCCGCTGCGGCAGTGGCATCGACCGACTATAGGAAACTGTTCGCTGAGCAGGTCCGTACTGACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGAGaa  <  2:937545/140‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                           
TGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGAGAA  >  NC_000913/125245‑125526

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