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,3870CT100.0% 147.8 / ‑6.7 52Q791* (CAG→TAG) aceEpyruvate dehydrogenase E1 component
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (22/28);  total (23/29)
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.

AGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATT‑‑GTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATC‑GCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGA  >  NC_000913/125247‑125523
                                                                                                                                               |                                                                                                                                        
gggtCCTGGCGAAAGATTTCTGCGTAGGTTCTGGCGTTTATAGCGTGACCTCCGTCACCGAGCTGGCGCGTGATGGTCAGGATT‑‑GTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATC‑GCTTAg                                                                                                                                        <  2:627255/140‑1 (MQ=255)
  aTCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACATCCTTAACGGAGC‑GG‑AAAAGCACACCAGATTTACGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATC‑GCTTAGGt                                                                                                                                      >  1:1048748/1‑143 (MQ=255)
          gAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATT‑‑GTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATC‑GCTTAGGTGATGAacg                                                                                                                              <  1:742648/143‑1 (MQ=255)
           acggactACGGCGTCGGGTCTGGCGGTTAGGGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATT‑‑GTGAACGCTGGAACATGCTGCCCCCGCTGGAAACTCCGCGCGTTCCGTATATC‑GCTTAGGTGATGAacga                                                                                                                             <  1:101619/137‑1 (MQ=255)
           aaaGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACATCCTTCACCGAGCTGGCGCGTGATGGTCAGGATT‑‑GTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTGCCGTACAGC‑GCTCGGGTGATGAacga                                                                                                                             >  2:234242/1‑143 (MQ=255)
           aaaGATTACGGCGTAGGTGCTGACGTTTATAGCGGGACCTCCTGCNCCGGGCTGGCGCGTGATGGTCAGGATT‑‑GTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATC‑GCTTAGGTGATGAacga                                                                                                                             <  2:527934/143‑1 (MQ=255)
             atagtACCGCGCAGGTTGTGACGGTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATT‑‑GTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATC‑GCTTAGGTGATGAacgacg                                                                                                                           <  2:251117/139‑1 (MQ=255)
              gATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATT‑‑GTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATC‑GCTTAGGTGATGAACGACGc                                                                                                                          >  2:975748/1‑143 (MQ=255)
                  aCGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATT‑‑GTGAACGCTGGGACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATC‑GCTTAGGTGAGGAACGACGCACCg                                                                                                                      >  2:632588/1‑143 (MQ=255)
                    ggCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATT‑‑GTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATC‑GCTTAGGTGATGCACGACCCGCCGGc                                                                                                                    >  1:617529/1‑143 (MQ=255)
                       gTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATT‑‑GTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATC‑GCTTAGGTGATGAACGACGCTCCGGCAGt                                                                                                                 <  1:883301/143‑1 (MQ=255)
                        tAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATT‑‑GTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATC‑GCTTAGGTGATGAACGACGCTCCGGCAGTg                                                                                                                <  2:352281/143‑1 (MQ=255)
                        tAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATT‑‑GTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATC‑GCTTAGGGGATGAACGACGCTCCGGCAGTg                                                                                                                >  1:449522/1‑143 (MQ=255)
                        tAGGTTCTGACGTTTATAGCGTGACCTCCCTCACCGAGCTGGCGCGTGATCGTCACGATT‑‑GTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATC‑GCTTAGGTGATGAACGACGCTCCGGCAGTg                                                                                                                <  1:270826/143‑1 (MQ=255)
                            ttCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATT‑‑GTGAACGCTGGAACGTGCTGGGCCGGCTGGAAACCCCGCGCGTTCCGTATATC‑GCTTCGGGGGTGAACGCCGCCCCGGCCGGGGcgt                                                                                                            >  2:913368/1‑141 (MQ=255)
                            ttCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATT‑‑GTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATC‑GCTTAGGTGATGAACGACGCTCCGGCAGTGGCAt                                                                                                            >  2:630337/1‑143 (MQ=255)
                             tCTGACGTTTATGGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATT‑‑GTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCACGTTCCGTATATC‑GCTTAGGTGATGAACGACGCTCCGGCAGTGGCATc                                                                                                           >  2:426489/1‑143 (MQ=255)
                                     gtataGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATT‑‑GTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATC‑GCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACt                                                                                                   <  2:1036357/142‑1 (MQ=255)
                                        taGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATT‑‑GTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATC‑GCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACtata                                                                                                >  2:123178/1‑143 (MQ=255)
                                          gCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATT‑‑GTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTAGATC‑GCTTAGGTGCTGAACGACGCTCCGGCAGTGGCAACTACAGAcgctata                                                                                              >  1:29131/1‑137 (MQ=255)
                                                    ctTCTCCGGGCTGGCGCGTGATGGTCAGGATT‑‑GTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATC‑GCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCg                                                                                    <  2:915817/143‑1 (MQ=255)
                                                    ctTCACCGAGCTGGCGCGTGATGGTCAGGATT‑‑GTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATC‑GCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCg                                                                                    >  1:393400/1‑143 (MQ=255)
                                                    ctTCACCGAGCTGGCGCGTGATGGTCAGGATT‑‑GTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATC‑GCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCg                                                                                    <  1:975748/143‑1 (MQ=255)
                                                      tCACCGAGCTGGCGCGTGATGGTCAGGATT‑‑GTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATC‑GCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCt                                                                                  >  1:1029174/1‑143 (MQ=255)
                                                      tCACCGAGCCGGCGCGTGATGGTCAGGATT‑‑GTGAACGCTGGCATGTTCTGATGAGACTTCAAACTACGCGATCTCCGTATATC‑GCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCt                                                                                  <  2:1024575/143‑1 (MQ=255)
                                                       cACCGAGCTGGCGCGTGATGGTCAGGATT‑‑GTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATC‑GCTTCGGTGATGAACGACGCCCCGGCAGTGGCATCTACCGACTATATGAACCTGttgcctg                                                                                 >  2:732514/1‑138 (MQ=255)
                                                       cACCGAGCTGGCGCGTGATGGTCAGGATT‑‑GTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATC‑GCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTg                                                                                 <  2:712914/143‑1 (MQ=255)
                                                             gCTGGCGCGTGATGGTCAGGATT‑‑GTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATC‑GCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAgg                                                                           <  1:630337/143‑1 (MQ=255)
                                                                ggCGCGTGATGGTCAGGATT‑‑GTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATC‑GCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTcc                                                                        <  2:13306/143‑1 (MQ=255)
                                                                 gcgcGTGATGGTCAGGATT‑‑GTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATC‑GCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGCGCAGGTCCg                                                                       >  1:8118/1‑143 (MQ=255)
                                                                   gcgTGATGGTCAGGATT‑‑GTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATC‑GCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGCCTATATGAAACTGGTGGCCGACCAGGTCCGTa                                                                     >  2:964632/1‑143 (MQ=255)
                                                                       ggtGGTCAGGATT‑‑GTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATC‑GCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTa                                                                 <  1:744888/141‑1 (MQ=255)
                                                                       gATGGTGAGGAGT‑‑GTGAGCGCTGGAACATGCTGCTCCCGCTGGAACCTCCGCGCGTTCCGTATATC‑GCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTa                                                                 <  2:788314/143‑1 (MQ=255)
                                                                            tCAGGATT‑‑GTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATC‑GCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTAc                                                            <  1:852204/143‑1 (MQ=255)
                                                                            tCAGGATT‑‑GTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATC‑GCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGATCCGTCCTTACGCAc                                                            >  2:334742/1‑143 (MQ=255)
                                                                             cAGGATG‑‑GTGAACGCTGTGTCATGGGGCACCCGCTGGAAACTCCGCGCGTTCCGTATATC‑GCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTAcc                                                           <  2:12117/143‑1 (MQ=255)
                                                                                  tt‑‑GTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATC‑GCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTg                                                      >  1:564104/1‑143 (MQ=255)
                                                                                               agcaCATGCTGCACCCGCTAGACGCTGCGCGCGTGCCGTACGGCGGCTC‑GGCGACGAACGACGCTCCGGGAGCGGCATCTACCGACTACATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACcgc                                           <  1:717670/140‑1 (MQ=255)
                                                                                                gAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATC‑GCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTCCTTACGTACCGGCGGCCGACTACcgcg                                          >  2:331196/1‑143 (MQ=255)
                                                                                                 aaCATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATC‑GCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGt                                         <  1:154237/143‑1 (MQ=255)
                                                                                                      gctgcCCCCCCTGGAAACTCCGCGCGTTCCGTATATC‑GCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTgg                                    <  1:62224/143‑1 (MQ=255)
                                                                                                      gctgcACCCGCTGGAAACTCCGCGCGTTCCGTATATC‑GCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTCACGACTACCGCGTACTgg                                    >  1:917949/1‑143 (MQ=255)
                                                                                                        tgcACCCGCTGGAAACTCCGCGCGTTCCGTATATC‑GCTTAGGTGATGAACGACGCTACGGCAGTGGCGTCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACCAACGCACCGGCTGACGAAGACCGCGCGCTGGGt                                  >  2:454856/1‑143 (MQ=255)
                                                                                                          cACCCGCTGGAAACTCCGCGCGTTCCGTATATC‑GCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTAc                                <  1:211759/143‑1 (MQ=255)
                                                                                                            cccGCTGGAAACTCCGCGCGTTCCGTATATC‑GCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACGGTTCGCTGAGCAGGTCCGTGCTTACGTTCCGACTGACAACTACCGCGTACTGGGCAcag                              >  1:37474/1‑141 (MQ=255)
                                                                                                                cTGGAAACTCCGCGCGTTCCGTATATC‑GCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATgg                          <  2:393400/143‑1 (MQ=255)
                                                                                                                   ggaaCTCCGCGCGTTCCGTATATC‑GCTTAGGTGATGATCGACGCTCCGGCAGTGGCATCTACCGACTATATAAAACTGTGTGCGGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCtt                       <  2:692741/141‑1 (MQ=255)
                                                                                                                   gAACCTCCGCGGGTTCCGTGTATC‑GCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGGGCAGGGCCGCACTTACGTTCCGGCTGACGACTACCGCGTACTGGGTACTGATGGCtt                       <  2:696933/143‑1 (MQ=255)
                                                                                                                          cgcgcgTTCCGTATATC‑GCGTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGt                <  2:253369/143‑1 (MQ=255)
                                                                                                                               gTTCCGTATATC‑GCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGa           <  1:291044/143‑1 (MQ=255)
                                                                                                                                 tCCGTATATC‑GCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACa         <  1:426489/143‑1 (MQ=255)
                                                                                                                                        atC‑GCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCTGACGTACGTACCGGCTGACGCCTACCGCGTACTGGGTACGGCGGGCTTCGGTCGTTCCGACAGCCGTga  >  2:907712/1‑143 (MQ=255)
                                                                                                                                               |                                                                                                                                        
AGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATT‑‑GTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATC‑GCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGA  >  NC_000913/125247‑125523

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: ‑