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.89.7% 100.0 / 1.7 29coding (2069/2664 nt)aceEpyruvate dehydrogenase E1 component
Reads supporting (aligned to +/- strand):  ref base T (1/0);  major base . (14/12);  minor base A (0/1);  total (16/13)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTAT  >  NC_000913/124954‑125220
                                                                                                                                   |                                                                                                                                       
ccACATTCAGTCGCTGACTATCCCGACCTGTATCTCTTACGACCCGGCTTACACTTAAGAAATTGCTATCATAATACATAACGGACTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGaa                                                                                                                             >  1:143322/1‑143 (MQ=255)
 cacaTTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTAATACACACTACGCTGAAc                                                                                                                            >  1:472985/1‑143 (MQ=255)
  acaTTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTCCTGTCCTCATACATAACGAACTGGAGCGTATGTACGGTGAAAAAAAAGAGAACGATAACTACTACACACTACGCTCaaag                                                                                                                           >  1:218680/1‑141 (MQ=255)
             cTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGAAGGTCTGGAGCGAATGTACGGAGAAAAACAAGAGAACGTTTACTACTACACACTACGCCAAACAAAAACAACcaa                                                                                                                >  2:425400/1‑142 (MQ=255)
                    tCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCCTGCATGACGGTCTGAAGCCCATGTACGGTGAAAAACAACAGAAAGTGTAATAGTTCCAACTAAGATAAACGAAAACTACCACATGccag                                                                                                         >  1:122849/1‑141 (MQ=255)
                                tctcTTACGGCCCGGCGTACGCGTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAg                                                                                             <  2:390354/143‑1 (MQ=255)
                                           ccGGCTTACCCTTACGAAGTTGCTGTCATAATGCATGACTGTCTGTAGCGTATGAACGGTGAAAAAGAAGAGAACGTCTACTACTACACACTACGCTCAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGaa                                                                                  <  1:139649/143‑1 (MQ=255)
                                            cGGCTTAGGCTGACGAAGTTGCTGTCATCATGAATGACGGGCTGGAGCGGATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAg                                                                                 <  1:329025/143‑1 (MQ=255)
                                            cGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAg                                                                                 <  1:438572/143‑1 (MQ=255)
                                                       tacgAAGTTGCTGTCATCATGCATGGCGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTaaa                                                                      <  2:329181/143‑1 (MQ=255)
                                                       tacgAAGTTGCTGTCATCATGCATGACGTTTTGGAGGGTATGTACGGTGAAAAACAAGAGAACGTTTACTAATACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTaaa                                                                      <  1:55287/143‑1 (MQ=255)
                                                          gAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGt                                                                   <  1:419642/143‑1 (MQ=255)
                                                           aaGTTGGTGTCATCATGCATGACGGTCTAGAGCTTATGTACGGTGAAATACAAGGTGACGTTTACTACTACACGCTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTa                                                                  <  2:370241/143‑1 (MQ=255)
                                                              ttGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATACGTAAAGatatag                                                               >  2:226652/1‑141 (MQ=255)
                                                                gCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGACGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTAc                                                             >  1:102044/1‑143 (MQ=255)
                                                                  tGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGGTGAGGAAGGTATCCGTAAAGGTATCTACaa                                                           >  2:461229/1‑143 (MQ=255)
                                                                   gTCAGCATGCATGGTGGTCTTGAGCGTGTGTGTGGTGAAAAACAAGATTATTGTTACTTCTACACACGAAGCTGAACGAGAACTATCACATGGCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACaaa                                                          <  2:370046/143‑1 (MQ=255)
                                                                            catgACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACt                                                 >  1:527878/1‑143 (MQ=255)
                                                                                        gAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTCTCCGTAACGGTATCTACAAACTCGAAACTAGTGAAGGTAGc                                     >  2:217424/1‑143 (MQ=255)
                                                                                        gAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAggacga                                     >  1:172242/1‑139 (MQ=255)
                                                                                          gCGTATGTACGGTGAAAAACAAGAGAGCGTTTACTACTACACGCTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCaa                                   <  2:363170/143‑1 (MQ=255)
                                                                                                 tccGGTGAAAAACAAGAGAACGTTTACTGCTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATTCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTaaa                            <  2:3388/141‑1 (MQ=255)
                                                                                                 tACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGCAGGTGCTGCGGAAGGTATCCGTAAAGGTATCTAAAAACTCGAAAATATTGAAGGTAGAAAAGGTaaa                            >  1:407803/1‑143 (MQ=255)
                                                                                                  aCGGTTAAAAACAAAAAAATGTTAACTAATAAACAATATACTGAACGAAAACTACCAAATGCCGGCAATGCAGGAAGGTGATGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAg                           <  1:48177/143‑1 (MQ=255)
                                                                                                    ggTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAAGTGAAGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGtt                         >  1:165778/1‑143 (MQ=255)
                                                                                                           aaCAAGAGAACATCTACTATTATTTATTACGATGAACGAAAACTACCACATGCCGACAATGACGGAAGGTCCTGAGGAAGGTATCCGTACAGGTATCTACAAACTAGAGACTATTGAAGGTAGCACAGGTAAAGTTCAgctcc                  >  2:223206/1‑141 (MQ=255)
                                                                                                              aaGAGAACGTTTACTACTATACTCTACGCTGAACGTAAACTGTCTTTTTTCTACAAAACAGGACGGTGCTATATAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCg               <  1:548699/143‑1 (MQ=255)
                                                                                                                 agaACGTTTACTACTACACACTACGCTGAACGAAAACTACAACATGCCGACAATGCCGCAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACCCGAAACGATTGCAGCTACCAAAGGTAAAATTCAGCTGATCGGCt            >  2:14091/1‑143 (MQ=255)
                                                                                                                           ctactacACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAAGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGCTTCTAt  >  1:455781/1‑143 (MQ=255)
                                                                                                                                   |                                                                                                                                       
CCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTAT  >  NC_000913/124954‑125220

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: