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

GTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATC  >  NC_000913/124947‑125221
                                                                                                                                          |                                                                                                                                        
ctCACAGCCACATTCAGCCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTa                                                                                                                                     <  2:259283/142‑1 (MQ=255)
 tCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTCGGTACGGTGAAAAACAAGAGAACGTTTGCTACTACACACTAc                                                                                                                                    >  2:362069/1‑143 (MQ=255)
     aGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGAGTGACCGCCCGCAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTg                                                                                                                                >  1:714894/1‑143 (MQ=255)
       ccACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGGCCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACGCACTACGCTgca                                                                                                                              >  1:967368/1‑141 (MQ=255)
             tCAGTCGCTGTCTGTCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGaaaa                                                                                                                        <  2:99739/143‑1 (MQ=255)
                gTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTa                                                                                                                     >  1:602024/1‑143 (MQ=255)
                   gcggACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACca                                                                                                                  <  1:598670/140‑1 (MQ=255)
                    cTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATAATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACcac                                                                                                                 >  1:224576/1‑143 (MQ=255)
                      gACTATCCGGAACTGGATCTCTTCCGACCCGGCCTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACAt                                                                                                               <  2:762592/143‑1 (MQ=255)
                           tCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCgg                                                                                                          >  2:218764/1‑143 (MQ=255)
                            cccGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGc                                                                                                         >  1:235079/1‑143 (MQ=255)
                              cGAACTGCATCTCATACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCaa                                                                                                       <  1:451278/143‑1 (MQ=255)
                                        cgctTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTAGCGCGTATGTACGATGAAACAGAAGAGCACGTTTACGACTACACACTACGCAGAACAAAGACTACCACATGCCGGCAATGCCGGAAgg                                                                                             >  2:147074/3‑143 (MQ=255)
                                           ttctgACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTTAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGc                                                                                          <  2:851649/139‑1 (MQ=255)
                                           ttACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCCGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGc                                                                                          <  1:379057/143‑1 (MQ=255)
                                                   cgccttgcggttacgAAGTGGCGGTCATGGTGCATGACGGGCTGGAGCGTATGGACGGTGAAAAAGAAGAGAGCGTTGACTACGACACGCTACGCTGAACGAAAACTACCAGATGCCGGCAATGCCGGAAGGTGCTGAGGAAg                                                                                  <  2:416477/140‑1 (MQ=255)
                                                       ttacgcttacgAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAACCTACAATATGGCTGTAATGCCGGAGGGTGCTGGGGCAGGGAt                                                                              >  1:662762/1‑143 (MQ=255)
                                                         acgcttacgAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATcc                                                                            >  1:660539/1‑143 (MQ=255)
                                                           gcttacgAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGCGATGAGGAAGGTATCCGt                                                                          >  2:681466/1‑143 (MQ=255)
                                                             ttacgAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTaa                                                                        >  1:319276/1‑143 (MQ=255)
                                                                 gAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTAAGGAAGGTATCCGTAAAGGt                                                                    >  1:840167/1‑143 (MQ=255)
                                                                 gAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGAGCTGAGGAAGGTATCCGTAAAGGt                                                                    >  1:849400/1‑143 (MQ=255)
                                                                   aGTTGCTGTCATCATGCATGACGGGCCGGAGCGTATGTACGGTGAAACACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTAt                                                                  <  1:158556/143‑1 (MQ=255)
                                                                      tGCTGTCAGCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTa                                                               <  1:847356/143‑1 (MQ=255)
                                                                       gCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAACGGTATCTAc                                                              >  2:957706/1‑143 (MQ=255)
                                                                          gTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGGCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACaaa                                                           <  2:17005/143‑1 (MQ=255)
                                                                            catcatGCATGACGGCCGGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACt                                                         <  1:628919/143‑1 (MQ=255)
                                                                                atgcatgACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTGTCCGTATCGGTATCTACAAATTCGaa                                                     >  1:233297/1‑143 (MQ=255)
                                                                                  gcatgACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAAAGCCGGAAGGAGCTGAGGAAGGGATCCGTAAAAGAATCTACAAAccgggaac                                                   >  2:808972/1‑136 (MQ=255)
                                                                                   tattgcGGTCTGGAGCGTATGCACGGTGAAAAACAAGAGAACGTTTACTACTACACACTGCGCTGAACGAAAACGACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACt                                                  <  2:406090/138‑1 (MQ=255)
                                                                                      gACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCCGCAATGCAGAAAGGTGCTGAGGAAGGTATCCGTAAAGGTCTCTACAAACTCGACACtgcc                                               >  2:866430/1‑140 (MQ=255)
                                                                                           tCTGGAGCGTACGTACGGTGAAAAACAAGAGAACGCTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGGATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAgg                                          <  2:801593/143‑1 (MQ=255)
                                                                                               gAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGc                                      >  1:155048/1‑143 (MQ=255)
                                                                                                 gCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCaa                                    >  2:938146/1‑143 (MQ=255)
                                                                                                   gTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAg                                  <  1:769514/143‑1 (MQ=255)
                                                                                                    tATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGACGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAGGGGTATCTACAAACGCGAAACTATGGAAGGGAGCAAAgg                                 >  2:89631/1‑143 (MQ=255)
                                                                                                           cgTGAAAAACAAGAGAACTTTTACTTCTATTCTCTTCTCTGAACAAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGtt                          <  1:410676/142‑1 (MQ=255)
                                                                                                                aaaaCAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACGCGAAACTGTTGAAGGTAGCAGAGGTAAAGTTGAgct                     >  1:422774/1‑143 (MQ=255)
                                                                                                                  aaCAAGAGAACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGCGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAgccgc                   >  2:681204/1‑143 (MQ=255)
                                                                                                                    ctagcgaACGTTGAGTAATACTCACTAAGGTTACCGAAAACCACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGCTATCCGTAAATTTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTc                 <  2:25838/138‑1 (MQ=255)
                                                                                                                         gaACGTTTACTACTACACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGCAGGTGCTGCGGAAGGTATACGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGAAAAGGTAAAGTTCAGCTGCTCGGCTc            >  1:460771/1‑143 (MQ=255)
                                                                                                                            cGTTTACTACTACACACTACAGAGAACGATAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCgg         >  2:859316/1‑143 (MQ=255)
                                                                                                                                   tactacACACTACGGTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCCAtg  >  1:700470/1‑142 (MQ=255)
                                                                                                                                          |                                                                                                                                        
GTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATC  >  NC_000913/124947‑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: