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,3870CT95.5% 120.8 / NA 45Q791* (CAG→TAG) aceEpyruvate dehydrogenase E1 component
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (15/27);  total (17/28)
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

AGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGAGA  >  NC_000913/125247‑125525
                                                                                                                                            |                                                                                                                                          
agaTCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTCAg                                                                                                                                          >  2:905407/1‑143 (MQ=255)
    ccTGGCGAAAGGTTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGtga                                                                                                                                      <  2:1033413/143‑1 (MQ=255)
    ccTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGtga                                                                                                                                      <  1:519163/143‑1 (MQ=255)
       ggCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGGGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGtgatga                                                                                                                                   >  1:88286/1‑143 (MQ=255)
         cGAAAGGTGACGGCGTAGGTCGGGTGGTTTATAGCGTGTGCTCCCCTAGCGGGCTGGCGCGTGATGTTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCTTATATCGCTTAGGTGATGAac                                                                                                                                 <  2:1152731/143‑1 (MQ=255)
                ttACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTAAGGTGATGAACGACGGTc                                                                                                                          >  1:459880/1‑143 (MQ=255)
                       gTAGGTGGTGACGTGGATCGCGTGACCTCCTGCGCCGGGCTGGCGCGTGACGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAGACGCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGt                                                                                                                   <  2:566493/143‑1 (MQ=255)
                        tAGGTTCTGACGTTGATGGCGCGACCTCCTTCACCGAGCTGGCGCGTGATTGTCAGGATTGTGAACGCTGTAATATGCTGCCCCCGCTGGAAACGCCGCGCGTTCCGTATATCTTTTAGGTGATGAACGACGCTCCGGCAGTg                                                                                                                  <  2:1051759/143‑1 (MQ=255)
                                gACGTTTTTAGCGTGACCTCCTTCACTGGGCTGGCCCGTGATGGTGAGGATTGTGACCGCTGGATCATGCGGCACCCGCTGGAAACTCCGCGCGTTCTGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTAc                                                                                                          <  2:975845/143‑1 (MQ=255)
                                           cGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGa                                                                                               <  1:1161878/143‑1 (MQ=255)
                                                   ctttCACTGCGCTGGAGCGTGAGGGGCAGGATTGGGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTc                                                                                       <  1:395411/141‑1 (MQ=255)
                                                   cctTCACCGAGCTGGCGCGTGAGGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTc                                                                                       <  1:1140494/143‑1 (MQ=255)
                                                    ctTCTCCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCg                                                                                      <  2:230284/143‑1 (MQ=255)
                                                    ctTCACCGGGCGGGCGCGGGATGGTCGGGATTGGGAGCGCTGGAAGATGCTGCACCCGCGGGAAGCTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCg                                                                                      <  1:435919/143‑1 (MQ=255)
                                                      tCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCt                                                                                    >  1:384322/1‑143 (MQ=255)
                                                              cTGCCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGt                                                                            <  2:490657/143‑1 (MQ=255)
                                                                           gTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATACCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGACACTGTTCGCTGCGCAGGTCCGTACTTACGTa                                                               >  1:3215/1‑143 (MQ=255)
                                                                            tCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTAc                                                              <  1:593691/143‑1 (MQ=255)
                                                                             cAGGACTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTAcc                                                             <  1:337807/143‑1 (MQ=255)
                                                                              aGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGCACATACGTACCg                                                            >  2:519612/1‑143 (MQ=255)
                                                                               gcaTTTTGGACGCTCGAGCATGCGTCACCCGCTGGATACTCCTCGCGTTCCGTACAGCGCTTAGGTGATGAACGACGCTCCGGCAGTGGGCTCTACAGACTCTGTGAAACTGTTCGCTGAGCCGGTACGTAAGTACGGAcctg                                                           <  2:203955/141‑3 (MQ=255)
                                                                                gATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTAATGACGTAGCGGc                                                          >  1:114911/1‑143 (MQ=255)
                                                                                  ctgtgAACGCTGGAACATGCGGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTg                                                        <  1:894442/142‑1 (MQ=255)
                                                                                       aaCGCTGGCACAGGCTTCGCCCGACGTAACCTCCGCGCGTTCCGTGTATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATAGGACACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTgacgac                                                   <  2:453598/143‑1 (MQ=255)
                                                                                       aaCGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTggcgac                                                   >  1:883705/1‑143 (MQ=255)
                                                                                           cTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTAcc                                               >  1:288678/1‑143 (MQ=255)
                                                                                             ggaaacagctgcACCCGGTGGAAACCCCGGGCGTTGCGTACATCGCTTAGGTGATGAACGACGCTCCGGCAGCGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACcgc                                             <  1:1202418/136‑1 (MQ=255)
                                                                                              gAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACcgcg                                            >  1:480438/1‑143 (MQ=255)
                                                                                                aCATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCAGCTAACGACTATATGAAACTGTTCGCTGAGCAGGTGCGTACTTACGTATCGGCAGAAGAGTACcgagga                                          >  2:670234/1‑139 (MQ=255)
                                                                                                 cATGCTGCACCCGCGGGAACCGCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTAc                                         <  1:1081066/143‑1 (MQ=255)
                                                                                                   tgctgtCCCCGCTGGAAACCGCGCGCGTTCCGTATATCGCGTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTg                                       <  1:68375/143‑1 (MQ=255)
                                                                                                    gctgcACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTgg                                      <  2:384322/143‑1 (MQ=255)
                                                                                                       gcACCCGCTGGAACCTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTa                                   <  1:337327/143‑1 (MQ=255)
                                                                                                        cACCCGCTGGAAACCCCGCGCGTGACGTATATCGCTTAGGTAATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTAc                                  <  2:329807/143‑1 (MQ=255)
                                                                                                           tggCTGGAAACGCCGCTCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGa                               <  2:879175/141‑1 (MQ=255)
                                                                                                             gCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATg                             >  1:24475/1‑143 (MQ=255)
                                                                                                             gCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATg                             >  1:506316/1‑143 (MQ=255)
                                                                                                                 gAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCtt                         <  1:905407/143‑1 (MQ=255)
                                                                                                                    gcTCCGCGCGTTCCGTATCTCGCTCTGGTGATGAACGACGCTCCGGCAGTGGCGTCTACCGACTATATGAGACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCgg                      <  2:683561/142‑1 (MQ=255)
                                                                                                                         gcgcgTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTCTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCCCGTACTGGGTACTGATGGCTTCGGGCGtt                 >  1:933471/1‑143 (MQ=255)
                                                                                                                          cgcgTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGCCTTCGGTCGTTc                >  2:869640/1‑143 (MQ=255)
                                                                                                                                    atatCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGGACTGATGGCTTCGGTCGTTCCGACAGCCGt      >  1:1032331/1‑143 (MQ=255)
                                                                                                                                    atatCGCTTAGGTGAGGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGt      <  1:527950/143‑1 (MQ=255)
                                                                                                                                      atCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTga    <  1:857036/143‑1 (MQ=255)
                                                                                                                                        cGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCCGACGACTACCGCGTACTGGGGACTGATGGCGTCGGCCGTCCCGACAGCCGTgagc  >  2:948917/1‑142 (MQ=255)
                                                                                                                                            |                                                                                                                                          
AGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGAGA  >  NC_000913/125247‑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: