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.93.5% 111.0 / 5.9 33coding (2069/2664 nt)aceEpyruvate dehydrogenase E1 component
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base . (16/13);  total (16/17)
Fisher's exact test for biased strand distribution p-value = 2.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTT  >  NC_000913/124964‑125216
                                                                                                                         |                                                                                                                                   
tCGGTGACTATCCAGAACTGTATCTTTTAAGACTCGGCTTACGTTTATGAAGTTGCTGTTATCATGTATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTAc                                                                                                               <  1:21871/143‑1 (MQ=255)
             cGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCggaca                                                                                                  >  2:267358/1‑140 (MQ=255)
             cGAACTGTATCTCTTACGACCCGGCTTACGCATACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGCACGTTTACTACTACACACTACGCTGAACGCAAAATAACAGAAGCAGGCaa                                                                                                  >  2:322713/1‑143 (MQ=255)
               aaCTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCTTGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATg                                                                                                <  1:459504/143‑1 (MQ=255)
                 cTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGCACGTATGTACGGTGAAAAACAAGAGAACGCTTACTACTACACAATACGCTGAACGGCCGCTACCAAATGCCCGCGACGcc                                                                                              >  2:49058/1‑143 (MQ=255)
                    tATCTCTTACGACCCGGCTTACTTTTACTAAGTTGCTGTCATCATGCATGACGGTCTGGAGAGTATGTACGGTGAAAAACAAGAGAACGTTGACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGa                                                                                           <  1:801952/143‑1 (MQ=255)
                    tATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGa                                                                                           >  1:651471/1‑143 (MQ=255)
                     aTCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGATCGTATGTACTTTTATTTTCTAGTTAACGTTTACTACTACACACGACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGaa                                                                                          <  1:666795/143‑1 (MQ=255)
                      tctcTTATTACCTGGCTTACGCTGAAGAAGTGGCTGTCAGCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAg                                                                                         <  1:320221/143‑1 (MQ=255)
                                cccGGCTTACGCTTACGAAGTTGCTGTCATCATTCATTACGTTCTTTATTTTATGTACGGTGAAAAACAAGAAAACGTGAACTACAACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGa                                                                               <  1:26998/143‑1 (MQ=255)
                                                   gTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGCTCAACAACAAGAGTACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTGAAGGTATc                                                            <  2:484505/143‑1 (MQ=255)
                                                        tGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAACGCAGGAAGGCGCCGAGGGAGGCGTCCGTGAAGGTATCcacac                                                       >  2:366152/1‑142 (MQ=255)
                                                         gTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGAGCAGAGAAAAGTATAAGTAAAGGCAACTACaaa                                                      >  1:725238/1‑143 (MQ=255)
                                                             tcatGCCTGTCGGTCTGGAGCGTATGTACTGTGACATACAAGAGAACGTTGACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCg                                                  <  1:567873/143‑1 (MQ=255)
                                                                tgcatgGCGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACGCGCTCTGCTGTACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGaaa                                               <  2:236877/143‑1 (MQ=255)
                                                                tgcatgACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACTCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGaaa                                               <  1:483033/143‑1 (MQ=255)
                                                                    tgACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTTAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTAt                                           <  1:267358/143‑1 (MQ=255)
                                                                         gTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGCAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAg                                      <  1:416013/143‑1 (MQ=255)
                                                                         gTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAg                                      <  2:651471/143‑1 (MQ=255)
                                                                                 cGTATTTACGGTGAAATACAAGAGATTGTTTACTTCTTTTTTTTTCTCTGAACTAAAACTACCACGAGCAGGCAATGCCGGAAGGTGCGGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCaaa                              <  1:55036/143‑1 (MQ=255)
                                                                                         cGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAACACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGAATCCGTAAAGGTATCTACAAACACGAAACTATTGAAGGTAGCAAAGGTAAAGt                      >  1:92087/1‑143 (MQ=255)
                                                                                          ggTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACACGAAACGATGGAAGGTAGCAAAGGAAAAGtt                     >  2:245215/1‑143 (MQ=255)
                                                                                          ggTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGATGAGAAAGGTATCCGGAAAGGTATCTACAAACTCGAAACTACTGAAGGTAGCAAAAGCAAAGtt                     >  2:164670/1‑143 (MQ=255)
                                                                                            tGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAGACTATTGAAGGACGCAAAGGTAAAGtcaa                   >  2:331205/1‑140 (MQ=255)
                                                                                              aaaaaCAAGAGAACGTTTACTACTACACACTACGCTGCACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATAGAAGGTAGCAAAGGTAAAGTTAAgc                 >  1:470359/1‑143 (MQ=255)
                                                                                              aaaaaCAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGCATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAgc                 >  1:589834/1‑143 (MQ=255)
                                                                                              aaaaaCAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCAAATGCCGGCAATGCCGGAAGGTCCTGAGGAAGATATCCGTAAAGGTATCTACCAACTCGAAACTCTTGCAGGTAGCAAACGGAAAGTTCAgc                 >  2:800092/1‑143 (MQ=255)
                                                                                               aaaaCAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGCGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACAATAGAAGGTAGCAAAGGAAAAGTTAAgct                >  1:669342/1‑143 (MQ=255)
                                                                                                     agagAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCTGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCgg          <  2:357735/143‑1 (MQ=255)
                                                                                                           cGTTTACTACTACACACTACGATGAACGAAAACTACAACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAAATATTGAAGGTAGCACAAGTACAGTTCAGCTGCTCGACTACgg    >  2:373839/1‑143 (MQ=255)
                                                                                                            gTTTCCTACTACCCACTCAGTTGAGCGAAAACTACCACCGGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGt   <  2:176019/143‑1 (MQ=255)
                                                                                                             tttACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGCAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACCAACTCGAAACTAGCGAAGGTAGCAAAGGTAAAGTTCACCTGCTCGGCTCCGGtt  >  2:533908/1‑143 (MQ=255)
                                                                                                             gttACTAATACACATTACGTAGAAAGAAAACTAACACAGGCCGGCAAGGGCGGAAGGGGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGtt  <  1:534470/142‑1 (MQ=255)
                                                                                                                         |                                                                                                                                   
TCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTT  >  NC_000913/124964‑125216

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