Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 125,387 C→T Q791* (CAG→TAG)  aceE → pyruvate dehydrogenase E1 component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913125,3870CT97.4% 108.4 / NA 40Q791* (CAG→TAG) aceEpyruvate dehydrogenase E1 component
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (20/18);  total (22/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

AAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGAGA  >  NC_000913/125258‑125525
                                                                                                                                 |                                                                                                                                          
tgaGATTGCGGCGTAGGTGCTGACGTGTATAGCGTGACCGCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAacga                                                                                                                               <  2:369122/141‑1 (MQ=255)
  aGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACATCCTTCACCGACCTCGCGCGGGTGGGTCGGGAAAGAGAACGCCGGAACATGCTGCACCCGCCGGAAACTCCGCGCGGCCCGTATAGCGCTCAGGCGATCAacgagc                                                                                                                             >  2:555118/1‑141 (MQ=255)
     ttACGGCGTAGGTTCTGACGTTTATAGCGTGAACTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCGGCGCGTTCCGTATATCGCGTAGGTGATGAACGACGCTc                                                                                                                          >  1:11935/1‑143 (MQ=255)
       aCGGCGTAGGTTCTGACGTTTATGGCGTGACCTCCTTCACCGAGCTGGCGCGGGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCg                                                                                                                        <  2:758752/143‑1 (MQ=255)
        cGGCGTAGGTTCTGACGGTTAGAGCGTGACCTCCGTCACCGAGCTGGCGCGCGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGCAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCgg                                                                                                                       <  2:56733/143‑1 (MQ=255)
         ggCGTAGGTTCTGACGTTTATAGCGTGACATCCTGCACCGAGCTAGCGCGTGATGGTCAGGATTGTGAAAGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGGGTATCGCTTGGGTGATGAACGGCGCCCCCGc                                                                                                                      >  1:719349/1‑143 (MQ=255)
                     gACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTAc                                                                                                          >  1:645953/1‑143 (MQ=255)
                            ataGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACtct                                                                                                   >  1:164824/1‑141 (MQ=255)
                            agaGCGTGCCCTCCCGCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAATCTCCTCTCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACtat                                                                                                   <  1:701076/141‑1 (MQ=255)
                                      ctcctTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGt                                                                                         >  1:263677/1‑143 (MQ=255)
                                      ctcctTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCCCCGGCAGTGGCGCCTACCGACTATATGAAACTGt                                                                                         >  2:112553/1‑143 (MQ=255)
                                      ctcctTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCCACGACGGTGGCATATACCGACTATATGAAACTga                                                                                         >  2:46501/1‑142 (MQ=255)
                                       tcctTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGGATATCGCTTAGGTGATGAACGACACTCCGGCAGTGGCATCTACCCACTATATGAAACTGtt                                                                                        >  1:698485/1‑143 (MQ=255)
                                        cctTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTc                                                                                       <  1:79296/143‑1 (MQ=255)
                                        cctCCAGCGAGCGGACGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGAGGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTc                                                                                       <  1:752196/143‑1 (MQ=255)
                                               cGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGc                                                                                <  1:203534/143‑1 (MQ=255)
                                               cGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGGGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGc                                                                                >  2:739569/1‑143 (MQ=255)
                                                 tgCTGGCGCGTGATGGTCAGGATTGGGGGCGCTGGAACATGCTGCACCCGCTGGAAACGCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAg                                                                              <  2:152151/142‑1 (MQ=255)
                                                   cTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCAGGCAGTAGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGt                                                                            >  2:566730/1‑143 (MQ=255)
                                                     cgcgcGTGATGGCCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTcc                                                                          <  1:540551/142‑1 (MQ=255)
                                                           tGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGCGCAGGTCCGTACtt                                                                    >  2:180478/1‑143 (MQ=255)
                                                           tGATGGGCAGCAGTGTGAACGCTGGAACAGGCCGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACtt                                                                    <  1:625572/143‑1 (MQ=255)
                                                           tGACGGTCCGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACtt                                                                    <  1:56794/143‑1 (MQ=255)
                                                                  cAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGtga                                                                                                                                      >  1:248647/1‑70 (MQ=255)
                                                                  cAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGtga                                                                                                                                      <  2:248647/70‑1 (MQ=255)
                                                                     gATTGTGAACGCTGGAACATGCTGCACCCGCAGGAAACTCCGCGCGTTCCGTATAACGATAAGCAGATGAGCGCCGCTCCGGCAGTGGCATCTACCGTCTATATAAAACTGTTCGCTGAGCAGGCCCGTACTAACGTAGCGGc                                                          >  2:514621/1‑143 (MQ=255)
                                                                            aaCGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTCCTTACGGACCGGCTgacgcc                                                   >  1:466044/1‑141 (MQ=255)
                                                                               gCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTAc                                                <  1:499812/143‑1 (MQ=255)
                                                                                          ctgcCCCCGCTGGAAACGCCGCGCGTGCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTggg                                     <  1:482966/143‑1 (MQ=255)
                                                                                           tgcACCCGCTGGAAACTCCGCGCGTTCCGTATATAGCGTAGGTGATGAACGACGCTCCGGCAGTGGCAACTACCGACTATCTGAGACTGTACGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGCCCACCGCGTACTAGGt                                    >  2:272306/1‑143 (MQ=255)
                                                                                                 ggCTGGAAACGCCGCGCGTGCCGTGTATCGCTTAGGTGATGAACGACGGTCCGGCAGTGGCATCTACCGACTAGATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGAt                              <  2:9080/142‑1 (MQ=255)
                                                                                                 cGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACCGTTCGCTGAACAGGTCCGTACTTACGTACCGGATGACGGCTACCGGGTACTGGGTACTGAt                              >  2:218500/1‑143 (MQ=255)
                                                                                                      gAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGCACTGATGGCtt                         >  2:709378/1‑143 (MQ=255)
                                                                                                             cgcgcgTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTCCGTACCGGCTGCCGACTACCGCGTACTGGGCACCGAGGGCTTCGGTCGt                  >  2:331545/1‑143 (MQ=255)
                                                                                                                   ttCCGTATATCGCTTAGGTGATGAACGACGCTCCGACAGTGGCAGCGGCCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTGGCGCGTGCTGGGGACTGATGGCTGCGGTGGGTCCGAc            >  2:55417/1‑143 (MQ=255)
                                                                                                                     ccGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAg          <  1:431696/143‑1 (MQ=255)
                                                                                                                     ccGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACCTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAg          <  1:657808/143‑1 (MQ=255)
                                                                                                                      cGTATACCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGc         <  1:509470/143‑1 (MQ=255)
                                                                                                                            tCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGCCAGCCGTgag   >  1:232310/1‑143 (MQ=255)
                                                                                                                             cGCTTAGGGGATGAACGACGCTCCGGCAGTGGCGTCTACCGATTATATGAAACTGTGCGCTGGGCCGGTCCGTACTTACGTACAGGCTGACGACTACCGCGTACTGGGTGCTGATGGCTTCGGGCGTTCCGACAGCCgaggga  >  2:417118/1‑138 (MQ=255)
                                                                                                                                 |                                                                                                                                          
AAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGAGA  >  NC_000913/125258‑125525

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