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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913124,2410CT94.7% 140.9 / ‑5.4 59Q409* (CAG→TAG) aceEpyruvate dehydrogenase E1 component
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (34/20);  total (36/23)
Fisher's exact test for biased strand distribution p-value = 5.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGGTCACGATCCGAAGAAAATCTACGCTGCATTCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACCAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACATCACCTTCCCGGAAGGTTCTGAAGAGCATACCTATCTGC  >  NC_000913/124104‑124373
                                                                                                                                         |                                                                                                                                    
gtggtCACGATCCGAAGAAAATCTACGCTGCATTCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAgggt                                                                                                                                 >  1:217069/1‑141 (MQ=255)
    tCACGATCCGAAGAAAATCTACGCTGCATTCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAATGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGa                                                                                                                             <  1:651100/143‑1 (MQ=255)
      aCGATCCGAAGAAAATCTACGCTGCATTCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTATTCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGGTAAGaaa                                                                                                                           >  2:1206780/1‑143 (MQ=255)
      aCGATCCGAAGAAAATCTACGCTGCATTCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGaaa                                                                                                                           >  2:721433/1‑143 (MQ=255)
       cgttCCGAAAAAAATCTACGCTGCATCCAATAAATCGCTGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGATTACGGCATGGGCGACGCGGCTGAATGTAAAAACATCGCGCACTAGGTTAAGaaaa                                                                                                                          <  2:375197/140‑1 (MQ=255)
       cGATCCGAAGAAAATCTACGCTGCAGTCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGaaaa                                                                                                                          <  1:1121539/143‑1 (MQ=255)
               agaaAATCTACGCTGCATTCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATg                                                                                                                  >  1:424741/1‑143 (MQ=255)
                  aaaTCTACGCTGCATTCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTACGTTAAGAAAATGAACATGGAc                                                                                                               >  2:323208/1‑143 (MQ=255)
                    aTCTACGCTGCATTCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCGCTAGGTTAAGAAAATGAACATGGACgg                                                                                                             >  1:330817/1‑143 (MQ=255)
                            tGCATTCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGta                                                                                                     >  2:352668/1‑142 (MQ=255)
                            tGCATTCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTc                                                                                                     >  1:538842/1‑143 (MQ=255)
                            tGCATTCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTc                                                                                                     >  2:399315/1‑143 (MQ=255)
                            tGCATTCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAACGAACATGGACGGTGTGCGTc                                                                                                     >  2:353308/1‑143 (MQ=255)
                                ttCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCatat                                                                                                 >  1:1217497/1‑143 (MQ=255)
                                   aaGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAGGAAAATGAACATGGAGGGTGTGCGTCATATcgg                                                                                              >  1:530749/1‑141 (MQ=255)
                                    aGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGCCataacagc                                                                                             >  2:676648/1‑138 (MQ=255)
                                          cgcAGGAAACCAAAGGCAAAGCGACAGTAATGCGAGCGCAGACCATTAAGGGTTACAGCACGAGCGCCTACGCTGAAGGTCAAAACATCGTGCACCAGGGAAAGAAGATGAAAAAGCACGACGTGCGGAATATCGGCGCCCga                                                                                       >  2:1193693/1‑142 (MQ=255)
                                           gcAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTCGGTTTCGAACTTGAACATGGACGGTGTGCTTCATATCCGCGACCGtt                                                                                      >  1:44220/1‑143 (MQ=255)
                                            taGGTAACCAAAGGAAAGGCGACAGTAATCCGTGCTCAGACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGCCACATGCAGCCACTTTGACGGTGTGCGTCATATCCGCGACCGttt                                                                                     <  1:1177524/142‑1 (MQ=255)
                                            cAGGGAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGttt                                                                                     <  2:486569/143‑1 (MQ=255)
                                              ggAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGCACATGGACGGTGTGCGTCATATCCGCGACCGTTTCa                                                                                   >  1:483285/1‑143 (MQ=255)
                                               gAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGCCGACGCGACTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTcca                                                                                  >  1:840699/1‑141 (MQ=255)
                                                aaaCCAACGGCAAATCGACAGTAATCCTGGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAAt                                                                                 <  2:760547/143‑1 (MQ=255)
                                                   ccAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATAGCCGCGACCGGTTCAAtgtg                                                                              >  1:1073207/1‑143 (MQ=255)
                                                    cAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGACGGTAAAAGCATCGCGCACTAGGGTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGc                                                                             >  1:1106474/1‑143 (MQ=255)
                                                    cAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGc                                                                             >  1:908497/1‑143 (MQ=255)
                                                    cAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCAGCGACCGCTTCAATGTGc                                                                             >  1:1183410/1‑143 (MQ=255)
                                                     aaaGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGcc                                                                            >  1:784690/1‑143 (MQ=255)
                                                        ggCAAAGCGACACTAATCCTTGCTCATACCATTCAAGGCTACGGCATGCGCGACGCGGCTCAAGGTAAAAACACCGCGCACCAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGTCCGTTTCAATGTGCCGgt                                                                         <  1:678243/143‑1 (MQ=255)
                                                           aaaGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCACATCAATGTTCCAGTGTc                                                                      >  1:766188/1‑143 (MQ=255)
                                                           aaaGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGGTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTc                                                                      >  1:939236/1‑143 (MQ=255)
                                                            gaccGACAGCGATCCTTGCTCAGACCGTTAAGGGTTCCGGCATGGGCGACGCGGCTGAAGGTAAGAACATCGCGCACTAGGTTAAGAGAATGAACATGGCCGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCt                                                                     <  2:534064/140‑1 (MQ=255)
                                                                  cAGTAATCCTTGCTCATACCCTTAAACGTCACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCCCTAGGTTACGAAAATGAACATGGACGGTTTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCa                                                               <  1:198583/143‑1 (MQ=255)
                                                                               tCATACCATTAAAGGTTACGGCCTGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGATAATGAACATGGACGGTGGGCGTCCTACCCGAGACCTCGTCCAGGTCCCGGGAGCTGCTGCAGAGATCGCAAAAc                                                  >  2:405851/1‑143 (MQ=255)
                                                                               tCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGCAGGTAAAAACATCGCGCACTAGGTTAAGAAAAGGAACATCGACGGTGTGAGTCATATCCGCGACCGTATCAATGTGCCGGGGTGTGATGCAGATATCGAAAGAc                                                  >  2:410987/1‑143 (MQ=255)
                                                                               tCATACCATGAAAGGTTACGGCATGGGGGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAAc                                                  <  2:222853/143‑1 (MQ=255)
                                                                                   aCCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATCTCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGcc                                              >  1:809930/1‑143 (MQ=255)
                                                                                   aCCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAGACTccc                                              >  1:218821/1‑143 (MQ=255)
                                                                                   aCCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGACGATATCCAAAAACTGcc                                              >  2:521751/1‑143 (MQ=255)
                                                                                    ccATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCg                                             <  2:45943/143‑1 (MQ=255)
                                                                                    ccATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCg                                             <  2:424741/143‑1 (MQ=255)
                                                                                       ttAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTAc                                          <  1:1120685/143‑1 (MQ=255)
                                                                                             gTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGGTATCGAGAAACTGCCGTCCATCAcc                                    >  1:228121/1‑143 (MQ=255)
                                                                                              ttACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGGTAAGAAAAGGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATTTGCCGGTGTCTGATGCAGATGTCGAAAACCTGCCGTACCTCCCCt                                   >  1:775789/1‑143 (MQ=255)
                                                                                                   gCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACATCACCTTCCCg                              >  1:1205573/1‑143 (MQ=255)
                                                                                                      tGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACATCACCTTCCCGGaa                           <  1:399315/143‑1 (MQ=255)
                                                                                                      cgggCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACATCACCTTCCCGGaa                           <  2:1153902/142‑1 (MQ=255)
                                                                                                           gACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACATCACCTTCCCGGAAGGTTc                      >  1:865378/1‑143 (MQ=255)
                                                                                                            aCGCGGTTGAAGGTAGAAACGTGGCGCACGAGGTTAAGAAAATGAACATGGACGGGGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAACACTGCCGTACATCACCTTCCCGGAAGGTTCt                     <  2:330856/143‑1 (MQ=255)
                                                                                                             ggcggcgGAAGGGAACAAAGTCGCGCAGCAGGGGAAGAAAATGACCATGGACGGTGTGCGTCATATCCGCGACCGTTTCAAGGGGCCGGTGTCGGATGCAGATATCGAAAAACCGCCGTACATCACCTTCCCGGAAGGTTCTg                    <  1:151836/142‑1 (MQ=255)
                                                                                                               cgGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGGAAAACTGCCGTACATCACCTTCCCGGAAGGTTCTGaa                  >  2:846326/1‑143 (MQ=255)
                                                                                                                 gCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACATCACCTTCCCGGAAGGTTCTGAaga                <  2:330817/143‑1 (MQ=255)
                                                                                                                   tGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTTCCGTACATCACCTTCCCGGATGGTTCTGacgtgc              >  1:24710/1‑138 (MQ=255)
                                                                                                                      aGGTAAACACATCGCGCACTAGGTTAAGAAAATGAACCTGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACATCACCTTCCCGGAAGGTTCTGAAGAGCATa           <  1:472074/143‑1 (MQ=255)
                                                                                                                       cttAAAAGCATCGCGAACTAGGTTAAGAAAAGGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACATCACCTTCCCGGAAGGTTCTGAAGAGCATAc          <  1:1089302/141‑1 (MQ=255)
                                                                                                                         tAAAAACATCGCGAACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCTGGTGTCTGATGCAGATATCGAAAAACTGCCGTACATCACCTTCCCGGAAGGTTCTGAAGAGCATACCt        <  1:173052/143‑1 (MQ=255)
                                                                                                                          aaaaaCATCGCGCACAAGGTTAAGAAAATGAACATGGACGGTGTGCGTAATATCCGCGACCGTTTCACTGTGCCGGTGTCTGATGCAGATATCGGAAAACTGCCGTGCATCACCTTCCCGCAAGGCTATGACGAGCagcgagc       >  2:20976/1‑137 (MQ=255)
                                                                                                                             aaCATCGCGCACTCGGTTAACAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTAAATCTCCGTCCCGGAAGGTTCTGAAGAGCATACCTATCt    <  2:380660/143‑1 (MQ=255)
                                                                                                                               ccgcgcgCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGGCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACATCACCTTCCCGGAAGGTTCTGAAGAGCATACCTATCTGc  <  2:851780/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                    
GTGGTCACGATCCGAAGAAAATCTACGCTGCATTCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACCAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACATCACCTTCCCGGAAGGTTCTGAAGAGCATACCTATCTGC  >  NC_000913/124104‑124373

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: