Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 125,085 Δ1 bp 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.95.7% 192.7 / 0.0 47coding (2069/2664 nt)aceEpyruvate dehydrogenase E1 component
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base . (17/28);  minor base G (1/1);  total (18/29)
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.98e-01

TCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATC  >  NC_000913/124948‑125221
                                                                                                                                         |                                                                                                                                        
tCACAGCCACATTCAGTCGCTGAGTATCCCGAACGGCACCGCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTAc                                                                                                                                    <  1:313133/143‑1 (MQ=255)
 cacaGCTACATTTAGTCGCTGACTATCCCGAACTGTGTCTCGTACGACCCGGCTGACGCGTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACg                                                                                                                                   <  1:385419/143‑1 (MQ=255)
      ccACATTCAGTCGCTTACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGTTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGaa                                                                                                                              <  2:921285/143‑1 (MQ=255)
      ccACATTCAGTCGCTGACTATCCTGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGaa                                                                                                                              <  2:81934/143‑1 (MQ=255)
       cacaTTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAAc                                                                                                                             >  1:879119/1‑143 (MQ=255)
          aTTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGaa                                                                                                                          <  2:717039/143‑1 (MQ=255)
             cAGTCGCTGACTATCCCGAACTGTATCTCTTACGCCCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAAc                                                                                                                       <  1:31918/143‑1 (MQ=255)
                 ctcTGACTATCCCGAACCGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGTCGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTCCGCTGAACGAAAACTAcc                                                                                                                   <  1:315909/141‑1 (MQ=255)
                 cGCTGACCATCCACAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTAcc                                                                                                                   <  1:373443/143‑1 (MQ=255)
                    tGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGACCGAAAACTACcaca                                                                                                                >  2:142892/1‑143 (MQ=255)
                       ctcgccagAACTTTATCTCGTACTACCCGGCTTACGCTTACGAAGTTGCGGTCATCATGCATGACGGGCTGGAGCGTATGTACTGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGc                                                                                                             <  2:800575/136‑1 (MQ=255)
                          tCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCgg                                                                                                          <  1:66282/143‑1 (MQ=255)
                            ccGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCa                                                                                                        >  1:639428/1‑143 (MQ=255)
                                    tATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGACATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGa                                                                                                <  1:219300/143‑1 (MQ=255)
                                              gACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCAGGCAATGCCGGAAGGTGCTGAg                                                                                      >  2:428103/1‑143 (MQ=255)
                                               aCCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAgg                                                                                     >  2:34450/1‑143 (MQ=255)
                                                cccGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGGACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGa                                                                                    >  1:779711/1‑143 (MQ=255)
                                                 ccGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTg                                                                                           >  1:242839/1‑135 (MQ=255)
                                                 ccGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTg                                                                                           <  2:242839/135‑1 (MQ=255)
                                                          gcttacgAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGCAAAACAAGAGAACGTATACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAGGGAATCCga                                                                          >  1:451761/1‑142 (MQ=255)
                                                                    ttGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCGGAGGAAGGTATCCGTAAAGGTATCt                                                                >  1:611732/1‑143 (MQ=255)
                                                                       cTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACa                                                             <  2:879119/143‑1 (MQ=255)
                                                                         ggcatcatCGATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGGTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACaaa                                                           <  2:815832/141‑1 (MQ=255)
                                                                          tcatcaTGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACGCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAAc                                                          >  2:590205/1‑143 (MQ=255)
                                                                          ggatggggcaggACGGGCTGGAGCGTATGGACGGTGAAAAACAAGAGAACGTTGACTACTACACACTGCGCTGAACGAAACCTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAAc                                                          <  2:610744/136‑1 (MQ=255)
                                                                            agcatgcatgACGGTCTGGAGGGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTc                                                        <  1:417366/141‑1 (MQ=255)
                                                                              catccatgACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGa                                                      <  1:544991/143‑1 (MQ=255)
                                                                                    tgACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGCAGGTGCTCAGGAAGGTCTCCGTAAAGGTACCTACAAGCTCGAACCTAt                                                >  1:186854/1‑143 (MQ=255)
                                                                                           cTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTGACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGt                                         <  2:585687/143‑1 (MQ=255)
                                                                                            tGGAGCTCATGTACGGTGCAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTa                                        <  2:520230/143‑1 (MQ=255)
                                                                                            tGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTg                                              <  1:693814/137‑1 (MQ=255)
                                                                                            tGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTg                                              >  2:693814/1‑137 (MQ=255)
                                                                                            tGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACACGAAACTATTGAAGGTa                                        >  2:576561/1‑143 (MQ=255)
                                                                                              gAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGc                                      >  2:778892/1‑143 (MQ=255)
                                                                                               ggCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTGCGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCa                                     <  2:305574/142‑1 (MQ=255)
                                                                                                 cGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCaaa                                   <  1:283205/143‑1 (MQ=255)
                                                                                                 cGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCaaa                                   <  1:521571/143‑1 (MQ=255)
                                                                                                  gTATGTACGAGGAAAAACAAGAGAACGTGTACGGGTGGGCACTACGCTCAACGAAAACTGCCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAg                                  <  2:190979/143‑1 (MQ=255)
                                                                                                   tATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTCACGGTATCTACAAACTCGAACCGATGGAAGGTAGCagagc                                 >  1:325327/1‑142 (MQ=255)
                                                                                                       tACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTaaa                             <  1:589469/143‑1 (MQ=255)
                                                                                                        actgTGAGAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAg                            <  1:118714/140‑1 (MQ=255)
                                                                                                               aaaaCAAGAGAACGTTTACTACTACAATCTACGCTAAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTAGCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAgct                     <  1:248067/143‑1 (MQ=255)
                                                                                                                 aaCAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAgctgc                   >  2:143509/1‑143 (MQ=255)
                                                                                                                   cAAGAGAACGTTTCCTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTc                 <  1:229676/143‑1 (MQ=255)
                                                                                                                   cAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTCCAGCTGCTc                 <  1:34450/143‑1 (MQ=255)
                                                                                                                    aaGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGCACAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCCGCTGCACg                >  2:500551/1‑143 (MQ=255)
                                                                                                                                  tactacACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGTTTCTATc  >  1:257564/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 < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: