Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 122,580 T→C 5.9% A133A (GCT→GCC aceF → pyruvate dehydrogenase, dihydrolipoyltransacetylase component E2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273122,5800TC 5.9% 628.2 / 11.2 206A133A (GCT→GCCaceFpyruvate dehydrogenase, dihydrolipoyltransacetylase component E2
Reads supporting (aligned to +/- strand):  ref base T (107/78);  new base C (4/8);  total (111/95)
Fisher's exact test for biased strand distribution p-value = 1.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 206 total aligned reads displayed.

GGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGC  >  CP009273/122514‑122646
                                                                  |                                                                  
 gATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCATGGTGAAAGTTGGCGAAAAAGTTGCAGCTGa                                                                  >  1:31491975/1‑68 (MQ=18)
  atatCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGa                                                                  >  1:3623630/1‑67 (MQ=255)
  atatCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGa                                                                  >  1:21081213/1‑67 (MQ=255)
  atatCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGa                                                                  >  1:33145261/1‑67 (MQ=255)
  atatCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGa                                                                  >  1:27187681/1‑67 (MQ=255)
  atatCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGAGAAAGTTGCCGCTGa                                                                  >  1:37024055/1‑67 (MQ=18)
   tatCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAAc                                                                >  1:28713128/1‑68 (MQ=255)
    atCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACa                                                               <  1:21563055/68‑1 (MQ=255)
    atCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACa                                                               >  1:15241560/1‑68 (MQ=255)
      cGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGt                                                             <  1:7287623/68‑1 (MQ=255)
       ggCTGCGACGAAGTTGAAGTGACCGATATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTc                                                            <  1:14819889/68‑1 (MQ=255)
       ggCAGCGACGAACTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAAAAGTc                                                            <  1:35892568/68‑1 (MQ=255)
        gCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTc                                                            >  1:18214240/1‑67 (MQ=255)
        gCAGCGACGAAGTTGAAGTGACAGAAGTCCTGGTGAAAGTTGGCGAAAAAGTTGCCGCTGAACAGTc                                                            >  1:22646236/1‑67 (MQ=14)
          tgtgacgaAGTTGAACTGACCGATATCCTGGTGAAAGTTGGCGACAAAGTTGAAGCTGAACAGTCGc                                                          <  1:30110019/64‑1 (MQ=17)
          aGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTc                                                            >  1:18612695/1‑65 (MQ=255)
          aGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTc                                                            >  1:25362472/1‑65 (MQ=255)
          aGCGACGAAGTTGAAGTGACCGAAATCCTGGCGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTc                                                            >  1:24338047/1‑65 (MQ=255)
           gCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTg                                                        >  1:20541741/1‑68 (MQ=255)
            cgatgaAGTTGAACTCACCGACATCCTGGTGAAAGTTGGCGAAAAAGTTGAAGCTGAACAGTCGCTGa                                                       <  1:19511932/68‑1 (MQ=11)
             gatgaAGTTGAACTCACCGATATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGa                                                       <  1:35568199/64‑1 (MQ=14)
             gatgaAGTTGAACTCACCGACATCCTGGTCAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGa                                                       <  1:33668948/64‑1 (MQ=1)
             gacgaAGTTGAAGTGACCGAAATCCTGGGGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCt                                                         <  1:34741236/65‑1 (MQ=255)
             gacgaAGTTGAACTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGa                                                       <  1:34422065/67‑1 (MQ=255)
              acgaAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAATCTGAACAGTCGCTGATc                                                     >  1:7795829/1‑68 (MQ=255)
               tgaAGTTGAACTGACCGACATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATc                                                     <  1:13416356/66‑1 (MQ=17)
               tgaAGTTGAACTCACCGACATCCTGGTCAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATc                                                     <  1:12451725/66‑1 (MQ=11)
               cgaAGTTGAAGTGACCGAGATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATc                                                     <  1:31933409/67‑1 (MQ=255)
               cgaAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATc                                                     <  1:30456888/67‑1 (MQ=255)
                 atgTTGAAGTGACCGAAATGCTGGTGAAAGTTGGCGATAAAGTTGCAGCTGAACAGTCTCTGATCAcc                                                  <  1:11339461/66‑1 (MQ=255)
                  aGTTGAATTGACCGATATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCAcc                                                  <  1:32376783/67‑1 (MQ=255)
                  aGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCAcc                                                  <  1:37015427/67‑1 (MQ=255)
                  aGTTGAACTGACCGATATCCTGGTCAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCAcc                                                  <  1:32443353/67‑1 (MQ=14)
                        aGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGaa                                            >  1:4299392/1‑67 (MQ=255)
                        aGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAg                                           <  1:6930951/68‑1 (MQ=255)
                        aGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAg                                           >  1:7340217/1‑68 (MQ=255)
                        aGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAg                                           >  1:13596498/1‑68 (MQ=255)
                        aGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAg                                           >  1:29128728/1‑68 (MQ=255)
                        aGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAg                                           >  1:14939065/1‑68 (MQ=255)
                        aGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAg                                           >  1:33290996/1‑68 (MQ=255)
                        aGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAg                                           >  1:19821577/1‑68 (MQ=255)
                        aGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAg                                           >  1:17825253/1‑68 (MQ=255)
                        aGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAg                                           >  1:4793272/1‑68 (MQ=255)
                        aGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAg                                           >  1:4485381/1‑68 (MQ=255)
                        aGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAg                                           <  1:16610107/68‑1 (MQ=255)
                        aGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCACTGATCACCGTa                                               >  1:31242749/1‑64 (MQ=255)
                          tGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGc                                         >  1:2687040/1‑68 (MQ=255)
                          tGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCACTGATCAc                                                   >  1:7217885/1‑58 (MQ=255)
                           gACCGAAATCCTGGTGAAAGTGGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCg                                        <  1:1447339/68‑1 (MQ=255)
                              cGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACa                                     <  1:29825033/68‑1 (MQ=255)
                                 taTCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAgg                                  <  1:2444329/67‑1 (MQ=255)
                                 aaTCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAgg                                  <  1:17136791/68‑1 (MQ=255)
                                 aaTACTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAgg                                  <  1:24767043/68‑1 (MQ=255)
                                  aTCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTa                                               >  1:30897709/1‑54 (MQ=255)
                                  aTCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGc                                 <  1:8334489/68‑1 (MQ=255)
                                   tCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCt                                >  1:163290/1‑68 (MQ=255)
                                   tCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCt                                >  1:30074420/1‑68 (MQ=255)
                                   tCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCt                                <  1:5120735/68‑1 (MQ=255)
                                   tCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCt                                >  1:9023418/1‑68 (MQ=255)
                                   tCCTGGTGAAAGTGGGCGATAAAGTTGAAGCGGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCt                                <  1:32609570/68‑1 (MQ=255)
                                    ccTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCt                                >  1:9777377/1‑67 (MQ=255)
                                    ccTGGTCAAAGTTGGCGATAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCtt                               <  1:7507772/68‑1 (MQ=255)
                                     cTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGGCGCTGATCACCGTAGAAGGCGACAAGGCTTc                              <  1:9308790/68‑1 (MQ=255)
                                     cTGGTCAAAGTGGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTc                              <  1:2555407/68‑1 (MQ=255)
                                      tGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCt                             <  1:13951423/68‑1 (MQ=255)
                                      tGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCt                             <  1:9652075/68‑1 (MQ=255)
                                      tGGTCAAAGTGGGCGAAAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCt                             <  1:6502706/68‑1 (MQ=11)
                                       ggTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTa                            >  1:1060456/1‑68 (MQ=255)
                                         tGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATg                          >  1:23430053/1‑68 (MQ=255)
                                         tGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATg                          >  1:20865011/1‑68 (MQ=255)
                                          gAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATgg                         >  1:35026754/1‑68 (MQ=255)
                                          gAAAGTTGGCGATAAAGTTGAAGCAGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATg                          <  1:1678343/67‑1 (MQ=255)
                                             aGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCtt                               >  1:25178841/1‑59 (MQ=255)
                                              gTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGt                     >  1:15351573/1‑68 (MQ=255)
                                              gTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGt                     >  1:4966472/1‑68 (MQ=255)
                                              gTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGt                     >  1:18835795/1‑68 (MQ=255)
                                              gTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGt                     >  1:7862811/1‑68 (MQ=255)
                                              gTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGt                     >  1:22814253/1‑68 (MQ=255)
                                                gggCGAGAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGtt                    <  1:23993174/66‑1 (MQ=255)
                                                gggCGACAAAGTTGAAGCTGAACAGTCGCTGATCCCCGTAGAAGGCGACAAGGCTTCTATGGAAGtt                    <  1:12909042/66‑1 (MQ=14)
                                                    gATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGc               <  1:31497830/68‑1 (MQ=255)
                                                      tAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTc             <  1:13195809/68‑1 (MQ=255)
                                                       aaaGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTcc            <  1:19960922/68‑1 (MQ=255)
                                                       aaaGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTcc            <  1:36564388/68‑1 (MQ=255)
                                                       aaaGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTcc            >  1:36990432/1‑68 (MQ=255)
                                                        aaGTTGAAGCGGAACAGTCGCTGATCACCGTAGAAGGCGGCAAGGCTTCTATGGAAGTTCCGGCTCCg           <  1:33231817/68‑1 (MQ=14)
                                                        aaGTTGAAGCGGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCg           <  1:24111778/68‑1 (MQ=255)
                                                          gTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGtt         >  1:34657053/1‑68 (MQ=255)
                                                            tGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGttt        >  1:26045751/1‑67 (MQ=255)
                                                            tGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGtt         <  1:1001391/66‑1 (MQ=255)
                                                            tGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTg       <  1:36790811/68‑1 (MQ=255)
                                                            tGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGtt         <  1:13148672/66‑1 (MQ=255)
                                                            tGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTg       <  1:15704362/68‑1 (MQ=255)
                                                            tGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTg       <  1:23897690/68‑1 (MQ=255)
                                                            tGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTg       <  1:33099990/68‑1 (MQ=255)
                                                            tGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCGTCTATGGAAGTTCCGGCTCCGttt        >  1:16486487/1‑67 (MQ=255)
                                                                gCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGcag    >  1:11141367/1‑65 (MQ=255)
                                                                gCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTg    >  1:20755504/1‑67 (MQ=255)
                                                                gCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCGCCGTTTGCTg    >  1:36080421/1‑67 (MQ=255)
                                                                 cTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGc  >  1:6868212/1‑68 (MQ=255)
                                                                  |                                                                  
GGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGC  >  CP009273/122514‑122646

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: