Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 125,085 2 bp→CA coding (2069‑2070/2664 nt) aceE → pyruvate dehydrogenase E1 component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913125,0850TC98.0% 146.2 / NA 50I690T (ATC→ACA) aceEpyruvate dehydrogenase E1 component
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (27/22);  minor base G (0/1);  total (27/23)
Fisher's exact test for biased strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_000913125,0860CA98.0% 137.4 / NA 50I690T (ATC→ACAaceEpyruvate dehydrogenase E1 component
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (26/22);  total (27/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATC  >  NC_000913/124948‑125221
                                                                                                                                         ||                                                                                                                                       
tCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTa                                                                                                                                     >  1:257681/1‑143 (MQ=255)
    aGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTAAGAAGTGGCTGTCCTCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGcg                                                                                                                                 >  1:368541/1‑142 (MQ=255)
      ccACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGa                                                                                                                               <  1:807265/143‑1 (MQ=255)
       cacaTTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCAACATGCATGACGGACTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTCTAATACTAACCAACGACGCTGaa                                                                                                                              >  2:637066/1‑143 (MQ=255)
        acaTTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAACACAAGAGAACGTTTACTCCTACACAACTACGCTGacc                                                                                                                             >  2:415231/1‑141 (MQ=255)
                tCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAAcaa                                                                                                                     >  2:764312/1‑141 (MQ=255)
                    tGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGCCTGGAGCGTATGTACGGTGAAAGACAAGAGAACGTTTACTACTACACCACTCCGCTGAACGAAGCCTACcac                                                                                                                 >  1:479699/1‑143 (MQ=255)
                          tCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGCGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTAAACAAGTACGCTGAACGAAAACTACGACATGCCg                                                                                                           >  2:162462/1‑143 (MQ=255)
                             cGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCa                                                                                                        <  1:433849/143‑1 (MQ=255)
                                   gTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCg                                                                                                  <  2:476538/143‑1 (MQ=255)
                                      tatctTACGACCCGGCTTACGCTTACGAAGTTGCTCTCATCATGCATGACGGTCTGAAGCGTATGTACCGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACACGCCGGCAATGCCGGaa                                                                                               <  2:999894/141‑1 (MQ=255)
                                       ctctTACGACCCTGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAg                                                                                              <  1:331208/143‑1 (MQ=255)
                                       ctctTACGACCCGGCTTACGCTTCCGAACTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTCCGGTGCAAAACAAGAGCACGTTTACTACTACACAACTACGCTGCCCGAAACCTCCCCCACGCCGCCAATGCCggccg                                                                                              >  1:1142517/1‑140 (MQ=255)
                                              gACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGa                                                                                       <  1:596426/143‑1 (MQ=255)
                                                  cGGCTTACGCGTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGaa                                                                                   <  2:72450/143‑1 (MQ=255)
                                                    gcttacgcttacGAAGTTGCTGTCAGCATGCATGACGGTCTGGAGCGGATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACGACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAgg                                                                                 <  1:1012229/143‑1 (MQ=255)
                                                      ttacgcttacgAAGTTGCTGTCATCATGCATGACGGACTGGAGAGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGtc                                                                               >  2:87528/1‑142 (MQ=255)
                                                         cgcttacgAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATcc                                                                            <  1:617177/143‑1 (MQ=255)
                                                          gcttacgAAATTGCTATCATAATGCATGACGATCTCGAGCGTAAGTACGGTGAAAAACAAGAGAACGTTTAGTACTGCACAACTAGGCTCAACGAAAACTAACACATGCCGGCAATGCCGGAAGGGGCTGAGGGACGTATcca                                                                           >  2:177322/1‑142 (MQ=255)
                                                               cgAAGTTGCTGTCATCATGCATGGCGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTTCCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAg                                                                      <  1:404060/143‑1 (MQ=255)
                                                                gAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTCTCCGTAAAgg                                                                     >  1:147376/1‑143 (MQ=255)
                                                                          gcatcatGCATGACGGTCTGGAGCGTATGTACGGTGAAGCGACGGAATGCGTGTGCTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACaaa                                                           <  2:1069750/142‑1 (MQ=255)
                                                                            atcatGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCAGGAAAGTCCAGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACACACt                                                         >  1:520407/1‑143 (MQ=255)
                                                                             tcacGCATGACGGGCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTc                                                        <  1:123827/143‑1 (MQ=255)
                                                                               atgcatgACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGCAGGTATCCGTAAAGGTATCTACAAACTCGa                                                      >  1:603036/1‑143 (MQ=255)
                                                                                 gcaggACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGaaa                                                    <  1:921491/143‑1 (MQ=255)
                                                                                   atgACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACt                                                  >  1:235548/1‑143 (MQ=255)
                                                                                     gACGGTCAGGGGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGCAGGGGCTGAGGAAGGTATCCGTAACGGTATCTACAAACTCGCAACTAt                                                >  2:613554/1‑143 (MQ=255)
                                                                                      aCGGTCTGGAGCGTACGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTAtt                                               <  2:598900/143‑1 (MQ=255)
                                                                                       cGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTg                                              >  1:409702/1‑143 (MQ=255)
                                                                                       cGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTAAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTg                                              >  2:139652/1‑143 (MQ=255)
                                                                                       cGGTCGTGAGCGTGTGTAGGGTGAAAAACAGGAGAGCGTTTACTACTGCACAACTGCGCTGAAGGGACACTCGCCCGCGCCGGCGATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTg                                              <  2:815248/143‑1 (MQ=255)
                                                                                            tGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCCGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAggg                                         >  2:475559/1‑142 (MQ=255)
                                                                                              gAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAg                                       >  1:783509/1‑143 (MQ=255)
                                                                                                gCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGCCAACGCCGGAAGGTGCTGAGGAAGGTATCCGCAAAGGTAGCTACAACACCGAAACTATCGCAGCTAGCa                                     >  2:1222591/1‑143 (MQ=255)
                                                                                                      gTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTa                               >  2:859652/1‑143 (MQ=255)
                                                                                                        aCGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACCAAAACCACCACATGCCGGCAACGCCGCAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTaaa                             <  1:87528/143‑1 (MQ=255)
                                                                                                          ggTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGt                           >  1:549676/1‑143 (MQ=255)
                                                                                                          ggTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGt                           >  1:1013100/1‑143 (MQ=255)
                                                                                                            tggaaaaCAAGAGAACGTGTACGACTACACAACTACGCTGGACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAGCTATTGAAGGTAGCAAAGGTAAAGTTc                         <  2:939263/140‑1 (MQ=255)
                                                                                                              aaaaaCAAGCGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAg                       <  2:498667/143‑1 (MQ=255)
                                                                                                               aacaCAAGAGAACGTTTACTAGTACAGAACGACGCGGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGGGCCGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAgc                      <  1:303545/140‑1 (MQ=255)
                                                                                                                 aaCAAGAGGACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAgctg                    <  1:202749/143‑1 (MQ=255)
                                                                                                                 aaCAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATGGAAGGTAGCAAAGGTAAAGttgcggcg                    >  2:729216/1‑137 (MQ=255)
                                                                                                                    aaGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACACACTCGAAAATATTGAAGGTAGCAAAGGTAAAGTTCAGATGCTc                 >  1:733440/1‑143 (MQ=255)
                                                                                                                    aaGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACCATTGACAGTAGCAAAGATATAGCTCAgctgcgc                 >  2:125333/1‑141 (MQ=255)
                                                                                                                           cGTTTACTACTACACAGCGACGCTGAACGTAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCg          <  2:1224058/143‑1 (MQ=255)
                                                                                                                              ttACTGCTACACTGCTACGCGGAACGAAACCTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGtt       <  2:918103/143‑1 (MQ=255)
                                                                                                                                  tactacACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTAt   >  2:517468/1‑143 (MQ=255)
                                                                                                                                   actacACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATc  <  1:764312/143‑1 (MQ=255)
                                                                                                                                         ||                                                                                                                                       
TCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATC  >  NC_000913/124948‑125221

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