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.97.5% 163.0 / ‑2.1 41coding (2069/2664 nt)aceEpyruvate dehydrogenase E1 component
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base . (20/19);  total (21/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTT  >  NC_000913/124957‑125216
                                                                                                                                |                                                                                                                                   
caTTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGa                                                                                                                      >  1:548763/1‑143 (MQ=255)
  ttCAGTCGCTGACCATCCCGAACTGCATCCCTTACGAACCAGATTACTCTGACGGAGTTGCTGTCATCCTGCAGAATGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGaaa                                                                                                                    <  2:80000/143‑1 (MQ=255)
        cGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGACCGAAAGCGAcc                                                                                                              >  1:94361/1‑143 (MQ=255)
         gCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACca                                                                                                             <  1:239766/143‑1 (MQ=255)
                 tCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCgg                                                                                                     >  1:524419/1‑143 (MQ=255)
                       tctgggtctcTTTCTACCCGGCTTATGCTTTCGATTTTTCTGTCATCATGCATGTCGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGc                                                                                               <  2:285826/137‑1 (MQ=255)
                        cTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTATGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATACCGACAATGcc                                                                                              >  2:444080/1‑143 (MQ=255)
                         tGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTAGGCTCGACGAAAACTACCACATGCCGGCAAGGcag                                                                                             >  2:623429/1‑141 (MQ=255)
                             tctcTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTCTGTACGGTGCAAAACAAGCGAACGTTTACTACTACCCGCTCCGCTGAACGAAAAATACCACATGCCGGCAATGCCGGAAg                                                                                         >  1:119847/1‑143 (MQ=255)
                             gctcgTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAg                                                                                         <  1:356372/142‑1 (MQ=255)
                                    ccaCCCGGCTTACGCTGACGAAGTCGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGa                                                                                  <  2:167946/141‑1 (MQ=255)
                                       cccGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGACAATGACGGAAGAAGCTGAgaa                                                                               >  1:70630/1‑141 (MQ=255)
                                        gcGGCGTACGGTGACGAAGGTGCTGTCATCAGGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGaa                                                                              <  2:471675/142‑1 (MQ=255)
                                         cGGCTTGCGCTGAAGAAGTTGCGGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAg                                                                             <  2:9546/143‑1 (MQ=255)
                                         cGGCTGACGCGGACGAAGTTGCTGTCAGCATGCATGGCGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAg                                                                             <  2:28000/143‑1 (MQ=255)
                                          ggCTTACGCGTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAgg                                                                            <  1:22002/143‑1 (MQ=255)
                                                    tacgAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTaaa                                                                  <  1:720273/143‑1 (MQ=255)
                                                    tacgAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTaaa                                                                  <  1:632630/143‑1 (MQ=255)
                                                     acgAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAg                                                                 <  2:524419/143‑1 (MQ=255)
                                                       gAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGt                                                               >  2:16973/1‑143 (MQ=255)
                                                                gTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGGGCTGAGGAAGGTATCCGTAAAGGTATCTACaga                                                      >  2:213260/1‑141 (MQ=255)
                                                                   atctttcatcACTGTCTTGATCGTATGTACGGTGAAAAACACAAGCACGTTTATTACTATACATTACGCTTAACGCTAACTACCACATGCCGGCAATGTCTGAAGGTGCTGAGGAATGTATCCTTTAAGGTATCTACAAACTc                                                   <  2:178984/137‑1 (MQ=255)
                                                                      atgcatgACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCACTGCCGGACGGTGATTAAAAAGAAATCAGTAAAGGTATCTAAAAACTCGaa                                                >  1:299579/1‑143 (MQ=255)
                                                                         catgACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAAcc                                             >  2:487420/1‑142 (MQ=255)
                                                                             aCGGTCTGGAGCGTATGTAAGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGGGATGAGGAAGGAAACCGTAAAGGTATCTACAAAGTCGAAAcacgtg                                         >  2:206486/1‑138 (MQ=255)
                                                                               ggTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGaa                                       <  2:301555/143‑1 (MQ=255)
                                                                                   tGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACAACAACCCGGCCATGCCGGAAGGAGCAGAGGAAGGTATCCGTAAAGGTATCTGCAAACTCGAAAACATTGAAGGTa                                   >  2:94908/1‑143 (MQ=255)
                                                                                    ggAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAg                                  <  1:152967/143‑1 (MQ=255)
                                                                                     gAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGc                                 >  2:693093/1‑143 (MQ=255)
                                                                                       gCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTAGTGAGGGGAGCaa                               >  1:90818/1‑143 (MQ=255)
                                                                                         gTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAg                             <  1:355924/143‑1 (MQ=255)
                                                                                              tACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTaaa                        <  2:14570/143‑1 (MQ=255)
                                                                                              tACGGTGAAAAACAAGAGAACGTTTACTACTACTATTTTTTTTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTaaa                        <  1:131641/143‑1 (MQ=255)
                                                                                                cGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGt                      >  1:703694/1‑143 (MQ=255)
                                                                                                  gTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAACGTAGCAAAGGTAAAgcta                    >  2:414865/1‑140 (MQ=255)
                                                                                                     aaaaaCAAGAGAACGTTTACTACTAGACACTACGCTGAACGACAGCTACCACATCCCAGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAgc                 <  1:302095/143‑1 (MQ=255)
                                                                                                       aaaCAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCGATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGCAGCCGGGGTAAAGCTCAgctg               >  2:658479/1‑143 (MQ=255)
                                                                                                           aaGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGATGGTAGCAAAGGTAAAGTTCAGCTGCTCg           >  2:165719/1‑143 (MQ=255)
                                                                                                            agagAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCgg          <  1:537019/143‑1 (MQ=255)
                                                                                                                aaCGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTAGCGAAGGTAGCAAAGGGCAAGTGCAGCTGCTCGGcccc      >  1:24663/1‑143 (MQ=255)
                                                                                                                    tttACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTCCAAACTCGAAACTATTGAAGGCAGCAAAGGCAAAGTTCAGCTGCTCGACTCCGGtt  >  2:61971/1‑143 (MQ=255)
                                                                                                                                |                                                                                                                                   
CATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTT  >  NC_000913/124957‑125216

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: