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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273122,2710CT 5.9% 2520.2 / 71.6 793A30A (GCC→GCTaceFpyruvate dehydrogenase, dihydrolipoyltransacetylase component E2
Reads supporting (aligned to +/- strand):  ref base C (241/478);  new base T (2/44);  total (267/526)
Fisher's exact test for biased strand distribution p-value = 5.02e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 821 total aligned reads displayed.

GGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGT  >  CP009273/122205‑122337
                                                                  |                                                                  
  aCATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGCAGCAGa                                                                  >  1:46184267/1‑67 (MQ=255)
  aCATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGa                                                                  >  1:10287894/1‑67 (MQ=255)
  aCATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGa                                                                  >  1:33570319/1‑67 (MQ=255)
  aCATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGa                                                                  >  1:30081542/1‑67 (MQ=255)
  aCATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGa                                                                  >  1:42269789/1‑67 (MQ=255)
  aCATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGa                                                                  <  1:31948360/67‑1 (MQ=255)
  aCATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGa                                                                  >  1:1620348/1‑67 (MQ=255)
  aCATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGa                                                                  >  1:28325269/1‑67 (MQ=255)
  aCATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGa                                                                  >  1:17440180/1‑67 (MQ=255)
  aCATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCAGa                                                                  >  1:12370430/1‑67 (MQ=255)
    aTCGGGGCTGATGAAGTTGAATTCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACa                                                               >  1:6094426/1‑68 (MQ=255)
    aTCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGCAGCTGAACa                                                               >  1:39008355/1‑68 (MQ=255)
    aTCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACa                                                               >  1:24639713/1‑68 (MQ=255)
    aTCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACa                                                               <  1:30891379/68‑1 (MQ=255)
     tCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAg                                                              <  1:5378513/68‑1 (MQ=255)
     tCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCAGAACAg                                                              >  1:8679310/1‑68 (MQ=255)
      cGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGt                                                             <  1:6977220/68‑1 (MQ=255)
      cGGGGCGGATGAAGTTGAAATCACCGAGAGCCTGGTCAAAGTGGGCGCCAAAGTTGAAGCCGAACAGt                                                             <  1:16037657/68‑1 (MQ=255)
       ggggCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTc                                                            >  1:22357554/1‑68 (MQ=255)
       ggggCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTc                                                            <  1:2170752/68‑1 (MQ=255)
       ggggCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTc                                                            <  1:45785229/68‑1 (MQ=255)
       ggggCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTc                                                            <  1:25350189/68‑1 (MQ=255)
       ggggCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTc                                                            <  1:20354502/68‑1 (MQ=255)
       ggggCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTc                                                            <  1:74716/68‑1 (MQ=255)
       ggggCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTc                                                            <  1:13333610/68‑1 (MQ=255)
       ggggCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTc                                                            <  1:1293264/68‑1 (MQ=255)
       ggggCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTc                                                            <  1:37740723/68‑1 (MQ=255)
       cgggCTGTTGAAGTTGAATTCTCCGAGATCCTGGTCAAAGTGGGGGACAAAGTTGAATCCGAACAGTc                                                            <  1:21216003/67‑1 (MQ=255)
        gggCTGATGAAGTTGAAATCACCGAGATGCTGGTCCAAGTGGGCGACAAAGTTGAAGCACAACAGTc                                                            >  1:30597846/1‑67 (MQ=37)
        gggCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTc                                                            >  1:16710579/1‑67 (MQ=255)
          gCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGCAGCCGAACAGTc                                                            >  1:46872424/1‑65 (MQ=255)
          gCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTc                                                            >  1:20793962/1‑65 (MQ=255)
          gCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTc                                                            >  1:28756567/1‑65 (MQ=255)
          gCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGc                                                          <  1:27010771/67‑1 (MQ=255)
          gCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGc                                                          <  1:3636888/67‑1 (MQ=255)
          gCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGc                                                          <  1:46489988/67‑1 (MQ=255)
            tgatgaAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGa                                                       <  1:33840217/68‑1 (MQ=18)
            tgatgaAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTg                                                        <  1:49246062/67‑1 (MQ=255)
            tgatgaAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGa                                                       <  1:43527688/68‑1 (MQ=255)
               tgaAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGCACAGTCGCTGATc                                                     >  1:44496794/1‑67 (MQ=255)
               tgaAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATc                                                     <  1:19540457/67‑1 (MQ=255)
               tgaAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATc                                                     <  1:11648752/67‑1 (MQ=255)
               tgaAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATc                                                     <  1:915094/67‑1 (MQ=255)
               tgaAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCAGAACAGTCGCTGATc                                                     >  1:25939107/1‑67 (MQ=255)
               tgaAGTTGAAATCACCGACATCCTGGTCAAAGTTGGCGACAAAGTTGAAGCTGAACAGTCGCTGATCa                                                    <  1:14436434/68‑1 (MQ=14)
                gaAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCAc                                                   <  1:32468271/68‑1 (MQ=255)
                  aGTTGAATTCACCGAGATGCTGGTCAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCAcc                                                  <  1:9566093/67‑1 (MQ=11)
                  aGTTGAAATCACCGATATCCTGGTCAAAGTGGGCGAAAAAGTTGAAGCCGAACAGTCGCTGATCACCg                                                 <  1:23079333/68‑1 (MQ=255)
                  aGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCAcc                                                  <  1:24961982/67‑1 (MQ=255)
                  aGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCAcc                                                  <  1:38564856/67‑1 (MQ=255)
                  aGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCAcc                                                  <  1:17956386/67‑1 (MQ=255)
                  aGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCAcc                                                  <  1:18091970/67‑1 (MQ=255)
                  aGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCAcc                                                  <  1:24157363/67‑1 (MQ=255)
                  aGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCg                                                 <  1:37363891/68‑1 (MQ=255)
                    ttGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTa                                               >  1:8203953/1‑68 (MQ=255)
                     tGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGt                                                <  1:43077450/66‑1 (MQ=255)
                         aTCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAg                                           >  1:29825751/1‑67 (MQ=255)
                          tCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGc                                         >  1:10745724/1‑68 (MQ=255)
                           cACCGATATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCTGAACAGTc                                                            <  1:19934772/48‑1 (MQ=17)
                            cccGAGAGCCAGGTCAAAGTGGGCGACAAAGTAGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGa                                       <  1:348042/67‑1 (MQ=18)
                            aCCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGa                                       <  1:15651131/68‑1 (MQ=255)
                             ccGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGaa                                            >  1:47451102/1‑62 (MQ=255)
                               gagaTCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACaa                                    <  1:41262866/68‑1 (MQ=255)
                                agaTCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACaag                                   >  1:23662993/1‑67 (MQ=255)
                                  aTCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACaa                                    >  1:31164122/1‑65 (MQ=255)
                                   tCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAg                                  <  1:41024816/66‑1 (MQ=255)
                                   tCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAg                                  <  1:50388935/66‑1 (MQ=255)
                                   tCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGcc                                <  1:26362337/68‑1 (MQ=255)
                                   tCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGcc                                <  1:51095205/68‑1 (MQ=255)
                                   tCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGcc                                <  1:18815407/68‑1 (MQ=255)
                                   tCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGcc                                <  1:39219544/68‑1 (MQ=255)
                                   tCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGcc                                <  1:27752986/68‑1 (MQ=255)
                                    ccTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGcg                                >  1:35684270/1‑66 (MQ=255)
                                    ccTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGct                                >  1:39603861/1‑66 (MQ=255)
                                    ccTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCct                               <  1:49641860/68‑1 (MQ=255)
                                    ccTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCct                               <  1:44943440/68‑1 (MQ=255)
                                    ccTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCct                               <  1:39941900/68‑1 (MQ=255)
                                     cTGGTCAAAGTGGGCGACAAAGTTGAACCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCctc                              <  1:35224490/68‑1 (MQ=255)
                                      tGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCctct                             <  1:44260210/68‑1 (MQ=255)
                                        gTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTAt                           >  1:17186047/1‑68 (MQ=255)
                                         tCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCTTCTATg                          >  1:43856493/1‑68 (MQ=255)
                                         tCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTAt                           <  1:11326651/67‑1 (MQ=255)
                                         tCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATg                          >  1:47832504/1‑68 (MQ=255)
                                          cAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCGTCTATgg                         >  1:32751876/1‑68 (MQ=255)
                                          cAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATg                          <  1:6557410/67‑1 (MQ=255)
                                          cAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATg                          <  1:47226833/67‑1 (MQ=255)
                                              gTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGt                     >  1:27301473/1‑68 (MQ=255)
                                                gggCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGt                     <  1:2913446/66‑1 (MQ=255)
                                                    gACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGtc               <  1:138629/68‑1 (MQ=255)
                                                      cAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGtctc             <  1:36837617/68‑1 (MQ=255)
                                                        aaGTTGAAGCCGAACAGTCGCTGATCACCCTAGAAGGCGACAAAGAGTCTATGGAAGTTCAGTCTcc            >  1:42562359/1‑67 (MQ=14)
                                                          gTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGctt                               <  1:22149550/46‑3 (MQ=255)
                                                            tGAAGCCGAACAGTCGCTGTTCACCGTAGAAGGCGACAAAGCCTCTATGTAAGTTCCGTCTCCGCAgg       <  1:48771145/68‑1 (MQ=18)
                                                            tGAAGCCGAACAGTCGCTGATCACCGTAGACGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAgg       <  1:33288730/68‑1 (MQ=255)
                                                            tGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAgg       <  1:22714246/68‑1 (MQ=255)
                                                            tGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAgg       <  1:18319917/68‑1 (MQ=255)
                                                            tGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAg        >  1:14818270/1‑67 (MQ=255)
                                                                  cGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCggg   <  1:41724319/66‑1 (MQ=255)
                                                                  cGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCggg   <  1:34422822/66‑1 (MQ=255)
                                                                  cGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTATATGGAAGTTCCGTCTCCGCAGGCggg   <  1:45370461/66‑1 (MQ=255)
                                                                  |                                                                  
GGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGT  >  CP009273/122205‑122337

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: