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,3870CT99.7% 111.1 / ‑7.3 50Q791* (CAG→TAG) aceEpyruvate dehydrogenase E1 component
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base T (27/20);  minor base G (2/0);  total (29/21)
Fisher's exact test for biased strand distribution p-value = 5.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGT  >  NC_000913/125250‑125521
                                                                                                                                         |                                                                                                                                       
tCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGtg                                                                                                                                    >  1:447441/1‑143 (MQ=255)
    ggCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCAGCGAGCTGGCGCGTGAGGGTAAGGCTTGTGAACGCTGGAACATGCTGCACGCGCTGGAAACTTCGCGCGTTCAGTATGTGGCTTAGGagatga                                                                                                                                >  2:433120/1‑143 (MQ=255)
       gAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCCACCGCTCGACTGCCGGCGCGTCCCGTATATCGCTTAGGttatcgacg                                                                                                                             >  1:529531/1‑138 (MQ=255)
       gAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAacg                                                                                                                             >  1:421768/1‑143 (MQ=255)
       gAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACAGGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAacc                                                                                                                             >  2:486196/1‑142 (MQ=255)
       gAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGAAGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAacg                                                                                                                             >  2:313491/1‑143 (MQ=255)
             ttACGGCGTAGGTCCTGACGTTTATAACGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCTCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTc                                                                                                                      <  1:239881/143‑1 (MQ=255)
                 ggCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTGCCGTATATCGCTTAGGTGATGAACGACGCTCCGGc                                                                                                                  >  1:542700/1‑143 (MQ=255)
                       ggTTCTGGCGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGc                                                                                                            <  2:399344/143‑1 (MQ=255)
                       ggTTCTGACGTTGATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGc                                                                                                            <  2:484821/143‑1 (MQ=255)
                              aCGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATCTCGCTTAGGTGCTGAACGACGCCCCGGCCGTGGCCTCTAcc                                                                                                     >  2:931002/1‑143 (MQ=255)
                              aCGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTAcc                                                                                                     >  1:5468/1‑143 (MQ=255)
                                        cGTGCCCCCCCTCACCGAGCTGGCCGCTGAGGGTCCGCGCGGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGa                                                                                           <  1:172213/143‑1 (MQ=255)
                                        cGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGCTTGTGAACGCTGGAACATGCTGCACCCGCTGGAACCTCCGCGCGTTCCGTATATCGCTTAGGTGGTGAACGACGCTCCGGCCGTGGCATCTCCCGACGATATGa                                                                                           >  1:217216/1‑143 (MQ=255)
                                        cGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCGTCTACCGACTATATGa                                                                                           >  1:444724/1‑143 (MQ=255)
                                          tGATCTCCTTCACCGAGCTGGGGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGaaa                                                                                         <  2:76809/143‑1 (MQ=255)
                                          tGACCGCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGaaa                                                                                         <  1:694622/143‑1 (MQ=255)
                                             cctcctTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAGACTCCGGGCGTTCGGTATGTCGCTGAGGTGGTGAAGGAAGCTCCGGCGGGGGCAGCTACGGCGTATAGGAGACTg                                                                                      >  1:5969/1‑143 (MQ=255)
                                               tcctTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGtt                                                                                    <  1:376254/143‑1 (MQ=255)
                                                gctTCACCGAGCTGGCGCGTGAGGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTc                                                                                   <  1:178393/142‑1 (MQ=255)
                                                  ttCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGGGATGGACGACGCGCCGGCAGTGGCAGCTACCGACTATATGAAACTGTgcgc                                                                                 >  1:533005/1‑143 (MQ=255)
                                                    cACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTg                                                                               >  1:858232/1‑143 (MQ=255)
                                                      ccGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTggg                                                                             >  1:795986/1‑141 (MQ=255)
                                                          gctctggcgTGATGGTCAGGATTGCGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCCTCTACCGACGATATGAAACTGTTCGCTGAGCAgg                                                                         <  2:335329/137‑1 (MQ=255)
                                                           cTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATCTCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACGATATCAACCTGATCCCTGAGCAGGt                                                                        >  2:869641/1‑143 (MQ=255)
                                                                  gTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGATCAGCAGGTCCGtttt                                                                 >  1:928617/1‑140 (MQ=255)
                                                                    gATGGTCAGGATAGTGAACGATGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTGCGTATATCGCGTAGGTGATGAACAACGCTCCGGCAGTGGCAGACACCGACTATATGAAACTGGTCGCGTAGCGGGGCCGTACtga                                                               >  2:341523/1‑141 (MQ=255)
                                                                     atggggaGGATTGTGAAGGCGGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTAc                                                              <  1:908143/137‑1 (MQ=255)
                                                                        gTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTa                                                           >  2:736068/1‑143 (MQ=255)
                                                                        gTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTGAGGGGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACCGCACCCTGAGCAGGTCCGTACTTACGTa                                                           >  2:367239/1‑143 (MQ=255)
                                                                        gTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCATAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTa                                                           >  2:341342/1‑143 (MQ=255)
                                                                         tCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCgtacgtacgtac                                                          >  1:697804/1‑143 (MQ=255)
                                                                         gcAGGATTGTGAACGCTGGAACATGCTGCACCCGGTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTAc                                                          <  2:284856/142‑1 (MQ=255)
                                                                          gaGGGTTGTGAACGCTGGAAGATGCTGCACGCGCTGGCAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTAcc                                                         <  2:350851/142‑1 (MQ=255)
                                                                           cggAGTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCg                                                        <  1:495700/142‑1 (MQ=255)
                                                                            ggATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTGCTTGCGTACCgg                                                       >  1:668191/1‑143 (MQ=255)
                                                                             gATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAGACTGTTCGCTGGGCAGGTCCGTACTTACGTACCGGc                                                      >  1:449164/1‑143 (MQ=255)
                                                                               ttGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTGTATTGCCTAGGTGCTGACAGCGCAGCCATCACTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACGTACGTACCGGccg                                                    >  1:690797/1‑141 (MQ=255)
                                                                                   gAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCATAGGTGGTGAACGACGCTCCGGCAGTGGCATGTACCGACTATATGAAGCTGTTCGCGGAGCAGGTCCGTACATACGTAGCGGCTgcaga                                                >  2:857603/1‑139 (MQ=255)
                                                                                      cGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCCCCGGTAGTGGCATCTACCGACTATAGGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTCGCGACTa                                             >  2:120367/1‑143 (MQ=255)
                                                                                             aCATGCTGCACCCGCGGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGCGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTa                                      <  1:478122/143‑1 (MQ=255)
                                                                                                 gctgcACCCGCTGGAAACTCCGCGCGTTCCGTGTATCGCTTAGGTGATGAACGACGCTTCGGAAATGGCATCTACCGACTATATGAAACTGTTCGCTGAGCGGGTCCGTACTTACGTACGGGCTGTCGACTACCGCGTACTgg                                  >  1:117301/1‑143 (MQ=255)
                                                                                                         cGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGAt                          <  1:569649/143‑1 (MQ=255)
                                                                                                            tGGAAACTCCGCGCGTTGCGTATATCGCTCAGGTCATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGc                       <  2:661180/143‑1 (MQ=255)
                                                                                                                 aCTCCGCGCGTTCCGTATATCGCTTAGGTGATGGATGGCGCTCCGTCAGGGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCgg                  <  2:412559/143‑1 (MQ=255)
                                                                                                                      gcgcgTTCTGTAGATCGCTTAGGTGATGACCGTCGTTCCGGCAGTGTCATCTACGGTTTTTCTGATACTGTTCTCTTAGCTGGGCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGtt             <  2:873167/143‑1 (MQ=255)
                                                                                                                        gcgTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTcc           <  1:915596/143‑1 (MQ=255)
                                                                                                                              cGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGc     <  1:736068/143‑1 (MQ=255)
                                                                                                                               gTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGATCGCTAAGCGGTTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTTGTTCCGACAGct    >  1:352534/1‑142 (MQ=255)
                                                                                                                                 gtatCGCTTAGGCGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTAGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGCTGGCAGCGGGCGTTCCGACAGCCGt  <  2:573580/142‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                       
TCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGT  >  NC_000913/125250‑125521

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: