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,0850TC94.1% 87.8 / ‑6.2 34I690T (ATC→ACA) aceEpyruvate dehydrogenase E1 component
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (15/17);  total (16/18)
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
*NC_000913125,0860CA97.0% 79.3 / ‑4.7 34I690T (ATC→ACAaceEpyruvate dehydrogenase E1 component
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (14/18);  total (16/18)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.69e-01

TCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATC  >  NC_000913/124948‑125221
                                                                                                                                         ||                                                                                                                                       
tCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTGACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTa                                                                                                                                     >  1:658443/1‑143 (MQ=255)
       cacaTTCAGTCGCGTATTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGaa                                                                                                                              <  2:270603/143‑1 (MQ=255)
          aTTCAGTCGCTGACTATCCCGAACTGTATCTAATCGGCTGCAATGAACGCTTACGAAGTTGCTGTCACCATGTATTCCCGTCCGCTCCGTAGTTGCGGTTAATATCATTATCACGTTTACTACTACACAACTACGCGTAACGa                                                                                                                           >  2:439235/1‑143 (MQ=255)
                  gCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTAcc                                                                                                                   >  1:272359/1‑143 (MQ=255)
                             cGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATCCCCGCa                                                                                                        >  2:650483/1‑143 (MQ=255)
                                              gACCCGGCTTAAGCGTTCAAATTGGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGa                                                                                       <  2:470273/143‑1 (MQ=255)
                                                cccGGCTTACGCTGACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAACGCCGGAAGGTGCCGAgg                                                                                     >  2:678590/1‑143 (MQ=255)
                                                  cGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGaa                                                                                   <  2:642820/143‑1 (MQ=255)
                                                  cGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGGATGTACGGAGAAAAACAAGAGAACGCTTACTACTACACGACTACGCTGAACGAAAAGTACCACATGCCAGCCATGCAGGAAGgcgcaacggaa                                                                                   >  2:359798/1‑136 (MQ=255)
                                                            ttacgAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGACGGTACCCgca                                                                         >  2:27361/1‑141 (MQ=255)
                                                               cgAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAAAACGTTCACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAg                                                                      <  1:369797/143‑1 (MQ=255)
                                                               cgAAGTTGCTGGCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTGCGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAg                                                                      <  2:138266/143‑1 (MQ=255)
                                                                gAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTGCACAACTACGCTGAACGAAAACTACCACATGCCGGCGACGCCGGAAGGAGCTGAGGAAGGCATCCGCaaacg                                                                     >  2:521677/1‑141 (MQ=255)
                                                                          tcatcaTGCATGACGGTCTGGAGCGTATGGACCGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGCCGGAGATCAGCAAGGTATCCGTAAAGGTATCTACaaa                                                           <  2:441662/143‑1 (MQ=255)
                                                                                        ggTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGTAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGa                                             <  1:391302/143‑1 (MQ=255)
                                                                                        ggTCTGGACCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGa                                             <  1:43584/143‑1 (MQ=255)
                                                                                         gTCTGGGGCGTAGGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGaa                                            <  2:505308/143‑1 (MQ=255)
                                                                                          tCTGGATCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAg                                           <  2:199121/143‑1 (MQ=255)
                                                                                           cTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAggg                                          >  1:402396/1‑143 (MQ=255)
                                                                                            tGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGt                                         >  1:245073/1‑143 (MQ=255)
                                                                                                     tGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACGACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGt                                <  2:92674/143‑1 (MQ=255)
                                                                                                       tACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGAGCTGAGGACGGTATCCGTAAAGGTATCTAAAAAGTCGACACTATTGAAGGTCGCAACGGTaa                              >  2:474295/1‑143 (MQ=255)
                                                                                                        aCGGTTAAAACGAAGAGAACGCGTACTACTACACAACTACGCTGAAGGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTaaa                             <  2:206500/143‑1 (MQ=255)
                                                                                                        aCGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTaaa                             <  1:439235/143‑1 (MQ=255)
                                                                                                        aCGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTaaa                             <  2:84955/143‑1 (MQ=255)
                                                                                                        aCGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACAAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAGAGGTATCTACAAACTCCAAACAATTGAAGGGAACAAAAGTaaa                             >  2:444996/1‑143 (MQ=255)
                                                                                                            tGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGCGGAAGCTATCCGCAAAGGTATCTACAAACTCGACACTATTCAAGGTAGCAAAGGtcagctcc                         >  2:193475/1‑136 (MQ=255)
                                                                                                            tGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGCAAAGGTATCTACAAACACGAAACAATTGAAGGTAGCAAAGGTAAACTTc                         >  2:644113/1‑143 (MQ=255)
                                                                                                                   cAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAgctgct                  >  2:298612/1‑143 (MQ=255)
                                                                                                                     agagAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGTCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCg                >  2:367577/1‑143 (MQ=255)
                                                                                                                      gagaACGTATACAACTACATAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCgg               <  1:96925/143‑1 (MQ=255)
                                                                                                                           cGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCg          <  1:542814/143‑1 (MQ=255)
                                                                                                                           cGTTTACTACTACACAACTACGCCGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCg          <  1:309806/143‑1 (MQ=255)
                                                                                                                                   actacACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATc  <  1:566265/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: