Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 121,435 T→C 7.2% G644G (GGT→GGC)  aceE → pyruvate dehydrogenase, decarboxylase component E1, thiamine triphosphate‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273121,4350TC 7.2% 153.3 / 3.6 56G644G (GGT→GGC) aceEpyruvate dehydrogenase, decarboxylase component E1, thiamine triphosphate‑binding
Reads supporting (aligned to +/- strand):  ref base T (30/22);  new base C (2/2);  total (32/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.95e-01

CGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACATCACTACG  >  CP009273/121287‑121579
                                                                                                                                                    |                                                                                                                                                
cGTGCTACATCTATTACTCGATGTTCGGCGTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGt                                                                                                                                                  <  2:739600/149‑1 (MQ=255)
  ttCTACATCTATTCCCCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTca                                                                                                                                                <  1:485730/149‑1 (MQ=255)
   tCTACATCTGTTACTCGATGTTCGGCTGCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTcac                                                                                                                                               <  1:152275/149‑1 (MQ=255)
    ccactgctacttcTCGAGGTTCGGCTTCCAGCGTATTGGCTATGTGCGCTGGGCGGCTGGCGACGAGCGAGCGCGTGGCGTCATGATCGGCGTTACTTCCGGTCGGACCACCCTGAACGGGGAATGTGGTCAGCACGAAGATGTTcaca                                                                                                                                              <  1:44359/137‑1 (MQ=255)
      aCATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCCCAGc                                                                                                                                            >  1:372059/1‑149 (MQ=255)
      aCATCTATTACTCGAGGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCATCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCAGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGc                                                                                                                                            >  1:549205/1‑149 (MQ=255)
       cATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCCCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCCCGAAGATGGTCACAGcc                                                                                                                                           >  1:158301/1‑149 (MQ=255)
          cTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCcaca                                                                                                                                        <  2:571522/149‑1 (MQ=255)
           tATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACAt                                                                                                                                       <  2:325957/149‑1 (MQ=255)
            aTTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACAtt                                                                                                                                      >  1:173356/1‑149 (MQ=255)
              tacttcgtGTTCGGCTTCTAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGTTCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCa                                                                                                                                    <  1:67769/142‑1 (MQ=255)
               ccgcgagggtCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAg                                                                                                                                   <  2:187356/140‑1 (MQ=255)
               aCTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGCCACAGCCACATTCAg                                                                                                                                   >  2:310108/1‑149 (MQ=255)
               aCTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGCCACAGCCACATTCAg                                                                                                                                   >  2:85431/1‑149 (MQ=255)
                          ggCTTCCAGCGTATTGGCGCTCTGTGCTGGGCGGCTGGCGCCCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTAt                                                                                                                        <  2:615831/149‑1 (MQ=255)
                           gCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTacc                                                                                                                       >  1:261315/1‑147 (MQ=255)
                           gCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATc                                                                                                                       >  2:57822/1‑149 (MQ=255)
                                 aGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAAc                                                                                                                 <  2:571529/149‑1 (MQ=255)
                                  gctgaTTGGCGATCTGGGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACt                                                                                                                <  2:34361/145‑1 (MQ=255)
                                                             tGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGgcttacgctt                                                                                     >  1:638087/1‑149 (MQ=255)
                                                                   ccAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAg                                                                               >  2:167363/1‑149 (MQ=255)
                                                                    cAGCAAGCTCGTTGCGTCCTGATCGGCGGTACTTCCGGCCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGt                                                                              <  1:698924/149‑1 (MQ=255)
                                                                     aGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCGGAACTGTATCTATTACGACCCGGCTGACGCTTACGAAGtt                                                                             >  1:678489/1‑149 (MQ=255)
                                                                      gCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGACGATGGTCACAGCCCCATTCTGTCGGTGACTATCCCGAACTGGATCTCTTACGTCCCGGCTTACGCTGACGAAGTTg                                                                            <  1:602523/149‑1 (MQ=255)
                                                                      gCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACAGTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGGTTACGCTTACGAAGTTg                                                                            <  1:642041/149‑1 (MQ=255)
                                                                          gcgcGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCATGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGAGCCGGCTTACGCTTACGAAGTTGCTGt                                                                        >  2:196239/1‑149 (MQ=255)
                                                                           cgcgTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGtc                                                                       <  1:665097/149‑1 (MQ=255)
                                                                            gcgTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGCCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGtca                                                                      <  1:310108/149‑1 (MQ=255)
                                                                            gcgTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGCCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGtca                                                                      <  1:85431/149‑1 (MQ=255)
                                                                                gtcTTCCTGCGTGGCGGTTGTTCCGGTCCTACCACCCTGTACGGCGAAGGTCTGCAGCACGAAGATGGTCAACGCCACATTTAGTCGCTGACTATTCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTcatcat                                                                  <  1:506605/147‑1 (MQ=255)
                                                                                ggCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTcatcat                                                                  <  1:66866/149‑1 (MQ=255)
                                                                                       tGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGAc                                                           >  2:739458/1‑149 (MQ=255)
                                                                                          tCGGCGGTACTTCCGGTCGCACCACCCTGACCGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGt                                                        <  1:57442/149‑1 (MQ=255)
                                                                                              cggTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTgg                                                    >  1:763741/1‑149 (MQ=255)
                                                                                              cggTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGAACCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTgg                                                    >  2:671438/1‑149 (MQ=255)
                                                                                               ggTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCGTACGCTTACGAAGTTGCTGTCATCATGCATGACGGGCGGGa                                                   >  2:134976/1‑149 (MQ=255)
                                                                                                 tACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGc                                                 <  2:173356/149‑1 (MQ=255)
                                                                                                    ttCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTa                                              >  1:148357/1‑149 (MQ=255)
                                                                                                    ttCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACCCTTAGGAAGTTGCTGTCATCATGGATGACGGTGTGGAGCGTa                                              >  1:151340/1‑149 (MQ=255)
                                                                                                      ccGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGACGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCAACCTGCATGACGCTCTGGAGCGTATg                                            >  1:783158/1‑149 (MQ=255)
                                                                                                         gTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTAc                                         >  1:190884/1‑149 (MQ=255)
                                                                                                         gTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTAAGtcc                                         >  1:733992/1‑147 (MQ=255)
                                                                                                          tCGTCCCACCCTGCTGCGCGACGGTCTGCGGCTAGGAGTTGGTCAGAGCTGCATTCAGGCGCTGACTATCCCGAACTGTGTGTGTTACGACCAGGCTTACGCTTACGAAGTCGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACg                                        <  1:439160/149‑1 (MQ=255)
                                                                                                              accaccCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGa                                    >  1:20268/1‑149 (MQ=255)
                                                                                                                      gAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTAGCTCGTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGCGCGTATGTACGGTGAAACACAag                            >  1:6826/1‑149 (MQ=255)
                                                                                                                       aaCGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGGACGGTGAAAAACAaga                           <  2:549205/149‑1 (MQ=255)
                                                                                                                            cGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACg                      <  1:609886/149‑1 (MQ=255)
                                                                                                                            cGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACg                      >  1:536665/1‑149 (MQ=255)
                                                                                                                                    tGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTtactac              >  2:412671/1‑149 (MQ=255)
                                                                                                                                     gcagcaCGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTtactact             >  1:766669/1‑149 (MQ=255)
                                                                                                                                     gcaccaCTAAGATGGTCACCGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTtactact             <  2:586312/149‑1 (MQ=255)
                                                                                                                                      cagcaCGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCACCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTtaatacta            >  1:574803/1‑149 (MQ=255)
                                                                                                                                       agcaCGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGGATGTACGGTGAAAAACAAGAGAACGTTtactactac           >  2:291975/1‑149 (MQ=255)
                                                                                                                                       agcaCGAAGATGGTCACAGCCACATTCAGTCGCTGAATATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGGGAAAAACAAGAGAACGTTtaatactac           >  1:794492/1‑149 (MQ=255)
                                                                                                                                         caCGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTGCGGTGAAAAACAAGAGAACGTTtactactacct         >  1:645860/1‑147 (MQ=255)
                                                                                                                                             aaGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACATCACt     >  1:770493/1‑149 (MQ=255)
                                                                                                                                                aTGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACATCACTACg  >  1:537276/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
CGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACATCACTACG  >  CP009273/121287‑121579

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