Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 125,085 Δ1 bp 100% coding (2069/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,0850T.96.2% 215.8 / ‑1.8 52coding (2069/2664 nt)aceEpyruvate dehydrogenase E1 component
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base . (23/27);  minor base A (1/0);  total (25/27)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-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.

TCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATC  >  NC_000913/124948‑125221
                                                                                                                                         |                                                                                                                                        
tCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTAc                                                                                                                                    >  1:120135/1‑143 (MQ=255)
  acaGCCACATTCAGTCGCTGACTATCCCGAACTGTAGCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGc                                                                                                                                  <  1:711832/143‑1 (MQ=255)
    aGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTAGGTACGGGGAAAAGGAAGAGAACGTTTACTACTA‑AACACTAGGCTg                                                                                                                                >  2:250516/1‑143 (MQ=255)
           ttCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTAC‑GCACTACGCTGAACGaaa                                                                                                                         >  1:84253/1‑143 (MQ=255)
                tCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGGCGGTCCGGAGCGGATGTCCGGTGAAAAACACGAGAACGTTTACTACTACA‑CACTACGCTGAACGaaaccgac                                                                                                                    >  2:622444/1‑138 (MQ=255)
                            ccGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGcc                                                                                                        >  1:450862/1‑142 (MQ=255)
                             cGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCaa                                                                                                       >  1:141287/1‑143 (MQ=255)
                                    tATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCAGGa                                                                                                >  2:494954/1‑143 (MQ=255)
                                    tATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTACGTACGGTGCAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGACAACTACCACATGCAGGCAACGCCggc                                                                                                >  1:8985/1‑142 (MQ=255)
                                      gccctTACGGCGCGGCTTACGCGTACGAAGTTGCGGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGCGAACGTTTACTACTACA‑CGCTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAg                                                                                              <  2:362884/140‑1 (MQ=255)
                                           tACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCt                                                                                         <  2:133942/143‑1 (MQ=255)
                                            aCGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTATGGCGCGCATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTg                                                                                        >  1:483525/1‑143 (MQ=255)
                                              gACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTgcg                                                                                      >  1:249880/1‑141 (MQ=255)
                                              gACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAg                                                                                      >  2:459450/1‑143 (MQ=255)
                                                 ccGGCTTACGCTTACGAAGTTGCTGTCATGATGCATGACGGTCTGGAGCGTATGTACGGCGAAACACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGaa                                                                                   <  2:683072/143‑1 (MQ=255)
                                                  cGGCTTACGCTTACGAAGTTGCTGTCCTCATTCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTCCCACA‑CACTACGCTGAACGAAACCTACCACACGCCGGCAATGCCGGAAGGTCCTGAGGAAg                                                                                  >  2:514065/1‑143 (MQ=255)
                                                  cGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAg                                                                                  <  1:726737/143‑1 (MQ=255)
                                                       tccgcgaacgAGGTGGCTGTCATCATGCATGACGGTCTGGATCGTATGGACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATc                                                                             <  2:115735/136‑1 (MQ=255)
                                                                gAAGTTGGTATCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGt                                                                    <  1:605413/143‑1 (MQ=255)
                                                                  aGTTGCTGTCATCAGGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTAt                                                                  <  2:225254/143‑1 (MQ=255)
                                                                   gtagCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATc                                                                 <  1:131091/140‑1 (MQ=255)
                                                                    ttGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGACGGTGCTGAGGAAGGTATCCGTAAAGGTATCt                                                                >  2:431972/1‑143 (MQ=255)
                                                                    ttGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCt                                                                <  1:718859/143‑1 (MQ=255)
                                                                    ttGCTGTCATCATGCATGACGGTCTGGAGCGTAGGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCt                                                                <  1:172083/143‑1 (MQ=255)
                                                                       cTGTCATCATGCATGACGGGCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCGGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACa                                                             <  1:780395/143‑1 (MQ=255)
                                                                        tGTCATCATGCATGACGGTCTGGCGCGTATGTACGGTGAAAAACAAGCGAACGTTTACTACTACA‑CACTACGCTAAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACaa                                                            >  1:25357/1‑143 (MQ=255)
                                                                         gTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACaaa                                                           <  1:779911/143‑1 (MQ=255)
                                                                            ctcatGCATGACGGTCTCGAGCGTATGTACCGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTc                                                        <  2:308208/142‑1 (MQ=255)
                                                                            atcatGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAACCTc                                                        >  1:72091/1‑143 (MQ=255)
                                                                                   atgACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTCAAGGTGTCTACAAACTAGAAACTa                                                 >  1:305591/1‑143 (MQ=255)
                                                                                       tcgTCTGGAGCGTATGGACGTTGAAAAACACGAGAACCTGTGCTACTGCA‑CGCTACGCTGAACGAAAACTGCCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGa                                             <  2:452456/141‑1 (MQ=255)
                                                                                        ggTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGaa                                            >  2:685025/1‑143 (MQ=255)
                                                                                          tCTGGTGCGTATGTACGGTGAAAAACCAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAgg                                          <  2:84253/143‑1 (MQ=255)
                                                                                              gAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGc                                      <  2:120135/143‑1 (MQ=255)
                                                                                              gAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTAACCGTAAAGGTATCTGCAAACTCGAAACTATTGAAGGTAGc                                      >  2:221575/1‑143 (MQ=255)
                                                                                                 tgTATGTACGGGGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCAGTCATTGGCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATTTATTAACTCGAAACTATTGAAGGTAGCaaa                                   <  2:686261/142‑1 (MQ=255)
                                                                                                  gTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAg                                  <  1:685025/143‑1 (MQ=255)
                                                                                                  gTATGTACGGGGAAAAACAAGAGAACGTTTACTACTACA‑CACCACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAg                                  <  2:42177/143‑1 (MQ=255)
                                                                                                   tATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAGTGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCACAgg                                 >  1:784655/1‑143 (MQ=255)
                                                                                                    aTGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTCTTGGAGGTCGCAATggg                                >  1:266913/1‑142 (MQ=255)
                                                                                                       tACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAGACTCGAAAGCATTGAAGGTGGCAGAGGTaga                             >  1:455843/1‑141 (MQ=255)
                                                                                                            tGAAACGCAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCa                        <  1:474371/143‑1 (MQ=255)
                                                                                                             gAAAAACAAGAGCACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCAGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAg                       <  2:436872/143‑1 (MQ=255)
                                                                                                             gAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAg                       <  1:459450/143‑1 (MQ=255)
                                                                                                               aaaaCAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAgct                     <  1:761779/143‑1 (MQ=255)
                                                                                                               aaaaCAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAACGCCGGACGGTGCTGAGGAAGGTATCCGTACAGGTATCTACAAACTCGAAACTACTGAAGGGAGCGAAGGGAACGttaacct                     >  2:627215/1‑138 (MQ=255)
                                                                                                                 aaCAAGAGAACGTTTACTACTACA‑CACTAGGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAACGGTAAAGttgacgtgc                   >  1:258684/1‑136 (MQ=255)
                                                                                                                  aCAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAgctgct                  <  1:535826/143‑1 (MQ=255)
                                                                                                                  aCAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGAAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAgctgct                  >  1:317473/1‑143 (MQ=255)
                                                                                                                            gTTTACTAGCACA‑CCCTACGCTGAACGAAAACTGCAACAGGCCGGCAATGCCGGAAGGTGCTGAGTAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTACTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGt        <  2:826934/143‑1 (MQ=255)
                                                                                                                             tttACTACCACA‑CACGACGCCGAACGAAACATACCCCATGCCGGCAATGCCGGAAGGTGCCGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGtt       <  1:416669/143‑1 (MQ=255)
                                                                                                                                  tactacA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCACCGGTTCTATc  >  1:643567/1‑143 (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: ‑