Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 122,613 G→A 5.4% K144K (AAG→AAA aceF → pyruvate dehydrogenase, dihydrolipoyltransacetylase component E2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273122,6130GA 5.4% 912.9 / 19.0 300K144K (AAG→AAAaceFpyruvate dehydrogenase, dihydrolipoyltransacetylase component E2
Reads supporting (aligned to +/- strand):  ref base G (155/119);  new base A (12/4);  total (173/127)
Fisher's exact test for biased strand distribution p-value = 1.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 311 total aligned reads displayed.

AATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACA  >  CP009273/122547‑122680
                                                                  |                                                                   
caTCCTGGTGAACGTGCGCGACAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAgg                                                                    <  1:27038147/67‑1 (MQ=11)
aaTCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAgg                                                                    <  1:18619292/68‑1 (MQ=255)
aaTCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAg                                                                     <  1:25225077/67‑1 (MQ=255)
  tCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGCACAGTCGCTGATCACCGTAGAAGGCGACAAGGCt                                                                  >  1:27083816/1‑68 (MQ=255)
  tCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCTCTGATCACCGTAGAAGGCGACAAAGCt                                                                  <  1:23771782/68‑1 (MQ=255)
  tCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCt                                                                  <  1:9075684/68‑1 (MQ=255)
  tCCTGGTGAAAGTGGGCGATAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCt                                                                  <  1:15536363/68‑1 (MQ=255)
   ccTGGTGAGAGTTGGCGAAAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCtt                                                                 <  1:2487049/68‑1 (MQ=255)
     tGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCt                                                               <  1:4684139/68‑1 (MQ=255)
               tGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGtt                                                      <  1:25661593/67‑1 (MQ=255)
                ggCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTcc                                                    >  1:17836545/1‑68 (MQ=255)
                   gATAAAGTTGAAGCTGAACAGTCTCTGATCACCGTAGAAGGCGACAAAGCtt                                                                 <  1:13502620/52‑1 (MQ=14)
                      aaaGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTcc                                              <  1:10225077/68‑1 (MQ=255)
                         gTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCGCCGtt                                           >  1:28002379/1‑68 (MQ=255)
                           tGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTg                                         <  1:7848063/68‑1 (MQ=255)
                           tGAAGCAGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTg                                         <  1:17052832/68‑1 (MQ=255)
                               gCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTg                                      >  1:20154990/1‑67 (MQ=255)
                               gCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCGCCGTTTGcag                                      >  1:11040484/1‑65 (MQ=9)
                                   aaCAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCAcc                                 >  1:245056/1‑68 (MQ=255)
                                    aCAGTCGCTGATCACCGTAGACGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCAcc                                 <  1:28919151/67‑1 (MQ=255)
                                    aCAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCAcc                                 <  1:24249608/67‑1 (MQ=255)
                                    aCAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCAcc                                 <  1:22543580/67‑1 (MQ=255)
                                    aCAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCAcc                                 <  1:24578561/67‑1 (MQ=255)
                                    aCAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCAcc                                 <  1:11332236/67‑1 (MQ=255)
                                    aCAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCAcc                                 <  1:27287575/67‑1 (MQ=255)
                                    aCAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCAGCg                                >  1:17296680/1‑68 (MQ=255)
                                    aCAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCg                                >  1:21087210/1‑68 (MQ=255)
                                    aCAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCg                                >  1:8312761/1‑68 (MQ=255)
                                    aCAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCg                                >  1:26835964/1‑68 (MQ=255)
                                    aCAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCg                                >  1:28535682/1‑68 (MQ=255)
                                      aGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGGCGTg                              >  1:27617433/1‑68 (MQ=255)
                                         cGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGt                               <  1:10506996/64‑1 (MQ=255)
                                         cGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGt                               <  1:19450081/64‑1 (MQ=255)
                                         cGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGt                               <  1:16543084/64‑1 (MQ=255)
                                         cGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGaaa                           >  1:7648533/1‑68 (MQ=255)
                                         cGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGaaa                           >  1:10558156/1‑68 (MQ=255)
                                         cGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCAGGCGGCGTGaag                           >  1:21521120/1‑67 (MQ=11)
                                          gCTGATCACCGTAGAAGGCGACAACGCTTCTATGGAAGTTCCGGCTCCGGTTGCGGGGAACGTTAAag                          >  1:21746163/1‑68 (MQ=14)
                                            tGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTg                                         <  1:12093683/51‑1 (MQ=255)
                                            tGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAagag                        <  1:10812781/68‑1 (MQ=255)
                                            tGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAagag                        <  1:8815788/68‑1 (MQ=255)
                                            tGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAagag                        <  1:16432711/68‑1 (MQ=255)
                                            tGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAaga                         >  1:26088730/1‑67 (MQ=255)
                                              aTCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGCCGTGAATgaga                       >  1:5525365/1‑67 (MQ=255)
                                              aTCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCCCCGTGAAAgaga                       >  1:6838777/1‑67 (MQ=255)
                                              aTCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAgaga                       >  1:26424133/1‑67 (MQ=255)
                                              aTCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAgaga                       >  1:13086753/1‑67 (MQ=255)
                                              aTCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAgaga                       >  1:8581816/1‑67 (MQ=255)
                                              aTCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAgaga                       >  1:5393322/1‑67 (MQ=255)
                                              aTCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCGCCGTTTGCTGGCAACGTGAAAgaga                       >  1:20445243/1‑67 (MQ=255)
                                              aTCACCGTAGAAGGCGACAACGCGTCTATGGAAGTTCCGTCTCCGAAGGCTGGCACCGTGAAAgaga                       >  1:3889714/1‑67 (MQ=11)
                                              aTCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGGCTCCGTTTGCTGGCAACGTGAAAgaga                       >  1:14058494/1‑67 (MQ=17)
                                               tCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATc                     >  1:22761376/1‑68 (MQ=255)
                                               tCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCAGGCAACGTGAAGGAGATc                     >  1:11420064/1‑68 (MQ=14)
                                               tCACCGTAGAAGGCGACAAAGCATCTATGGAAGTTCCGGCGCCGAATGCTGGCTACGTGAAAGAGAtg                     >  1:16859780/1‑67 (MQ=11)
                                                cACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGtt                                           <  1:5765940/45‑1 (MQ=255)
                                                cACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTg                                         <  1:2083820/47‑1 (MQ=255)
                                                cACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCAcc                                 <  1:13258951/55‑1 (MQ=255)
                                                cACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCa                    >  1:3788371/1‑68 (MQ=255)
                                                 aCCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCCCCGTGAAAGAGATCaa                   >  1:14841518/1‑68 (MQ=255)
                                                 aCCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATGaa                   >  1:2434452/1‑68 (MQ=255)
                                                  ccttagtagGCTACTAGTCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCaaa                  <  1:15184452/65‑1 (MQ=25)
                                                  ccGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCaaa                  <  1:4517832/68‑1 (MQ=255)
                                                   cGTAGAATGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGaa                            <  1:1627906/57‑1 (MQ=255)
                                                    gTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCCCCGTGAATGATATCAAAGt                >  1:8755589/1‑68 (MQ=18)
                                                    gTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCCCCGTGAAAGAGATCAAAGt                >  1:1001390/1‑68 (MQ=255)
                                                    gTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGt                >  1:10188002/1‑68 (MQ=255)
                                                    gTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGt                >  1:7374524/1‑68 (MQ=255)
                                                    gTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGt                >  1:15639776/1‑68 (MQ=255)
                                                    gTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTATGCTGGCACCGTGAAAGAGATCAAAGt                >  1:14573737/1‑68 (MQ=255)
                                                    gTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCGCCGTTTGCTGGCGCCGTGAAGGAGCTCAAAGt                >  1:518710/1‑68 (MQ=11)
                                                    gTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGTTTGCTGGCAACGTGAAAGAGATCAAAGt                >  1:23078572/1‑68 (MQ=14)
                                                      aGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCCGCGTGAAAGATATCAAAGtca              >  1:11678040/1‑66 (MQ=18)
                                                      aGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAATGAGATCAAAGTGa              >  1:4031460/1‑68 (MQ=255)
                                                         aGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCCCCGTGAAAGAGATCAAAGTGAACg           >  1:2845525/1‑68 (MQ=255)
                                                         aGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAGGAGATGAAAGTCAACg           >  1:2134951/1‑68 (MQ=17)
                                                         aGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGCTGAAAGTCAACg           >  1:1205273/1‑68 (MQ=17)
                                                         aGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACg           >  1:9447612/1‑68 (MQ=255)
                                                         aGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACg           >  1:696442/1‑68 (MQ=255)
                                                         aGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACg           >  1:28078680/1‑68 (MQ=255)
                                                         aGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACg           >  1:21834293/1‑68 (MQ=255)
                                                         aGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACg           >  1:2937836/1‑68 (MQ=255)
                                                         aGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCAACGTGAAAGATATCAAAGTGAACg           >  1:5841666/1‑68 (MQ=255)
                                                         aGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGTTTGCTGGCACCGTTAAAGAGATCAAAGTGaaa            >  1:25750286/1‑66 (MQ=14)
                                                           gCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCGACGTGAATGAGCTCAAAGTGAACGTg         >  1:13754190/1‑68 (MQ=17)
                                                           gCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTg         >  1:15235700/1‑68 (MQ=255)
                                                           gCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTg         >  1:24037919/1‑68 (MQ=255)
                                                            cGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGt          <  1:3395748/66‑1 (MQ=255)
                                                            cGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGt          <  1:1823489/66‑1 (MQ=255)
                                                              aCAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGt      >  1:19779760/1‑68 (MQ=255)
                                                               cAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGt      <  1:322617/67‑1 (MQ=255)
                                                               cAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGt      <  1:15935522/67‑1 (MQ=255)
                                                               cAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGt      <  1:6103186/67‑1 (MQ=255)
                                                               cAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGt      <  1:14888465/67‑1 (MQ=255)
                                                               cAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGt      <  1:17499005/67‑1 (MQ=255)
                                                               cAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGt      <  1:7928884/67‑1 (MQ=255)
                                                               cAACGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGt      <  1:15563813/67‑1 (MQ=255)
                                                                aaGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGa    <  1:12628830/68‑1 (MQ=255)
                                                                 aGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGAc   <  1:14998155/68‑1 (MQ=255)
                                                                  ggCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACa  >  1:12337546/1‑68 (MQ=255)
                                                                  |                                                                   
AATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACA  >  CP009273/122547‑122680

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: