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,2410CT98.0% 142.3 / NA 49Q409* (CAG→TAG) aceEpyruvate dehydrogenase E1 component
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (18/30);  total (19/30)
Fisher's exact test for biased strand distribution p-value = 3.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGATCCGAAGAAAATCTACGCTGCATTCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACCAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACATCACCTTCCCGGAAGGTTCTGAAGAGCATACCTATCTGCACGC  >  NC_000913/124111‑124377
                                                                                                                                  |                                                                                                                                        
cgctCCGAAGCAAATCGACGCTGCATTCAAGACAGCGCAGGAAACCAAGGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGaaaa                                                                                                                              <  1:340003/140‑1 (MQ=255)
 gATCCCAAGAAAATCCACGCGGCATTCAAGAAAGCGCAGGAGACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAAt                                                                                                                             <  2:451824/143‑1 (MQ=255)
  aTCCGAAGAAAATCTACGCTGCATTCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATg                                                                                                                            >  1:318915/1‑143 (MQ=255)
    ctgaagaaAATTTACGCTGCATTCAAGAAAGCGCAGGAAACCAAAGGCTAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGaa                                                                                                                          <  2:372284/141‑1 (MQ=255)
     cGAAGAAAATCTACGCTGCATTCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAAc                                                                                                                         >  1:72243/1‑143 (MQ=255)
       aagaaAATCTACGCTGCATGCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACAt                                                                                                                       <  1:158438/143‑1 (MQ=255)
       aacaaaaTCTCCGCTGCATGCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACAt                                                                                                                       <  2:430618/140‑1 (MQ=255)
        agaaAATCTACGCTGCATTCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAGATGAACATg                                                                                                                      >  1:237714/1‑143 (MQ=255)
         gaaAATCTACGCTGCATTCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATCAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATgg                                                                                                                     <  2:453879/143‑1 (MQ=255)
             aTCTACGCTGCATTCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACgg                                                                                                                 >  1:904888/1‑143 (MQ=255)
               cTACGCTGCATTCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGgtg                                                                                                               <  1:388049/143‑1 (MQ=255)
                tACGCTGCATTCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGgtgt                                                                                                              <  1:789129/143‑1 (MQ=255)
                 aCGCTGCATTCAAGAAAGCGCAGGCAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCGGCAGGGAACAACATCGCCCACTCGGCTAAGAAAACGACCATGGACGGcgtg                                                                                                             >  2:180102/1‑143 (MQ=255)
                   gCTGCATTCAATACAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCg                                                                                                           <  2:638208/143‑1 (MQ=255)
                            aaGAAAGCGCAGGAAACCAACGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCcg                                                                                                  <  1:249940/143‑1 (MQ=255)
                                      aGGAAACCAAAGGCAAAGCGACCGCAATCCTTGCTCATACCATTAAAGGTAACGGCATGGCCGACGCGGCTGACGGTACAAACATCGCGCACCAGGTTAAGGCGGTGAACAAAGACGGTCTGCGTCATATCCGCGGGCGTTTc                                                                                        >  2:206052/1‑143 (MQ=255)
                                        gAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCaa                                                                                      >  1:775378/1‑143 (MQ=255)
                                        gAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACAGTTTCaa                                                                                      >  1:776072/1‑143 (MQ=255)
                                        gAAACCAAAGGCAAAGCGACAGGAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCaa                                                                                      >  1:697698/1‑143 (MQ=255)
                                         acaCCAAAGGCAAAGCGACAGTAATCCTTGCTCAGACCAGGAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAAGGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAAt                                                                                     <  1:958509/141‑1 (MQ=255)
                                           aCCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAAtgt                                                                                   <  1:785219/143‑1 (MQ=255)
                                               aaggcaaggcGACAGTAACCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCg                                                                               <  1:354930/143‑1 (MQ=255)
                                               aaGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCg                                                                               >  1:39958/1‑143 (MQ=255)
                                                aGGCAAAGCGACAGTAACCCTTGCTCATCCCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCgg                                                                              <  1:711476/143‑1 (MQ=255)
                                                 ggCAGAGCGACAGTAATCCGTGCGCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACTGCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGgt                                                                             <  2:264076/143‑1 (MQ=255)
                                                   cAAAGGCACAGTAATCCTTGCTCATACCATTAAAGGGTACGGCATGGGCGACGCGGCTGAAGGTAAAAACGTCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGAGCGTCATATCCGCGACCGTTTCAATGTGCCGgtgt                                                                           <  2:876471/143‑1 (MQ=255)
                                                   cAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGgtgt                                                                           <  2:39958/143‑1 (MQ=255)
                                                   cAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTCCAATGTGCCGgtgt                                                                           >  1:328727/1‑143 (MQ=255)
                                                      agggACAGTAATCCTTGCTCACACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTg                                                                        <  1:588604/140‑1 (MQ=255)
                                                        cGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAAGGTGCCGGTGTCTGAt                                                                      >  1:744071/1‑143 (MQ=255)
                                                          aCAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGc                                                                    >  2:24487/1‑143 (MQ=255)
                                                           cAGTAATACTGGCTCATACCATAAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCa                                                                   <  1:73141/143‑1 (MQ=255)
                                                                     tGCTCATACCATTAACGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGaaa                                                         <  2:560006/143‑1 (MQ=255)
                                                                     ggcgcATACCATTACAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATACCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGaaa                                                         <  1:698035/139‑1 (MQ=255)
                                                                     cgCTCATACCATTACAGGTTCCAGCATGGGCGCCGCGGCTGCAGGTACAGACATCGTGCAATTACTTAGCTAATTGTACATGGACGGTGGGCGTCCTACCCGAGCCCATTCAAATGTGTCGGTGTCTGATGCAGATATCGaaa                                                         <  2:300388/142‑1 (MQ=255)
                                                                           tACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGc                                                   >  2:68084/1‑143 (MQ=255)
                                                                                ttAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATGCGCGACCGGTTCAATGTGCCGGTGTCTGATGCGGATACCGACAACCCGCCGGAc                                              >  2:448137/1‑143 (MQ=255)
                                                                                  acgggTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACAt                                            <  2:582364/140‑1 (MQ=255)
                                                                                   aaGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACATc                                           >  1:85076/1‑143 (MQ=255)
                                                                                      gTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGGCTGGTGCAGATAGCGAAAAACTGCCGCACaccacc                                        >  1:93932/1‑143 (MQ=255)
                                                                                          cGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACATCACCTTcc                                    <  1:24487/143‑1 (MQ=255)
                                                                                           ggCATAGGCGACGAGACTGAAGGTAAAAGCATCGCGTACTAGGTGCACACAATAAACATGTGCGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAACAACTGCCGTACATCACCTTccc                                   <  2:172033/143‑1 (MQ=255)
                                                                                               tGGGCAACGCGGCTGAAGTTAAAAACAGCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACATCACCTTCCCGGaa                               <  2:145567/143‑1 (MQ=255)
                                                                                                       gcgGCTGAAGGTAAGAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACATGACCTTCCCGGGAGGGTCTGa                       >  1:917377/1‑143 (MQ=255)
                                                                                                              aaGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGAGTGCGTCATATCCGCGACCGTTTCAATGTTCCGGTGTCTGATGCGGATATCGAAAAACCGCCATCCAACACCTCCCGGGAAGGGTCTGAAGAGCac                >  2:292926/1‑142 (MQ=255)
                                                                                                                  tAAAAACATCGCGCACTAACTAAAGAAAATGAACATGGACGGTGTGCGTCAAATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAATAAATGTAGTACATCACCTTCCCGGAAGGTTCTGAAGAGCATACCt            <  2:712680/143‑1 (MQ=255)
                                                                                                                    aaaaCATCGCGCACTCGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACATCACCTTCCCGGAAGGTTCTGAAGAGCATACCTAt          <  2:318915/143‑1 (MQ=255)
                                                                                                                        cATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACATCACCTTCCCGGAAGGTTCTGAAGAGTATACCTATCTGc      <  1:813759/143‑1 (MQ=255)
                                                                                                                            gtgcACTAGGTTAAGAAAATGAGCATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACATCACCTTCCCGGAAGGTTCTGAAGAGCATACCTATCTGCACGc  <  2:728955/141‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                        
CGATCCGAAGAAAATCTACGCTGCATTCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACCAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACATCACCTTCCCGGAAGGTTCTGAAGAGCATACCTATCTGCACGC  >  NC_000913/124111‑124377

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: