Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 125,387 C→T 100% 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,3870CT93.7% 105.4 / ‑4.5 46Q791* (CAG→TAG) aceEpyruvate dehydrogenase E1 component
Reads supporting (aligned to +/- strand):  ref base C (1/1);  major base T (20/19);  minor base G (2/2);  total (24/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.

GATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGAGAAC  >  NC_000913/125248‑125527
                                                                                                                                           |                                                                                                                                            
gaTCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAgg                                                                                                                                           <  1:895592/143‑1 (MQ=255)
     tGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAAATACGCGCGTTCCGGATATCCCTTAGGggatg                                                                                                                                      >  2:187049/1‑143 (MQ=255)
      ggCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTAAGGATTGTGAACGCTGGAACATGCTGCACCCGCCGGAAACTCCGCGCGTTCCGGATATCACGCAGGtgaggg                                                                                                                                     >  2:107766/1‑140 (MQ=255)
       gCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTGAGGCGATGaa                                                                                                                                    >  2:637917/1‑143 (MQ=255)
        cGAAAGATTACGGCGGAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAac                                                                                                                                   <  1:633190/143‑1 (MQ=255)
                  cGGCGTAGGTTCTGACGGTTCTAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCgg                                                                                                                         <  2:248594/143‑1 (MQ=255)
                       tAGGTTCTGACGTTTATAGCGTGACCTCCTTCCCCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTg                                                                                                                    <  1:783539/143‑1 (MQ=255)
                          gTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGCGGCa                                                                                                                 >  2:137116/1‑143 (MQ=255)
                                  gTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGa                                                                                                         <  1:731902/143‑1 (MQ=255)
                                        aGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTatat                                                                                                   <  1:114878/143‑1 (MQ=255)
                                          cGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGa                                                                                                 >  1:880061/1‑143 (MQ=255)
                                            gcacttgctTCACCGGGCGGGCGCGTGATGGCCAGGATTGTGAACGCTGGAACATGCTGCCCCCGCTGTAAACTCCGCGCGTTCCGTCTATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGaaa                                                                                               <  2:771284/136‑1 (MQ=255)
                                              aCCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACACGCTGGAAACTCCGCGCGTTCCGTAGATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCGTCTACCGACTCTATGAAACt                                                                                             >  1:703457/1‑143 (MQ=255)
                                                 tcctTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGACCGCCGCTCCGTCATGGGCATCTACCGACTATATGAAACTGtt                                                                                          >  1:620690/1‑143 (MQ=255)
                                                   ctTCCCCGAGCTGGAGCGTGACGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCg                                                                                        <  2:220808/143‑1 (MQ=255)
                                                   ctTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAGCTGTTCg                                                                                        >  1:368685/1‑143 (MQ=255)
                                                   ctTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTAGATGAAACTGTTCg                                                                                        >  2:301336/1‑143 (MQ=255)
                                                         cGAGGTGGCGCGGGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGc                                                                                  <  1:127908/143‑1 (MQ=255)
                                                             cttgggcgTGCAGGCCAGGAGTGTGGACGCTGGAACATGCTGCACCCGCTGGAACCTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGt                                                                              <  1:861633/138‑1 (MQ=255)
                                                               ggCGCGTGATGCTCAGGATTGCGAACGCTGGAACATGCTGCACCCGCTGGAAACGCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTcc                                                                            <  1:610289/143‑1 (MQ=255)
                                                                 cgcgTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGt                                                                          <  1:770748/143‑1 (MQ=255)
                                                                      gATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTCGGCGATTAACGTCGCTCCGGCAGACGCAACTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGGACttg                                                                     >  1:468730/1‑142 (MQ=255)
                                                                           tCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTCCGTAc                                                                >  1:390632/1‑143 (MQ=255)
                                                                                  tgtgAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCTGTATATCGCTTAGGTGATGAACGACGTTCCGGCAGTGGTAGCTACCGACTATATGAAACGGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTga                                                         <  1:742489/143‑1 (MQ=255)
                                                                                   gtgAACGCTGGAACATGCTGCACGCGCTGGAAACTCCGCGCGTTCCGTATATCGCATAGGTGATGACCGACGCTCCGGCAGTGGCATCCACCGCCTATATGAAACTGTTCGCTGCGCAGGTCCGTACTTACGTACCTGCAgac                                                        >  2:438521/1‑143 (MQ=255)
                                                                                    tgAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGAGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAAATGTTCGCTGAGCAGGTCCGTACTTACGGACCGGGTgacg                                                       >  1:99866/1‑143 (MQ=255)
                                                                                     gAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAAGGACGGTCAGGCAGTGGCATCTACCGACTATTTGAAACGGTGCGGTGAGCAGGTCCGTAGTTACGTACCGGCTgcgga                                                      >  1:170846/1‑139 (MQ=255)
                                                                                        cGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTa                                                   >  1:387822/1‑143 (MQ=255)
                                                                                         gCTGGAACAATCTCCACCCGCTTGAAACTCTCCGCGTCCCGTATATCGCTTAGGTGACGAACGACGCTCAGGCAGTGGCATCCACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTAc                                                  <  2:53547/143‑1 (MQ=255)
                                                                                             gAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACcgcg                                              >  2:470011/1‑143 (MQ=255)
                                                                                             gAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTAAGGTGCTGAACGACGCCGCGGCAGGGGCATGTACCGGCTATGTGGAACTGTTCGCTGGGCGGGGCCGCACTTACGTACCGGCCGACGACCACcgcg                                              >  2:628264/1‑143 (MQ=255)
                                                                                             gAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCCGAGGTGATGAACGACGCTCCGGCCGTGGCATCTACCGACTATCTGAACCTGTTCGCCAAGCCGGTCCGTACTAACCTACCGGACCACGACTGCcgcg                                              >  2:789558/1‑143 (MQ=255)
                                                                                                  ctctgaccccGCTCGAAGCGCCGTGCGGTCCGTACATCGCTGATGTGATGAACGACGCTCCGGCAGTGGTATCTACCGACTATATGAAGCTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTg                                         <  2:224501/136‑1 (MQ=255)
                                                                                                   gctccACCCGGTGGAAGCTCCGTGCGTCCCTTACATCGCTTAGGTGTTGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAATCTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTgg                                        <  2:423078/143‑1 (MQ=255)
                                                                                                       cACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTAc                                    >  2:315503/1‑143 (MQ=255)
                                                                                                       cACCCGCGGGAACCTCCGCGCGTTCCGTATATCGCGTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTAc                                    <  1:187755/143‑1 (MQ=255)
                                                                                                                     tCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTGCTGGGTACTGATGGCTTCGGTc                      >  1:461573/1‑143 (MQ=255)
                                                                                                                       cgcgcgTTCCGTATATCGCTTAGGTGATGACCGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGt                    <  1:137116/143‑1 (MQ=255)
                                                                                                                           cgTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCg                <  1:637917/143‑1 (MQ=255)
                                                                                                                                cGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGCAACTGTTCGCTGAGCGATTCCGTATTTAGGTACGCGCTGACGGGAACCGCGTACTGGGTACTGATGGCTTCGGTCGTTGCGACAgg           >  1:716136/1‑142 (MQ=255)
                                                                                                                                 gTAGATCGCGGAGGTCATGAACGACGATGCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGcc          <  1:597403/143‑1 (MQ=255)
                                                                                                                                    tatCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTg       <  2:628073/143‑1 (MQ=255)
                                                                                                                                    tatCGCTTAGGTGATGAACGACGCTCCGGCAGCGGCATCCACCGACTGTATGAAACTGTTCGCCGAGCGGGTGCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACGGACGGCTTCGGGCGTCCCGACAGCGGTg       >  2:243393/1‑143 (MQ=255)
                                                                                                                                       cGCTTAGGTGATGAACGACGCTGCGGCAGTGGCATCGACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACGGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGCGAGGGGTgaga    >  1:505513/1‑143 (MQ=255)
                                                                                                                                       cGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTAGTGATTATGTGAGACTGTGCGCGGAGAGGGTCAGTTTGCACCCACCGGCTGACGACTACCGCGTACTGGCCACTGACGGCTTCCGTCGTTCCGACAGCCGCgaga    >  1:673431/1‑143 (MQ=255)
                                                                                                                                         cTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGAGAAc  <  1:301336/143‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                            
GATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGAGAAC  >  NC_000913/125248‑125527

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: