Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 120,914 A→C T471P (ACC→CCC)  aceE → pyruvate dehydrogenase (acetyl‑transferring), homodimeric type

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273120,9140AC75.4% 126.1 / 42.1 61T471P (ACC→CCC) aceEpyruvate dehydrogenase (acetyl‑transferring), homodimeric type
Reads supporting (aligned to +/- strand):  ref base A (8/7);  new base C (22/24);  total (30/31)
Fisher's exact test for biased strand distribution p-value = 7.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

TGTCTGATGCAGATATCGAAAAACTGCCGTACATCACCTTCCCGGAAGGTTCTGAAGAGCATACCTATCTGCACGCTCAGCGTCAGAAACTGCACGGTTATCTGCCAAGCCGTCAGCCGAACTTCACCGAGAAGCTTGAGCTGCCGAGCCTGCAAGACTTCGGCGCGCTGTTGGAAGAGCAGAGCAAAGAGATCTCTACCACTATCGCTTTCGTTCGTGCTCTGAACGTGATGCTGAAGAACAAGTCGATCAAAGA  >  NZ_CP009273/120789‑121044
                                                                                                                             |                                                                                                                                  
tgtCTGATGCAGATATCGAAAAACTGCCGTACATCACCTTCCCGGAAGGTTCTGAAGAGCATACCTATCTGCACGCTCAGCGTCAGAAACTGCACGGTTATCTGCCAAGCCGTCAGCCGAACTTCCCCGAGAAGCtt                                                                                                                         >  1:735889/1‑137 (MQ=255)
    tGATGCAGATATCGAAAAACTGCCGTACATCACCTTCCCGGAAGGTTCTGAAGAGCATACCTATCTGCACGCTCAGCGTCAGAAACTGCACGGTTATCTGCCAAGCCGTCAGCCGAACTTCCCCGAGAAGCTTGAGCTg                                                                                                                   >  2:765345/1‑139 (MQ=255)
      aTGCAGATATCGAAAAACTGCCGTACATCACCTTCCCGGAAGGTTCTGAAGAGCATACCTATCTGCACGCTCAGCGTCAGAAACTGCACGGTTATCTGCCAAGCCGTCAGCCGAACTTCCCCGAGAAGCTTGAGCTGcc                                                                                                                 >  1:241640/1‑139 (MQ=255)
      aTGCAGATATCGAAAAACTGCCGTACATCACCTTCCCGGAAGGTTCTGAAGAGCATACCTATCTGCACGCTCAGCGTCAGAAACTGCACGGTTATCTGCCAAGCCGTCAGCCGAACTTCCCCGAGAAGCTTGAGCTGcc                                                                                                                 >  2:766354/1‑139 (MQ=255)
       tGCAGATATCGAAAAACTGCCGTACATCACCTTCCCGGAAGGTTCTGAAGAGCATACCTATCTGCACGCTCAGCGTCAGAAACTGCACGGTTATCTGCCAAGCCGTCAGCCGAACTTCACCGAGAAGCTTGAGCTGcc                                                                                                                 >  1:584358/1‑138 (MQ=255)
            atatCGAAAAACTGCCGTACATCACCTTCCCGGCAGGTTCTGAAGAGCATACCTATCTGCACGCTCAGCGTCAGAAACTGCACGGTTATCTGACAAGCCGTCAGCCGCACTTCACCGAGAAGCTTGAGCTGCCGAGCCt                                                                                                           >  2:40238/1‑139 (MQ=255)
            atatCGAAAAACTGCCGTACATCACCTTCCCGGAAGGTTCTGAAGAGCATACCTATCTGCACGCTCAGCGTCAGAAACTGCACGGTTATCTGCCAAGCCGTCAGCCGAACTTCCCCGAGAAGCTTGAGCTGCCGAGCCt                                                                                                           >  1:202809/1‑139 (MQ=255)
            atatCGAAAAACTGCCGTACATCACCTTCCCGGAAGGTTCTGAAGAGCATACCTATCTGCACGCTCAGCGTCAGAAACTGCACGGTTATCTGCCAAGCCGTCAGCCGAACTTCCCCGAGAAGCTTGAGCTGCCGAGCCt                                                                                                           >  2:42751/1‑139 (MQ=255)
             tatCGAAAAACTGCCGTACATCACCTTCCCGGAAGGTTCTGAAGAGCATACCTATCTGCACGCTCAGCGTCAGAAACTGCACGGTTATCTGCCAAGCCGTCAGCCGAACTTCCCCGAGAAGCTTGAGCTGCCGAGCCTg                                                                                                          <  2:735889/139‑1 (MQ=255)
                           cgtccgtCACCTTCCCGGAAGGTTCTGAAGAGCATCCCTATCTGCACGCTCAGCGTCAGAAACTGCACGGTTATCTGCCAAGCCGTCAGCCGAACTTCCCCGAGAAGCTTGAGCTGCCGAGCCTGCAAGACTTCGgcgc                                                                                            <  2:207038/133‑1 (MQ=255)
                           cGTACATCACCTTCCCGGAAGGTTCTGAAGAGCATACCTATCTGCACGCTCAGCGTCAGAAACTGCACGGTTATCTGCCAAGCCGTCAGCCGAACTTCACCGAGAAGCTTGAGCTGCCGAGCCTGCAAGACTTCGgcgc                                                                                            <  2:378150/139‑1 (MQ=255)
                            gTACATCACCTTCCCGGAAGGTTCTGAAGAGCATACCTATCTGCACGCTCAGCGTCAGAAACTGCACGGTTATCTACCAAGCCGTCAGCCGAACTTCACCGAGAAGCTTGAGCTGCCGAGCCTGCAAGACTTCGgcgcg                                                                                           <  1:2946/139‑1 (MQ=255)
                            gTACAGCACCTTCCCGGAAGGTTCTGAAGAGCATACCTATCTGCACGCTCAGCGTCAGAAACTGCACGGTTATCTGCCAAGCCGTCAGCCGAACTTCCCCGAGAAGCTTGAGCTGCCGAGCCTGCAAGACTGCGgcgcg                                                                                           <  2:26371/139‑1 (MQ=255)
                                  cACCTTCCCGGAAGGTTCTGAAGAGCATACCTATCTGCACGCTCAGCGTCAGAAACTGCACGGTTATCTGCCAAGCCGTCAGCCGAACTTCCCCGAGAAGCTTGAGCTGCCGAGCCTGCAAGACTTCGGCGCGCTGTTg                                                                                     >  2:351158/1‑139 (MQ=255)
                                   aCCTTCCCGGAAGGTTCTGAAGAGCATACCTATCTGCACGCTCAGCGTCAGAAACTGCACGGTTATCTGCCAAGCCGTCAGCCGAACTTCCCCGAGAAGCTTGAGCTGCCGAGCCTGCAAGACTTCGGCGCGCTGTTgg                                                                                    <  1:766354/139‑1 (MQ=255)
                                    ccTTCCCGGAAGGTTCTGAAGAGCATACCTATCTGCACGCTCAGCGTCAGAAACTGCACGGTTATCTGCCAAGCCGTCAGCCGAACTTCCCCGAGAAGCTTGAGCTGCCGAGCCTGCAAGACTTCGGCGCGCt                                                                                         <  1:841595/133‑1 (MQ=255)
                                    ccTTCCCGGAAGGTTCTGAAGAGCATACCTATCTGCACGCTCAGCGTCAGAAACTGCACGGTTATCTGCCAAGCCGTCAGCCGAACTTCCCCGAGAAGCTTGAGCTGCCGAGCCTGCAAGACTTCGGCGCGCt                                                                                         >  2:841595/1‑133 (MQ=255)
                                       tCCCGGAAGGTTCTGAAGAGCATACCTATCTGCACGCTCAGCGTCAGAAACTGCACGGTTATCTGCCAAGCCGTCAGCCGAACTTCCCCGAGAAGCTTGAGCTGCCGAGCCTGCAAGACTTCGgcgcg                                                                                           >  1:418160/1‑128 (MQ=255)
                                       tCCCGGAAGGTTCTGAAGAGCATACCTATCTGCACGCTCAGCGTCAGAAACTGCACGGTTATCTGCCAAGCCGTCAGCCGAACTTCCCCGAGAAGCTTGAGCTGCCGAGCCTGCAAGACTTCGgcgcg                                                                                           <  2:418160/128‑1 (MQ=255)
                                        cccGGAAGGTTCTGAAGAGCATACCTATCTGCACGCTCAGCGTCAGAAACTGCACGGTTATCTGCCAAGCCGTCAGCCGAACTTCCCCGAGAAGCTTGAGCTGCCGAGCCTGCAAGACTTCGGCGCGCTGTTGGAagag                                                                               >  1:493400/1‑139 (MQ=255)
                                         ccGGAAGGTTCTGAAGAGCATACCTATCTGCACGCTCAGCGTCAGAAACTGCACGGTTATCTGCCAAGCCGTCAGCCGAACTTCACCGAGAAGCTTGAGCTGCCGAGCCTGCAAGACTTCGGCGCGCTGTTGGAagagc                                                                              <  1:510212/139‑1 (MQ=255)
                                         ccGGAAGGTTCTGAAGAGCATACCTACCGTCACGCTCAGCGTCAGAAACTGCACGGTTATCCGCCAAGCCGTCAGCCGAACTTCCCCGAGAAGCTTGAGCTGCCGAGCCTGCAAGACTTCGGCGCGCTGTTGGAagagc                                                                              <  2:219508/139‑1 (MQ=255)
                                                  tCTGAAGAGCATACCTATCTGCACGCTCAGCGTCAGAAACTGCACGGTTATCTGCCAAGCCGTCAGCCGAACTTCCCCGAGAAGCTTGAGCTGCCGAGCCTGCAAGACTTCGGCGCGCTGTTGGAAgagcagagca                                                                        >  2:733540/1‑136 (MQ=255)
                                                  tCTGAAGAGCATACCTATCTGCACGCTCAGCGTCAGAAACTGCACGGTTATCTGCCAAGCCGTCAGCCGAACTTCACCGAGAAGCTTGAGCTGCCGAGCCTGCAAGACTTCGGCGCGCTGTTGGAagagc                                                                              >  2:452650/1‑130 (MQ=255)
                                                  tCTGAAGAGCATACCTATCTGCACGCTCAGCGTCAGAAACTGCACGGTTATCTGCCAAGCCGTCAGCCGAACTTCACCGAGAAGCTTGAGCTGCCGAGCCTGCAAGACTTCGGCGCGCTGTTGGAagagc                                                                              <  1:452650/130‑1 (MQ=255)
                                                      aaGAGCATACCTATCTGCACGCTCAGCGTCAGAAACTGCACGGTTATCTGCCAAGCCGTCAGCCGAACTTCCCCGAGAAGCTTGAGCTGCCGAGCCTGCAAGACTTCGGCGCGCTGTTGGAAGAGCAGAGCAAAGAGAt                                                                 <  1:493257/139‑1 (MQ=255)
                                                      aaGAGCATACCTATCTGCACGCTCAGCGTCAGAAACTGCACGGTTATCTGCCAAGCCGTCAGCCGAACTTCACCGAGAAGCTTGAGCTGCCGAGCCTGCAAGACTTCGGCGCGCTGTTGGAAGAGCAGAGCAAagag                                                                   <  2:584358/137‑1 (MQ=255)
                                                      aaGAGCATACCTATCTGCACGCTCAGCGTCAGAAACTGCACGGTTATCTGCCAAGCCGTCAGCCGAACTTCACCGAGAAGCTTGAGCTGCCGAGCCTGCAAGACTTCGGCGCGCTGTTGGAAGAGCAGAGCAAagag                                                                   <  1:493246/137‑1 (MQ=255)
                                                      aaGAGCATACCTATCTGCACGCGCAGCGTCAGAAACTGCACGGTTATCTGCCAAGCCGTCAGCCGAACTTCCCCGAGAAGCTTGAGCTGCCGAGCCTGCAAGACTTCGGCGCGCTGTTGGAAGAGCAGAGCAAagag                                                                   <  2:22752/137‑1 (MQ=255)
                                                           cATACCTATCTGCACGCTCAGCGTCAGAAACTGCACGGTTATCTGCCAAGCCGTCAGCCGAACTTCCCCGAGAAGCTTGAGCTGCCGAGCCTGCAAGACTTCGgcgc                                                                                            <  2:531498/107‑1 (MQ=255)
                                                           cATACCTATCTGCACGCTCAGCGTCAGAAACTGCACGGTTATCTGCCAAGCCGTCAGCCGAACTTCCCCGAGAAGCTTGAGCTGCCGAGCCTGCAAGACTTCGgcgc                                                                                            >  1:531498/1‑107 (MQ=255)
                                                           cATACCTATCTGCACGCTCAGCGTCAGAAACTGCACGGTTATCTGCCAAGCCGTCAGCCGAACTTCCCCGAGAAGCTTGAGCTGCCGAGCCTGCAAGACTTCGGCGCGCTGt                                                                                       >  1:82890/1‑112 (MQ=255)
                                                           cATACCTATCTGCACGCTCAGCGTCAGAAACTGCACGGTTATCTGCCAAGCCGTCAGCCGAACTTCCCCGAGAAGCTTGAGCTGCCGAGCCTGCAAGACTTCGGCGCGCTGt                                                                                       <  2:82890/112‑1 (MQ=255)
                                                           cATACCTATCTGCACGCTCAGCGTCAGAAACTGCACGGTTATCTGCCAAGCCGTCAGCCGAACTTCCCCGAGAAGCTTGAGCTGCCGAGCCTGCAAGACTTCGGCGCGCTGTTGGAAGAGCAGAGCAAAGAGATCTCTa                                                            >  2:344129/1‑139 (MQ=255)
                                                           cATACCTATCTGCACGCTCAGCGTCAGAAACTGCACGGTTATCTGCCAAGCCGTCAGCCGAACTTCACCGAGAAGCTTGAGCTGCCGAGCCTGCAAGACTTCGGCGCGCTGTTGGAAGAGCAGAGCAAAGAGATCTCTa                                                            >  2:460631/1‑139 (MQ=255)
                                                           cATACCTATCTGCACGCCCAGCGTCAGAAACTGCACGGTTATCTGCCAAGCCGTCAGCCGAACTTCCCCGAGAAGCTTGAGCTGCCGAGCCTGCAAGACTTCGGCGCGCTGTTGGAAGAGCAGAGCAAAGAGATCTCTa                                                            >  2:821088/1‑139 (MQ=255)
                                                             tACCTATCTGCACGCTCAGCGTCAGAAACTGCACGGTTATCTGCCAAGCCGTCAGCCGAACTTCCCCGAGAAGCTTGAGCTGCCGAGCCTGCAAGACTTCGGCGCGCTGTTGGAAGAGCAGAGCAAAGAGATCTCTAc                                                           >  2:324659/1‑138 (MQ=255)
                                                                   tCTGCACGCTCAGCGTCAGAAACTGCACGGTTATCTGCCAAGCCGTCAGCCGAACTTCCCCGAGAAGCTTGAGCTGCCGAGCCTGCAAGACTTCGGCGCGCTGTTGGAAGAGCAGAGCAAagag                                                                   <  1:91107/124‑1 (MQ=255)
                                                                   tCTGCACGCTCAGCGTCAGAAACTGCACGGTTATCTGCCAAGCCGTCAGCCGAACTTCCCCGAGAAGCTTGAGCTGCCGAGCCTGCAAGACTTCGGCGCGCTGTTGGAAGAGCAGAGCAAagag                                                                   >  2:91107/1‑124 (MQ=255)
                                                                   tCTGCACGCTCAGCGTCAGAAACTGCACGGTTATCTGCCAAGCCGTCAGCCGAACTTCCCCGAGAAGCTTGAGCTGCCGAGCCTGCAAGACTTCGGCGCGCTGTTGGAAGAGCAGAGCAAAGAGATCTCTACCACTATc                                                    >  2:10016/1‑139 (MQ=255)
                                                                   tCTGCACGCTCAGCGTCAGAAACTGCACGGTTATCTGCCAAGCCGTCAGCCGAACTTCACCGAGAAGCTTGAGCTGCCGAGCCTGCAAGACTTCGGCGCGCTGTTGGAAGAGCAGAGCAAAGAGATCTCTACCACTATc                                                    >  2:530557/1‑139 (MQ=255)
                                                                   tCTGCACGCTCAGCGTCAGAAACTGCACGGTTATCTGCCAAGCCGTCAGCCGAACTTCACCGAGAAGCTTGAGCTGCCGAGCCTGCAAGACTTCGGCGCGCTGTTGGAAGAGCAGAGCAAAGAGATCTCTACCACTATc                                                    >  2:358600/1‑139 (MQ=255)
                                                                     tGCACGCTCAGCGTCAGAAACTGCACGGTTATCTGCCAAGCCGTCAGCCGAACTTCCCCGAGAAGCTTGAGCTGCCGAGCCTGCAAGACTTCGGCGCGCTGTTGGAAGAGCAGAGCAAAGAGATCTCTACCACTATCGc                                                  >  1:117486/1‑139 (MQ=255)
                                                                              aGCGTCAGAAACTGCACGGTTATCTGCCAAGCCGTCAGCCGAACTTCCCCGAGAAGCTTGAGCTGCCGAGCCTGCAAGACTTCGGCGCGCTGTTGGAAGAGCAGAGCAAaga                                                                    >  1:133306/1‑112 (MQ=255)
                                                                              aGCGTCAGAAACTGCACGGTTATCTGCCAAGCCGTCAGCCGAACTTCCCCGAGAAGCTTGAGCTGCCGAGCCTGCAAGACTTCGGCGCGCTGTTGGAAGAGCAGAGCAAaga                                                                    <  2:133306/112‑1 (MQ=255)
                                                                                    aGAAACTGCACGGTTATCTGCCAAGCCGTCAGCCGAACTTCCCCGAGAAGCTTGAGCTGCCGAGCCTGCAAGACTTCGGCGCGCTGTTGGAAGAGCAGAGCAAAGAGATCTCTACCACTATCGCTTTCGTTCGTGctct                                   <  1:593892/139‑1 (MQ=255)
                                                                                         cTGCACGGTTATCTGCCAAGCCGTCAGCCGAACTTCACCGAGAAGCTTGAGCTGCCGAGCCTGCAAGACTTCGGCGCGCTGTTGGAAGAGCAGAGCAAAGAGATCTCTACCACTATCGCTTTCGTTCGTGCTCTGAACg                              <  1:460631/139‑1 (MQ=255)
                                                                                            cACGGTTATCTGCCAAGCCGTCAGCCGAACTACCCCGAGAAGCTTGAGCTGCCGAGCCTg                                                                                                          <  2:699930/60‑1 (MQ=255)
                                                                                            cACGGTTATCTGCCAAGCCGTCAGCCGAACTACCCCGAGAAGCTTGAGCTGCCGAGCCTg                                                                                                          >  1:699930/1‑60 (MQ=255)
                                                                                                gTTATCTGCCAAGCCGTCAGCCGAACTTCCCCGAGAAGCTTGAGCTGCCGAGCCTGCAAGACTTCGGCGCGCTGTTGGAAGAGCAGAGCAAAGAGATCTCTACCACTATCGCTTTCGTTCGTGCTCTGAACGTGATGCt                       <  2:148599/139‑1 (MQ=255)
                                                                                                        ccAAGCCGTCAGCCGAACTTCCCCGAGAAGCTTGAGCTGCCGAGCCTGCAAGACTTCGGCGCGCTGTTGGAAGAGCAGAGCAAAGAGATCTCTACCACTATCGCTTTCGTTCGTGCTCTGAACGTGATGCTgaagaa                 <  1:765345/137‑1 (MQ=255)
                                                                                                          aaGCCGTCAGCCGAACTTCCCCGAGAAGCTTGAGCTGCCGAGCCTGCAAGACTTCGGCGCGCTGTTGGAAGAGCAGAGCAAAGAGATCTCTACCACTATCGCTTTCGTTCGTGCTCTGAACGTGATGCTGAAGAACAAg             <  1:324659/139‑1 (MQ=255)
                                                                                                           aGCCGTCAGCCGAACTTCCCCGAGAAGCTTGAGCTGCCGAGCCTGCAAGACTTCGGCGCGCTGTTGGAAGAGCAGAGCAAAGAGATCTCTACCACTATCGCTTTCGTTCGTGCTCTGAACGTGATGCTGAAGAACAAg             <  1:196530/138‑1 (MQ=255)
                                                                                                           aGCCGTCAGCCGAACTTCCCCGAGAAGCTTGAGCTGCCGAGCCTGCAAGACTTCGGCGCGCTGTTGGAAGAGCAGAGCAAAGAGATCTCTACCACTATCGCTTTCGTTCGTGCTCTGAACGTGATGCTGAAGAACAAg             >  2:196530/1‑138 (MQ=255)
                                                                                                           aGCCGTCAGCCGAACTTCCCCGAGAAGCTTGAGCTGCCGAGCCTGCAAGACTTCGGCGCGCTGTTGGAAGAGCAGAGCAAAGAGATCTCTACCACTATCGCTTTCGTTCGTGCTCTGAACGTGATGCTGAAGAACAAGt            <  2:117486/139‑1 (MQ=255)
                                                                                                           aGCCGTCAGCCGAACTTCCCCGAGAAGCTTGAGCTGCCGAGCCTGCAAGACTTCGGCGCGCTGTTGGAAGAGCAGAGCAAAGAGATCTCTACCACTATCGCTTTCGTTCGTGCTCTGAACGTGATGCTGAAGAACAAGt            <  1:54972/139‑1 (MQ=255)
                                                                                                           aGCCGTCAGCCGAACTTCACCGAGAAGCTTGAGCTGCCGAGCCTGCAAGACTTCGGCGCGCTGTTGGAAGAGCAGAGCAAAGAGATCTCTACCACTATCGCTTTCGTTCGTGCTCTGAACGTGATGCTGAAGAACAAGt            >  2:614786/1‑139 (MQ=255)
                                                                                                                   gCCGAACTTCCCCGAGAAGCTTGAGCTGCCGAGCCTGCAAGACTTCGGCGCGCTGTTGGAAGAGCAGAGCAAAGAGATCTCTACCACTATCGCTTTCGTTCGTGCTCTGAACGTGATGCTGAAGAACAAGTCGATCaaa    <  2:371155/139‑1 (MQ=255)
                                                                                                                     cGAACTTCCCCGAGAAGCTTGAGCTGCCGAGCCTGCAAGACTTCGGCGCGCTGTTGGAAGAGCAGAGCAAAGAGATCTCTACCACTATCGCTTTCGTTCGTGCTCTGAACGTGATGCTGAAGAACAAGTCGATCAAAGa  <  1:344129/139‑1 (MQ=255)
                                                                                                                     cGAACTTCCCCGAGAAGCTTGAGCTGCCGAGCCTGCAAGACTTCGGCGCGCTGTTGGAAGAGCAGAGCAAAGAGATCTCTACCACTATCGCTTTCGTTCGTGCTCTGAACGTGATGCTGAAGAACAAGTCGATCAAAGa  <  1:153450/139‑1 (MQ=255)
                                                                                                                       aaCTTCACCGAGAAGCTTGAGCTGCCGAGCCTGCAAGACTTCGGCGCGCTGTTGGAAGAGCAGAGCAAAGAGATCTCTACCACTATCGCTTTCGTTCGTGCTCTGAACGTGATGCTGAAGAACAAGTCGATCa      >  2:414989/1‑133 (MQ=255)
                                                                                                                             |                                                                                                                                  
TGTCTGATGCAGATATCGAAAAACTGCCGTACATCACCTTCCCGGAAGGTTCTGAAGAGCATACCTATCTGCACGCTCAGCGTCAGAAACTGCACGGTTATCTGCCAAGCCGTCAGCCGAACTTCACCGAGAAGCTTGAGCTGCCGAGCCTGCAAGACTTCGGCGCGCTGTTGGAAGAGCAGAGCAAAGAGATCTCTACCACTATCGCTTTCGTTCGTGCTCTGAACGTGATGCTGAAGAACAAGTCGATCAAAGA  >  NZ_CP009273/120789‑121044

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: