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,3870CT94.1% 121.6 / ‑5.8 51Q791* (CAG→TAG) aceEpyruvate dehydrogenase E1 component
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base T (20/28);  total (22/29)
Fisher's exact test for biased strand distribution p-value = 5.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGA  >  NC_000913/125249‑125523
                                                                                                                                           |                                                                                                                                        
aTCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGACACTCCGCGCGCTCCGTCTATCGCTTAGGt                                                                                                                                      >  2:318228/1‑143 (MQ=255)
       cGAAAGTTTACGGCGTAGGTTCTGTCGTTTATAGCGTGACTTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAac                                                                                                                               <  2:30009/143‑1 (MQ=255)
       cGAAAGAGTGCGGCGTAGGTTCTGACGTTTATAGCGTCACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAac                                                                                                                               <  1:724483/143‑1 (MQ=255)
           agcgtACGGCGTAGGTGCTGACGTTTATAGCGTGACCTCCTGCACCGAGCGGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGTTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAacgacg                                                                                                                           <  2:569408/139‑1 (MQ=255)
                  ggCGTGGGTTCTGACGTTTATAGCGTGACCTCCTGCACCGAGCTGGCGCGTGAGGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGc                                                                                                                    <  2:178166/143‑1 (MQ=255)
                      tAGGTTCTGACGTTTATAGCGTGCCCTACTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTg                                                                                                                <  2:446563/143‑1 (MQ=255)
                      ggggTTCTGACGTTTATCGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGCATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTATGTGATGAACGACGCTCCGGCAGTg                                                                                                                <  2:72382/141‑1 (MQ=255)
                                  tttATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGAc                                                                                                    >  1:894475/1‑143 (MQ=255)
                                    tataGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACtg                                                                                                  >  1:903174/1‑142 (MQ=255)
                                    tataGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATGGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTGTATCGCTTACGTGATGAACGACGCTCCGGCCGTGGCATCTGCCGCCta                                                                                                  >  1:229021/1‑143 (MQ=255)
                                    tataGCGCGGCCCCCTGCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACta                                                                                                  <  2:804352/143‑1 (MQ=255)
                                      taGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACtata                                                                                                <  1:351795/143‑1 (MQ=255)
                                       aGCTTGACCTCCTTCACCGAGCTGGCGCGTGCTGCTCAGGATTGTGAACGCTGGAACATGCTGCAACCGCTTGAAACTCCCCGCGGGCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTatat                                                                                               <  1:661470/143‑1 (MQ=255)
                                         cGTGACCTGCGTCACCGATCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGa                                                                                             <  2:202404/143‑1 (MQ=255)
                                                tcctTCACCGAGCTGGCGGGTGATGGTCAGGATTGTGAACGCTGTAACATGGTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGGGATGAACGACGCTCCGGACGCGGGCTATACCGACACTGTGAAACTggt                                                                                      >  2:570844/1‑141 (MQ=255)
                                                tcctTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAGCTGCGGGCGTTCCCCATCAACACTCCACTTGGTTCGTCGCTCCGGCAGTGGCATCTACCGACTAGATGAAACCGtt                                                                                      >  2:190208/1‑143 (MQ=255)
                                                tcctTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACACTCCGGCAGTGGCATCTACCGACTATATGAAACTGtt                                                                                      >  2:632773/1‑143 (MQ=255)
                                                  ctccgccGGGCTGGTGCGTGATGCTCAGCACTCTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTACCGTCTATCGCGTAGGTGATGAACGACGCTCCGGCAGGGGCATCTACCGACTATATGAAACTGTTCg                                                                                    <  2:331295/138‑1 (MQ=255)
                                                  ctTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGTAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCg                                                                                    <  1:454260/143‑1 (MQ=255)
                                                  ctTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTAGATGAAACTGGTCg                                                                                    >  1:802814/1‑143 (MQ=255)
                                                                ggcggcgtGGTCAGGGTTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGt                                                                      <  2:667910/136‑1 (MQ=255)
                                                                 gcggtacggTCAGGAGTGTGAGGGCTGCAACATGCTGCACCCGCTGGAAACTCCGCGCGTTGCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTa                                                                     <  1:268757/136‑1 (MQ=255)
                                                                 gcgTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTGCCGTATAGCGCGTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTGGGCGGGGCAGGTCCGTa                                                                     >  1:559677/1‑143 (MQ=255)
                                                                 gcgGGATGGTCGGGAGTGTGAACGCTGGAACCTGCTGCACCCTCGGGAAACTCCGCGCGTTCCGGATATCGCTTGGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTa                                                                     <  2:121408/143‑1 (MQ=255)
                                                                              gATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGc                                                        >  1:84207/1‑143 (MQ=255)
                                                                              gATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGGACTTGCGTACCGGc                                                        >  1:858173/1‑143 (MQ=255)
                                                                                 tgtgAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGCGCAGGTCCGTACTTACGTACCGGCTgg                                                     >  2:605374/1‑142 (MQ=255)
                                                                                 tgtgAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCGCAGGTGATGAACGACGCTCCGCCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGCGCAGGTCCGTACTTACGTACCGGCTga                                                     >  2:836519/1‑143 (MQ=255)
                                                                                  gtgAACGCTGTAACATGCTGCACCCGGTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTgac                                                    <  2:712149/143‑1 (MQ=255)
                                                                                  gtgAACGCTGGAGCATGCTGCACCCGCGGGGGGCGGAGCCCCTCCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTgac                                                    <  1:376174/143‑1 (MQ=255)
                                                                                   tgAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTCCTTACGTACCAGCTgacg                                                   >  2:649839/1‑143 (MQ=255)
                                                                                    gagcgcgGGCACATGCTGCACCCGCGGGAAACGCAGCGCGTGCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTgacga                                                  <  1:16412/140‑1 (MQ=255)
                                                                                     aaCGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTgacgac                                                 >  1:312932/1‑143 (MQ=255)
                                                                                          tGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCGGGTCCGTGCTTGCGTACCGGCTGCCGACGACcg                                            >  1:685670/1‑143 (MQ=255)
                                                                                           ggAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAGCTGTTCGCTGAGCAGGTCCGTACTTACGTACAGACTAACGGCTACcgc                                           >  2:872005/1‑143 (MQ=255)
                                                                                             aaCATGCTGCCCCCGCTGGAAACTCCGCGCGTTCCAGACATCTCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGt                                         <  1:375767/143‑1 (MQ=255)
                                                                                                         cctcTGGAACCTCCGCGCGTGCGGTATATCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGa                             <  2:766739/140‑1 (MQ=255)
                                                                                                         ccGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGCTCCGTACTCACGTACCGGCTGACGACTACCGCGTACTGGGTACTGa                             >  2:614058/1‑143 (MQ=255)
                                                                                                         ccGCTGGAAACTCCGCGCGTTCCGTATAGCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATAGGAAACTGGTCGCTGAGCAGGTGCGTCATTACGTACCGGCTGAAGACTACCGCGTACTGGGTACTGa                             >  2:614054/1‑143 (MQ=255)
                                                                                                          cGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAACCTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGAt                            <  1:488661/143‑1 (MQ=255)
                                                                                                             tGGAAACTCCGCGCGTTCCGTGTATCGATTCGGTGAGGAACGACGCTCCGGCAGTGGGATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGc                         <  1:700301/143‑1 (MQ=255)
                                                                                                             tGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGAGGAACGACGCTCCGGCAGTGGCATCTCCCGACTATATGAAACTGTTCGCTGAGCAGGTAGGTACTTACGTACCGGCTAGCGACCACCGCGCAATGGCTACTGATgac                         >  2:668811/1‑141 (MQ=255)
                                                                                                               gAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGAt                            <  2:375936/138‑1 (MQ=255)
                                                                                                               gAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGAt                            >  1:375936/1‑138 (MQ=255)
                                                                                                                        cgcgTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTc              <  1:435361/143‑1 (MQ=255)
                                                                                                                         gcgTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTTGGTACTGATGGCTTCCGTCgaccc             >  2:369935/1‑139 (MQ=255)
                                                                                                                           tatctgTATATCACTTAGTTGATAAACGACGCTCCTGCAGTGGCTTTTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGa           <  2:241265/138‑1 (MQ=255)
                                                                                                                           gTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGa           <  1:499339/143‑1 (MQ=255)
                                                                                                                              cggtttatCGCTTAGGTGATGAACGTCGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAg        <  2:431846/138‑1 (MQ=255)
                                                                                                                                     atCGCTTAGGTGATGACCGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTga  <  1:315754/143‑1 (MQ=255)
                                                                                                                                     atCGCTTAGGTGAGGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTga  <  1:28119/143‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                        
ATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGA  >  NC_000913/125249‑125523

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